<div dir="ltr"><div>As an example. I wish to perform a non-periodic geometry relaxation of the [N(CH3)4]+ molecule. I have been unable to find any documentation for a simple example.</div><div><br></div><div>I have used the following input file to relax the standard N(CH3)3 molecule.</div><div><br></div><div>    &GLOBAL</div><div>      PROJECT CP2K</div><div>      RUN_TYPE GEO_OPT</div><div>      PRINT_LEVEL LOW</div><div>      EXTENDED_FFT_LENGTHS TRUE</div><div>    &END GLOBAL</div><div>    &FORCE_EVAL</div><div>      METHOD QUICKSTEP</div><div>      STRESS_TENSOR ANALYTICAL</div><div>      &SUBSYS</div><div>        &KIND C</div><div>          BASIS_SET DZVP-GTH-PBE</div><div>          POTENTIAL GTH-PBE-q4</div><div>        &END KIND</div><div>        &KIND N</div><div>          BASIS_SET DZVP-GTH-PBE</div><div>          POTENTIAL GTH-PBE-q5</div><div>        &END KIND</div><div>        &KIND H</div><div>          BASIS_SET DZVP-GTH-PBE</div><div>          POTENTIAL GTH-PBE-q1</div><div>        &END KIND</div><div>        &CELL</div><div>          ABC 10.0 10.0 10.0</div><div>          PERIODIC NONE</div><div>        &END CELL</div><div>        &TOPOLOGY</div><div>          &CENTER_COORDINATES TRUE</div><div>          &END CENTER_COORDINATES</div><div>          COORD_FILE_NAME NC3H9.xyz</div><div>          COORD_FILE_FORMAT XYZ</div><div>        &END TOPOLOGY</div><div>      &END SUBSYS</div><div>      &DFT</div><div>        BASIS_SET_FILE_NAME CP2K/Basis_Set/BASIS_SET</div><div>        POTENTIAL_FILE_NAME CP2K/Basis_Set/GTH_POTENTIALS</div><div>        &POISSON</div><div>          PERIODIC NONE</div><div>          POISSON_SOLVER WAVELET</div><div>        &END POISSON</div><div>        &QS</div><div>          METHOD GPW</div><div>          EPS_DEFAULT 1.0E-8</div><div>          EPS_PGF_ORB 1.0E-8</div><div>        &END QS</div><div>        &MGRID</div><div>          CUTOFF 350</div><div>        &END MGRID</div><div>        &XC</div><div>          &XC_FUNCTIONAL PBE</div><div>          &END XC_FUNCTIONAL</div><div>        &END XC</div><div>        &SCF</div><div>          SCF_GUESS ATOMIC</div><div>          EPS_SCF 1.0E-6</div><div>          MAX_SCF 60</div><div>          &OT</div><div>            MINIMIZER CG</div><div>            PRECONDITIONER FULL_SINGLE_INVERSE</div><div>          &END OT</div><div>          &OUTER_SCF</div><div>            EPS_SCF 1.0E-6</div><div>            MAX_SCF 10</div><div>          &END OUTER_SCF</div><div>        &END SCF</div><div>      &END DFT</div><div>      &PRINT</div><div>        &FORCES ON</div><div>        &END FORCES</div><div>      &END PRINT</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>    &END MOTION</div><div><br></div><div>What would be the best way to amend the input file and the optimal algorithm to use?</div><div>Thank you in advance.</div></div>