/Tutorials_CP2K

Tutorials on CP2K calculations

Primary LanguagePythonCreative Commons Zero v1.0 UniversalCC0-1.0

CP2K tutorials for electronic structure calculations

This repository contains input examples for calculations associated with CP2k package, including the NonAdiabatic Molecular Dynamics (NA-MD) calculations via CP2K/Libra interface.

1. General instructions

This repository summarizes inputs for different types of electronic structure calculations. The same files are also available in other places as well (project_cp2k_libra, project_perovskite_crystal_symmetry) for use in nonadiabatic dynamics. However, here we provide detailed instructions on how to use CP2K and how the functionality and timings change with different inputs.

Here we use CP2K v6.1 compiled with Intel parallel studio 2019. For TD-DFT with hybrid functionals, we will use CP2K v7.1 which is compiled with GCC-8.3 compiler. This is because in lower versions of CP2K the TD-DFT calculation does not converge for hybrid functional due to a problem with converging ADMM calculations. It is also worth noting that the computed results are done using Intel(R) Xeon(R) E5-2699 v4 @ 2.20GHz CPUs.

For CP2K installation we recommend the use of ./install_too_chain.sh in the tools/toolchain folder of CP2K. For compilation with Intel parallel studio you can use the instructions given in XConfigure website.

Here, we use a 2D perovskite of (BA)2PbI4 as a test and which its cif file is available from the crystallography website. We will use the unit cell with 156 atoms and a 2x2x2 supercell with 1248 atoms and check the performance of CP2K calculations functionality.

2. Types of calculations

  • geometry preparation
  • energy (single-point) calculations
  • convergence with respect to basis set parameters
  • geometry and cell optimization
  • single-point calculations at the TD-DFT level
  • single-point calculations with hybrid density functionals
  • single-point calculations of huge systems
  • molecular dynamics
  • time-overlap calculations using cube files
  • nonadiabatic coupling calculations

The legacy folder contains a number of recent tutorials by Mohammad Shakiba and Brendan Smith, but they are yet to be organized and revised. These tutorials contain detailed and very useful instructions which will eventually be migrated to the main folder here.

We highly welcome improving the functionality of the input files in this repository. So, please feel free to share your inputs and timings with us if you used these inputs.

3. CP2K Installation

The instructions on how the CP2K is compiled with Intel Parallel Studio can be found in this link.