<br>I have recently downloaded CP2K on my pc.
I want to perform a geometry optimization of an isolated Fe (CO) 5 molecule contained in a cubic cell with side 12 angstroms.
I created this input:<br><div><div>&GLOBAL</div><div> PROJECT FeCO5</div><div> RUN_TYPE GEO_OPT</div><div> PRINT_LEVEL LOW</div><div>&END GLOBAL</div><div>&FORCE_EVAL</div><div> METHOD QS</div><div> &SUBSYS</div><div> &CELL</div><div> ABC 12.0 12.0 12.0</div><div> &END CELL</div><div> &COORD</div><div> Fe 0.0000 0.0000 0.0000<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> C 0.0000 0.0000 1.8283<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> C 0.0000 1.8202 0.0000<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> C 1.5763 -0.9101 0.0000<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> C -1.5763 -0.9101 0.0000<span style="white-space:pre"> </span></div><div> C 0.0000 0.0000 -1.8283<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> O 0.0000 0.0000 2.9665<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> O 0.0000 2.9626 0.0000<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> O 2.5657 -1.4813 0.0000<span style="white-space:pre"> </span> <span style="white-space:pre"> </span></div><div> O -2.5657 -1.4813 0.0000<span style="white-space:pre"> </span> </div><div> O 0.0000 0.0000 -2.9665<span style="white-space:pre"> </span></div><div> &END COORD</div><div> &KIND Fe</div><div> BASIS_SET DZVP-GTH-PADE</div><div> POTENTIAL GTH-PADE-q8</div><div> &END KIND</div><div> &KIND C</div><div> BASIS_SET DZVP-GTH-PADE</div><div> POTENTIAL GTH-PADE-q4</div><div> &END KIND</div><div> &KIND O</div><div> BASIS_SET DZVP-GTH-PADE</div><div> POTENTIAL GTH-PADE-q6</div><div> &END KIND</div><div> &END SUBSYS</div><div> &DFT</div><div> BASIS_SET_FILE_NAME ./BASIS_SET</div><div> POTENTIAL_FILE_NAME ./POTENTIAL</div><div> &QS</div><div> EPS_DEFAULT 1.0E-7</div><div> &END QS</div><div> &MGRID</div><div> CUTOFF 400</div><div> NGRIDS 4</div><div> REL_CUTOFF 60</div><div> &END MGRID</div><div> &SCF</div><div> SCF_GUESS ATOMIC</div><div> EPS_SCF 1.0E-05</div><div> MAX_SCF 200</div><div> &DIAGONALIZATION T</div><div> ALGORITHM STANDARD</div><div> &END DIAGONALIZATION</div><div> &MIXING T</div><div> ALPHA 0.5</div><div> METHOD PULAY_MIXING</div><div> NPULAY 5</div><div> &END MIXING</div><div> &PRINT</div><div> &RESTART OFF</div><div> &END RESTART</div><div> &END PRINT</div><div> &END SCF</div><div> &XC</div><div> &XC_FUNCTIONAL PADE</div><div> &END XC_FUNCTIONAL</div><div> &END XC</div><div> &END DFT</div><div>&END FORCE_EVAL</div><div>&MOTION</div><div> &GEO_OPT</div><div> TYPE MINIMIZATION</div><div> MAX_DR 1.0E-03</div><div> MAX_FORCE 1.0E-03</div><div> RMS_DR 1.0E-03 </div><div> RMS_FORCE 1.0E-03</div><div> MAX_ITER 200</div><div> OPTIMIZER CG</div><div> &CG</div><div> MAX_STEEP_STEPS 0</div><div> RESTART_LIMIT 9.0E-01</div><div> &END CG</div><div> &END GEO_OPT</div><div> &CONSTRAINT</div><div> &FIXED_ATOMS</div><div> COMPONENTS_TO_FIX XYZ</div><div> LIST 1</div><div> &END FIXED_ATOMS</div><div> &END CONSTRAINT</div><div>&END MOTION </div></div><div> </div><div><br>I wanted to know if this input into CP2K logic makes sense.
I conclude by asking if I can run the calculation with the command "cp2k.popt -i FeCO5.inp -oFeCO5.out &". </div><div>Thank you.<br></div><div><br><br></div>