|
CP2K 2.4 (Revision 12889)
|
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' |
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)
| 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)
This method would not work for a non periodic system It should be generalized like the calculation of Hartree
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:| 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.
1.7.3