20221018

•Use Spartan or WebMO to optimize the structure of the following molecule with B3LYP  theory with 6-31+G(d). Notice that different multiplicity might cause different structures, compare the energies of the structures too.

•Check and correct the data in our homework these weeks.

•Plot the eletron density of He and He- and use MP2 theory with d-aug-cc-pVTZ and aug-cc-pVTZ  basis sets.

• B3LYP/6-31+G(d)-square planar

   -2941.32293 hartree

• B3LYP/6-31+G(d)-tetrahedral

   -2941.33094 hartree

• 如何繪製電子密度圖

Step1:先利用 MP2/aug-cc-pvtz 計算He原子能量

Step2:將計算完後產生的chk 檔做成formchk檔

Step3:建立xyz.dat

需要自己去計算要掃的平面以及範圍，第一行:-1代表讀取此文字檔後面三點為起始座標，只掃xy平面

第二、三、四行:分別為xyz，第一個數字:number of point， (-)表單位bohr，(+)表單位Angstrom，每一向量掃描一次在yz方向取60個點

Step3:將formchk做成cube檔

(ref:https://gaussian.com/cubegen/)

Step4:以./ 執行rcubes檔，生成*.dat (含有座標以及電子密度的文字檔)

2:計算時所使用記憶體 density:計算電子密度指令 He--daug.fchk He--daug.cube:formchk檔做成cube檔 -1<xyz.dat:為(60)^2個點

Step5:將複製*dat內座標及電子密度至excel繪製出電子密度圖

• Analysis