CP2K 2.4 (Revision 12889)
Functions | Variables

qs_linres_nmr_epr_common_utils Namespace Reference

given the response wavefunctions obtained by the application of the (rxp), p, and ((dk-dl)xp) operators, here the current density vector (jx, jy, jz) is computed for the 3 directions of the magnetic field (Bx, By, Bz) More...

Functions

subroutine, public mult_G_ov_G2_grid (cell, pw_pool, rho_gspace, funcG_times_rho, idir, my_chi, error)
 Given the current density on the PW grid in reciprcal space (obtained by FFT), calculate the integral {r}[ ((r-r') x j(r))/|r-r'|^3 ] = Bind(r') which in reciprcal space reads (for G/=0) i G/|G|^2 x J(G)

Variables

CHARACTER(len=*), parameter,
private 
moduleN = 'qs_linres_nmr_epr_common_utils'

Detailed Description

given the response wavefunctions obtained by the application of the (rxp), p, and ((dk-dl)xp) operators, here the current density vector (jx, jy, jz) is computed for the 3 directions of the magnetic field (Bx, By, Bz)

History
created 02-2006 [MI]
Author:
MI

Function Documentation

subroutine,public qs_linres_nmr_epr_common_utils::mult_G_ov_G2_grid ( TYPE(cell_type),pointer  cell,
TYPE(pw_pool_type),pointer  pw_pool,
TYPE(pw_p_type),pointer  rho_gspace,
TYPE(pw_p_type)  funcG_times_rho,
INTEGER,intent(in)  idir,
REAL(dp),intent(in)  my_chi,
TYPE(cp_error_type),intent(inout)  error 
)

Given the current density on the PW grid in reciprcal space (obtained by FFT), calculate the integral {r}[ ((r-r') x j(r))/|r-r'|^3 ] = Bind(r') which in reciprcal space reads (for G/=0) i G/|G|^2 x J(G)

Note:
The G=0 component is not comnputed here, but can be evaluated through the susceptibility and added to the shift in a second time

This method would not work for a non periodic system It should be generalized like the calculation of Hartree

Author:
MI

Definition at line 58 of file qs_linres_nmr_epr_common_utils.f90.

References pw_types::COMPLEXDATA1D, mathconstants::gaussi, pw_pool_types::pw_pool_create_pw(), pw_methods::pw_transfer(), and timings::timeset().

Referenced by qs_linres_epr_ownutils::epr_ind_magnetic_field(), and qs_linres_nmr_shift::nmr_shift().

Here is the call graph for this function:

Here is the caller graph for this function:


Variable Documentation

CHARACTER(len=*),parameter,private qs_linres_nmr_epr_common_utils::moduleN = 'qs_linres_nmr_epr_common_utils'

Definition at line 40 of file qs_linres_nmr_epr_common_utils.f90.