第一篇:DP優(yōu)化 學習小結
DP優(yōu)化 學習小結
2010-08-27 15:55
這兩天簡單學習了一些動態(tài)規(guī)劃的優(yōu)化技巧,感覺自己的數(shù)學實在太弱了,雖然大概知道怎么用,但很多證明自己理解的還不是很透。
關于用單調隊列優(yōu)化的動態(tài)規(guī)劃以前整理過一些,可以看這里:【單調隊列 學習小結】
四邊形不等式優(yōu)化
主要是對于dp(i,j)的決策點 s(i,j),通過單調性證明 s(i,jpo[p2].y)/(po[p1].xsum[arr[head]]);po[i].x = a * sum[i];
po[i].y = dp[i] + a * sum[i] * sum[i]po[p2].y)/(po[p1].xf[p2]f[p2]f[arr[head]]f[arr[head]]sum[p2]p2);
}
inlinelonglong SS(int p1, int p2){
return num[p1 +1]sum[p1] + p1 * num[p1 +1]sum[p2] + p2 * num[p2 +1]);
}
inlinebool rate(int p1, int p2, int p3){
if(GG(p1, p2)* SS(p2, p3)>= GG(p2, p3)* SS(p1, p2))returnfalse;elsereturntrue;
}
int main(int argc, char** argv){
while(scanf(“%d%d”, &n, &m)!= EOF){
for(int i =1;i <= n;i++)num[i] = next_int();
sort(num +1, num + n +1);
sum[0] =0;
for(int i =1;i <= n;i++)sum[i] = sum[i-1] + num[i];dp[0] =0;
int head =0, rear =0;
arr[++rear] =0;
for(int i = m;i <= n;i++){
for(int k = head +1;k <= rear;k++){
if(cal(i, arr[k])<= cal(i, arr[head]))head +=1;elsebreak;
}
dp[i] = cal(i, arr[head]);
if(im +1;
while(head +2<= rear){
if(rate(arr[rear-2], arr[rear-1], arr[rear])==false){arr[rear-1] = arr[rear];
rear-=1;
} elsebreak;
}
}
cout << dp[n] << endl;
}
return(EXIT_SUCCESS);
}
第二篇:《IT職業(yè)行為優(yōu)化》學習小結
《IT職業(yè)行為優(yōu)化》學習小結
世界在不斷地向前發(fā)展,日積月累,一直到現(xiàn)在,在這日新月異的社會里,IT行業(yè)也在以驚人的速度在發(fā)展,并引領著時代的潮流。
這學期的IT課程的時間雖然不長,但是老師又給我們輸入了許多新的元素,我們的思想意識也慢慢的成熟了許多,老師的講解讓我們懂得如何才能才能成為一名合格的IT人員。《IT職業(yè)行為優(yōu)化》這門 課它讓我學到了很多以前沒學過的東西,讓我了解到很多前人在職業(yè)生涯的親身經(jīng)歷,覺得受益匪淺。它讓我在以后的職業(yè)生涯中增加了自信的籌碼,在以后的職業(yè)中多了一份職業(yè)技巧。下面是我通過學習體會出的感想:
一、我們應該有團結合作的精神
通過實驗“同舟共濟”,讓我明白了團結是職場必不可少的職業(yè)素養(yǎng)。為了培養(yǎng)團隊精神,我們應該在課后積極的參加各種社會團體的工作。與他人分工合作、分享成果、互助互惠的過程中,我們可以體會團隊精神的重要性。
有首歌唱的好“團結就是力量”。團隊合作往往能激發(fā)出團體不可思議的潛力,集體協(xié)作干出的成果往往能超過成員個人業(yè)績的總和。正所謂“同心山成玉,協(xié)力土變金?!币粋€團體,如果組織渙散,人心浮動,人人自行其是,甚至搞“窩里斗”,何來生機與活力?又何談干事創(chuàng)業(yè)?在一個缺乏凝聚力的環(huán)境里,個人再有雄心壯志,再有聰明才智,也不可能得到充分發(fā)揮!只有嚴密有序的集體組織和高效的團隊協(xié)作,懂得團結協(xié)作克服重重困難,甚至創(chuàng)造奇跡。當今社會,隨著知識經(jīng)濟時代的到來,各種知識、技術不斷推陳出新,競爭日趨緊張激烈,社會需求越來越多樣 化,使人們在工作學習中所面臨的情況和環(huán)境極其復雜。在很多情況下,單靠個人能力已很難完全處理各種錯綜復雜的問題并采取切實高效的行動。所有這些都需要 人們組成團體,并要求組織成員之間進一步相互依賴、相互關聯(lián)、共同合作,建立合作團隊來解決錯綜復雜的問題,并進行必要的行動協(xié)調,開發(fā)團隊應變能力和持 續(xù)的創(chuàng)新能力,依靠團隊合作的力量創(chuàng)造奇跡。
二、要懂得表達技巧和理解含義
通過撕紙實驗,我了解了表達與理解的重要性。一句話可以產(chǎn)生很多種不同的效果,只有把話表達清楚,才能夠讓別人理解。
語言表達能力是現(xiàn)代大學生必備的基本素質之一。在和別人說話的過程中用詞必須準確,不能產(chǎn)生歧義。我們可以多聽、多讀、多說、多寫的方法逐漸的糾正我們的錯誤。
由于經(jīng)濟的迅速發(fā)展,人們之間的交往密切,語言表達能力必須加強,這是作為現(xiàn)代大學生的必備之物。我們不僅要有新的思想和見解,還要在別人面前充分的表達出來。
三、要敢于面對挫折
任何一個活著的人,都會有或大或小的困難,但不同的人對困難的理解程度不同,有人認為遇到困難,解決困難才有進步。社會也才能得到發(fā)展,而一些消極悲觀的人把困難當成挫折,遇到挫折采取逃避、放棄。我們想成功的人,應該正確對待挫折。人生沒有坦途,必須經(jīng)歷風浪,才能駛向成功的彼岸。失敗是成功之母,沒有搏斗,那便沒有成功。挫折,它不是絆腳石,它卻是成功的助跑器,所以用一顆平常的心去對待挫折,只要對自己有信心,那么自己就肯定充滿力量。不要過分的顧及別人,過分的去注意別人對自己的看法。要以無所謂、平常的心態(tài)去對待任何挫折。
四、要注意語言的優(yōu)化
言談舉止是一個人精神面貌的 體現(xiàn),要開朗、熱情,讓人感覺隨和親切,平易近人,容易接觸。
在語言方面,交談總的要求是:文明、禮貌、準確。語言是交談的載體,交談者對它應當高度重視,精心斟酌,這是不言而喻的。在語言文明方面,我們要做到:不講粗話、臟話、黑話、葷話、怪話和氣話。在語言禮貌方面,我們應該使用一些禮貌用語,比如說:您好、請、謝謝、對不起、再見等等,這些詞的使用時博得他人好感與體諒的最為簡單易行的做法。在語言準確方面,要注意以下的問題:
1、發(fā)音要準確;
2、語速要適度;
3、口氣要謙和;
4、內容要簡明;
5、土話要少用;
6、外話要慎用。
只有使用文明、禮貌、準確的語言和他人交談,才能促進彼此的溝通。
五、要有清晰的責任感
在我們的生活中,有些事情我們可以不去做,但責任要求我們去做,甚至責任要求我們完成一些以我們的能力很難完成的事情。如果你做到了,得到的不僅僅是心理上的坦蕩和安然,你的精神和責任意識會感染別人,然后別人會因為你的影響而更有責任感。
在一個企業(yè)中,并不是所有的員工都能對自己的工作有強烈的責任感,但是如果他周圍的同事、整個企業(yè)環(huán)境都充滿一種責任的氛圍,那么他也會深受感染,進而能夠承擔起自己的責任。因為,他發(fā)現(xiàn),承擔責任并不是件很困難和痛苦的事情;相反,擔當起責任會給他一種驕傲的感覺,因為他在這個企業(yè)中同樣是重要的、不可或缺的。與其逃避責任,不如勇敢地承擔起來,說不定你的勇敢會成為你成功的契機。
雖然多數(shù)人是能夠明白自己的責任的,但是堅持不懈地承擔下去卻不是很多人能夠做到的。有些時候,堅持就是一種責任,尤其是在最困難的時候。
具有清晰的責任感是現(xiàn)代大學生的必備行為,是成為IT職業(yè)人員必備要素。沒有清晰的責任感的人不能算完人,他將不能在社會上好好的生存,也不可能被人錄用。所以我們應從此刻做起,好好的培養(yǎng)自己的責任心。
IT職業(yè)行為優(yōu)化課對我來說,對我的影響極深,除了以上以上五點以外,還有很多,我就不一一說明了。只要我們記住要想成為一個合格的復合型的IT人才,不僅要把書本學的知識連在一起使用,還要不斷的從生活汲取經(jīng)驗就行了。
第三篇:優(yōu)化小結
優(yōu)化第一步:
確定分子構型,可以根據(jù)對分子的了解通過GVIEW和CHEM3D等軟件來構建,但更多是通過實驗數(shù)據(jù)來構建(如根據(jù)晶體軟件獲得高斯直角坐標輸入文件,軟件可在大話西游上下載,用GVIEW可生成Z-矩陣高斯輸入文件),需要注意的是分子的原子的序號是由輸入原子的順序或構建原子的順序決定來實現(xiàn)的,所以為實現(xiàn)對稱性輸入,一定要保證第一個輸入的原子是對稱中心,這樣可以提高運算速度。我算的分子比較大,一直未曾嘗試過,希望作過這方面工作的朋友能補全它。以下是從本論壇,大話西游及宏劍公司上下載的帖子。將鍵長相近的,如
B12 1.08589 B13 1.08581 B14 1.08544 鍵角相近的,如
A6 119.66589 A7 120.46585 A8 119.36016 二面角相近的 如
D10-179.82816 D11-179.71092 都改為一致,聽說這樣可以減少變量,提高計算效率,是嗎? 在第一步和在以后取某些鍵長鍵角相等,感覺是一樣的。
只是在第一步就設為相等,除非有實驗上的證據(jù),不然就是純粹的憑經(jīng)驗了。在前面計算的基礎上,如果你比較信賴前面的計算,那么設為相等,倒還有些依據(jù)。但是,設為相等,總是冒些風險的。對于沒有對稱性的體系,應該是沒有絕對的相等的。
或許可以這么試試: 先PM3,再B3LYP/6-31G.(其中的某些鍵長鍵角設為相等),再B3LYP/6-31G(放開人為設定的那些鍵長鍵角相等的約束)。
比如鍵長,鍵角,還有是否成鍵的問題,Gview看起來就是不精確,不過基本上沒問題,要是限制它們也許就有很大的問題,能量上一般會有差異,有時還比較大
如果要減少優(yōu)化參數(shù),不是僅僅將相似的參數(shù)改為一致,而是要根據(jù)對稱性,采用相同 的參數(shù)。例如對苯分子分子指定部分如下:
C
C 1 B1
C 2 B2 1 A1
C 3 B3 2 A2 1 D1
C 4 B4 3 A3 2 D2
C 1 B5 2 A4 3 D3
H 1 B6 2 A5 3 D4
H 2 B7 1 A6 6 D5
H 3 B8 2 A7 1 D6
H 4 B9 3 A8 2 D7
H 5 B10 4 A9 3 D8
H 6 B11 1 A10 2 D9
B1 1.395160
B2 1.394712
B3 1.395427
B4 1.394825
B5 1.394829
B6 1.099610
B7 1.099655
B8 1.099680
B9 1.099680
B10 1.099761
B11 1.099604
A1 120.008632
A2 119.994165
A3 119.993992
A4 119.998457
A5 119.997223
A6 119.980770
A7 120.012795
A8 119.981142
A9 120.011343
A10 120.007997
D1-0.056843
D2 0.034114
D3 0.032348
D4-179.972926
D5 179.953248
D6 179.961852
D7-179.996436
D8-179.999514
D9 179.989175
參數(shù)很多,但是通過對稱性原則,并且采用亞原子可以將參數(shù)減少為:
X
X 1 B0
C 1 B1 2 A1
C 1 B1 2 A1 3 D1
C 1 B1 2 A1 4 D1
C 1 B1 2 A1 5 D1
C 1 B1 2 A1 6 D1
C 1 B1 2 A1 7 D1
H 1 B2 2 A1 8 D1
H 1 B2 2 A1 3 D1
H 1 B2 2 A1 4 D1
H 1 B2 2 A1 5 D1
H 1 B2 2 A1 6 D1
H 1 B2 2 A1 7 D1
B0 1.0
B1 1.2
B2 2.2
A1 90.0
D1 60.0
對于這兩個工作,所用的時間為57s和36s,對稱性為C01和D6H,明顯后者要遠遠優(yōu)于前者。
好,我真的很想知道這句話“你手工把坐標寫好輸入,輸入為C2對稱性”是什么意思?老是見到版面上有人說對輸入施加對稱性限制,就是不明白是怎么回事。您能詳細解釋一下嗎?有例子就更好了。麻煩了謝謝!
比如,你要輸入CO2,那么你這么輸入 C 0.0 0.0 0.0 O 0.0 0.0 1.203212313123 O 0.0 0.0 1.203212313123 寫得很精確,那就是對稱性輸入了。這樣g98可自動判定為D∞h 如果你輸入的差不多,比如: C 0.0 0.0 0.0 O 0.0 0.0 1.203212313123 O 0.0 0.0 1.204212313123 精度不夠,那么g98就只能判定為 C∞h了。
再給你舉個例子,比如甲烷的Td點群的輸入坐標,考慮對稱性,把C放在直角坐標的原點,四個氫放在四個角的頂點。C 0.0.0.H 1.1.1.H-1.-1.1.H 1.-1.-1.H-1.-1 1.這樣出來的結果點群判斷就會是正確的,當然這里的1我是隨便取的一個值,只是來說明這個問題而矣。
謝謝二位指教,也就是說第一個輸入的原子應該是對稱中心,其它輸入的原子應保持數(shù)值上的對應性。那對于用Gview或3D生成的分子構型,也可以根據(jù)這個原則將處于對稱位置的原子的坐標改為一致了?如果分子很大(如我這次算的C14H14N4S20)有三個環(huán),整體性的對稱輸入就有些麻煩,是不是可以先生成分子構型,然后利用對稱性對它的笛卡兒坐標或矩陣坐標進行修改,將相近的坐標的數(shù)據(jù)改為一致?
“第一個輸入的原子應該是對稱中心,其它輸入的原子應保持數(shù)值上的對應性”沒錯,做到這點就很方便了,“對于用Gview或3D生成的分子構型,也可以根據(jù)這個原則將處于對稱位置的原子的坐標改為一致”那還不如用手動輸入,大體系很難辦了,如果體系太大,不用保持嚴格的對稱性吧.有機分子,原子多了以后構型就沒有什么對稱性了。
但對于有些系統(tǒng),比如 C60, 你當然得按照它的實際的對稱性來輸入。如果你希望你的計算按照某種對稱性來做,那就一定要手工輸入嚴格的 滿足對稱性的坐標.用軟件生成的坐標,不是嚴格滿足對稱性,即使 用symm=loose,有時也不一定能判斷出來。
各位大蝦是用什么軟件建模的呀,我是用SGI里sybyl軟件進行構像搜索后取最低能量構像導入高斯程序的;有的是用chem3D直接建模,這樣子有沒有什么問題啊,應該怎么做,謝謝
通常也就是這些方法,建議還是認真的核實分子結構的是否合理,然后再試著增大迭代圈數(shù),降低收斂判據(jù).用什么軟件建模并不是問題,對于我們量化計算,如vliant所說,核實分子的合理構型才是最重要的,包 括對個原子間的成鍵分析。就算不同的軟件,也不可能保證給出的構型合理,因此我們一般都提倡手動輸入坐標,就是為了加強對構型的理解,只有在這個基礎上使用軟件建模才能盡可能保證不出錯。一個合理的計算結果不能寄希望于軟件。
我們都知道常用的分子坐標輸入有: 內坐標 直角坐標 混合坐標
分子構型輸入方法的不同會對計算過程產(chǎn)生影響么?或者說,對于哪些體系,用哪種方法更好優(yōu)化? “混合坐標”怎么講?
我通常用內坐標,它能很好的控制對稱性,且時間要短。
“混合坐標”是直角坐標和內坐標的混合使用,高斯手冊上有講這部分內容。有看過講哪種體系使用哪種方法更好優(yōu)化,好像這都是憑經(jīng)驗的結果。
在“免費資源”論壇中有本《計算化學》,其中第八章(67-72頁)好像提到你所關心的問題,你可以下來看一下,或許有所幫助。
如果你使用缺省的優(yōu)化方案且初始幾何構型完全相同,那么,不論你使用何種坐標輸入方式,優(yōu)化結果應該是一樣的。簡言之,優(yōu)化結果只取決于優(yōu)化方案和初始幾何構型,與坐標錄入方式無甚關聯(lián)。
于對稱性限制的優(yōu)化來說,顯然以內坐標輸入,且選用opt=z-matrix可以減少變量,使得優(yōu)化效率提高。免費資源中《gauss入門》有介紹。
您冗余坐標是指有過剩的內坐標。冗余內坐標算法可以用于結構優(yōu)化。
優(yōu)化第二步:
對于大分子,一般要先預優(yōu)化(如用PM3或AM1),有時是用部分優(yōu)化,然后才用能達到實驗精度方法和基組進行優(yōu)化,就我所了解大多數(shù)情況下,6-31G(D)基組就差不多。
不好意思,我是新手,剛剛接觸guassian,最近要算兩個較大有機分子中原子的靜電荷急相關能量和幾何性質,分子中已經(jīng)超過了80個原子(包括氫),用chemdraw畫出圖形,在chem3d中生成的input文件,可是,軟件報錯——2070~~ 好象說是連接(last link)錯誤,怎么回事,哪里出錯了,各位指點一下~ 怎么,先用pm3預優(yōu)化就可以了?? 我的意思是,初試構象對計算量有這么大的影響么??
我以前也遇到過相似情況,同一個分子,相同的計算條件,用gview畫出的模型進行計算報2070的錯,用chemdraw畫的模型,chem3d生成的input文件進行計算就OK,完成可以計算完成。我還以為你畫分子式的時候出什么問題了呢~~
只想優(yōu)化H原子,其他重原子不做優(yōu)化。怎樣進行分子坐標輸入,下面的坐標輸入方法對嗎?就是僅H原子的參數(shù)設置成Vairable.%chk=sample #p HF/6-31G* opt
Title line 0 1 cu 7.577854 5.926989 4.935050 cu 9.091644 8.489386 4.644861 cu 10.605857 5.766733 4.988313 p 10.432217 3.600341 5.598077 p 12.551342 6.888289 4.775263 p 7.164457 9.584280 4.209578 p 11.125158 9.457819 4.661391 p 7.503791 3.715690 5.278502 p 5.670551 7.057133 4.652208 s 9.100712 6.979012 6.564150 s 9.088921 6.492596 3.113104 n 8.884097 3.048614 5.950712 n 5.827205 8.602518 4.003874 n 12.406788 8.511755 5.225239 H 1 B14 2 A13 7 D12 H 1 B15 2 A14 7 D13 H 1 B16 2 A15 7 D14 H 1 B17 2 A16 7 D15 H 2 B18 1 A17 3 D16 H 2 B19 1 A18 3 D17 H 2 B20 1 A19 3 D18 H 2 B21 1 A20 3 D19 H 3 B22 1 A21 2 D20 H 3 B23 1 A22 2 D21 H 3 B24 1 A23 2 D22 H 3 B25 1 A24 2 D23 H 1 B26 2 A25 7 D24 H 2 B27 1 A26 9 D25 B14 2.990841 B15 2.990841 B16 2.985833 B17 2.985833 B18 3.009872 B19 3.009872 B20 2.998051 B21 2.998051 B22 3.015712 B23 3.015712 B24 3.002042 B25 3.002042 B26 3.489189 B27 3.477113 A13 161.290094 A14 146.871331 A15 108.517371 A16 96.114237 A17 106.549047 A18 100.057030 A19 155.415250 A20 145.871488 A21 163.425916 A22 138.991013 A23 105.076579 A24 99.387878 A25 66.120736 A26 65.475681 D12 64.540783 D13-53.098055 D14-163.347541 D15 152.573174 D16 152.164746 D17-162.671553 D18-56.766147 D19 42.781046 D20 46.122424 D21-46.534404 D22 150.210575 D23-164.670929 D24-78.866142 D25 109.837455 你如果要這么做的話,必須加opt=Z-matrix,這樣的話它就只優(yōu)化后面的H,前面的直角坐標就不會優(yōu)化了。你只用了opt的話,那就是全優(yōu)化了 你也可以通過molredunt來實現(xiàn)部分優(yōu)化 固定原子的方法有許多,最簡單的是這樣的: #rhf/sto-3g opt nosymm test the two C atoms are frozen 0,1 8 0-1.000000 0.000000 0.000000 6-1 0.000000 1.000000 0.000000 6-1 0.000000-1.000000 0.000000 8 0 1.000000 0.000000 0.000000 我不知道你的那種方法行不行,但我這個肯定可以。前面我都貼出來了:)這可是我的心血噢。========= opt=z-matrix的部分優(yōu)化。%Chk=reac-hf3 #p RHF/STO-3G opt(z-matrix,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 x2 y2 z2 H 0 x3 y3 z3 H 0 x4 y4 z4 O-1 2.16296 0.18089-0.15503 H 0 x6 y6 z6 H 0 x7 y7 z7 H 0 x8 y8 z8 x2=-2.2088 y2=0.39705 z2=0.46255 x3=-2.14589 y3=-0.9279 z3=-0.4899 x4=-1.88451 y4=-1.06386 z4=1.11629 x6=1.36256 y6=0.28882 z6=-0.67401 x7=2.32973 y7=0.98241 z7=0.34634 x8=2.90378 y8=0.01009 z8=-0.74115 *********************************************** opt=modredundant的部分優(yōu)化。
分子坐標部分可以為具體的數(shù)值,而不用為參數(shù)形式。%Chk=reac-hf3 #p RHF/STO-3G opt(modredundant,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 H 0 H 0 O-1 2.16296 0.18089-0.15503 H 0 H 0 H 0 1 5 F
感謝各位大俠的幫助!
看了之后有個不明白的地方,還請大俠不吝賜教!在opt=modredundant的部分優(yōu)化時最后一行的 1 5 F 是什么意思?
就是將第一和第五號原子的距離固定
問一個菜鳥的問題:N 后面的-1和H后面的0是什么意思?是像手冊上所說的把坐標增加所需要的值嗎? 見笑了!
-1指的是固定,0指的是優(yōu)化,前面helpme貼的輸入方法也是這樣的
部分優(yōu)化輸入格式之二
opt=z-matrix的部分優(yōu)化。%Chk=reac-hf3 #p RHF/STO-3G opt(z-matrix,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 x2 y2 z2 H 0 x3 y3 z3 H 0 x4 y4 z4 O-1 2.16296 0.18089-0.15503 H 0 x6 y6 z6 H 0 x7 y7 z7 H 0 x8 y8 z8 x2=-2.2088 y2=0.39705 z2=0.46255 x3=-2.14589 y3=-0.9279 z3=-0.4899 x4=-1.88451 y4=-1.06386 z4=1.11629 x6=1.36256 y6=0.28882 z6=-0.67401 x7=2.32973 y7=0.98241 z7=0.34634 x8=2.90378 y8=0.01009 z8=-0.74115 *********************************************** opt=modredundant的部分優(yōu)化。
分子坐標部分可以為具體的數(shù)值,而不用為參數(shù)形式。%Chk=reac-hf3 #p RHF/STO-3G opt(modredundant,maxcyc=500)freq=noraman scf(maxcyc=200)reactant fopt 1 1 N-1-1.75389-0.48105 0.31424 H 0 H 0 H 0 O-1 2.16296 0.18089-0.15503 H 0 H 0 H 0 1 5 F
凍結坐標的目的是什么了?我雖然會做但不太清楚Gaussian為什么會有這個功能.還有最近做了一個體系.有一個鍵如果不凍結的話.就優(yōu)化不出那個反應中間體.那個鍵會斷.我這樣凍結鍵優(yōu)化出來的結構計算頻率也無虛頻.它能看作為一個minima嗎?這個結構可以作為我反應機理的中間體嗎?我用溶劑化優(yōu)化結構可以在不凍結坐標的情況得到一個和這個結構很接近的結構.凍結坐標有時是有用的,舉幾個例子:
1,如優(yōu)化溶質分子和nH2O的supermolecule,周圍的水分子,可能結構變化不大,這時就可以固定水的鍵長或鍵角,減少工作量。
2,如用cluster模擬大塊固體,這時也許就要按照晶格參數(shù)固定原子鍵長和鍵角。3,又如不固定NH3-H-H2O中N-O鍵長是得不到質子傳遞勢壘。
4,還有在優(yōu)化穩(wěn)定構型或過渡態(tài)得時候也有時用到部分優(yōu)化得到接近得結果,再放開優(yōu)化,這樣逐步得到你想要得結果。
凍結坐標對一些大的體系,超分子體系好象是很有用的.好象做生物分子的反應時氨基酸殘基必須先固定.但我這個體系如果不固定,只要一放開,結構就跑了會優(yōu)化到其它構型上去.就算stepsize設得很小也沒有用.優(yōu)化第三步:
優(yōu)化結果分析,一般要進行頻率計算,無虛頻,表示是極小值;有虛頻,表示是過渡態(tài)。如果出現(xiàn)虛頻,可根據(jù)虛頻的簡正模式進行調整以獲得極小值,我還沒有碰到這種情況,希望各位朋友能給些詳細的資料。
如前信所說,鍵長鍵角有差別是很正常的,因為氣相分子和晶體內分子的結構不可能完 全吻合,只可能近似相同。
另外,對不同元素的限制主要看基組是否能計算該類元素,例如你所采用的方法中用了
6-31G(D)基組,該基組支持H-Kr的元素,即可以計算Br原子但是不能夠計算I原子。如果你所
用的基組不支持你所計算的某個原子可以采用擴展基組來進行計算,或者降低基組的精度。
我用Gaussian優(yōu)化一個分子結構,由于目前實驗上沒有辦法得到純凈的樣品來實驗驗證,那么我怎么判斷優(yōu)化得到的結構是正確的?
一些主要的鍵長和鍵角的實驗值或者別人計算的值是否和與你計算的分子類似? 這個判斷不好做。
所謂正確不正確,如果計算出來的頻率沒有虛頻,就表示你找到了一個局域穩(wěn)定點。要確定是不是全局穩(wěn)定點很難,只有通過別的方式比如HOMO-LUMO能級差等信息看你算出來態(tài)是不是穩(wěn)定。
也覺得高級的方法和大基組的結果更可信,謝謝各位高手的指點!再問一個問題,如果出現(xiàn)虛頻如何修改得到分子穩(wěn)定的結構
用GV看一下,根椐產(chǎn)生虛頻的振動模式調整結構
消去虛頻有什么技巧?我有一個計算結果,虛頻有三個,都只有-30~-20左右,如果結構改動過大,一不小心,虛頻就會更大。頭痛ing~~~
那你就把結構朝振動的方向改小一點試一下 你
確實據(jù)我了解,消除虛頻的主要方法是改變構型。
其次在計算上還可以嘗試:nosymm;加大循環(huán)次數(shù);提高收斂度;iop(1/8=1)等。其實我自己也是有這個問題,對于勢能面很平,較小的虛頻很難消除。首先,1在優(yōu)化時采用Scf(tight)的選項,增加收斂的標準。再去計算頻率。如果 還有虛頻,參見下一步。
2.對稱性的影響,很多情況下的虛頻是由于分子本身的對稱性造成的。這樣,在優(yōu)化時,如果必要,要將對稱性降低,還有,輸入文件有時是用內坐標。建議如果有虛頻的話,將內坐標改成直角坐標優(yōu)化。
3.如果上述方法還有虛頻,看一下虛頻,找到強度較大的,將在頻率中 產(chǎn)生的原子的振動坐標加到相應的輸入文件中。這樣,重新計算。直到 虛頻沒有。實際上,如果分子柔性較大,很難找到最低點,這是電子結構計算的問題,這種情況下,需要動力學的 東西,用構象搜尋的辦法解決。如:模擬退火,最陡下降法,淬火法等。將得到的能量最低的構象做 一般的電子結構計算,這樣,應當沒有問題。
不要講你還沒有得到最穩(wěn)定的結構,那么,是你的分子有問題,要么計算錯了,要么就是游離出現(xiàn)代計算的范疇 zixia上大蝦的回復: 將內坐標改成直角坐標優(yōu)化”這一方法沒有多大作用,我試過,還是一樣的結果 另外如果是根據(jù)振動頻率的模式來調整分子結構的話,還要注意調整完后的分子結構 對應的能量是不是比未調整之前有所降低,如果沒有,說明還是不行。
alwens兄所講的第三點,經(jīng)過我的猜測和與alwens的探討,終于知道應該怎么做: 分子結構輸入用直角坐標,將這些產(chǎn)生虛頻的振動坐標直接與原始坐標相加。alwens兄進一步提到: 這個方法有時候不太好。
不一定將坐標全部加上。有時候可以選擇加1/2等。主要的目的是去掉由分子對稱所形成的虛頻。
象上面的那個計算,完全可以再利用scf(tight)的選項消掉1-2個虛頻。強度很小。
又多學了一招,自我感覺很幸福中??^_^ 怎么把內坐標改成直角坐標啊 看你的輸出文件
Full point group CS NOp 2 Largest Abelian subgroup CS NOp 2 Largest concise Abelian subgroup CS NOp 2 Standard orientation:-------Center Atomic Atomic Coordinates(Angstroms)Number Number Type X Y Z-------1 46 0-0.242506 1.028743 2.204311 2 46 0-0.476137 2.619863 0.000000 3 46 0-0.242506 1.028743-2.204311 4 46 0-0.242506-1.618028-1.382208 5 46 0-0.242506-1.618028 1.382208 6 46 0-1.725240 0.349678 0.000000
這一部分就是,如果沒有,看 Input orientation:也行。
在opt的選項中加入“tight”或“saddle=0”,若是計算過度態(tài),應是“saddle=1”。但是,也不一定能消除虛頻!最好,將結構重新做一下,注意分子的構型,如CH3的三個H與相鄰基團是交叉,還是重疊!請教虛頻
請問如何知道虛頻的結構是什么,是不是頻率計算負值都可以認為是虛頻,那我用gaussianview 看有三個負值,接下來我想做過渡態(tài),用IRC,那我應該怎么弄呢?謝謝指導 你的輸出文件中這部分就是看是否有虛頻的
A“ A' A” Frequencies---65.9781 56.7053 59.7893 Red.masses--105.9032 105.9032 105.9032 Frc consts--0.2716 0.2006 0.2231 IR Inten--4.2436 0.0151 1.4107 Raman Activ--0.0000 0.0000 0.0000 Depolar--0.0000 0.0000 0.0000 Atom AN X Y Z X Y Z X Y Z 1 46-0.49 0.13-0.15-0.41-0.16 0.08-0.10 0.15 0.22 2 46 0.00 0.00 0.16 0.64-0.01 0.00 0.00 0.00-0.18 3 46 0.49-0.13-0.15-0.41-0.16-0.08 0.10-0.15 0.22 4 46-0.22 0.05 0.08 0.05-0.06 0.09-0.26-0.19 0.00 5 46 0.22-0.05 0.08 0.05-0.06-0.09 0.26 0.19 0.00 6 46 0.00 0.00 0.44-0.06 0.23 0.00 0.00 0.00-0.31 7 46 0.00 0.00-0.23-0.11 0.01 0.00 0.00 0.00 0.52 8 46 0.00 0.00-0.24 0.23 0.20 0.00 0.00 0.00-0.48 Frequencies 那一行數(shù)值為負,即有虛頻。有三個負值,說明找的不是過渡態(tài)。但是如果另外兩個虛頻很小,而另外一個虛頻很大,而且它的振動模式對應于你所計算的反應模式(比如反就前后某一鍵的拉長,縮短等),在計算所用的方法充許校正范圍內,也還可以認為是找到過渡態(tài)。但是一般要求能夠把另外兩個虛頻去掉。方法是要據(jù)振動模式,調整結構,再優(yōu)化。
般認為,過渡態(tài)僅有一個虛頻,中間體沒有虛頻!建議仔細閱讀《Exploring Chemistry with Electronic Structure Methos》和《G98W程序使用入門》,后者稍微簡單一點,廈大的免費資源里有!
有一個虛頻是過渡態(tài)的必要條件;有幾個虛頻肯定不行了,有時候有一個虛頻也不一定就找準了,做IRC成功就可靠了。另:復合物也會有小的虛頻,只要不大于100,也可認為是一極小值。
沒錯,有些體系可能存在several or many local minima,那么找到的具有一個虛頻的過渡態(tài)可能連接兩個loccal minima,而非連接反應物和產(chǎn)物,因此必須做IRC來判斷.關于虛頻的物理意義
自恨理論功底淺薄,所以特地在此向諸位高手求教:
在尋找過渡態(tài)計算中opt=ts, 需要計算頻率,如果有一個虛頻就說明優(yōu)化所得的構型是個 過渡態(tài)。有時在做部分限制優(yōu)化opt=minimal時,也能得到一個虛頻(一般發(fā)生在frozen 部分的化學鍵上)。我想問一下,虛頻的數(shù)值大小本身具有什么物理意義。后一種情況是否也可以認定為過渡態(tài)。
我的理解:
虛頻就是在這個振動方向上力常數(shù)是負的
過渡態(tài)之所以有一個虛頻,是因為它是一個鞍點,鞍點上有一個負的梯度的方向 分子沿這個方向振動時將轉化為反應物或產(chǎn)物。就象從山上掉下來,受到的力可以 認為是負的。
但并不是有一個虛頻就是過渡態(tài),還要在這個虛頻的振動方向上分別指向反應物和 產(chǎn)物
有一個虛頻并不一定是你想要的過渡態(tài),你要看虛頻的振動方向是不是對應著反應的方向。從概念上說,虛頻就是化學反應方向上的簡正坐標的二次偏導小于零。
另外,根據(jù)我的經(jīng)驗,這個虛頻的負值應該大一些,如果這個值太小,那也不足信。不知對不對,呵呵!
優(yōu)化輸出信息:
這說明優(yōu)化失敗了,沒有收斂,你可以嘗試:
1。降低收斂標準(算一般的有機分子不要緊,如果是過渡金屬,那就不要了)2。增大循環(huán)次數(shù) 3。換基組和算法 4。改變初始猜測
是有機分子阿,呵呵,用的是AM1,還是太大嗎hsi
這種不收斂的情況可能是你最初用軟件建模的時候構型不合理,仔細考慮一下你的分子的構型是否合理,是否還存在其它的構型,單純地降低收斂標準意義不大,不過你可以先嘗試一下增大循環(huán)次數(shù),最多scf(maxcyc=300),再不收斂的話就不是循環(huán)的問題,是構型不對。
S能不能說說,π軌道和σ軌道的組成上有什么差別?這方面我比較欠缺。ure.You can directly analyze the components of the molecule orbitals from the population in the *.log file.我在版面上說過了,當初的討論還真是熱烈:)你可以翻看一下
我這里簡述幾句:
π軌道主要是由npx,npy,nd-1,nd+1組成σ軌道主要是由ns,npz(鍵軸方向),nd0
還有,nd+2,nd-2主要是組成δ軌道。
opt 是可以restart,前面的中間文件是少不了的。1, %chk=...#..opt Stopped 2, %chk=...#..opt(restart)guess=read geom=allcheck-----------------對弈freq的restart,是要采用數(shù)值計算方法才可以restart。(以下我省略chk和方法的輸入,自己加上)1, freq(numerical)2, freq(restart,numerical)guess=read
OPT=restart不能進行,why 因為死機OPT被中斷,我采用OPT=restart繼續(xù)進行,但是出現(xiàn)以下情況:
Restoring state from the checkpoint file “Mn7O14V11O2-PCOPT1-2.chk”.FileIO operation on non-existent file.FileIO: IOper= 2 IFilNo(1)=-997 Len= 120004 IPos= 0 Q= 8710320
dumping /fiocom/, unit = 1 NFiles = 16 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 667136 FType=2 FMxFil=10000 Number 0 0 0 0 0 508 522 598 Base 144896 149504 142848 22016 444416 143360 143872 172032 End 146944 172032 143360 142336 667136 143375 144572 172034 End1 146944 172032 143360 142336 667136 143872 144896 172544 Wr Pntr 143360 149504 142848 22016 199168 143360 143872 172032 Rd Pntr 143360 149504 142948 142020 199168 143360 143872 172032 Number 634 991 992 994 995 996 998 999 Base 199168 147456 146944 20480 142336 21504 20992 172544 End 444168 149097 146949 20510 142346 21554 21042 198796 End1 444416 149504 147456 20992 142848 22016 21504 199168 Wr Pntr 199168 147456 146944 20480 142336 21504 20992 172544 Rd Pntr 199168 149097 146949 20510 142346 21554 21042 172545
dumping /fiocom/, unit = 2 NFiles = 5 SizExt = 0 WInBlk = 512 defal = F LstWrd = 267264 FType=2 FMxFil=10000 Number 0 508 522 634 998 Base 266245 20480 20545 21245 20495 End 267264 20495 21245 266245 20545 End1 267264 20495 21245 266245 20545 Wr Pntr 266245 20480 20545 21245 20495 Rd Pntr 266245 20480 20545 21245 20495
dumping /fiocom/, unit = 3 NFiles = 1 SizExt = 524288 WInBlk = 512 defal = T LstWrd = 66048 FType=2 FMxFil=10000 Number 0 Base 20480 End 66048 End1 66048 Wr Pntr 20480 Rd Pntr 20480 Error termination in NtrErr: NtrErr Called from FileIO.為什么我的file IO會不存在?是什么參數(shù)設漏了嗎?請各位大蝦支持,謝謝!
在你前一次的計算中,第一單點的計算未完成,所以,checkpoint文檔中尚無優(yōu)化過程的信息,無法restart.過渡態(tài)優(yōu)化過程中遇到一個這樣的問題
我在過渡態(tài)優(yōu)化過程中用這樣的方法 # hf/6-31g(d)opt=(ts,calcfc)freq當運行到L999時Link died請問可能是什么原因?是不是我的過渡態(tài)輸入有問題?另外,尋找過渡態(tài)有哪些技巧,請指教!
你的計算在設定的步數(shù)里面沒有找到過渡態(tài)。這有兩個可能:一是初始輸入的構型不好。二是循環(huán)次數(shù)不夠(可以maxcycle=大一點的數(shù))。具體還是等高手來指點吧。
對于決定收斂結構的標準,做了一些很小但很重要的變動。當力比截斷值小兩個數(shù)量級時(即,極限值的1/100),結構即被認為是收斂的,即便是位移比截斷值還大。這個變動有利于大分子的優(yōu)化,因為它們的最小值周圍可能有非常平緩的勢能曲面。(G98manual p106)所以你的opt有2個yes 頻率計算完成后還有一次優(yōu)化,可能這次就優(yōu)化到3個了
!做頻率分析的時候只需要用Freq就可以,不需要加Raman選項,對嗎? freq=noraman不作拉曼強度計算,可以提高10%-20%的計算速度.算一個分子的電離能,里面有I,CL 在優(yōu)化起正離子的時候老是沒有沒有結果 體系的能量趨勢不固定,有時大有時小,MAXCYCLE=80后出現(xiàn)2個YES,然后DIE LINK 999 但是COPY。OUT最后的構型再去優(yōu)化后 結構4個全是NO 現(xiàn)在還在算,但是還是NO 怎么辦? 要OPT=CALCALL嗎?有什么用? 另外,還有一個題外的問題 這樣寫對嗎?
OPT(MAXCYCLE=80,STEPSIZE=10)可以嗎?
增大優(yōu)化次數(shù),optcyc=N;N給大一些。opt=(maxcycle=N)也一樣
前面的老帖子已經(jīng)講到這些,好象在免費資源里的Gaussian98 tips.一部分引用如下:
First, whenever you encounter SCF convergence failures you need to look at the progress of the SCF by using #P in place of # which will print out the status at each iteration.It is useful to know if the SCF is converging slowly, oscillating, taking large jumps etc.Slow convergence can be a sign that a higher order converger like SCF=QC would be helpful.Oscillations can be a sign that a small HOMO/LUMO gap is causing a flip in the occupied orbitals and SCF=VShift=n is useful.Large jumps can indicate that a poor initial guess or a change in the electronic state from the initial guess has occured and normal convergence will eventually achieved but after a larger number of cycles and a better initial guess would help.Second, whenever possible it is helpful to start with a smaller basis set, i.e.STO-3G, LANL1MB, etc for single points.This is cheaper and so you can often generate a better initial guess which can be used with GUESS=READ and the larger basis set and without the more expensive options.In your case the SCF converges initially quite smoothly and then bounces around without making progress.This makes SCF=QC a likely approach to clean up this problem, at least at the initial point.You might want to look at adding IOP(5/13=1)which allows the SCF to print out the MO's even in the event that convergence was not reached.This should only be used with single points, i.e.remove OPT.As to optimization constraints.To freeze atomic coordinates you should turn on NoSymm.This method of freezing structures causes problems for redundant internals which are being looked into but NoSymm fixes the most serious ones.下面是血紅素的化學合成的模型,由于分子有點大,又加上有個Fe,所以優(yōu)化非常困難,請求斑竹支 高招:
(1)這是一個高度對稱的分子,如何利用加虛原子,提高計算速度?
(2)優(yōu)化Fe,我采用了LANL2DZ,請問有沒有更好的基組?怎么手工加基組?(3)由于變量太多,您認為我通過怎樣的辦法,減少變量來提高計算速度? 非常感謝??!做出文章了,我一定掛你的大名,呵呵。
S O 1 B1 Fe 2 B2 1 A1 N 3 B3 2 A2 1 D1 N 3 B4 2 A3 1 D2 N 3 B5 2 A4 1 D3 N 3 B6 2 A5 1 D4 C 4 B7 3 A6 2 D5 C 8 B8 4 A7 3 D6 C 9 B9 8 A8 4 D7 C 4 B10 3 A9 2 D8 C 5 B11 3 A10 2 D9 C 12 B12 5 A11 3 D10 C 13 B13 12 A12 5 D11 C 5 B14 3 A13 2 D12 C 6 B15 3 A14 2 D13 C 16 B16 6 A15 3 D14 C 17 B17 16 A16 6 D15 C 6 B18 3 A17 2 D16 C 7 B19 3 A18 2 D17 C 20 B20 7 A19 3 D18 C 21 B21 20 A20 7 D19 C 7 B22 3 A21 2 D20 C 23 B23 7 A22 3 D21 H 24 B24 23 A23 7 D22 C 11 B25 4 A24 3 D23 H 26 B26 11 A25 4 D24 C 16 B27 6 A26 3 D25 H 28 B28 16 A27 6 D26 C 20 B29 7 A28 3 D27 H 30 B30 20 A29 7 D28 H 1 B31 3 A30 2 D29 H 9 B32 8 A31 4 D30 H 10 B33 9 A32 8 D31 H 13 B34 12 A33 5 D32 H 14 B35 13 A34 12 D33 H 17 B36 16 A35 6 D34 H 18 B37 17 A36 16 D35 H 21 B38 20 A37 7 D36 H 22 B39 21 A38 20 D37 B1 4.21561119 B2 1.64053835 B3 2.01170682 B4 2.00223534 B5 2.01140630 B6 2.01948030 B7 1.38502525 B8 1.45027869 B9 1.36206816 B10 1.38792662 B11 1.38546948 B12 1.45027336 B13 1.36191526 B14 1.38377173 B15 1.38899517 B16 1.45046091 B17 1.36116892 B18 1.38527197 B19 1.39574494 B20 1.44813170 B21 1.36193008 B22 1.39696571 B23 1.37779525 B24 1.08374598 B25 1.38063627 B26 1.08377153 B27 1.37911143 B28 1.08372124 B29 1.37753366 B30 1.08361344 B31 1.37864549 B32 1.07802759 B33 1.07801256 B34 1.07797961 B35 1.07795952 B36 1.07798544 B37 1.07794319 B38 1.07803016 B39 1.07806743 A1 6.31961440 A2 94.07195309 A3 95.30512753 A4 94.38669042 A5 94.58713162 A6 127.10303272 A7 109.89556522 A8 107.10994453 A9 126.54307980 A10 126.81263629 A11 109.85724161 A12 107.06302128 A13 126.91593763 A14 126.68048694 A15 109.77898998 A16 107.13028054 A17 126.98352930 A18 126.67933314 A19 109.65642922 A20 107.43086078 A21 126.85650291 A22 125.36309437 A23 117.60829677 A24 125.56378588 A25 117.60334742 A26 125.52442366 A27 117.63863540 A28 125.46104554 A29 117.64536447 A30 98.48331793 A31 124.49018481 A32 128.37464062 A33 124.51087771 A34 128.41640662 A35 124.47708155 A36 128.39197501 A37 124.32863598 A38 128.24305473 D1-108.51711460 D2 161.17665430 D3 70.83221148 D4-18.82295908 D5 85.79136868 D6-174.21655652 D7 0.09638172 D8-87.07481057 D9 86.40576040 D10-175.60372334 D11-0.05103700 D12-88.28171080 D13 88.94769251 D14-174.64817315 D15 0.35121394 D16-84.34232126 D17 84.15970765 D18-172.91552513 D19 0.72935930 D20-85.63442533 D21-5.66240036 D22 179.49688145 D23-6.19734957 D24-179.54995066 D25 4.93855843 D26 179.77279426 D27 6.47288831 D28-179.56429342 D29-89.13791474 D30-179.79660459 D31-179.88612542 D32-179.83142492 D33-179.77466949 D34-179.60732447 D35-179.84766166 D36-179.39904616 D37 179.72163849
我也在做大分子的計算,希望能交流交流。我的電子郵件:chsun@imr.ac.cn 我的經(jīng)驗是先做PM3優(yōu)化,然后再用精度高的方法和基組。
你的前三個原子,不是一條直線嗎?如果是直線,那對稱性就可以是C4v了。如果不是直線,那至少要做成cs對稱性吧?把前三個分子所在的平面作為對稱面。你現(xiàn)在的對稱性是沒有——雖然看上去很對稱。Fe用贗勢基組不知道會怎么樣,你可以看看別人的文章。
做優(yōu)化的時候,可以固定一些外層的原子,只讓和中間的原子相連的動。
大約有3M 另外,F(xiàn)e、O、S分子是在一條直線上么?如果不是,你可以將吡啶環(huán)上的參數(shù)先減少了,然后再加Fe、O、S,肯定能夠減少不少參數(shù),但是不知道由沒有實際意義了
1)我談一下我做過的體系加虛原子的經(jīng)驗,希望對你有幫助:一般情況下,虛原子加在分子的對稱中心,或鍵軸Cn軸上,如要加第二個,可放在垂直于Cn軸的方面上。(2)你可以到基組庫去尋找基組,可能要比單單的Lanl2Dz效果好
(3)你的體系高度對稱,肯定有不少的CH鍵長鍵角是一樣的(或者差別很?。宦稍O為一個參數(shù),讓它同步優(yōu)化。
PS:你的輸入文件的參數(shù)也太多了點
我看了一下你的結構,你說是高度對稱?我想問你一下這是你手建的嗎?從你的坐標看好象沒有對稱性。32號H怎么位置這么怪!我認為你的結果最高對稱性是 C(4V),最好手建一下為好(如不考慮外場作用)。
這種體系的計算現(xiàn)在非常熱門也非常有難度,上次我們聽的以色列教授介紹的也是這方面的內容,這種體系以我們的水平可能很難處理好。你可能要找到做這方面的專家指導才會有很好的效果,否則很難分析清楚計算出的結果。
上次以色列教授沙松過來我們這邊就是介紹的有關P450的計算,雖然我沒有聽懂什么,不過以我自己計算含金屬的小分子經(jīng)驗來看,關于這種體系很難處理的部分就在于如果控制中心金屬原子鐵和周圍配體的配位情況,以及了解中心金屬原子的電子占據(jù),如果你僅是靠軟件產(chǎn)生體系的輸入坐標,可能計算很難達到你想要的結果。這些坐標的產(chǎn)生可以要靠自己的經(jīng)驗來步步設置,所以就像你所說的太細節(jié)的東西是很難說清楚的,因為他們在計算過程中擁有比你更多的經(jīng)驗。只能給你這些參考性建議了,沒有什么實質性的幫助。
我的輸入文件是:
%chk=zux3
#p hf/lanl2dz opt symm=loose scf=(maxcycle=200)test
zux3
0 1
0.000000 0.000000 1.463625 6-2.059783 0.151849 1.234523 8-3.189332 0.235187 0.995024 6 0.153438 2.073727 1.558845 8 0.238770 3.225066 1.644723 6 0.000000 0.000000 3.544137 8 0.000000 0.000000 4.703699 6 0.000000 0.000000-0.664692 6 0.000000 0.000000-1.891052 6 0.000000 0.000000-3.251669 7-1.186271-0.125670-3.989741 6-2.427060-0.359632-3.218024 6-1.381218 0.710276-5.201852 1-2.263789-1.158724-2.497766 1-2.749787 0.537294-2.683916 1-3.209142-0.660105-3.913517 1-0.432645 1.124763-5.529600 1-1.810550 0.110340-6.005013 1-2.062473 1.534158-4.980011 6 2.059783-0.151849 1.234523 8 3.189332-0.235187 0.995024 6-0.153438-2.073727 1.558845 8-0.238770-3.225066 1.644723 7 1.186271 0.125670-3.989741 6 2.427060 0.359632-3.218024 1 2.263789 1.158724-2.497766 1 2.749787-0.537294-2.683916 1 3.209142 0.660105-3.913517 6 1.381218-0.710276-5.201852 1 0.432645-1.124763-5.529600 1 1.810550-0.110340-6.005013 1 2.062473-1.534158-4.980011
分子的點群是c2,但是優(yōu)化第二輪是就出錯了,Omega: Change in point group or standard orientation.Error termination via Lnk1e in d:g98wl202.exe.如果我不用關鍵詞symm=loose,優(yōu)化時分子點群就變成了c1。
請高手指教,這是何原因呢?若我想保持分子仍是c2,改如何處理?
3x!換內坐標吧
opt=Z-matrix 你用內坐標可能會好一點吧?而且也能保證你的對稱性的。[em12][em12] Angle between quadratic step and forces= 29.72 degrees.Linear search not attempted--first point.Iteration 1 RMS(Cart)= 0.00036483 RMS(Int)= 0.00000015 Iteration 2 RMS(Cart)= 0.00000015 RMS(Int)= 0.00000006 Variable Old X-DE/DX Delta X Delta X Delta X New X(Linear)(Quad)(Total)R1 2.05001-0.00006 0.00000-0.00017-0.00017 2.04984 R2 2.83834 0.00000 0.00000 0.00009 0.00009 2.83844 R3 2.05240-0.00007 0.00000-0.00023-0.00023 2.05218 R4 2.05240-0.00007 0.00000-0.00023-0.00023 2.05218 R5 2.03466-0.00003 0.00000-0.00007-0.00007 2.03459 R6 2.47511 0.00010 0.00000 0.00015 0.00015 2.47526 R7 2.02743-0.00009 0.00000-0.00024-0.00024 2.02719 R8 2.51854 0.00015 0.00000 0.00038 0.00038 2.51892 A1 1.94852-0.00009 0.00000-0.00069-0.00069 1.94783 A2 1.88382 0.00007 0.00000 0.00042 0.00042 1.88424 A3 1.88382 0.00007 0.00000 0.00042 0.00042 1.88424 A4 1.93636-0.00006 0.00000-0.00034-0.00034 1.93602 A5 1.93636-0.00006 0.00000-0.00034-0.00034 1.93602 A6 1.87204 0.00008 0.00000 0.00063 0.00063 1.87267 A7 2.06277 0.00007 0.00000 0.00058 0.00058 2.06335 A8 2.15661 0.00001 0.00000 0.00012 0.00012 2.15674 A9 2.06380-0.00008 0.00000-0.00071-0.00071 2.06310 A10 2.19767 0.00004 0.00000 0.00041 0.00041 2.19808 A11 2.13077-0.00008 0.00000-0.00050-0.00050 2.13027 A12 1.95474 0.00004 0.00000 0.00009 0.00009 1.95483 D1 3.14159 0.00000 0.00000 0.00000 0.00000 3.14159 D2 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 D3 1.03955 0.00001 0.00000 0.00017 0.00017 1.03972 D4-2.10204 0.00001 0.00000 0.00017 0.00017-2.10187 D5-1.03955-0.00001 0.00000-0.00017-0.00017-1.03972 D6 2.10204-0.00001 0.00000-0.00017-0.00017 2.10187 D7 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 D8 3.14159 0.00000 0.00000 0.00000 0.00000 3.14159 D9 3.14159 0.00000 0.00000 0.00000 0.00000 3.14159 D10 0.00000 0.00000 0.00000 0.00000 0.00000 0.00000 Item value Threshold Converged? Maximum Force 0.000148 0.000450 YES RMS Force 0.000059 0.000300 YES Maximum Displacement 0.000912 0.001800 YES RMS Displacement 0.000365 0.001200 YES Optimization completed.--Stationary point found.詳細解釋一下好嗎,這一段幾乎全看不懂
Iteration 1 RMS(Cart)= 0.00036483 RMS(Int)= 0.00000015 Iteration 2 RMS(Cart)= 0.00000015 RMS(Int)= 0.00000006 這應該是按笛卡爾坐標和內坐標計算的均方誤差,比較的標準不清楚!
Variable列是輸入的變量;Old X列是上一次的變量值;-DE/DX是能量對坐標的一階導數(shù),是判斷駐點的依據(jù);Delta X是與上次相比變量的變化;New X是新的變量數(shù)值!
SCF Done: E(RHF)=-113.866331170 A.U.after 12 cycles Convg = 0.1529D-08-V/T = 2.0020 S**2 = 0.0000 NROrb= 34 NOA= 8 NOB= 8 NVA= 26 NVB= 26 最后一行是什么東東?
總軌道數(shù)34,Alpha占居軌道8個,Beta占居軌道8個,Alpha虛軌道26個,Beta虛軌道26個!NROrb=number of orb NOA= numberof alpha occupied orb NVA=number of virtual alpha orb
結構優(yōu)化的時候有沒有可能出現(xiàn)這種情況:
即構型在兩個或幾個結構循環(huán)變化,總也不能達到終點 可能出現(xiàn)么?如果這樣的話該如何處理?非線形搜索? 多謝!
“構型在兩個或幾個結構循環(huán)變化,總也不能達到終點,可能出現(xiàn)么?”這個問題是你實際遇到的,還是憑空想出來的?幾個結構在勢能面上對應著局部最小,自洽場的結果是會到達一個局部最小值,而不會出現(xiàn)在幾個個結構循環(huán)變化。我覺得你想說的是收斂振蕩的問題吧?
他指的可能是結構優(yōu)化振蕩,我曾經(jīng)用過iop(1/8)里的一些選項,對控制結構優(yōu)化有一定幫助
我現(xiàn)在遇到的問題是這樣的,優(yōu)化了很多圈,但每次的能量總是在很小的范圍內起伏
只能在Force上達到兩個YES,而Displacement上總是NO,我覺得就象是在幾個結構中振蕩一樣,這是所謂的優(yōu)化振蕩么?
iop(1/8)里有設定優(yōu)化步長的功能:
0 DXMAXT = 0.1 BOHR OR RADIAN(L103, Estm or UnitFC).= 0.3 Bohr or Radian(L103, Read or CalcFC).= 0.2 Bohr or Radian(L105).= 0.3 Bohr or Radian(L113, L114).N DXMAXT = 0.01 * N 哪里有關于iop(1/8)的詳細內容呢??
在免費資源版有IOP的pdf文件
謝supi!使用了iop(1/8=1)結構終于優(yōu)化完了,但這說明什么呢?體系的能量對于結構的微小變化很敏感?還是說勢能面比較平坦?
版上的今天剛貼的Guassian tips這篇貼子中最后有提到收斂遇到的一種情況就是“taking large jumps”.對于優(yōu)化默認采用的步長為0.3Bohr對你優(yōu)化的結構來說太大,所以你看到結果能量總是在很小的范圍內起伏.因此你使用IOP(1/8=1)設置步長為0.01Bohr,而達到收斂的最終結果,說明這個結構的勢能面是起伏不平的,對結構的變化非常敏感.不過,這種情況也像是oscillating,建議用SCF(Vshift=300)試一下,看能不能收斂到同一結果?
試了好幾次,總覺得使用二次收斂方法scf=qc,對收斂沒有什么太大的效果,反而不如以前.而且對于oscillating采用SCF(Vshift=N)來說,由于DFT已經(jīng)使用N=200,即使增大N的值對收斂的結果也無明顯效果,比不上減小優(yōu)化步長.但是對于減少步長來說,機時消耗太大,雖然在能量上比前者來得低,但是仍可能達不到收斂的預期效果,因此這些都說明了一個問題:The initial guess isn't very good.關于如果處理這些問題,《計算化學》上給出了一些更好的建議。
看了關于學術腐敗的討論,很慚愧,學習高斯已經(jīng)有好幾個月,還停留在Opt階段,而且還有不少問題,先問兩個吧,請莫見笑。
1.我發(fā)現(xiàn)構型優(yōu)化時,Opt(tight)和只有Opt時,機時差別若干倍,但結果好像差的不是很多,為什么? 2.對于大分子,可以先用小基組優(yōu)化,而后再用更高級基組優(yōu)化,具體寫輸入文件時,小基組時的分子說明部分用初始分子構型的數(shù)據(jù),而更高基組優(yōu)化時的分子說明部分應該用小基組的計算結果,怎么寫呢?
。opt=tight采用了更為嚴格的收斂標準,使用選項Opt=Tight的優(yōu)化計算比使用默認截止點的優(yōu)化計算要多算幾步。所用的時間當然也大得多了。要注意我們使用的能量的原子單位是個很大的數(shù)值1au=627.51kcal/mol,也許用兩種方法算出來的能量只差了0。05個單位,可換算一下就有30多的kcal/mol!2。你可以先用小基組算,然后用前面的chk文件來作大基組的計算。在控制行加入geom=check,guess=read,就不用分子構型輸入了
結構優(yōu)化不收斂,一般有什么OPT=ReadFc, opt=CalcFc, opt=calcAll,或者設置maxcycle,opt=restart等等,但有沒有人知道那種最為有效?具體方法怎樣?當然有示例是最好的了!呵呵:)談談我的想法: 我感覺結構不收斂可能與你的初始猜測有關,也許先固定某些鍵長鍵角再優(yōu)化,事后再放開優(yōu)化會比較容易收斂.采用低級算法先計算只是其中的方法之一,你為何一直要用這種方法?有什么依據(jù)? 如果你順利解決了你的計算,那就談談你的經(jīng)驗如何?
我覺得采用OPT=CALCFC并不一定得到好的結果,有時候收不收斂實際和我們處理的化學問題相關,而不僅僅是用計算上的技巧進行解決.比如兩個原子間電子占據(jù)沒有sigma成鍵,但是存在強的sigma非鍵作用,那么它們之間的斥力很強,計算就很難收斂形成穩(wěn)定的體系.因此在處理收斂的問題時,我們應該多分析不收斂的問題實質,是電子的什么運動狀態(tài)引起的,而非看到不收斂就想利用一些key來解決,我自己也試過一些key,包括opt=calcfc,allcfc,NoDIIS,QC等等(包括一些強制性的手段),并不是很有效,因為它們都是從技術的角度來解決問題.如果一個體系能形成穩(wěn)定的體系,那么它不依靠這些key也是能解決的.反而如果一體系不收斂只是利用了這些技術就收斂了,那么這種結果如何可信?前面的貼子中,有些技巧諸如調整原子間的距離和鍵角,這些手段才是比較有效的,因為它直接影響了原子間的電子運動狀態(tài),即波函數(shù),我們有時候初始猜測不正確就是因為程序給出的電子的占據(jù)狀態(tài)錯誤而直接影響了原子間的成鍵,這時就需要對電子占據(jù)的軌道進行調整,Guess=alter,這對于過渡金屬的計算來說特別要注意,因為G98總不能給出正確的初始猜測.謝謝諸位意見,收益非淺。我起初一直認為我的初始模型是正確的,但后來我試過4個模型,才發(fā)現(xiàn)初始結構非常重要。開始三個是MnO2理想的晶體結構,通過改變鍵長希望結構優(yōu)化收斂,但是失敗。我又重新定義采用MnO2的真實晶胞參數(shù),現(xiàn)在結果有些眉目,能收斂只是沒有正常結束,可能是我吸附分子的坐標出了點問題,出現(xiàn)
Error on Z-matrix card number 23 angle Alpha is outside the valid range of 0 to 180.所以我深切體會到樓上師兄的話:調整原子間的距離和鍵角,這些手段才是比較有效的。我現(xiàn)在才發(fā)覺,G98真是一個費內存,費時間的東西,很多計算要的是經(jīng)驗,你計算的垃圾多了,你才知道你需要什么!這個過程真夠悲慘。
第四篇:總結數(shù)位DP算法
數(shù)位dp是一種計數(shù)用的dp,一般就是要統(tǒng)計一個區(qū)間[le,ri]內滿足一些條件數(shù)的個數(shù)。比如,[1,10000] 中統(tǒng)計不含有4的數(shù)。
所謂數(shù)位dp,字面意思就是在數(shù)位上進行dp咯。就是對數(shù)字每一位每一位遞推
此類題目最基本的暴力方法:
1.for(int i=le;i<=ri;i++)
2.if(Check(i))ans++;
而數(shù)位DP就是從最低(高)位起,一位一位的放數(shù)字,然后記憶化一下,累加一下
有兩種方法,一是遞推,二是記憶化搜索
一,記憶化搜索:
思路來自: 數(shù)位dp總結之從入門到模板 假設題目要求是不含有62的數(shù)
狀態(tài)定義:d[pos][pre] 表示當前枚舉到pos位置,且pos+1位的數(shù)字是pre,此時滿足題意的數(shù)字的個數(shù)(也即是pre==6時,pos該位置不能放2)還要個數(shù)組a[i]保存第i位的數(shù)字,如213,a[0]=3,注意是從右往左數(shù)
有個問題是枚舉第pos位數(shù)時,此位置放數(shù)字的范圍要判斷一下,比如題目給出在[1,894] 枚舉的時候要判斷是否在894以內
比如,213,第一位放了2,那么第二位就只能放0~1,所以模板中用了個limit判斷pos前的幾位數(shù)字是否與n一樣,true的話只能枚舉0~a[pos],false就是0~9,不然比題目要求的213大了
還有個問題是前導0的問題,假如枚舉5位數(shù),你放的時候前2位都是00,那數(shù)字不變成3位了嘛,所以需要個lead保存前幾位是否都是0,當然這是看題意的,有時候題目不要求,可以直接省去
好了,看模板:
1.typedef long long ll;2.int a[20];
3.ll dp[20][state];//不同題目狀態(tài)不同
4.ll dfs(int pos,/*state變量*/,bool lead/*前導零*/,bool limit/*數(shù)位上界變量*/)//不是每個題都要判斷前導零
5.{
6.//遞歸邊界,既然是按位枚舉,最低位是0,那么pos==-1說明這個數(shù)我枚舉完了
7.if(pos==-1)return 1;/*這里一般返回1,表示你枚舉的這個數(shù)是合法的,那么這里就需要你在枚舉時必須每一位都要滿足題目條件,也就是說當前枚舉到pos位,一定要保證前面已經(jīng)枚舉的數(shù)位是合法的。不過具體題目不同或者寫法不同的話不一定要返回1 */ 8.//第二個就是記憶化(在此前可能不同題目還能有一些剪枝)
9.if(!limit &&!lead && dp[pos][state]!=-1)return dp[pos][state];10./*常規(guī)寫法都是在沒有限制的條件記憶化,這里與下面記錄狀態(tài)是對應,具體為什么是有條件的記憶化后面會講*/
11.int up=limit?a[pos]:9;//根據(jù)limit判斷枚舉的上界up;這個的例子前面用213講過了
12.ll ans=0;13.//開始計數(shù)
14.for(int i=0;i<=up;i++)//枚舉,然后把不同情況的個數(shù)加到ans就可以了
15.{
16.if()...17.else if()...18.ans+=dfs(pos-1,/*狀態(tài)轉移*/,lead && i==0,limit && i==a[pos])//最后兩個變量傳參都是這樣寫的
19./*這里還算比較靈活,不過做幾個題就覺得這里也是套路了
20.大概就是說,我當前數(shù)位枚舉的數(shù)是i,然后根據(jù)題目的約束條件分類討論
21.去計算不同情況下的個數(shù),還有要根據(jù)state變量來保證i的合法性,比如題目
22.要求數(shù)位上不能有62連續(xù)出現(xiàn),那么就是state就是要保存前一位pre,然后分類,23.前一位如果是6那么這意味就不能是2,這里一定要保存枚舉的這個數(shù)是合法*/
24.}
25.//計算完,記錄狀態(tài)
26.if(!limit &&!lead)dp[pos][state]=ans;
27./*這里對應上面的記憶化,在一定條件下時記錄,保證一致性,當然如果約束條件不需要考慮lead,這里就是lead就完全不用考慮了*/
28.return ans;29.}
30.ll solve(ll x)31.{
32.int pos=0;
33.while(x)//把數(shù)位都分解出來
34.{
35.a[pos++]=x%10;//個人老是喜歡編號為[0,pos),看不慣的就按自己習慣來,反正注意數(shù)位邊界就行
36.x/=10;37.}
38.return dfs(pos-1/*從最高位開始枚舉*/,/*一系列狀態(tài) */,true,true);//剛開始最高位都是有限制并且有前導零的,顯然比最高位還要高的一位視為0嘛
39.}
40.int main()41.{
42.ll le,ri;
43.while(~scanf(“%lld%lld”,&le,&ri))44.{
45.//初始化dp數(shù)組為-1,這里還有更加優(yōu)美的優(yōu)化,后面講 46.printf(“%lldn”,solve(ri)-solve(le-1));47.} 48.}
注意:
那個if(!limit &&!lead &&dp[pos][state]!=-1)return dp[pos][state];limit 的數(shù)字必須要枚舉,不能直接返回,每次都要算
雖然這會導致重復,但這可以解決狀態(tài)沖突,而且重復計算的數(shù)字也很少 舉例如下:
題目:不能出現(xiàn)連續(xù)的11(11、112、211都是不合法的)那么我們開始枚舉:
要枚舉3位數(shù),已經(jīng)枚舉了兩位01_,要枚舉最后一位,此時狀態(tài)為d[0][1] 即:在枚舉個位,且前一位為1,那么顯然得出d[0][1]=9 開始新的一輪枚舉,枚舉到11_,此時狀態(tài)也是d[0][1] 因為已經(jīng)有9這個值了,所以返回了,但很明顯答案是0,是錯的 當然可以多開一維防止狀態(tài)沖突
可以看看數(shù)位DP模板題: HDU 2089 不要62 數(shù)位DP.二,遞推方法
思路來自:初探數(shù)位dp
狀態(tài)定義:d[i][j] 有i位數(shù)字,且第一位為j,在 0~j-1 + 000....999的符合題意的個數(shù),如 d[4][3] 就是在 3000~3999 的符合題意的個數(shù)
還要個數(shù)組a[i]保存第i位的數(shù)字,如213,a[1]=3,注意是從右往左數(shù)(下面是從1開始數(shù)起了)
這樣狀態(tài)定義的能更加方便,可以預處理,因為當一個數(shù)字的第一位比題目要求的第一位小后,后面的幾位能000..~999..如4269,如果第一位枚舉 3 _ _ _,那么后三位可以任取
模板如下:
1.for(int i=1;i<=7;i++)//枚舉位數(shù)
2.{
3.for(int j=0;j<10;j++)//枚舉第i位可能出現(xiàn)的數(shù)
4.{
5.for(int k=0;k<10;k++)//枚舉第i-1位可能出現(xiàn)的數(shù)
6.{
7.if(j!=4&&!(j==6&&k==2))//符合題意的條件
8.dp[i][j] += dp[i-1][k];9.} 10.} 11.}
以HDU 2089,解釋怎么算出答案(不含4,62的數(shù)字)
1.#include
2.#include
4.#include
5.using namespace std;6.int d[10][10],digit[10];
7.//d[i][j] 表示有i位數(shù)字,且第一位是j的數(shù)字的 滿足題意的數(shù)量
8.void init()9.{
10.d[0][0]=1;
11.for(int i=1;i<=7;i++)12.for(int j=0;j<=9;j++)13.for(int k=0;k<=9;k++)14.if(j!=4&&!(j==6&&k==2))15.d[i][j]+=d[i-1][k];16.}
17.int solve(int x)// [0,x)
18.{
19.int len=0;20.while(x){
21.digit[++len]=x%10;22.x/=10;23.}
24.digit[len+1]=0;25.int ans=0;
26.for(int i=len;i>=1;i--){
27.for(int j=0;j 28.if(j!=4&&!(j==2&&digit[i+1]==6))29.ans+=d[i][j];30.31.if(digit[i]==4||(digit[i+1]==6&&digit[i]==2))32.break;33.} 34.return ans;35.} 36.int main(int argc, char const *argv[])37.{ 38.int n,m;39.init(); 40.while(cin>>n>>m,n+m)41.cout< 42.return 0;43.} 假設一個數(shù)3229 得出 0000~0999 的個數(shù) 1000~1999 的個數(shù) 2000~2999 的個數(shù) 000~099 的個數(shù) 100~199 的個數(shù) 00~99 的個數(shù) 10~19 的個數(shù) 0~8 的個數(shù) 累加就是答案了 所以該區(qū)間是[0,n)是取不到的n的,注意計算的時候要加一個1 下面是一些題目: HDU 2089 不要62和4 HDU 3555 含49的數(shù) HDU 3652 含13且可以被13整除 codeforces 55d A 一個數(shù)字可以被它所有非零數(shù)整除的個數(shù) POJ 3252 Round Numbers HDU 4734 F(x)HDU 3709 Balanced Number HYSBZ 1799 self 同類分布 URAL 1057 Amount of Degrees * HDU 4507 吉哥系列故事——恨7不成妻 * 總結: 可能要用到的數(shù)位DP的題目類型: 1~10^18,求某區(qū)間(很大),有特定要求的數(shù)字的個數(shù) 如求mod,求和,可以整除各位數(shù),不出現(xiàn)某些數(shù)...框架: int DFS(intpos,......)//DFS一位一位放數(shù)字,求出答案,函數(shù)的參數(shù)保存題目要求的狀態(tài) int solve(int n)//把n一位一位拆分,求出[1,n] 的符合要求的值 難點:定義好狀態(tài)! 1.dp狀態(tài)要找好,不要出現(xiàn)狀態(tài)重疊現(xiàn)象,注意前導0有沒有影響 2.題目有求和sum,可能會很大,但可以轉化為保存sum對一個數(shù)求mod的值 3.有時候dp狀態(tài)定義不好可能要求每次DFS都要memset一下,換換思路想想通用的狀態(tài)定義,如sum從加法改為減法 《IT職業(yè)行為優(yōu)化》課程學習總結 大一開始學習IT職業(yè)行為優(yōu)化,在剛開始的學習中我就對這門課程的學習方式很感興趣。以小組為學習單位,靠小組內的同學互相協(xié)作來完成學習。在學習中多的不是作業(yè),多的而是小組成員間的討論以及對各項有用資料的收集。為的是完美的完成小組的任務。對此我的收獲頗豐,這不但鍛煉了我的團隊意識同時也讓我明白了一種好的職業(yè)行為可以決定你在那條職業(yè)道路上能走多遠。 在第一堂課的學習中我便感受到了這門課程的重要性,對于IT職業(yè)道德的理解以及職業(yè)行為的優(yōu)化有著很現(xiàn)實的看法。在大一機房維護組的招新中我有幸的進入機房工作,但是在機房的衛(wèi)生維護中感到很是頭疼。每天到機房要面對的就是同學們留下的零食垃圾和早飯,餐巾紙也是隨處可見,這無疑是加重了我們機房維護人員的工作量。在機房門口都有垃圾桶,可每次看見滿桌的雜務和空蕩蕩的垃圾桶,這種悲劇性對比讓我心里很不是滋味。這種行為特別需要優(yōu)化。因為人們生活在職業(yè)活動之中,不斷受到職業(yè)活動和職業(yè)道德的影響,又通過自己把它傳遞到另一個職業(yè)或者許許多多職業(yè)中去,從而引起多種多樣的連鎖反應。所以一種良好的職業(yè)道德是多么的重要。 從學習中我明白了如何養(yǎng)成良好的職業(yè)道德行為。第一要在日常生活中培養(yǎng),第二要在專業(yè)學習中訓練,第三在社會實踐中體驗,第四在自我修養(yǎng)中提高,第五在職業(yè)活動中強化。這五步是我在學習后要去經(jīng)歷的。 IT行業(yè)是一個日新月異的行業(yè),每時每刻到在產(chǎn)生新的元素.我們將來做這個行業(yè)就要有IT人員的基本素養(yǎng),要有規(guī)范化,標準化的代碼編寫習慣;需求理解能力;要具備復用性,模塊化的思維能力;要有經(jīng)常測試的習慣;要具備學習和總結的能力;同樣,高水平的交流能力也非常重要;最后最終要的是做一件事必須要有耐心,有毅力.當一個熱病能夠意識到自己的責任時,他又在完善自己的道路上邁進了一大步 即將步入社會,我們需要順應時代要求,把握歷史機遇,全面了解行業(yè)、了解企業(yè),主動去感受企業(yè)文化、適應社會規(guī)范,為將來成功就業(yè)鋪好路基。為此,我們需要加強職業(yè)素質培養(yǎng)與訓練,用自己的聰明才智和良好的職業(yè)素養(yǎng),去把握職業(yè)發(fā)展的機遇,為自己今后的職業(yè)生涯開拓寬廣而通暢的道路。 在學習這門課時,老師讓我們做了幾個實驗來讓我們體會將來的路.同時我們將實驗寫成報告來總結經(jīng)驗,讓我學習到了很多.通過“插手實驗”,讓我懂得了手與腦的分工合作,這樣不僅能提高工作效率還能加強工作質量,手與腦分工的明確,會讓自己更加的聰穎.通過實驗“感受榮譽”,讓我了解到的榮譽給予人的一種無與倫比的激勵,一個人獲得容易能激勵一個人的發(fā)展,然而,利弊應相當,同樣也會使人驕傲,所以有一個正確的榮譽觀對于我們的職業(yè)生涯非常重要.我們要正確的對待榮譽.注重道德社會價值的整合,創(chuàng)造形成大學生正確榮譽觀的社會氛圍。亞里士多德說過:“德性是一種習慣”,而習慣的養(yǎng)成是環(huán)境熏陶的結果。當前,社會價值取向的多元化和彼 此間的沖突是造成大學生榮譽感淡化的一個重要原因。無視甚至是抹殺這種多元性顯然是不現(xiàn)實的,這就要求我們必須積極營造有利于形成大學生正確榮譽觀的社會氛圍。即堅持一定的道德導向,通過法律和教育等手段使人們知道“應該”做什么,“不應該”做什么;怎樣做是“善”的,怎樣做是“惡”的,形成一個正確的善惡、榮辱觀。 這門課程讓我知道了什么叫做職業(yè)化,所謂職業(yè)化就是一種工作狀態(tài)的標準化、規(guī)范化、制度化,即在合適的時間、合適的地點,用合適的方式,說合適的話,做合適的事。 它讓我知道了從事一種職業(yè)就要有一個好的職業(yè)心理素質,職業(yè)心理素質是職業(yè)素質的一種,職業(yè)素質是指勞動者對社會職業(yè)了解與適應能力的一種綜合體現(xiàn),其主要表現(xiàn)在職業(yè)興趣、職業(yè)能力及職業(yè)情況等方面。第五篇:IT職業(yè)行為優(yōu)化小結