Skip to main content

Questions tagged [cp2k]

For questions about (or related to) the software CP2K.

Filter by
Sorted by
Tagged with
1 vote
0 answers
41 views

How to Converge the CUTOFF and REL_CUTOFF

I'm doing the "How to Converge the CUTOFF and REL_CUTOFF" tutorial, but I can't get the results after cutoff_run.sh using .ssmp . The calculation cashes and end with the message below: ...
Inna's user avatar
  • 19
5 votes
0 answers
208 views

Help regarding Geometry Optimization of a bimetallic surface

I am trying to perform a geometry optimization (GEO_OPT) of a FCC 111 Pt$_{0.5}$Au$_{0.5}$ bimetallic surface (10x10x3) generated from ASE, with a lattice parameter of 3.9 $\overset{\circ}{\mathrm{A}}$...
ARKA PRAVA SARKAR's user avatar
4 votes
0 answers
29 views

Query regarding individual polarization calculation using CP2K

I am seeking assistance in calculating the individual polarization matrix from CP2K. While utilizing the Polar input under FORCE_EVAL -> PROPERTIES -> LINRES -> POLAR, I observed that the ...
Simantini Paul's user avatar
3 votes
0 answers
46 views

Set coordination number cutoff in CP2K's D3 dispersion

The D3 dispersion correction has 2 distance cutoff values. One for the calculation of coordination numbers, one for the energy contribution. In other implementations, such as vasp and DFTB+ (tblite), ...
Fuzzy's user avatar
  • 928
3 votes
0 answers
33 views

Normal coordinate analysis in TRAVIS

I have run a few AIMD calculations of some impurities trapped in a solid using CP2K and TRAVIS to simulate the resulting infrared spectrum. While the infrared spectrum was simulated well, there is no ...
isolated matrix's user avatar
2 votes
0 answers
86 views

Is memory per core more import than # of memory channels for high performance VASP calculation?

Cross posted on StackOverflow. I ask this question cause recently my department want to purchase a HPC cluster mainly for first principle calculation. For now we have three choices: (1) Per node dual ...
Tieyuan Bian's user avatar
4 votes
1 answer
196 views

Can CP2K be used for slab calculations?

I am trying to map out potential energy of a slab surface + cation interaction, since VASP has been giving me very not-so-reliable results I wanted to know if CP2K have the ability to do such ...
Xdrake's user avatar
  • 541
1 vote
0 answers
68 views

How to calculate constraint force for angle constraint in CP2K NVT MD simulation? [closed]

I am running a CP2K NVT molecular dynamics simulation and want to apply an angle constraint between three atoms. In my previous simulations, I used a bond constraint between two atoms, and it was easy ...
Ashique Lal's user avatar
  • 1,621
1 vote
0 answers
62 views

Getting CP2K to print the COLVAR info [closed]

Is it possible to make CP2K to print the values of all the defined COLVARs? I did find the following subsection in COLVAR section ...
Ashique Lal's user avatar
  • 1,621
3 votes
1 answer
116 views

Model coordination complex using GROMACS or CP2K

Is it possible to compute a 3D model of a calcium EDTA complex and it's energy? If so, how should I go about doing it or where can I obtain relevant information regarding computation methods? I ...
beardeadclown's user avatar
6 votes
1 answer
458 views

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

Sometimes I have observed this phenomenon where my calculations' SCF iterations keep on fluctuating between two values, and not converging. I have observed it in all sort of systems, ranging from ...
ipcamit's user avatar
  • 545
8 votes
1 answer
326 views

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

I have been studying liquid systems (particularly deep eutectic solvents) using classical molecular dynamics (MD) simulations. The main drawback of classical MD is the arbitrariness of the potential ...
user35952's user avatar
  • 343
2 votes
0 answers
184 views

CP2K geometry optimization SCF takes 3000s for a step [closed]

I am new to CP2K and doing a geometry optimization for a cluster. At the beginning everything looks good, each SCF step takes about 13 seconds, the energy drops quickly and the convergence looks like ...
Jessie A's user avatar
8 votes
1 answer
191 views

Total force calculated with cp2k not zero

I performed DFT calculations for a box of 64 water molecules under Periodic Boundary Condition using CP2K. I found that the forces on each atom do not sum up to zero, as shown in the attached figure. ...
andy90's user avatar
  • 305
7 votes
3 answers
510 views

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

I want to calculate the Hessian matrix for molecular configurations that are not at the minimum of the Potential Energy Surface, using cp2k. But it seems that cp2k requires the configuration to be ...
andy90's user avatar
  • 305
8 votes
1 answer
476 views

How to run 2nd generation CPMD in CP2K?

I believe CP2K by default implements Born–Oppenheimer MD. I have seen many recent papers that use 2nd generation Car–Parrinello MD using CP2K. But I couldn't find any direct way of switching to 2nd-...
Ashique Lal's user avatar
  • 1,621
10 votes
0 answers
464 views

Setup for 4-point flexible water model [closed]

I tried to search online if there are examples or suggestions on how to set up a 4-point flexible water model, such as TIP4P/2005f, qTIP4P/f, TIP4P/$\epsilon$ FLEX, etc., but there's no clear example ...
mykd's user avatar
  • 2,342
7 votes
0 answers
84 views

Molecular dynamics to reproduce dispersion interactions [closed]

I'm trying to simulate the behavior of chromophore molecules incorporated into a polymer under an external electric field. I've tried to use classical molecular dynamics (GROMACS) but it's unable to ...
Roman's user avatar
  • 2,373
9 votes
1 answer
337 views

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

I'm doing a polarization calculation for the first time. To specify one of the input parameters (the "reference point") correctly, I have to figure out if my system is "charged." ...
NTS's user avatar
  • 759
5 votes
0 answers
47 views

Kinetic energy cut-off for Gaussian-type orbitals when applying periodic boundary conditions [duplicate]

I have a solid with PBC (periodic boundary conditions) and I want to perform a DFT simulation. I use GTO as the basis. In a paper I read, I can transform the GTO to a Bloch function (crystalline ...
Suslik's user avatar
  • 271
10 votes
1 answer
103 views

Reading complex potential into CP2K

Is it possible to read a complex potential into CP2K? I have generated a potential using FDTD at a given frequency, which results in a complex result. The potential is to be read in to a real time ...
Jason M Gray's user avatar
7 votes
0 answers
104 views

CP2K Spherical Cell [closed]

I'm trying to run a gas phase cluster calculation in CP2K. I want to create an outside force to constrain water molecules that are in the simulation within a radius of a central ion. I want to create ...
dtg67's user avatar
  • 71
7 votes
1 answer
623 views

cp2k conserved quantity changed when restart NVT MD?

The details are descripted elsewhere, but got no answer for days, so I post it here for help. As title said, I noticed a relatively big change in conserved quantity and don't know why. I've upload ...
jtyang's user avatar
  • 295
14 votes
1 answer
316 views

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

CP2K implements the Gaussian and Augmented Planewaves (GAPW) approach for all-electron calculations. My understanding is that the GAPW method involves using atom-centered Gaussian type orbitals to ...
Hayden S's user avatar
  • 940
12 votes
0 answers
155 views

Periodic polarizable QM/MM embedding [closed]

Are there any standalone open source software codes available for periodic polarizable QM/MM embedding MD simulations? In my knowledge, CP2K only has the option for electrostatic embedding and other ...
mykd's user avatar
  • 2,342
8 votes
1 answer
635 views

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

The structure was created using Stillinger Weber with LAMMPS. The command I used was the following: ...
cp2kuser's user avatar
15 votes
0 answers
290 views

CP2K vs BigDFT comparison [closed]

I typically run DFT calculations with 1000 to 5000 atoms using CP2K. This works fine but I'm interested in BigDFT also. Is there anyone here that has experience using BigDFT and is able to compare ...
MD Simulation's user avatar
12 votes
1 answer
629 views

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

I have been having trouble with the VIBRATIONAL_ANALYSIS routine. Namely, the routine keeps timing out. I am trying to find the vibrational modes of Pyridine in a non periodic system, as well as $\ce{...
Jason M Gray's user avatar
8 votes
1 answer
298 views

Reading in E-field with CP2K

I have figured out a way to input a custom electric field in a .cube format, but I am not sure I am doing it in the proper way. Currently my file contains information for the magnitude of the electric ...
Jason M Gray's user avatar
7 votes
1 answer
200 views

What are the default units for the core hamiltonian and the density AO matrices in CP2K?

I have a few questions about default units. I am trying to find the default units for the density matrix and the core hamiltonian. I am assuming they are in e/(bohr)$^3$ and Hartrees respectively. I ...
Jason M Gray's user avatar
19 votes
1 answer
608 views

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

I would be very grateful for some newbie-level advice from a thermodynamics guru. I ran NPT simulations on a particular system (in CP2K software) to get fluid densities for use in fluid dynamics ...
NTS's user avatar
  • 759
10 votes
2 answers
295 views

What are some methods for modeling bulk phase infrared spectra?

I am interested in predicting the infrared spectra of bulk phase molecules, it seems that AIMD (ab initio molecular dynamics) is the current best approach. I have found a tutorial on using CP2K to do ...
Cavenfish's user avatar
  • 1,222
8 votes
1 answer
414 views

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

It is possible to use different types of thermostats in CP2K by specifying them in the corresponding section of the input as mentioned here. But there is no mention of specifying the barostat type in ...
Rashid Rafeek's user avatar
15 votes
2 answers
1k views

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

I'm just now starting to use the program CP2K, and I am having trouble understanding the difference between the optimizer options (CG, BFGS, LBFGS). I have seen ...
Cavenfish's user avatar
  • 1,222