CP2K Open Source Molecular Dynamics

User Tools

Site Tools

Trace:
exercises:common:chg

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revisionPrevious revision
Next revision		Previous revision
exercises:common:chg [2022/10/18 15:01] – jglanexercises:common:chg [2023/12/04 13:49] (current) – fnunes
Line 6: Line 6:
  
 It requires three single-point calculations of system AB, A, and B to print the [[https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/E_DENSITY_CUBE.html|electron density]] or [[https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/TOT_DENSITY_CUBE.html|total density]]. It requires three single-point calculations of system AB, A, and B to print the [[https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/E_DENSITY_CUBE.html|electron density]] or [[https://manual.cp2k.org/cp2k-2022_1-branch/CP2K_INPUT/FORCE_EVAL/DFT/PRINT/TOT_DENSITY_CUBE.html|total density]].
 +
 +One can use [[https://www.cp2k.org/tools:cubecruncher | Cubecruncher]] to manipulate the CUBE files.
 +
 +<code>
 +cubecruncher.x -i AB.cube   -subtract A.cube -o AB_A.cube
 +cubecruncher.x -i AB_A.cube -subtract B.cube -o chg_dif.cube
 +</code>
 +
 +This type of analysis is also useful in the case of time-dependent calculations. In this case, the difference would be taken over time, with respect to some reference snapshot (e.g. ground state). 
 +
 +$\Delta \rho (t-t_0) = \rho(t) - \rho(t_0)$,
 +
 +where $\rho(t)$ is the electron density at time $t$, and $\rho(t_0)$ is the reference electron density.
 +
 +<code>
 +cubecruncher.x -i elec_dens_t.cube -subtract elec_dens_t0.cube -o diff_t-t0.cube
 +</code>
exercises/common/chg.1666105303.txt.gz · Last modified: by jglan