CP/MM Tutorial
CP/MM Tutorial
CP/MM Tutorial
Create successful ePaper yourself
Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.
NUMBER OF PLANE WAVES FOR WAVEFUNCTION CUTOFF: 34686<br />
NUMBER OF PLANE WAVES FOR DENSITY CUTOFF: 277389<br />
****************************************************************<br />
Here we see how <strong>CP</strong>MD generates the initial guess for the wavefunction<br />
optimization. In this case it uses a superposition of atomic wavefunctions using an<br />
(internal) minimal Slater basis:<br />
*** RINFORCE| SIZE OF THE PROGRAM IS 16052/ 127264 kBYTES ***<br />
*** FFTPRP| SIZE OF THE PROGRAM IS 19764/ 128656 kBYTES ***<br />
GENERATE ATOMIC BASIS SET<br />
C SLATER ORBITALS<br />
2S ALPHA= 1.6083 OCCUPATION= 2.00<br />
2P ALPHA= 1.5679 OCCUPATION= 2.00<br />
O SLATER ORBITALS<br />
2S ALPHA= 2.2458 OCCUPATION= 2.00<br />
2P ALPHA= 2.2266 OCCUPATION= 4.00<br />
H SLATER ORBITALS<br />
1S ALPHA= 1.0000 OCCUPATION= 1.00<br />
INITIALIZATION TIME:<br />
0.39 SECONDS<br />
*** GMOPTS| SIZE OF THE PROGRAM IS 24656/ 141780 kBYTES ***<br />
*** PHFAC| SIZE OF THE PROGRAM IS 24764/ 152468 kBYTES ***<br />
*** ATOMWF| SIZE OF THE PROGRAM IS 25784/ 154300 kBYTES ***<br />
ATRHO| CHARGE(R-SPACE): 24.000000 (G-SPACE): 24.000000<br />
****************************************************************<br />
* ATOMIC COORDINATES *<br />
****************************************************************<br />
1 C 0.000000 0.000000 0.000000<br />
2 C 2.583394 0.000000 -0.913367<br />
3 C -1.291698 -2.237285 -0.913367<br />
4 O 0.000000 0.000000 2.305466<br />
5 H 2.583394 0.000000 -2.971279<br />
6 H 3.553503 1.680278 -0.227396<br />
7 H 3.553503 -1.680278 -0.227396<br />
8 H -1.291698 -2.237285 -2.971279<br />
9 H -0.321588 -3.917563 -0.227396<br />
10 H -3.231915 -2.237285 -0.227396<br />
****************************************************************<br />
From the following statement in the output you can see that the Hessian has been<br />
initialized from a simple guess assuming a molecule with some specified bonds. This<br />
behavior can be controlled with the keyword HESSIAN. For bulk systems or<br />
complicated molecules, it may be better to start from a unit Hessian instead:<br />
INITIALIZE EMPIRICAL HESSIAN<br />
><br />
2 1 3 1 4 1 5 2 6 2<br />
7 2 8 3 9 3 10 3<br />
><br />
5 1 5 3 6 1 6 4 7 1<br />
7 3 7 4 8 1 8 2 9 1<br />
9 2 9 4 10 1 10 4<br />
TOTAL NUMBER OF MOLECULAR STRUCTURES: 1<br />
****************************************************************<br />
* ATOMIC COORDINATES *<br />
****************************************************************<br />
1 C 0.000000 0.000000 0.000000<br />
2 C 2.583394 0.000000 -0.913367<br />
3 C -1.291698 -2.237285 -0.913367<br />
4 O 0.000000 0.000000 2.305466<br />
5 H 2.583394 0.000000 -2.971279<br />
6 H 3.553503 1.680278 -0.227396<br />
7 H 3.553503 -1.680278 -0.227396<br />
8 H -1.291698 -2.237285 -2.971279<br />
9 H -0.321588 -3.917563 -0.227396<br />
10 H -3.231915 -2.237285 -0.227396