版塊導(dǎo)航: 正在加載中...

應(yīng)《網(wǎng)絡(luò)安全法》要求，自2017年10月1日起，未進(jìn)行實(shí)名認(rèn)證將不得使用互聯(lián)網(wǎng)跟帖服務(wù)。為保障您的帳號(hào)能夠正常使用，請(qǐng)盡快對(duì)帳號(hào)進(jìn)行手機(jī)號(hào)驗(yàn)證，感謝您的理解與支持！

24小時(shí)熱門版塊排行榜

北京石油化工學(xué)院2026年研究生招生接收調(diào)劑公告

返回列表

本帖產(chǎn)生 1 個(gè) QC強(qiáng)帖，點(diǎn)擊這里進(jìn)行查看

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

[交流] 【小卒個(gè)人文集】用NLO-Calculator來計(jì)算二階非線性光學(xué)響應(yīng)系數(shù)張量及其各向異性已有23人參與

用NLO-Calculator來計(jì)算二階非線性光學(xué)響應(yīng)系數(shù)張量及其各向異性
密度泛函小卒

通過Gaussian 09的CPDFT或CPHF計(jì)算超極化率，進(jìn)而再用NLO-Calculator來計(jì)算二階非線性光學(xué)響應(yīng)系數(shù)張量及其各向異性。這件事情已經(jīng)完成了。Gaussian 09是商業(yè)程序。NLO Calculator (Version 0.2)是我自己寫的程序，用的編程語言是Python 2.7。
可執(zhí)行程序可以在美國(guó)化學(xué)會(huì)網(wǎng)站下載。網(wǎng)址http://pubs.acs.org/doi/suppl/10.1021/jp3079293，程序很小，zip壓縮包的體積是4.9M，解壓直接能用，不用安裝。
論文是J. Phys. Chem. A, 2012, 116(41), pp 10249-10256。這個(gè)論文第一作者是我老婆，第二作者才是我。老婆指揮我干活兒，我只負(fù)責(zé)把老婆的要求翻譯成程序而已。
程序說明書也同時(shí)附在可執(zhí)行程序的文件夾里，解壓就能看到。還有一個(gè)例子，也在程序文件夾里
特別鳴謝中山大學(xué)的Denis，南開大學(xué)的coolrainbow！(排名不分先后，都感謝)
源碼公布如下：

CODE:

# This is a calculating program for second-order nonlinear

# optical parameters based on Gaussian 09 output file

# Dongdong Qi, NLO Calculator, Version 0.2

# University of Science and Technology Beijing (USTB)

# Beijing 100083, China

# E-mail: qdd-china@foxmail.com

import math

freq = {}

xxx = {}

xxy = {}

xyx = {}

xyy = {}

xxz = {}

xzx = {}

xyz = {}

xzy = {}

xzz = {}

yxx = {}

yyx = {}

yxy = {}

yyy = {}

yxz = {}

yzx = {}

yyz = {}

yzy = {}

yzz = {}

zxx = {}

zxy = {}

zyx = {}

zyy = {}

zxz = {}

zzx = {}

zyz = {}

zzy = {}

zzz = {}

x={}

y={}

z={}

tot={}

R1={}

R2={}

R3={}

R4={}

R5={}

R6={}

R7={}

R8={}

R9={}

R10={}

R11={}

beta2zzz={}

beta2xzz={}

betaHRS={}

DR={}

rho={}

betaJ1={}

betaJ3={}

phiJ1={}

phiJ3={}

vec={}

Bers={}

IHV={}

IHVjx={}

IHVjy={}

def Calc(file,line):

    # define global var

    global freq, xxx, xxy, xyx, xyy, xxz, xzx, xyz, xzy, xzz, yxx, yyx, yxy, \

           yyy, yxz, yzx, yyz, yzy, yzz, zxx, zxy, zyx, zyy, zxz, zzx, zyz, zzy, zzz

    global i, x, y, z, tot, R1, R2, R3, R4, R5, R6, R7, R8, R9, R10, R11, beta2zzz, \

           beta2xzz, betaHRS, DR, rho, betaJ1, betaJ3, phiJ3, phiJ1, ux, uy, uz, vec, \

           Bers, IHV, IHVjx, IHVjy

    # read log file

    freq[i] = float(line[49:60])

    file.readline()

    GL1 = file.readline() #xxx Line in Gausian output file

    GL2 = file.readline() #xxy/xyx Line in Gausian output file

    GL3 = file.readline() #xyy Line in Gausian output file

    GL4 = file.readline() #xxz/xzx Line in Gausian output file

    GL5 = file.readline() #xyz/xzy Line in Gausian output file

    GL6 = file.readline() #xzz Line in Gausian output file

    GL7 = file.readline() #yxx Line in Gausian output file

    GL8 = file.readline() #yyx/yxy Line in Gausian output file

    GL9 = file.readline() #yyy Line in Gausian output file

    GL10 = file.readline() #yxz/yzx Line in Gausian output file

    GL11 = file.readline() #yyz/yzy Line in Gausian output file

    GL12 = file.readline() #yzz Line in Gausian output file

    GL13 = file.readline() #zxx Line in Gausian output file

    GL14 = file.readline() #zxy/zyx Line in Gausian output file

    GL15 = file.readline() #zyy Line in Gausian output file

    GL16 = file.readline() #zxz/zzx Line in Gausian output file

    GL17 = file.readline() #zyz/zzy Line in Gausian output file

    GL18 = file.readline() #zzz Line in Gausian output file

    xxx[i] = -float(GL1[8:18])*math.pow(10,int(GL1[19:22]))

    xxy[i] = -float(GL2[8:18])*math.pow(10,int(GL2[19:22]))

    xyy[i] = -float(GL3[8:18])*math.pow(10,int(GL3[19:22]))

    xxz[i] = -float(GL4[8:18])*math.pow(10,int(GL4[19:22]))

    xyz[i] = -float(GL5[8:18])*math.pow(10,int(GL5[19:22]))

    xzz[i] = -float(GL6[8:18])*math.pow(10,int(GL6[19:22]))

    yxx[i] = -float(GL7[8:18])*math.pow(10,int(GL7[19:22]))

    yyx[i] = -float(GL8[8:18])*math.pow(10,int(GL8[19:22]))

    yyy[i] = -float(GL9[8:18])*math.pow(10,int(GL9[19:22]))

    yxz[i] = -float(GL10[8:18])*math.pow(10,int(GL10[19:22]))

    yyz[i] = -float(GL11[8:18])*math.pow(10,int(GL11[19:22]))

    yzz[i] = -float(GL12[8:18])*math.pow(10,int(GL12[19:22]))

    zxx[i] = -float(GL13[8:18])*math.pow(10,int(GL13[19:22]))

    zxy[i] = -float(GL14[8:18])*math.pow(10,int(GL14[19:22]))

    zyy[i] = -float(GL15[8:18])*math.pow(10,int(GL15[19:22]))

    zxz[i] = -float(GL16[8:18])*math.pow(10,int(GL16[19:22]))

    zyz[i] = -float(GL17[8:18])*math.pow(10,int(GL17[19:22]))

    zzz[i] = -float(GL18[8:18])*math.pow(10,int(GL18[19:22]))

    xyx[i] = -float(GL2[8:18])*math.pow(10,int(GL2[19:22]))

    xzx[i] = -float(GL4[8:18])*math.pow(10,int(GL4[19:22]))

    xzy[i] = -float(GL5[8:18])*math.pow(10,int(GL5[19:22]))

    yxy[i] = -float(GL8[8:18])*math.pow(10,int(GL8[19:22]))

    yzx[i] = -float(GL10[8:18])*math.pow(10,int(GL10[19:22]))

    yzy[i] = -float(GL11[8:18])*math.pow(10,int(GL11[19:22]))

    zyx[i] = -float(GL14[8:18])*math.pow(10,int(GL14[19:22]))

    zzx[i] = -float(GL16[8:18])*math.pow(10,int(GL16[19:22]))

    zzy[i] = -float(GL17[8:18])*math.pow(10,int(GL17[19:22]))

    # print Gaussian beta data

    print "\n============= Freq", i, "============================================="

    if i == 0 : print 'Freq:', freq[i], 'a.u.\nbeta( lambda = inf , Static)'

    else : print 'Freq:', freq[i], 'a.u. (', int(round(freq[i]*219474.6307)), \

         'cm-1 )\nbeta( lambda =', int(round(10000000.0/(freq[i]*219474.6307))), \

         'nm , Dynamic)'

    print "------------ Input -----------------------------------------------"

    print "beta(xxx) =", int(round(xxx[i])), "; beta(xxy) =", int(round(xxy[i])), \

          "; beta(xxz) =", int(round(xxz[i]))

    print "beta(xyx) =", int(round(xyx[i])), "; beta(xyy) =", int(round(xyy[i])), \

          "; beta(xyz) =", int(round(xyz[i]))

    print "beta(xzx) =", int(round(xzx[i])), "; beta(xzy) =", int(round(xzy[i])), \

          "; beta(xzz) =", int(round(xzz[i]))

    print "beta(yxx) =", int(round(yxx[i])), "; beta(yxy) =", int(round(yxy[i])), \

          "; beta(yxz) =", int(round(yxz[i]))

    print "beta(yyx) =", int(round(yyx[i])), "; beta(yyy) =", int(round(yyy[i])), \

          "; beta(yyz) =", int(round(yyz[i]))

    print "beta(yzx) =", int(round(yzx[i])), "; beta(yzy) =", int(round(yzy[i])), \

          "; beta(yzz) =", int(round(yzz[i]))

    print "beta(zxx) =", int(round(zxx[i])), "; beta(zxy) =", int(round(zxy[i])), \

          "; beta(zxz) =", int(round(zxz[i]))

    print "beta(zyx) =", int(round(zyx[i])), "; beta(zyy) =", int(round(zyy[i])), \

          "; beta(zyz) =", int(round(zyz[i]))

    print "beta(zzx) =", int(round(zzx[i])), "; beta(zzy) =", int(round(zzy[i])), \

          "; beta(zzz) =", int(round(zzz[i]))

    print "------------ Output ----------------------------------------------"

    # Start calculation

    # Step 1: beta(total)

    x[i] = ((xyy[i] + xzz[i] + yxy[i] + zxz[i] + yyx[i] + zzx[i]) / 3.0) +  xxx[i]

    y[i] = ((yxx[i] + yzz[i] + xyx[i] + zyz[i] + xxy[i] + zzy[i]) / 3.0) +  yyy[i]

    z[i] = ((zxx[i] + zyy[i] + xzx[i] + yzy[i] + xxz[i] + yyz[i]) / 3.0) +  zzz[i]

    print 'beta(x) =', int(round(x[i])), '; beta(y) =', int(round(y[i])), \

          '; beta(z) =', int(round(z[i]))

    tot[i] = math.sqrt(math.pow(x[i],2) + math.pow(y[i],2) + math.pow(z[i],2))

    print 'beta(total) =', int(round(tot[i]))

    # step 2:  and 

    R1[i] = (math.pow(xxx[i],2) + math.pow(yyy[i],2) + math.pow(zzz[i],2))

    R2[i] = (math.pow(xxy[i],2) + math.pow(xxz[i],2) + math.pow(yyx[i],2) + \

             math.pow(yyz[i],2) + math.pow(zzx[i],2) + math.pow(zzy[i],2))

    R3[i] = (xxx[i]*xyy[i] + xxx[i]*xzz[i] + yyy[i]*yxx[i] + yyy[i]*yzz[i] + \

             zzz[i]*zxx[i] + zzz[i]*zyy[i])

    R4[i] = (xyy[i]*yyx[i] + xzz[i]*zzx[i] + yxx[i]*xxy[i] + yzz[i]*zzy[i] + \

             zxx[i]*xxz[i] + zyy[i]*yyz[i])

    R5[i] = (xxx[i]*yyx[i] + xxx[i]*zzx[i] + yyy[i]*xxy[i] + yyy[i]*zzy[i] + \

             zzz[i]*xxz[i] + zzz[i]*yyz[i])

    R6[i] = (math.pow(yxx[i],2) + math.pow(zxx[i],2) + math.pow(xyy[i],2) + \

             math.pow(zyy[i],2) + math.pow(xzz[i],2) + math.pow(yzz[i],2))

    R7[i] = (xxy[i]*yzz[i] + xxz[i]*zyy[i] + yyx[i]*xzz[i] + yyz[i]*zxx[i] + \

             zzx[i]*xyy[i] + zzy[i]*yxx[i])

    R8[i] = (xyy[i]*xzz[i] + yxx[i]*yzz[i] + zxx[i]*zyy[i]) * 2.0

    R9[i] = (xxy[i]*zzy[i] + xxz[i]*yyz[i] + yyx[i]*zzx[i]) * 2.0

    R10[i] = (math.pow(xyz[i],2) + math.pow(xzy[i],2) + math.pow(yxz[i],2) + \

              math.pow(yzx[i],2) + math.pow(zyx[i],2) + math.pow(zxy[i],2))

    R11[i] = (xyz[i]*yxz[i] + xzy[i]*zxy[i] + yzx[i]*zyx[i]) * 2.0

    beta2zzz[i] = (R1[i]*15.0 + R2[i]*12.0 + R3[i]*6.0 + R4[i]*12.0 + \

                   R5[i]*12.0 + R6[i]*3.0 + R7[i]*4.0 + R8[i] + R9[i]*4.0 + \

                   R10[i]*2.0 + R11[i]*4.0)/105.0

    beta2xzz[i] = (R1[i]*3.0 + R2[i]*8.0 + R3[i]*4.0 - R4[i]*6.0 - \

                   R5[i]*6.0 + R6[i]*9.0 - R7[i]*2.0 + R8[i]*3.0 - \

                   R9[i]*2.0 + R10[i]*6.0 - R11[i]*2.0)/105.0

    betaHRS[i] = math.sqrt(beta2zzz[i] + beta2xzz[i])

    DR[i] = beta2zzz[i]/beta2xzz[i]

    if DR[i] < 1.5 :

        print "Be careful! DR is lower than 1.5."

    print " =", int(round(beta2zzz[i])), ";  =", \

          int(round(beta2xzz[i]))

    print "beta(HRS)(-2w;w,w) =", int(round(betaHRS[i])), "; DR =", DR[i]

    if DR[i] > 1.5 :

        rho[i] = math.sqrt(7.0*(9.0-DR[i])/(12.0*DR[i]-18.0))

        print "rho =", rho[i]

        betaJ1[i] = betaHRS[i]/math.sqrt((2.0/9.0)+(2.0*rho[i]*rho[i]/21.0))

        betaJ3[i] = rho[i]*betaJ1[i]

        phiJ3[i] = rho[i]/(1.0+rho[i])

        phiJ1[i] = 1.0-phiJ3[i]

        print "|beta(J=1)| =", betaJ1[i], "; |beta(J=3)| =", betaJ3[i]

        print "phi[beta(J=1)] =", phiJ1[i], "; phi[beta(J=3)] =", phiJ3[i]

    else :

        print "\n      ____ Attention ___________________________________\n"

        print "       If DR is lower than 1.5, the calculations of"

        print "       rho, |beta(J=1)|, |beta(J=3)|, phi[beta(J=1)],"

        print "       and phi[beta(J=3)] are unreliable. In addition,"

        print "       please check whether the half of the incident"

        print "       light freqency (lambda/2) is near to the"

        print "       absorption band, which is the direct reason of"

        print "       the resonance for SFG. For details, please"

        print "       check any NLO textbook."

        print "      __________________________________________________\n"

    # step 3: IHV list

    Bers[i] = 7.0 * beta2xzz[i] - beta2zzz[i]

    angle = -180

    while (angle < 181):

        IHV[i,angle] = (math.pow(math.cos(float(angle)*(2.0*math.pi)/360.0),4))*beta2xzz[i] + \

                       (math.pow(math.sin(float(angle)*(2.0*math.pi)/360.0),4))*beta2zzz[i] + \

                       (math.pow(math.cos(float(angle)*(2.0*math.pi)/360.0),2))*\

                       (math.pow(math.sin(float(angle)*(2.0*math.pi)/360.0),2))*Bers[i]

        IHVjx[i,angle] = IHV[i,angle] * math.cos(float(angle)*(2.0*math.pi)/360.0)

        IHVjy[i,angle] = IHV[i,angle] * math.sin(float(angle)*(2.0*math.pi)/360.0)

        angle = angle + 1

    print "If you need the I(2w,PsiV) values, please check them in I2wPsiV.txt."

    i=i+1

# Main Programe

print "       __________________ NLO Calculator _________________\n"

print "                       v 0.2 ( June 21, 2012 )\n"

print "       Theory: [1] J. Chem. Phys. 2012, 136, 024506."

print "               [2] Phys. Chem. Chem. Phys. 2008, 10, 6223.\n"

print "Cite:   Dongdong Qi, NLO Calculator, Version 0.2"

print "        University of Science and Technology Beijing (USTB)"

print "        Beijing 100083, China\nE-mail: qdd-china@foxmail.com"

filename = raw_input('Please input the name of Gaussian log/out file: ')

file = open(filename, 'r')

line = None

i = 0

freqN = input('How many special frequencies in this log/out file? ')

while (line != "EOF") and (i < freqN+1) :

    line = file.readline()

    if line.startswith(' Property number 6 -- Beta(w,w,-2w) frequency'):

        Calc(file,line)

file.close()

f = open('I2wPsiV.txt','w')

f.write("       __________________ NLO Calculator _________________\n\n")

f.write("                       v 0.2 ( June 21, 2012 )\n")

f.write("       Theory: [1] J. Chem. Phys. 2012, 136, 024506.\n")

f.write("               [2] Phys. Chem. Chem. Phys. 2008, 10, 6223.\n")

f.write("       ___________________________________________________\n")

f.write("\n___________________________________________________________________\n")

f.write("\nCite:   Dongdong Qi, NLO Calculator, Version 0.2")

f.write("        University of Science and Technology Beijing (USTB)")

f.write("        Beijing 100083, China\nE-mail: qdd-china@foxmail.com")

f.write("___________________________________________________________________\n\n")

i = -1

while (i < freqN):

    i = i+1

    f.write("Freq ")

    f.write(repr(i))

    f.write(": ")

    f.write(repr(freq[i]))

    f.write(" a.u.\nlambda = ")

    f.write(repr(int(round(freq[i]*219474.6307))))

    f.write(" cm-1 \n")

    if (i!=0): f.write(repr(int(round(10000000.0/(freq[i]*219474.6307)))))

    else: f.write("inf")

    f.write(" nm\n\n")

    if (i!=0): f.write("Dynamic beta(HRS)(-2w;w,w): ")

    else: f.write("Static beta(HRS)(-2w;w,w): ")

    f.write(repr(int(round(betaHRS[i]))))

    f.write("\n\n")

    anglex = -180

    while (anglex < 181):

        f.write("Psi = ")

        f.write(repr(anglex))

        f.write(" ; [Cartesian Coord.] I(2w,PsiV) = ")

        f.write(repr(int(round(IHV[i,anglex]))))

        f.write(" ; [Polar Coord.] I(2w,PsiV,x) = ")

        f.write(repr(int(round(IHVjx[i,anglex]))))

        f.write(" and I(2w,PsiV,y) = ")

        f.write(repr(int(round(IHVjy[i,anglex]))))

        if (anglex == 0): f.write(" ===== Horizontal [I(2w,HV)] =====")

        elif (anglex == 90): f.write(" ===== Vertical [I(2w,VV)] =====")

        elif (anglex == 45): f.write(" ===== Right-circular =====")

        elif (anglex == -45): f.write(" ===== Left-circular =====")

        f.write("\n")

        anglex = anglex + 1

    f.write("________________________________________\n\n")

f.close()

raw_input('Press  to exit.')

[ Last edited by yjcmwgk on 2013-1-19 at 18:09 ]

回復(fù)此樓

» 收錄本帖的淘帖專輯推薦

wien2k常見問題集錦	青椒不辣辣椒不青	vasp	量化計(jì)算
計(jì)算化學(xué)	wien2k-電磁	自己感覺好的東西

» 本帖已獲得的紅花（最新10朵）

gongxd325

popeye.cn

pyncell

» 猜你喜歡

2026湖北師范大學(xué)化學(xué)化工學(xué)院招收研究生調(diào)劑已經(jīng)有0人回復(fù)
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué)) 已經(jīng)有0人回復(fù)
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)? 已經(jīng)有160人回復(fù)
0703調(diào)劑已經(jīng)有2人回復(fù)
2026年博士研究生招生已經(jīng)有1人回復(fù)
304求調(diào)劑已經(jīng)有6人回復(fù)
武漢紡織大學(xué)全國(guó)重點(diǎn)實(shí)驗(yàn)室楊應(yīng)奎教授團(tuán)隊(duì)擬招17-19位碩博生！已經(jīng)有0人回復(fù)
材料與化工專碩接受跨專業(yè)調(diào)劑（重慶市屬高校）已經(jīng)有0人回復(fù)
武漢紡織大學(xué)全國(guó)重點(diǎn)實(shí)驗(yàn)室楊應(yīng)奎教授團(tuán)隊(duì)擬招17-19位碩博生！已經(jīng)有0人回復(fù)
上海工程技術(shù)大學(xué)生態(tài)與功能紡織品開發(fā)研究團(tuán)隊(duì)招收學(xué)碩專碩研究生已經(jīng)有0人回復(fù)

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

二階非線性光學(xué)效應(yīng)的表征已經(jīng)有6人回復(fù)
關(guān)于交換系數(shù)J的計(jì)算已經(jīng)有11人回復(fù)
使用Multiwfn分析Gaussian09的極化率、第一超極化率的輸出已經(jīng)有12人回復(fù)
求助，四維張量和三階張量是不是指同一個(gè)東西？已經(jīng)有9人回復(fù)
二階張量雙點(diǎn)乘三階張量已經(jīng)有13人回復(fù)
二階變系數(shù)差分方程已經(jīng)有11人回復(fù)
raw data的說明已經(jīng)有3人回復(fù)
關(guān)于線性系統(tǒng)的輸出響應(yīng)的問題已經(jīng)有6人回復(fù)
晶體的物理性與對(duì)稱性已經(jīng)有232人回復(fù)
求一個(gè)態(tài)求和（SOS）計(jì)算非線性光學(xué)系數(shù)的程序已經(jīng)有4人回復(fù)
【求助】拉曼極化率已經(jīng)有7人回復(fù)

1樓 2012-11-18 22:11:05

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

回帖支持 ( 顯示支持度最高的前 50 名 )

wuhen116

木蟲 (正式寫手)

應(yīng)助: 0 (幼兒園)
金幣: 6489.3
散金: 874
紅花: 4
帖子: 895
在線: 961.8小時(shí)
蟲號(hào): 674041
注冊(cè): 2008-12-13
性別: GG
專業(yè): 光譜分析

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

樓主很牛啊，樓主老婆（zhlj？）好像更牛啊，99年就有一作了（ACS），到現(xiàn)在都30幾篇文章了，大牛啊。

贊一下(1人)

回復(fù)此樓

15樓2012-11-20 10:00:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

引用回帖:

15樓: Originally posted by wuhen116 at 2012-11-20 10:00:38
樓主很牛啊，樓主老婆（zhlj？）好像更牛啊，99年就有一作了（ACS），到現(xiàn)在都30幾篇文章了，大牛啊。

那是重名……我和我老婆現(xiàn)在都是博三……99年那會(huì)兒我們還在接受義務(wù)教育

贊一下(1人)

回復(fù)此樓

» 本帖已獲得的紅花（最新10朵）

小范范1989

17樓2012-11-20 18:10:03

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

引用回帖:

23樓: Originally posted by 偉明 at 2012-12-04 21:54:08
我是想問你程序里的公式跟我們組一直用的這個(gè)是否一樣呢？還是這個(gè)系數(shù)選擇影響不大？

未命名.jpg

不一樣。我用的是動(dòng)態(tài)的解析的公式。你們這幾個(gè)公式是靜態(tài)的數(shù)值的公式�，F(xiàn)在國(guó)內(nèi)大部分確實(shí)還是在用靜態(tài)的數(shù)值的公式

這種解析的動(dòng)態(tài)的計(jì)算，現(xiàn)在貌似在國(guó)內(nèi)我暫時(shí)還是“獨(dú)此一家別無分店”

贊一下(1人)

回復(fù)此樓

24樓2012-12-05 08:49:55

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

普通回帖

ApolloYang

木蟲 (正式寫手)

應(yīng)助: 15 (小學(xué)生)
金幣: 2209.6
散金: 114
紅花: 4
帖子: 390
在線: 236.9小時(shí)
蟲號(hào): 582657
注冊(cè): 2008-07-25
性別: GG
專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)

樓主你很強(qiáng)！

回復(fù)此樓

2樓2012-11-18 22:17:38

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

諸葛神鷹

榮譽(yù)版主 (職業(yè)作家)

嶶笑殗葢蕜傷，瀋黙?釋①切

應(yīng)助: 98 (初中生)
貴賓: 11.445
金幣: 31032.6
散金: 2612
紅花: 30
沙發(fā): 7
帖子: 3689
在線: 889.1小時(shí)
蟲號(hào): 1815241
注冊(cè): 2012-05-14
專業(yè): 人格心理學(xué)
管轄: 環(huán)境

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

好牛啊，那天要計(jì)算么，找你教教怎么使用，呵呵

贊一下

回復(fù)此樓

3樓2012-11-18 22:19:20

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

引用回帖:

3樓: Originally posted by 諸葛神鷹 at 2012-11-18 22:19:20
好牛啊，那天要計(jì)算么，找你教教怎么使用，呵呵

根本不用教
你去ACS把程序下載下來
然后輸入文件名
然后回車
然后就等著好了
后面的事情程序包辦

贊一下

回復(fù)此樓

4樓2012-11-18 22:26:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

cenwanglai

榮譽(yù)版主 (知名作家)

老和山猥瑣派九段

應(yīng)助: 46 (小學(xué)生)
貴賓: 8.842
金幣: 7440.4
散金: 9654
紅花: 118
沙發(fā): 23
帖子: 5306
在線: 1961.4小時(shí)
蟲號(hào): 537452
注冊(cè): 2008-04-01
性別: GG
專業(yè): 理論和計(jì)算化學(xué)
管轄: 仿真模擬

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

不過，我不知道在說什么~

贊一下

回復(fù)此樓

5樓2012-11-18 23:15:26

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

gongxd325

至尊木蟲 (文壇精英)

應(yīng)助: 177 (高中生)
金幣: 37777.1
散金: 13187
紅花: 109
帖子: 34003
在線: 1844.8小時(shí)
蟲號(hào): 233962
注冊(cè): 2006-03-30
性別: MM
專業(yè): 理論和計(jì)算化學(xué)

送鮮花一朵

不錯(cuò)的夫妻檔

回復(fù)此樓

一個(gè)人做點(diǎn)好事并不難難的是一輩子做好事!!!

6樓2012-11-19 08:08:54

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yanrding

金蟲 (正式寫手)

應(yīng)助: 114 (高中生)
金幣: 192.9
散金: 601
紅花: 10
帖子: 594
在線: 207.5小時(shí)
蟲號(hào): 1000562
注冊(cè): 2010-04-19
性別: GG
專業(yè): 理論和計(jì)算化學(xué)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

mark一下。以后慢慢研究。

贊一下

回復(fù)此樓

華北電力大學(xué)（北京）數(shù)理學(xué)院，計(jì)算化學(xué)，凝聚態(tài)物理。歡迎合作！招碩士/博士研究生，招老師。

7樓2012-11-19 08:34:51

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

南飛探

禁蟲 (正式寫手)

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

本帖內(nèi)容被屏蔽

8樓2012-11-19 08:42:25

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

popeye.cn

榮譽(yù)版主 (職業(yè)作家)

聰聰水手

應(yīng)助: 3 (幼兒園)
貴賓: 2.172
金幣: 7204.4
散金: 3967
紅花: 26
沙發(fā): 47
帖子: 4831
在線: 570.5小時(shí)
蟲號(hào): 399505
注冊(cè): 2007-06-11
性別: GG
專業(yè): 電燈泡; 副業(yè):挖墻角
管轄: 教師之家

★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖

送鮮花一朵

我也要娶好多好多這樣的老婆.

贊一下

回復(fù)此樓

我趕著驢，你騎著馬，咱倆攜手闖天涯！-★--★--★-★-★--交頸千年，相隨萬里！--★-★-★-

9樓2012-11-19 08:48:04

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

yjcmwgk

禁蟲 (文壇精英)

密度泛函·小卒

QC強(qiáng)帖: 38
應(yīng)助: 165 (高中生)
貴賓: 16.999
金幣: 554.6
散金: 39880
紅花: 395
沙發(fā): 2
帖子: 10406
在線: 1528.1小時(shí)
蟲號(hào): 448368
注冊(cè): 2007-11-01
專業(yè): 理論和計(jì)算化學(xué)
管轄: 科研工具資源

引用回帖:

5樓: Originally posted by cenwanglai at 2012-11-18 23:15:26
不過，我不知道在說什么~

也許是太偏門了吧……偏門到只能發(fā)JPC-A的地步，在B刊和C刊里根本找不到對(duì)應(yīng)類型，只能選A刊的“General Theoretical Chemistry”類型

贊一下

回復(fù)此樓

10樓2012-11-19 09:31:01

已閱回復(fù)此樓關(guān)注TA 給TA發(fā)消息送TA紅花 TA的回帖

相關(guān)版塊跳轉(zhuǎn) 我要訂閱樓主 yjcmwgk 的主題更新

返回列表

普通表情龍兔虎貓高級(jí)回復(fù) (可上傳附件)

最具人氣熱帖推薦 [查看全部]		作者	回/看	最后發(fā)表

[考研] 求調(diào)劑 +4	QiMing7 2026-03-25	5/250	2026-03-29 21:10 by 唐沐兒
[考研] 調(diào)劑310 +12	溫柔的晚安 2026-03-25	13/650	2026-03-29 20:01 by 無際的草原
[考研] 化學(xué)工程085602 305分求調(diào)劑 +22	RichLi_ 2026-03-25	22/1100	2026-03-29 19:51 by 無際的草原
[考研] 299求調(diào)劑 +10	15188958825 2026-03-25	10/500	2026-03-29 17:51 by 王亮_大連醫(yī)科大
[考研] 085600，專業(yè)課化工原理，321分求調(diào)劑 +5	大饞小子 2026-03-28	5/250	2026-03-29 08:56 by qingfeng258
[考研] 本科雙非材料，跨考一志愿華電085801電氣，283求調(diào)劑，任何專業(yè)都可以 +6	芝士雪baoo 2026-03-28	8/400	2026-03-29 08:16 by 松花缸1201
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +5	0857求調(diào)劑. 2026-03-24	6/300	2026-03-28 22:04 by zhq0425
[考研] 本科新能源科學(xué)與工程，一志愿華理能動(dòng)285求調(diào)劑 +3	AZMK 2026-03-27	5/250	2026-03-28 16:19 by xxxsssccc
[考研] 265求調(diào)劑 +8	小木蟲085600 2026-03-27	8/400	2026-03-27 22:16 by 無際的草原
[考研] 07化學(xué)280分求調(diào)劑 +10	722865 2026-03-23	10/500	2026-03-27 15:51 by Plutoqq
[考研] 085600材料與化工調(diào)劑 +10	A-哆啦Z夢(mèng) 2026-03-23	16/800	2026-03-27 15:13 by caszguilin
[論文投稿] Journal of Mechanical Science and Technology +3	Russ_ss 2026-03-25	5/250	2026-03-27 10:49 by 陸小果畫大餅
[考研] 304材料求調(diào)劑 +4	鐘llll 2026-03-26	4/200	2026-03-27 03:42 by wxiongid
[考研] 調(diào)劑 +4	柚柚yoyo 2026-03-26	4/200	2026-03-26 20:43 by fmesaito
[考研] 289求調(diào)劑 +17	碩星赴 2026-03-23	17/850	2026-03-26 16:18 by 不吃魚的貓
[考研] 一志愿南京郵電大學(xué) 288分材料考研求調(diào)劑 +3	jl0720 2026-03-26	3/150	2026-03-26 13:39 by zzll406
[考研] 化學(xué)調(diào)劑一志愿上海交通大學(xué)336分-本科上海211 +4	小魚愛有機(jī) 2026-03-25	4/200	2026-03-26 10:19 by aa331100
[考研] 293求調(diào)劑 +7	加一一九 2026-03-24	7/350	2026-03-25 12:02 by userper
[考研] 材料考研調(diào)劑生 +3	黃粱一夢(mèng)千年 2026-03-24	3/150	2026-03-24 17:00 by barlinike
[考研] 328求調(diào)劑 +4	LHHL66 2026-03-23	4/200	2026-03-23 14:55 by lbsjt

24小時(shí)熱門版塊排行榜

[交流] 【小卒個(gè)人文集】用NLO-Calculator來計(jì)算二階非線性光學(xué)響應(yīng)系數(shù)張量及其各向異性 已有23人參與

» 收錄本帖的淘帖專輯推薦

» 本帖已獲得的紅花（最新10朵）

» 猜你喜歡

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:

» 本帖已獲得的紅花（最新10朵）

[交流] 【小卒個(gè)人文集】用NLO-Calculator來計(jì)算二階非線性光學(xué)響應(yīng)系數(shù)張量及其各向異性已有23人參與

» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助: