03.01.2014 Views

CP/MM Tutorial

CP/MM Tutorial

CP/MM Tutorial

SHOW MORE
SHOW LESS

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

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!