# Tag Info

Accepted

### What exactly is the "Optimizer" in CP2K geometry optimization?

The question says: I am having trouble understanding the difference between the optimizer options (CG, BFGS, LBFGS). These are three optimization algorithms that evidently are implemented in CP2K ...

### What exactly is the "Optimizer" in CP2K geometry optimization?

Optimizer refers to the algorithm used for finding the minimum of a function. In this specific case the optimizer is referring to "geometry optimization" which seeks to find the local minimum of the ...
• 8,017
Accepted

### Is it possible to use the Parrinello-Rahman barostat for NPT simulations in CP2K?

The barostat used by CP2K is the one developped by Glenn J. Martyna et al. J. Chem. Phys. 101, 4177 (1994). The Parrinello-Rahman barostat is not implemented in CP2K.
• 1,421
Accepted

### SCF energy keeps on fluctuating between two values, instead of converging, why?

The cause of the oscillations in the SCF energies is probably a "sloshing instability"; the two most common types are "charge sloshing" and "occupancy sloshing". Sloshing ...
• 8,052
Accepted

### Total force calculated with cp2k not zero

There are generally two main reasons why the total force can be nonzero (in any DFT code): Your forces are only correct when the calculation has converged fully to the ground state, and so in ...
• 8,052
Accepted

### Is basis set superposition error reduced when using the GAPW method?

Short answer: no. The idea of the GAPW method described in Theor. Chem. Acc. 103, 124 (1999) is simply to speed up the evaluation of the Coulomb and exchange-correlation contributions. Quoting from ...
• 19.7k
Accepted

### Is there any reason not to sum the kinetic and potential energy from an NPT simulation to get internal energy?

Energy in an NPT simulation is not conserved, but (once equilibrated), it will fluctuate around an average value, and that average value has meaning. That is the ensemble average for your NPT and is a ...
• 4,366

### Ab-initio molecular dynamics of liquid systems: best tools and resources

CP2K CP2K is an open source electronic structure and molecular dynamics software package to perform atomistic simulations of solid-state, liquid, molecular, and biological systems. It is especially ...
• 940
Accepted

### Calculate Hessian for configurations that are not at the minimum of the Potential Energy Surface

It is difficult to understand what you mean by "correct" answer, especially when we are discussing a "not-geometrically-optimized configuration" as you have called it, so I will ...
• 1,411

### How to run 2nd generation CPMD in CP2K?

A similar question has been asked on the cp2k mailing list. The relevant settings are the MOTION / MD / ...
• 2,669

### Calculate Hessian for configurations that are not at the minimum of the Potential Energy Surface

The Hessian at a non-stationary point is well-defined, since it is the second derivative of the energy wrt. nuclear coordinates. The transformation to vibrational normal coordinates, however, assumes ...
• 1,796

### Calculate Hessian for configurations that are not at the minimum of the Potential Energy Surface

You could in principle compute the entire Hessian for any configuration of a molecule. The problem is that we are typically not interested in the entire Hessian: we want the Hessian for the internal ...
• 15.6k
Accepted

### Geometric optimization taking a long time. Are my convergence criteria too tight in CP2K?

After a bit of hit and trial, I am able to run the pyridine sample, using CP2K v7.0 (git:d15a9c53f), however I think this can be run on any CP2K version beyond v2.5....
• 2,342

### Reading in E-field with CP2K

From looking at the very long wavelength you specified, I assume your intent is to apply a static field to the calculation? I read here that EFIELD is intended for ...
• 2,669
Accepted

### What are some methods for modeling bulk phase infrared spectra?

Vibrational maps and density of states In water AIMD community, to avoid excessively heavy calculations of molecular moments, researchers use velocity autocorrelation function (which converge much ...
• 2,342
Accepted

### cp2k conserved quantity changed when restart NVT MD?

In the input file you define &THERMOSTAT &CSVR TIMECON 9.9999999999999986E+01 &END CSVR &END THERMOSTAT So I guess you intended ...
• 1,421
Accepted

### Model coordination complex using GROMACS or CP2K

Here is how you can get the structure from PubChem, read the SMILES into Avogadro2, do a MMFF94 (classical) optimization, and then a single-point energy calculation using NWChem (RHF, 6-31G*). The ...

### Can CP2K be used for slab calculations?

CP2K employs Gaussian basis sets and can do 3D (crystals), 2D (slabs), 1D (wires) and 0D (atoms and molecules) periodic calculations naturally, without inducing fictitious interactions of periodic ...
• 19.7k

### What is a "charged system" in this specific context?

I have not used CP2K, but based on the documentation, this should refer to the charge of the whole system. So if you don't have a reason to believe your system is charged, you should treat it as ...
• 15.6k

### What are some methods for modeling bulk phase infrared spectra?

I'm not an expert in MD/AIMD software so I can really only discuss methods rather than an actual workflow/software that can carry it out. One approach that has been explored for accelerating AIMD and ...
• 15.6k

### Failure of CP2K density functional software to converge for a 64 Si atom amorphous structure

After modifying your file, I was able to get a converged energy of -250.37 hartrees, using CP2K v7.0 `(git:d15a9c53f), however I think this can be run on any CP2K version beyond v2.5. ...
• 2,342