(git:b6ef100)
Loading...
Searching...
No Matches
transport.F
Go to the documentation of this file.
1!--------------------------------------------------------------------------------------------------!
2! CP2K: A general program to perform molecular dynamics simulations !
3! Copyright 2000-2026 CP2K developers group <https://cp2k.org> !
4! !
5! SPDX-License-Identifier: GPL-2.0-or-later !
6!--------------------------------------------------------------------------------------------------!
7
8! **************************************************************************************************
9!> \brief routines for DFT+NEGF calculations (coupling with the quantum transport code OMEN)
10!> \par History
11!> 12.2012 created external_scf_method [Hossein Bani-Hashemian]
12!> 05.2013 created rotines to work with C-interoperable matrices [Hossein Bani-Hashemian]
13!> 07.2013 created transport_env routines [Hossein Bani-Hashemian]
14!> 11.2014 switch to CSR matrices [Hossein Bani-Hashemian]
15!> 12.2014 merged [Hossein Bani-Hashemian]
16!> 04.2016 added imaginary DM and current density cube [Hossein Bani-Hashemian]
17!> \author Mohammad Hossein Bani-Hashemian
18! **************************************************************************************************
19MODULE transport
20 USE iso_c_binding, ONLY: c_associated,&
21 c_bool,&
22 c_double,&
23 c_f_procpointer,&
24 c_int,&
25 c_loc,&
26 c_null_ptr,&
27 c_ptr
28 USE atomic_kind_types, ONLY: atomic_kind_type,&
29 get_atomic_kind
30 USE bibliography, ONLY: bruck2014,&
31 cite_reference
32 USE cp_control_types, ONLY: dft_control_type
33 USE cp_dbcsr_api, ONLY: &
34 dbcsr_convert_csr_to_dbcsr, dbcsr_convert_dbcsr_to_csr, dbcsr_copy, dbcsr_create, &
35 dbcsr_csr_create, dbcsr_csr_create_from_dbcsr, dbcsr_csr_dbcsr_blkrow_dist, &
36 dbcsr_csr_print_sparsity, dbcsr_csr_type, dbcsr_deallocate_matrix, dbcsr_desymmetrize, &
37 dbcsr_has_symmetry, dbcsr_set, dbcsr_type, dbcsr_type_no_symmetry
38 USE cp_dbcsr_cp2k_link, ONLY: cp_dbcsr_to_csr_screening
39 USE cp_log_handling, ONLY: cp_get_default_logger,&
40 cp_logger_get_default_unit_nr,&
41 cp_logger_type
42 USE cp_output_handling, ONLY: cp_p_file,&
43 cp_print_key_finished_output,&
44 cp_print_key_should_output,&
45 cp_print_key_unit_nr
46 USE cp_realspace_grid_cube, ONLY: cp_pw_to_cube
47 USE input_section_types, ONLY: section_get_ivals,&
48 section_vals_get,&
49 section_vals_get_subs_vals,&
50 section_vals_type,&
51 section_vals_val_get
52 USE kinds, ONLY: dp
53 USE particle_list_types, ONLY: particle_list_type
54 USE particle_methods, ONLY: get_particle_set
55 USE particle_types, ONLY: particle_type
56 USE physcon, ONLY: boltzmann,&
57 e_charge,&
58 evolt,&
59 h_bar
60 USE pw_env_types, ONLY: pw_env_get,&
61 pw_env_type
62 USE pw_methods, ONLY: pw_zero
63 USE pw_pool_types, ONLY: pw_pool_type
64 USE pw_types, ONLY: pw_c1d_gs_type,&
65 pw_r3d_rs_type
66 USE qs_environment_types, ONLY: get_qs_env,&
67 qs_environment_type,&
68 set_qs_env
69 USE qs_kind_types, ONLY: get_qs_kind,&
70 qs_kind_type
71 USE qs_ks_types, ONLY: qs_ks_env_type
72 USE qs_linres_current, ONLY: calculate_jrho_resp
73 USE qs_linres_types, ONLY: current_env_type
74 USE qs_rho_types, ONLY: qs_rho_get,&
75 qs_rho_type
76 USE qs_subsys_types, ONLY: qs_subsys_get,&
77 qs_subsys_type
78 USE transport_env_types, ONLY: cp2k_csr_interop_type,&
79 cp2k_transport_parameters,&
80 csr_interop_matrix_get_info,&
81 csr_interop_nullify,&
82 transport_env_type
83#include "./base/base_uses.f90"
84
85 IMPLICIT NONE
86 PRIVATE
87
88 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'transport'
89
90 PUBLIC :: transport_env_create, transport_initialize, external_scf_method
91 PUBLIC :: qs_scf_post_transport
92
93!> interface between C/C++ and FORTRAN
94 INTERFACE c_func_interface
95! **************************************************************************************************
96!> \brief C routine that takes the S and H matrices as input and outputs a P matrix
97!> \param cp2k_transport_params transport parameters read form a CP2K input file
98!> \param s_mat C-interoperable overlap matrix
99!> \param ks_mat C-interoperable Kohn-Sham matrix
100!> \param p_mat C-interoperable density matrix
101!> \param imagp_mat C-interoperable imaginary density matrix
102!> \author Mohammad Hossein Bani-Hashemian
103! **************************************************************************************************
104 SUBROUTINE c_scf_routine(cp2k_transport_params, s_mat, ks_mat, p_mat, imagp_mat) BIND(C)
105 IMPORT :: c_int, c_ptr, cp2k_csr_interop_type, cp2k_transport_parameters
106 TYPE(cp2k_transport_parameters), VALUE, INTENT(IN) :: cp2k_transport_params
107 TYPE(cp2k_csr_interop_type), VALUE, INTENT(IN) :: s_mat
108 TYPE(cp2k_csr_interop_type), VALUE, INTENT(IN) :: ks_mat
109 TYPE(cp2k_csr_interop_type), INTENT(INOUT) :: p_mat
110 TYPE(cp2k_csr_interop_type), INTENT(INOUT) :: imagp_mat
111 END SUBROUTINE c_scf_routine
112 END INTERFACE c_func_interface
113
114CONTAINS
115
116! **************************************************************************************************
117!> \brief creates the transport environment
118!> \param[inout] qs_env the qs_env containing the transport_env
119!> \author Mohammad Hossein Bani-Hashemian
120! **************************************************************************************************
121 SUBROUTINE transport_env_create(qs_env)
122 TYPE(qs_environment_type), POINTER :: qs_env
123
124 CHARACTER(LEN=*), PARAMETER :: routinen = 'transport_env_create'
125
126 INTEGER :: handle
127 TYPE(dft_control_type), POINTER :: dft_control
128 TYPE(section_vals_type), POINTER :: input
129 TYPE(transport_env_type), POINTER :: transport_env
130
131 CALL timeset(routinen, handle)
132
133 CALL get_qs_env(qs_env, &
134 transport_env=transport_env, &
135 dft_control=dft_control, &
136 input=input)
137
138 cpassert(.NOT. ASSOCIATED(transport_env))
139
140 ALLOCATE (transport_env)
141
142 CALL transport_init_read_input(input, transport_env)
143 CALL transport_set_contact_params(qs_env, transport_env)
144
145 CALL set_qs_env(qs_env, transport_env=transport_env)
146
147 CALL timestop(handle)
148
149 END SUBROUTINE transport_env_create
150
151! **************************************************************************************************
152!> \brief intitializes all fields of transport_env using the parameters read from
153!> the corresponding input section
154!> \param[inout] input the input file
155!> \param[inout] transport_env the transport_env to be initialized
156!> \author Mohammad Hossein Bani-Hashemian
157! **************************************************************************************************
158 SUBROUTINE transport_init_read_input(input, transport_env)
159 TYPE(section_vals_type), POINTER :: input
160 TYPE(transport_env_type), INTENT(INOUT) :: transport_env
161
162 CHARACTER(len=*), PARAMETER :: routinen = 'transport_init_read_input'
163
164 INTEGER :: contact_bandwidth, contact_injsign, &
165 contact_natoms, contact_start, handle, &
166 i, n_contacts, stride_contacts
167 INTEGER, DIMENSION(:), POINTER :: i_vals
168 LOGICAL :: contact_explicit, injecting_contact, &
169 obc_equilibrium, one_circle
170 TYPE(section_vals_type), POINTER :: beyn_section, contact_section, &
171 pexsi_section, transport_section
172
173 CALL timeset(routinen, handle)
174
175 transport_section => section_vals_get_subs_vals(input, "DFT%TRANSPORT")
176 contact_section => section_vals_get_subs_vals(transport_section, "CONTACT")
177 beyn_section => section_vals_get_subs_vals(transport_section, "BEYN")
178 pexsi_section => section_vals_get_subs_vals(transport_section, "PEXSI")
179 CALL section_vals_get(contact_section, explicit=contact_explicit, n_repetition=n_contacts)
180
181 NULLIFY (i_vals)
182! read from input
183 CALL section_vals_val_get(transport_section, "TRANSPORT_METHOD", i_val=transport_env%params%method)
184 CALL section_vals_val_get(transport_section, "INJECTION_METHOD", i_val=transport_env%params%injection_method)
185 CALL section_vals_val_get(transport_section, "REAL_AXIS_INTEGRATION_METHOD", &
186 i_val=transport_env%params%rlaxis_integration_method)
187 CALL section_vals_val_get(transport_section, "QT_FORMALISM", i_val=transport_env%params%qt_formalism)
188 CALL section_vals_val_get(transport_section, "LINEAR_SOLVER", i_val=transport_env%params%linear_solver)
189 CALL section_vals_val_get(transport_section, "MATRIX_INVERSION_METHOD", &
190 i_val=transport_env%params%matrixinv_method)
191 CALL section_vals_val_get(transport_section, "CONTACT_FILLING", i_val=transport_env%params%transport_neutral)
192 CALL section_vals_val_get(transport_section, "N_KPOINTS", i_val=transport_env%params%n_kpoint)
193 CALL section_vals_val_get(transport_section, "NUM_INTERVAL", i_val=transport_env%params%num_interval)
194 CALL section_vals_val_get(transport_section, "TASKS_PER_ENERGY_POINT", &
195 i_val=transport_env%params%tasks_per_energy_point)
196 CALL section_vals_val_get(transport_section, "TASKS_PER_POLE", i_val=transport_env%params%tasks_per_pole)
197 CALL section_vals_val_get(transport_section, "NUM_POLE", i_val=transport_env%params%num_pole)
198 CALL section_vals_val_get(transport_section, "GPUS_PER_POINT", i_val=transport_env%params%gpus_per_point)
199 CALL section_vals_val_get(transport_section, "N_POINTS_INV", i_val=transport_env%params%n_points_inv)
200 CALL section_vals_val_get(transport_section, "COLZERO_THRESHOLD", r_val=transport_env%params%colzero_threshold)
201 CALL section_vals_val_get(transport_section, "EPS_LIMIT", r_val=transport_env%params%eps_limit)
202 CALL section_vals_val_get(transport_section, "EPS_LIMIT_CC", r_val=transport_env%params%eps_limit_cc)
203 CALL section_vals_val_get(transport_section, "EPS_DECAY", r_val=transport_env%params%eps_decay)
204 CALL section_vals_val_get(transport_section, "EPS_SINGULARITY_CURVATURES", &
205 r_val=transport_env%params%eps_singularity_curvatures)
206 CALL section_vals_val_get(transport_section, "EPS_MU", r_val=transport_env%params%eps_mu)
207 CALL section_vals_val_get(transport_section, "EPS_EIGVAL_DEGEN", r_val=transport_env%params%eps_eigval_degen)
208 CALL section_vals_val_get(transport_section, "EPS_FERMI", r_val=transport_env%params%eps_fermi)
209 CALL section_vals_val_get(transport_section, "ENERGY_INTERVAL", r_val=transport_env%params%energy_interval)
210 CALL section_vals_val_get(transport_section, "MIN_INTERVAL", r_val=transport_env%params%min_interval)
211 CALL section_vals_val_get(transport_section, "TEMPERATURE", r_val=transport_env%params%temperature)
212 CALL section_vals_val_get(transport_section, "DENSITY_MIXING", r_val=transport_env%params%dens_mixing)
213 CALL section_vals_val_get(transport_section, "CSR_SCREENING", l_val=transport_env%csr_screening)
214
215 ! logical*1 to logical*4 , l_val is logical*1 and c_bool is equivalent to logical*4
216 CALL section_vals_val_get(transport_section, "OBC_EQUILIBRIUM", l_val=obc_equilibrium)
217 IF (obc_equilibrium) THEN
218 transport_env%params%obc_equilibrium = .true.
219 ELSE
220 transport_env%params%obc_equilibrium = .false.
221 END IF
222
223 CALL section_vals_val_get(transport_section, "CUTOUT", i_vals=i_vals)
224 transport_env%params%cutout = i_vals
225
226 CALL section_vals_val_get(beyn_section, "TASKS_PER_INTEGRATION_POINT", &
227 i_val=transport_env%params%tasks_per_integration_point)
228 CALL section_vals_val_get(beyn_section, "N_POINTS_BEYN", i_val=transport_env%params%n_points_beyn)
229 CALL section_vals_val_get(beyn_section, "N_RAND", r_val=transport_env%params%n_rand_beyn)
230 CALL section_vals_val_get(beyn_section, "N_RAND_CC", r_val=transport_env%params%n_rand_cc_beyn)
231 CALL section_vals_val_get(beyn_section, "SVD_CUTOFF", r_val=transport_env%params%svd_cutoff)
232 CALL section_vals_val_get(beyn_section, "ONE_CIRCLE", l_val=one_circle)
233 IF (one_circle) THEN
234 transport_env%params%ncrc_beyn = 1
235 ELSE
236 transport_env%params%ncrc_beyn = 2
237 END IF
238
239 CALL section_vals_val_get(pexsi_section, "ORDERING", i_val=transport_env%params%ordering)
240 CALL section_vals_val_get(pexsi_section, "ROW_ORDERING", i_val=transport_env%params%row_ordering)
241 CALL section_vals_val_get(pexsi_section, "VERBOSITY", i_val=transport_env%params%verbosity)
242 CALL section_vals_val_get(pexsi_section, "NP_SYMB_FACT", i_val=transport_env%params%pexsi_np_symb_fact)
243
244 IF (contact_explicit) THEN
245 transport_env%params%num_contacts = n_contacts
246 stride_contacts = 5
247 transport_env%params%stride_contacts = stride_contacts
248 ALLOCATE (transport_env%contacts_data(stride_contacts*n_contacts))
249
250 DO i = 1, n_contacts
251 CALL section_vals_val_get(contact_section, "BANDWIDTH", i_rep_section=i, i_val=contact_bandwidth)
252 CALL section_vals_val_get(contact_section, "START", i_rep_section=i, i_val=contact_start)
253 CALL section_vals_val_get(contact_section, "N_ATOMS", i_rep_section=i, i_val=contact_natoms)
254 CALL section_vals_val_get(contact_section, "INJECTION_SIGN", i_rep_section=i, i_val=contact_injsign)
255 CALL section_vals_val_get(contact_section, "INJECTING_CONTACT", i_rep_section=i, l_val=injecting_contact)
256
257 IF (contact_natoms <= 0) cpabort("Number of atoms in contact region needs to be defined.")
258
259 transport_env%contacts_data((i - 1)*stride_contacts + 1) = contact_bandwidth
260 transport_env%contacts_data((i - 1)*stride_contacts + 2) = contact_start - 1 ! C indexing
261 transport_env%contacts_data((i - 1)*stride_contacts + 3) = contact_natoms
262 transport_env%contacts_data((i - 1)*stride_contacts + 4) = contact_injsign
263
264 IF (injecting_contact) THEN
265 transport_env%contacts_data((i - 1)*stride_contacts + 5) = 1
266 ELSE
267 transport_env%contacts_data((i - 1)*stride_contacts + 5) = 0
268 END IF
269 END DO
270 transport_env%params%contacts_data = c_loc(transport_env%contacts_data(1))
271 ELSE
272 cpabort("No contact region is defined.")
273 END IF
274
275 CALL timestop(handle)
276
277 END SUBROUTINE transport_init_read_input
278
279! **************************************************************************************************
280!> \brief initializes the transport environment
281!> \param ks_env ...
282!> \param[inout] transport_env the transport env to be initialized
283!> \param[in] template_matrix template matrix to keep the sparsity of matrices fixed
284!> \author Mohammad Hossein Bani-Hashemian
285! **************************************************************************************************
286 SUBROUTINE transport_initialize(ks_env, transport_env, template_matrix)
287 TYPE(qs_ks_env_type), POINTER :: ks_env
288 TYPE(transport_env_type), INTENT(INOUT) :: transport_env
289 TYPE(dbcsr_type), INTENT(IN) :: template_matrix
290
291 CHARACTER(len=*), PARAMETER :: routinen = 'transport_initialize'
292
293 INTEGER :: handle, numnodes, unit_nr
294 TYPE(cp_logger_type), POINTER :: logger
295
296 CALL timeset(routinen, handle)
297
298 CALL cite_reference(bruck2014)
299
300 logger => cp_get_default_logger()
301 IF (logger%para_env%is_source()) THEN
302 unit_nr = cp_logger_get_default_unit_nr(logger, local=.true.)
303 ELSE
304 unit_nr = -1
305 END IF
306
307 numnodes = logger%para_env%num_pe
308
309 IF (dbcsr_has_symmetry(template_matrix)) THEN
310 CALL dbcsr_copy(transport_env%template_matrix_sym, template_matrix)
311 CALL dbcsr_desymmetrize(transport_env%template_matrix_sym, transport_env%template_matrix_nosym)
312 ELSE
313 CALL dbcsr_copy(transport_env%template_matrix_nosym, template_matrix)
314 CALL dbcsr_copy(transport_env%template_matrix_sym, template_matrix)
315 END IF
316
317 ALLOCATE (transport_env%dm_imag)
318 CALL dbcsr_create(transport_env%dm_imag, "imaginary DM", &
319 template=template_matrix, matrix_type=dbcsr_type_no_symmetry)
320 CALL dbcsr_set(transport_env%dm_imag, 0.0_dp)
321
322 CALL dbcsr_create(transport_env%csr_sparsity, "CSR sparsity", &
323 template=transport_env%template_matrix_sym)
324 CALL dbcsr_copy(transport_env%csr_sparsity, transport_env%template_matrix_sym)
325
326 CALL cp_dbcsr_to_csr_screening(ks_env, transport_env%csr_sparsity)
327
328 IF (.NOT. transport_env%csr_screening) CALL dbcsr_set(transport_env%csr_sparsity, 1.0_dp)
329 CALL dbcsr_csr_create_from_dbcsr(transport_env%template_matrix_nosym, &
330 transport_env%s_matrix, &
331 dbcsr_csr_dbcsr_blkrow_dist, &
332 csr_sparsity=transport_env%csr_sparsity, &
333 numnodes=numnodes)
334
335 CALL dbcsr_csr_print_sparsity(transport_env%s_matrix, unit_nr)
336
337 CALL dbcsr_convert_dbcsr_to_csr(transport_env%template_matrix_nosym, transport_env%s_matrix)
338
339 CALL dbcsr_csr_create(transport_env%ks_matrix, transport_env%s_matrix)
340 CALL dbcsr_csr_create(transport_env%p_matrix, transport_env%s_matrix)
341 CALL dbcsr_csr_create(transport_env%imagp_matrix, transport_env%s_matrix)
342
343 CALL timestop(handle)
344
345 END SUBROUTINE transport_initialize
346
347! **************************************************************************************************
348!> \brief SCF calcualtion with an externally evaluated density matrix
349!> \param[inout] transport_env transport environment
350!> \param[in] matrix_s DBCSR overlap matrix
351!> \param[in] matrix_ks DBCSR Kohn-Sham matrix
352!> \param[inout] matrix_p DBCSR density matrix
353!> \param[in] nelectron_spin number of electrons
354!> \param[in] natoms number of atoms
355!> \param[in] energy_diff scf energy difference
356!> \param[in] iscf the current scf iteration
357!> \param[in] extra_scf whether or not an extra scf step will be performed
358!> \par History
359!> 12.2012 created [Hossein Bani-Hashemian]
360!> 12.2014 revised [Hossein Bani-Hashemian]
361!> \author Mohammad Hossein Bani-Hashemian
362! **************************************************************************************************
363 SUBROUTINE external_scf_method(transport_env, matrix_s, matrix_ks, matrix_p, &
364 nelectron_spin, natoms, energy_diff, iscf, extra_scf)
365
366 TYPE(transport_env_type), INTENT(INOUT) :: transport_env
367 TYPE(dbcsr_type), INTENT(IN) :: matrix_s, matrix_ks
368 TYPE(dbcsr_type), INTENT(INOUT) :: matrix_p
369 INTEGER, INTENT(IN) :: nelectron_spin, natoms
370 REAL(dp), INTENT(IN) :: energy_diff
371 INTEGER, INTENT(IN) :: iscf
372 LOGICAL, INTENT(IN) :: extra_scf
373
374 CHARACTER(len=*), PARAMETER :: routinen = 'external_scf_method'
375
376 TYPE(cp2k_csr_interop_type) :: imagp_mat, ks_mat, p_mat, s_mat
377
378 PROCEDURE(c_scf_routine), POINTER :: c_method
379 INTEGER :: handle
380
381 CALL timeset(routinen, handle)
382
383 CALL c_f_procpointer(transport_env%ext_c_method_ptr, c_method)
384 IF (.NOT. c_associated(transport_env%ext_c_method_ptr)) THEN
385 CALL cp_abort(__location__, &
386 "MISSING C/C++ ROUTINE: The TRANSPORT section is meant to be used together with an external "// &
387 "program, e.g. the quantum transport code OMEN, that provides CP2K with a density matrix.")
388 END IF
389
390 transport_env%params%n_occ = nelectron_spin
391 transport_env%params%n_atoms = natoms
392 transport_env%params%energy_diff = energy_diff
393 transport_env%params%evoltfactor = evolt
394 transport_env%params%e_charge = e_charge
395 transport_env%params%boltzmann = boltzmann
396 transport_env%params%h_bar = h_bar
397 transport_env%params%iscf = iscf
398 transport_env%params%extra_scf = LOGICAL(extra_scf, C_BOOL) ! C_BOOL and Fortran logical don't have the same size
399
400 CALL csr_interop_nullify(s_mat)
401 CALL csr_interop_nullify(ks_mat)
402 CALL csr_interop_nullify(p_mat)
403 CALL csr_interop_nullify(imagp_mat)
404
405 CALL dbcsr_copy(transport_env%template_matrix_sym, matrix_s, keep_sparsity=.true.)
406 CALL convert_dbcsr_to_csr_interop(transport_env%template_matrix_sym, transport_env%s_matrix, s_mat)
407
408 CALL dbcsr_copy(transport_env%template_matrix_sym, matrix_ks, keep_sparsity=.true.)
409 CALL convert_dbcsr_to_csr_interop(transport_env%template_matrix_sym, transport_env%ks_matrix, ks_mat)
410
411 CALL dbcsr_copy(transport_env%template_matrix_sym, matrix_p, keep_sparsity=.true.)
412 CALL convert_dbcsr_to_csr_interop(transport_env%template_matrix_sym, transport_env%p_matrix, p_mat)
413
414 CALL dbcsr_copy(transport_env%template_matrix_sym, matrix_s, keep_sparsity=.true.)
415 CALL convert_dbcsr_to_csr_interop(transport_env%template_matrix_sym, transport_env%imagp_matrix, imagp_mat)
416
417 CALL c_method(transport_env%params, s_mat, ks_mat, p_mat, imagp_mat)
418
419 CALL convert_csr_interop_to_dbcsr(p_mat, transport_env%p_matrix, transport_env%template_matrix_nosym)
420 CALL dbcsr_copy(matrix_p, transport_env%template_matrix_nosym)
421
422 CALL convert_csr_interop_to_dbcsr(imagp_mat, transport_env%imagp_matrix, transport_env%template_matrix_nosym)
423 CALL dbcsr_copy(transport_env%dm_imag, transport_env%template_matrix_nosym)
424
425 CALL timestop(handle)
426
427 END SUBROUTINE external_scf_method
428
429! **************************************************************************************************
430!> \brief converts a DBCSR matrix to a C-interoperable CSR matrix
431!> \param[in] dbcsr_mat DBCSR matrix to be converted
432!> \param[inout] csr_mat auxiliary CSR matrix
433!> \param[inout] csr_interop_mat C-interoperable CSR matrix
434!> \author Mohammad Hossein Bani-Hashemian
435! **************************************************************************************************
436 SUBROUTINE convert_dbcsr_to_csr_interop(dbcsr_mat, csr_mat, csr_interop_mat)
437
438 TYPE(dbcsr_type), INTENT(IN) :: dbcsr_mat
439 TYPE(dbcsr_csr_type), INTENT(INOUT) :: csr_mat
440 TYPE(cp2k_csr_interop_type), INTENT(INOUT) :: csr_interop_mat
441
442 CHARACTER(LEN=*), PARAMETER :: routinen = 'convert_dbcsr_to_csr_interop'
443
444 INTEGER :: handle
445 INTEGER, ALLOCATABLE, DIMENSION(:) :: nrows_local_all, first_row_all
446 INTEGER(C_INT), DIMENSION(:), POINTER :: colind_local, rowptr_local, nzerow_local
447 REAL(c_double), DIMENSION(:), POINTER :: nzvals_local
448 TYPE(cp_logger_type), POINTER :: logger
449
450 CALL timeset(routinen, handle)
451
452 logger => cp_get_default_logger()
453
454! dbcsr to csr
455 CALL dbcsr_convert_dbcsr_to_csr(dbcsr_mat, csr_mat)
456
457! csr to csr_interop
458 rowptr_local => csr_mat%rowptr_local
459 colind_local => csr_mat%colind_local
460 nzerow_local => csr_mat%nzerow_local
461 nzvals_local => csr_mat%nzval_local%r_dp ! support real double percision for now
462
463 IF (SIZE(rowptr_local) == 0) THEN
464 csr_interop_mat%rowptr_local = c_null_ptr
465 ELSE
466 csr_interop_mat%rowptr_local = c_loc(rowptr_local(1))
467 END IF
468
469 IF (SIZE(colind_local) == 0) THEN
470 csr_interop_mat%colind_local = c_null_ptr
471 ELSE
472 csr_interop_mat%colind_local = c_loc(colind_local(1))
473 END IF
474
475 IF (SIZE(nzerow_local) == 0) THEN
476 csr_interop_mat%nzerow_local = c_null_ptr
477 ELSE
478 csr_interop_mat%nzerow_local = c_loc(nzerow_local(1))
479 END IF
480
481 IF (SIZE(nzvals_local) == 0) THEN
482 csr_interop_mat%nzvals_local = c_null_ptr
483 ELSE
484 csr_interop_mat%nzvals_local = c_loc(nzvals_local(1))
485 END IF
486
487 associate(mp_group => logger%para_env, mepos => logger%para_env%mepos, num_pe => logger%para_env%num_pe)
488 ALLOCATE (nrows_local_all(0:num_pe - 1), first_row_all(0:num_pe - 1))
489 CALL mp_group%allgather(csr_mat%nrows_local, nrows_local_all)
490 CALL cumsum_i(nrows_local_all, first_row_all)
491
492 IF (mepos == 0) THEN
493 csr_interop_mat%first_row = 0
494 ELSE
495 csr_interop_mat%first_row = first_row_all(mepos - 1)
496 END IF
497 END associate
498 csr_interop_mat%nrows_total = csr_mat%nrows_total
499 csr_interop_mat%ncols_total = csr_mat%ncols_total
500 csr_interop_mat%nze_local = csr_mat%nze_local
501 IF (csr_mat%nze_total > huge(csr_interop_mat%nze_total)) THEN
502 cpabort("overflow in nze")
503 END IF
504 csr_interop_mat%nze_total = int(csr_mat%nze_total, kind=kind(csr_interop_mat%nze_total))
505 csr_interop_mat%nrows_local = csr_mat%nrows_local
506 csr_interop_mat%data_type = csr_mat%nzval_local%data_type
507
508 CALL timestop(handle)
509
510 CONTAINS
511! **************************************************************************************************
512!> \brief cumulative sum of a 1d array of integers
513!> \param[in] arr input array
514!> \param[out] cumsum cumulative sum of the input array
515! **************************************************************************************************
516 SUBROUTINE cumsum_i(arr, cumsum)
517 INTEGER, DIMENSION(:), INTENT(IN) :: arr
518 INTEGER, DIMENSION(SIZE(arr)), INTENT(OUT) :: cumsum
519
520 INTEGER :: i
521
522 cumsum(1) = arr(1)
523 DO i = 2, SIZE(arr)
524 cumsum(i) = cumsum(i - 1) + arr(i)
525 END DO
526 END SUBROUTINE cumsum_i
527
528 END SUBROUTINE convert_dbcsr_to_csr_interop
529
530! **************************************************************************************************
531!> \brief converts a C-interoperable CSR matrix to a DBCSR matrix
532!> \param[in] csr_interop_mat C-interoperable CSR matrix
533!> \param[inout] csr_mat auxiliary CSR matrix
534!> \param[inout] dbcsr_mat DBCSR matrix
535!> \author Mohammad Hossein Bani-Hashemian
536! **************************************************************************************************
537 SUBROUTINE convert_csr_interop_to_dbcsr(csr_interop_mat, csr_mat, dbcsr_mat)
538
539 TYPE(cp2k_csr_interop_type), INTENT(IN) :: csr_interop_mat
540 TYPE(dbcsr_csr_type), INTENT(INOUT) :: csr_mat
541 TYPE(dbcsr_type), INTENT(INOUT) :: dbcsr_mat
542
543 CHARACTER(LEN=*), PARAMETER :: routinen = 'convert_csr_interop_to_dbcsr'
544
545 INTEGER :: data_type, handle, ncols_total, &
546 nrows_local, nrows_total, nze_local, &
547 nze_total
548 INTEGER, DIMENSION(:), POINTER :: colind_local, nzerow_local, rowptr_local
549 REAL(dp), DIMENSION(:), POINTER :: nzvals_local
550
551 CALL timeset(routinen, handle)
552
553! csr_interop to csr
554 CALL csr_interop_matrix_get_info(csr_interop_mat, &
555 nrows_total=nrows_total, ncols_total=ncols_total, nze_local=nze_local, &
556 nze_total=nze_total, nrows_local=nrows_local, data_type=data_type, &
557 rowptr_local=rowptr_local, colind_local=colind_local, &
558 nzerow_local=nzerow_local, nzvals_local=nzvals_local)
559
560 csr_mat%nrows_total = nrows_total
561 csr_mat%ncols_total = ncols_total
562 csr_mat%nze_local = nze_local
563 csr_mat%nze_total = nze_total
564 csr_mat%nrows_local = nrows_local
565 csr_mat%nzval_local%data_type = data_type
566
567 csr_mat%rowptr_local = rowptr_local
568 csr_mat%colind_local = colind_local
569 csr_mat%nzerow_local = nzerow_local
570 csr_mat%nzval_local%r_dp = nzvals_local
571
572! csr to dbcsr
573 CALL dbcsr_convert_csr_to_dbcsr(dbcsr_mat, csr_mat)
574
575 CALL timestop(handle)
576
577 END SUBROUTINE convert_csr_interop_to_dbcsr
578
579! **************************************************************************************************
580!> \brief extraxts zeff (effective nuclear charges per atom) and nsgf (the size
581!> of spherical Gaussian basis functions per atom) from qs_env and initializes
582!> the corresponding arrays in transport_env%params
583!> \param[in] qs_env qs environment
584!> \param[inout] transport_env transport environment
585!> \author Mohammad Hossein Bani-Hashemian
586! **************************************************************************************************
587 SUBROUTINE transport_set_contact_params(qs_env, transport_env)
588 TYPE(qs_environment_type), POINTER :: qs_env
589 TYPE(transport_env_type), INTENT(INOUT) :: transport_env
590
591 INTEGER :: i, iat, ikind, natom, nkind
592 INTEGER, DIMENSION(:), POINTER :: atom_list
593 REAL(kind=dp) :: zeff
594 TYPE(atomic_kind_type), DIMENSION(:), POINTER :: atomic_kind_set
595 TYPE(atomic_kind_type), POINTER :: atomic_kind
596 TYPE(particle_type), DIMENSION(:), POINTER :: particle_set
597 TYPE(qs_kind_type), DIMENSION(:), POINTER :: qs_kind_set
598
599 CALL get_qs_env(qs_env, nkind=nkind, natom=natom)
600 CALL get_qs_env(qs_env, particle_set=particle_set, &
601 qs_kind_set=qs_kind_set, &
602 atomic_kind_set=atomic_kind_set)
603
604 ALLOCATE (transport_env%nsgf(natom))
605 ALLOCATE (transport_env%zeff(natom))
606 CALL get_particle_set(particle_set, qs_kind_set, nsgf=transport_env%nsgf)
607
608 ! reference charges
609 DO ikind = 1, nkind
610 CALL get_qs_kind(qs_kind_set(ikind), zeff=zeff)
611 atomic_kind => atomic_kind_set(ikind)
612 CALL get_atomic_kind(atomic_kind, atom_list=atom_list)
613 DO iat = 1, SIZE(atom_list)
614 i = atom_list(iat)
615 transport_env%zeff(i) = zeff
616 END DO
617 END DO
618
619 IF (natom == 0) THEN
620 transport_env%params%nsgf = c_null_ptr
621 transport_env%params%zeff = c_null_ptr
622 ELSE
623 transport_env%params%nsgf = c_loc(transport_env%nsgf(1))
624 transport_env%params%zeff = c_loc(transport_env%zeff(1))
625 END IF
626
627 END SUBROUTINE transport_set_contact_params
628
629!**************************************************************************************************
630!> \brief evaluates current density using the imaginary part of the density matrix and prints it
631!> into cube files
632!> \param[in] input the input
633!> \param[in] qs_env qs environment
634!> \author Mohammad Hossein Bani-Hashemian
635! **************************************************************************************************
636 SUBROUTINE transport_current(input, qs_env)
637 TYPE(section_vals_type), POINTER :: input
638 TYPE(qs_environment_type), POINTER :: qs_env
639
640 CHARACTER(len=*), PARAMETER :: routinen = 'transport_current'
641
642 CHARACTER(len=14) :: ext
643 CHARACTER(len=2) :: sdir
644 INTEGER :: dir, handle, unit_nr
645 LOGICAL :: do_current_cube, do_transport, mpi_io
646 TYPE(cp_logger_type), POINTER :: logger
647 TYPE(current_env_type) :: current_env
648 TYPE(dbcsr_type), POINTER :: zero
649 TYPE(dft_control_type), POINTER :: dft_control
650 TYPE(particle_list_type), POINTER :: particles
651 TYPE(pw_c1d_gs_type) :: gs
652 TYPE(pw_c1d_gs_type), DIMENSION(:), POINTER :: rho_g
653 TYPE(pw_env_type), POINTER :: pw_env
654 TYPE(pw_pool_type), POINTER :: auxbas_pw_pool
655 TYPE(pw_r3d_rs_type) :: rs
656 TYPE(pw_r3d_rs_type), DIMENSION(:), POINTER :: rho_r
657 TYPE(qs_rho_type), POINTER :: rho
658 TYPE(qs_subsys_type), POINTER :: subsys
659 TYPE(section_vals_type), POINTER :: dft_section
660 TYPE(transport_env_type), POINTER :: transport_env
661
662 CALL timeset(routinen, handle)
663
664 logger => cp_get_default_logger()
665 dft_section => section_vals_get_subs_vals(input, "DFT")
666 CALL get_qs_env(qs_env=qs_env, &
667 subsys=subsys, &
668 pw_env=pw_env, &
669 transport_env=transport_env, &
670 do_transport=do_transport, &
671 dft_control=dft_control, &
672 rho=rho)
673 CALL qs_subsys_get(subsys, particles=particles)
674 CALL pw_env_get(pw_env, auxbas_pw_pool=auxbas_pw_pool)
675 CALL qs_rho_get(rho, rho_r=rho_r, rho_g=rho_g)
676
677 do_current_cube = btest(cp_print_key_should_output(logger%iter_info, input, &
678 "DFT%TRANSPORT%PRINT%CURRENT"), cp_p_file)
679
680 ! check if transport is active i.e. the imaginary dm has been actually passed to cp2k by the external code
681 IF (do_transport) THEN
682 IF (c_associated(transport_env%ext_c_method_ptr)) THEN
683
684 ! calculate_jrho_resp uses sab_all which is not associated in DFTB environment
685 IF (dft_control%qs_control%dftb) cpabort("Current is not available for DFTB.")
686 IF (dft_control%qs_control%xtb) cpabort("Current is not available for xTB.")
687
688 ! now do the current
689 IF (do_current_cube) THEN
690 current_env%gauge = -1
691 current_env%gauge_init = .false.
692
693 CALL auxbas_pw_pool%create_pw(rs)
694 CALL auxbas_pw_pool%create_pw(gs)
695
696 NULLIFY (zero)
697 ALLOCATE (zero)
698 CALL dbcsr_create(zero, template=transport_env%dm_imag)
699 CALL dbcsr_copy(zero, transport_env%dm_imag)
700 CALL dbcsr_set(zero, 0.0_dp)
701
702 DO dir = 1, 3
703 CALL pw_zero(rs)
704 CALL pw_zero(gs)
705 CALL calculate_jrho_resp(zero, transport_env%dm_imag, &
706 zero, zero, dir, dir, rs, gs, qs_env, current_env, &
707 retain_rsgrid=.true.)
708 SELECT CASE (dir)
709 CASE (1)
710 sdir = "-x"
711 CASE (2)
712 sdir = "-y"
713 CASE (3)
714 sdir = "-z"
715 END SELECT
716 ext = sdir//".cube"
717 mpi_io = .true.
718 unit_nr = cp_print_key_unit_nr(logger, dft_section, "TRANSPORT%PRINT%CURRENT", &
719 extension=ext, file_status="REPLACE", file_action="WRITE", &
720 log_filename=.false., mpi_io=mpi_io)
721 CALL cp_pw_to_cube(rs, unit_nr, "Transport current", particles=particles, &
722 stride=section_get_ivals(dft_section, "TRANSPORT%PRINT%CURRENT%STRIDE"), &
723 mpi_io=mpi_io)
724 CALL cp_print_key_finished_output(unit_nr, logger, dft_section, "TRANSPORT%PRINT%CURRENT", &
725 mpi_io=mpi_io)
726 END DO
727
728 CALL dbcsr_deallocate_matrix(zero)
729 CALL auxbas_pw_pool%give_back_pw(rs)
730 CALL auxbas_pw_pool%give_back_pw(gs)
731 END IF
732 END IF
733
734 END IF
735
736 CALL timestop(handle)
737
738 END SUBROUTINE transport_current
739
740!**************************************************************************************************
741!> \brief post scf calculations for transport
742!> \param[in] qs_env qs environment
743!> \author Mohammad Hossein Bani-Hashemian
744! **************************************************************************************************
745 SUBROUTINE qs_scf_post_transport(qs_env)
746 TYPE(qs_environment_type), POINTER :: qs_env
747
748 CHARACTER(len=*), PARAMETER :: routinen = 'qs_scf_post_transport'
749
750 INTEGER :: handle
751 TYPE(section_vals_type), POINTER :: input
752
753 CALL timeset(routinen, handle)
754
755 NULLIFY (input)
756
757 CALL get_qs_env(qs_env=qs_env, &
758 input=input)
759
760 CALL transport_current(input, qs_env)
761
762 CALL timestop(handle)
763
764 END SUBROUTINE qs_scf_post_transport
765
766END MODULE transport
cp_dbcsr_api::dbcsr_create
Definition cp_dbcsr_api.F:194
pw_methods::pw_zero
Definition pw_methods.F:83
atomic_kind_types
Define the atomic kind types and their sub types.
Definition atomic_kind_types.F:23
atomic_kind_types::get_atomic_kind
subroutine, public get_atomic_kind(atomic_kind, fist_potential, element_symbol, name, mass, kind_number, natom, atom_list, rcov, rvdw, z, qeff, apol, cpol, mm_radius, shell, shell_active, damping)
Get attributes of an atomic kind.
Definition atomic_kind_types.F:138
bibliography
collects all references to literature in CP2K as new algorithms / method are included from literature...
Definition bibliography.F:21
bibliography::bruck2014
integer, save, public bruck2014
Definition bibliography.F:36
cp_control_types
Defines control structures, which contain the parameters and the settings for the DFT-based calculati...
Definition cp_control_types.F:12
cp_dbcsr_api
Definition cp_dbcsr_api.F:8
cp_dbcsr_api::dbcsr_has_symmetry
logical function, public dbcsr_has_symmetry(matrix)
...
Definition cp_dbcsr_api.F:928
cp_dbcsr_api::dbcsr_convert_dbcsr_to_csr
subroutine, public dbcsr_convert_dbcsr_to_csr(dbcsr_mat, csr_mat)
...
Definition cp_dbcsr_api.F:352
cp_dbcsr_api::dbcsr_deallocate_matrix
subroutine, public dbcsr_deallocate_matrix(matrix)
...
Definition cp_dbcsr_api.F:235
cp_dbcsr_api::dbcsr_convert_csr_to_dbcsr
subroutine, public dbcsr_convert_csr_to_dbcsr(dbcsr_mat, csr_mat)
...
Definition cp_dbcsr_api.F:336
cp_dbcsr_api::dbcsr_desymmetrize
subroutine, public dbcsr_desymmetrize(matrix_a, matrix_b)
...
Definition cp_dbcsr_api.F:525
cp_dbcsr_api::dbcsr_copy
subroutine, public dbcsr_copy(matrix_b, matrix_a, name, keep_sparsity, keep_imaginary)
...
Definition cp_dbcsr_api.F:371
cp_dbcsr_api::dbcsr_csr_create_from_dbcsr
subroutine, public dbcsr_csr_create_from_dbcsr(dbcsr_mat, csr_mat, dist_format, csr_sparsity, numnodes)
...
Definition cp_dbcsr_api.F:469
cp_dbcsr_api::dbcsr_set
subroutine, public dbcsr_set(matrix, alpha)
...
Definition cp_dbcsr_api.F:1204
cp_dbcsr_cp2k_link
Routines that link DBCSR and CP2K concepts together.
Definition cp_dbcsr_cp2k_link.F:14
cp_dbcsr_cp2k_link::cp_dbcsr_to_csr_screening
subroutine, public cp_dbcsr_to_csr_screening(ks_env, csr_sparsity)
Apply distance screening to refine sparsity pattern of matrices in CSR format (using eps_pgf_orb)....
Definition cp_dbcsr_cp2k_link.F:560
cp_log_handling
various routines to log and control the output. The idea is that decisions about where to log should ...
Definition cp_log_handling.F:41
cp_log_handling::cp_logger_get_default_unit_nr
recursive integer function, public cp_logger_get_default_unit_nr(logger, local, skip_not_ionode)
asks the default unit number of the given logger. try to use cp_logger_get_unit_nr
Definition cp_log_handling.F:576
cp_log_handling::cp_get_default_logger
type(cp_logger_type) function, pointer, public cp_get_default_logger()
returns the default logger
Definition cp_log_handling.F:234
cp_output_handling
routines to handle the output, The idea is to remove the decision of wheter to output and what to out...
Definition cp_output_handling.F:25
cp_output_handling::cp_print_key_unit_nr
integer function, public cp_print_key_unit_nr(logger, basis_section, print_key_path, extension, middle_name, local, log_filename, ignore_should_output, file_form, file_position, file_action, file_status, do_backup, on_file, is_new_file, mpi_io, fout)
...
Definition cp_output_handling.F:853
cp_output_handling::cp_print_key_finished_output
subroutine, public cp_print_key_finished_output(unit_nr, logger, basis_section, print_key_path, local, ignore_should_output, on_file, mpi_io)
should be called after you finish working with a unit obtained with cp_print_key_unit_nr,...
Definition cp_output_handling.F:1073
cp_output_handling::cp_p_file
integer, parameter, public cp_p_file
Definition cp_output_handling.F:103
cp_output_handling::cp_print_key_should_output
integer function, public cp_print_key_should_output(iteration_info, basis_section, print_key_path, used_print_key, first_time)
returns what should be done with the given property if btest(res,cp_p_store) then the property should...
Definition cp_output_handling.F:345
cp_realspace_grid_cube
A wrapper around pw_to_cube() which accepts particle_list_type.
Definition cp_realspace_grid_cube.F:12
cp_realspace_grid_cube::cp_pw_to_cube
subroutine, public cp_pw_to_cube(pw, unit_nr, title, particles, zeff, stride, max_file_size_mb, zero_tails, silent, mpi_io)
...
Definition cp_realspace_grid_cube.F:48
input_section_types
objects that represent the structure of input sections and the data contained in an input section
Definition input_section_types.F:15
input_section_types::section_get_ivals
integer function, dimension(:), pointer, public section_get_ivals(section_vals, keyword_name)
...
Definition input_section_types.F:1001
input_section_types::section_vals_get_subs_vals
recursive type(section_vals_type) function, pointer, public section_vals_get_subs_vals(section_vals, subsection_name, i_rep_section, can_return_null)
returns the values of the requested subsection
Definition input_section_types.F:735
input_section_types::section_vals_get
subroutine, public section_vals_get(section_vals, ref_count, n_repetition, n_subs_vals_rep, section, explicit)
returns various attributes about the section_vals
Definition input_section_types.F:708
input_section_types::section_vals_val_get
subroutine, public section_vals_val_get(section_vals, keyword_name, i_rep_section, i_rep_val, n_rep_val, val, l_val, i_val, r_val, c_val, l_vals, i_vals, r_vals, c_vals, explicit)
returns the requested value
Definition input_section_types.F:1052
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
particle_list_types
represent a simple array based list of the given type
Definition particle_list_types.F:15
particle_methods
Define methods related to particle_type.
Definition particle_methods.F:14
particle_methods::get_particle_set
subroutine, public get_particle_set(particle_set, qs_kind_set, first_sgf, last_sgf, nsgf, nmao, basis, ncgf)
Get the components of a particle set.
Definition particle_methods.F:123
particle_types
Define the data structure for the particle information.
Definition particle_types.F:19
physcon
Definition of physical constants:
Definition physcon.F:68
physcon::boltzmann
real(kind=dp), parameter, public boltzmann
Definition physcon.F:129
physcon::e_charge
real(kind=dp), parameter, public e_charge
Definition physcon.F:106
physcon::h_bar
real(kind=dp), parameter, public h_bar
Definition physcon.F:103
physcon::evolt
real(kind=dp), parameter, public evolt
Definition physcon.F:183
pw_env_types
container for various plainwaves related things
Definition pw_env_types.F:14
pw_env_types::pw_env_get
subroutine, public pw_env_get(pw_env, pw_pools, cube_info, gridlevel_info, auxbas_pw_pool, auxbas_grid, auxbas_rs_desc, auxbas_rs_grid, rs_descs, rs_grids, xc_pw_pool, vdw_pw_pool, poisson_env, interp_section)
returns the various attributes of the pw env
Definition pw_env_types.F:113
pw_methods
Definition pw_methods.F:25
pw_pool_types
Manages a pool of grids (to be used for example as tmp objects), but can also be used to instantiate ...
Definition pw_pool_types.F:24
pw_types
Definition pw_types.F:24
qs_environment_types
Definition qs_environment_types.F:14
qs_environment_types::get_qs_env
subroutine, public get_qs_env(qs_env, atomic_kind_set, qs_kind_set, cell, super_cell, cell_ref, use_ref_cell, kpoints, dft_control, mos, sab_orb, sab_all, qmmm, qmmm_periodic, mimic, sac_ae, sac_ppl, sac_lri, sap_ppnl, sab_vdw, sab_scp, sap_oce, sab_lrc, sab_se, sab_xtbe, sab_tbe, sab_core, sab_xb, sab_xtb_pp, sab_xtb_nonbond, sab_almo, sab_kp, sab_kp_nosym, sab_cneo, particle_set, energy, force, matrix_h, matrix_h_im, matrix_ks, matrix_ks_im, matrix_vxc, run_rtp, rtp, matrix_h_kp, matrix_h_im_kp, matrix_ks_kp, matrix_ks_im_kp, matrix_vxc_kp, kinetic_kp, matrix_s_kp, matrix_w_kp, matrix_s_ri_aux_kp, matrix_s, matrix_s_ri_aux, matrix_w, matrix_p_mp2, matrix_p_mp2_admm, matrix_vhxc, rho, rho_xc, pw_env, ewald_env, ewald_pw, active_space, mpools, input, para_env, blacs_env, scf_control, rel_control, kinetic, qs_charges, vppl, xcint_weights, rho_core, rho_nlcc, rho_nlcc_g, ks_env, ks_qmmm_env, wf_history, scf_env, local_particles, local_molecules, distribution_2d, dbcsr_dist, molecule_kind_set, molecule_set, subsys, cp_subsys, oce, local_rho_set, rho_atom_set, task_list, task_list_soft, rho0_atom_set, rho0_mpole, rhoz_set, rhoz_cneo_set, ecoul_1c, rho0_s_rs, rho0_s_gs, rhoz_cneo_s_rs, rhoz_cneo_s_gs, do_kpoints, has_unit_metric, requires_mo_derivs, mo_derivs, mo_loc_history, nkind, natom, nelectron_total, nelectron_spin, efield, neighbor_list_id, linres_control, xas_env, virial, cp_ddapc_env, cp_ddapc_ewald, outer_scf_history, outer_scf_ihistory, x_data, et_coupling, dftb_potential, results, se_taper, se_store_int_env, se_nddo_mpole, se_nonbond_env, admm_env, lri_env, lri_density, exstate_env, ec_env, harris_env, dispersion_env, gcp_env, vee, rho_external, external_vxc, mask, mp2_env, bs_env, kg_env, wanniercentres, atprop, ls_scf_env, do_transport, transport_env, v_hartree_rspace, s_mstruct_changed, rho_changed, potential_changed, forces_up_to_date, mscfg_env, almo_scf_env, gradient_history, variable_history, embed_pot, spin_embed_pot, polar_env, mos_last_converged, eeq, rhs, do_rixs, tb_tblite)
Get the QUICKSTEP environment.
Definition qs_environment_types.F:539
qs_environment_types::set_qs_env
subroutine, public set_qs_env(qs_env, super_cell, mos, qmmm, qmmm_periodic, mimic, ewald_env, ewald_pw, mpools, rho_external, external_vxc, mask, scf_control, rel_control, qs_charges, ks_env, ks_qmmm_env, wf_history, scf_env, active_space, input, oce, rho_atom_set, rho0_atom_set, rho0_mpole, run_rtp, rtp, rhoz_set, rhoz_tot, ecoul_1c, has_unit_metric, requires_mo_derivs, mo_derivs, mo_loc_history, efield, rhoz_cneo_set, linres_control, xas_env, cp_ddapc_env, cp_ddapc_ewald, outer_scf_history, outer_scf_ihistory, x_data, et_coupling, dftb_potential, se_taper, se_store_int_env, se_nddo_mpole, se_nonbond_env, admm_env, ls_scf_env, do_transport, transport_env, lri_env, lri_density, exstate_env, ec_env, dispersion_env, harris_env, gcp_env, mp2_env, bs_env, kg_env, force, kpoints, wanniercentres, almo_scf_env, gradient_history, variable_history, embed_pot, spin_embed_pot, polar_env, mos_last_converged, eeq, rhs, do_rixs, tb_tblite)
Set the QUICKSTEP environment.
Definition qs_environment_types.F:1128
qs_kind_types
Define the quickstep kind type and their sub types.
Definition qs_kind_types.F:23
qs_kind_types::get_qs_kind
subroutine, public get_qs_kind(qs_kind, basis_set, basis_type, ncgf, nsgf, all_potential, tnadd_potential, gth_potential, sgp_potential, upf_potential, cneo_potential, se_parameter, dftb_parameter, xtb_parameter, dftb3_param, zatom, zeff, elec_conf, mao, lmax_dftb, alpha_core_charge, ccore_charge, core_charge, core_charge_radius, paw_proj_set, paw_atom, hard_radius, hard0_radius, max_rad_local, covalent_radius, vdw_radius, gpw_type_forced, harmonics, max_iso_not0, max_s_harm, grid_atom, ngrid_ang, ngrid_rad, lmax_rho0, dft_plus_u_atom, l_of_dft_plus_u, n_of_dft_plus_u, u_minus_j, hund_j, u_of_dft_plus_u, j_of_dft_plus_u, alpha_of_dft_plus_u, beta_of_dft_plus_u, j0_of_dft_plus_u, occupation_of_dft_plus_u, dispersion, bs_occupation, magnetization, no_optimize, addel, laddel, naddel, orbitals, max_scf, eps_scf, smear, u_ramping, u_minus_j_target, eps_u_ramping, proj_shell_charge, lr_atom, do_mtlr, u_j_loop, ao_coef, init_u_ramping_each_scf, reltmat, ghost, monovalent, floating, name, element_symbol, pao_basis_size, pao_model_file, pao_potentials, pao_descriptors, nelec)
Get attributes of an atomic kind.
Definition qs_kind_types.F:491
qs_ks_types
Definition qs_ks_types.F:15
qs_linres_current
given the response wavefunctions obtained by the application of the (rxp), p, and ((dk-dl)xp) operato...
Definition qs_linres_current.F:17
qs_linres_current::calculate_jrho_resp
subroutine, public calculate_jrho_resp(mat_d0, mat_jp, mat_jp_rii, mat_jp_riii, ib, idir, current_rs, current_gs, qs_env, current_env, soft_valid, retain_rsgrid)
Calculation of the idir component of the response current density in the presence of a constant magne...
Definition qs_linres_current.F:578
qs_linres_types
Type definitiona for linear response calculations.
Definition qs_linres_types.F:12
qs_rho_types
superstucture that hold various representations of the density and keeps track of which ones are vali...
Definition qs_rho_types.F:18
qs_rho_types::qs_rho_get
subroutine, public qs_rho_get(rho_struct, rho_ao, rho_ao_im, rho_ao_kp, rho_ao_im_kp, rho_r, drho_r, rho_g, drho_g, tau_r, tau_g, rho_r_valid, drho_r_valid, rho_g_valid, drho_g_valid, tau_r_valid, tau_g_valid, tot_rho_r, tot_rho_g, rho_r_sccs, soft_valid, complex_rho_ao)
returns info about the density described by this object. If some representation is not available an e...
Definition qs_rho_types.F:229
qs_subsys_types
types that represent a quickstep subsys
Definition qs_subsys_types.F:12
qs_subsys_types::qs_subsys_get
subroutine, public qs_subsys_get(subsys, atomic_kinds, atomic_kind_set, particles, particle_set, local_particles, molecules, molecule_set, molecule_kinds, molecule_kind_set, local_molecules, para_env, colvar_p, shell_particles, core_particles, gci, multipoles, natom, nparticle, ncore, nshell, nkind, atprop, virial, results, cell, cell_ref, use_ref_cell, energy, force, qs_kind_set, cp_subsys, nelectron_total, nelectron_spin)
...
Definition qs_subsys_types.F:137
transport_env_types
CP2K transport environment and related C-interoperable types.
Definition transport_env_types.F:17
transport_env_types::csr_interop_nullify
subroutine, public csr_interop_nullify(csr_interop_mat)
nullifies (and zeroizes) a C-interoperable CSR matrix
Definition transport_env_types.F:178
transport_env_types::csr_interop_matrix_get_info
subroutine, public csr_interop_matrix_get_info(csr_interop_mat, nrows_total, ncols_total, nze_local, nze_total, nrows_local, data_type, first_row, rowptr_local, colind_local, nzerow_local, nzvals_local)
gets the fields of a C-interoperable CSR matrix
Definition transport_env_types.F:222
transport
routines for DFT+NEGF calculations (coupling with the quantum transport code OMEN)
Definition transport.F:19
transport::transport_env_create
subroutine, public transport_env_create(qs_env)
creates the transport environment
Definition transport.F:122
transport::external_scf_method
subroutine, public external_scf_method(transport_env, matrix_s, matrix_ks, matrix_p, nelectron_spin, natoms, energy_diff, iscf, extra_scf)
SCF calcualtion with an externally evaluated density matrix.
Definition transport.F:365
transport::transport_initialize
subroutine, public transport_initialize(ks_env, transport_env, template_matrix)
initializes the transport environment
Definition transport.F:287
transport::qs_scf_post_transport
subroutine, public qs_scf_post_transport(qs_env)
post scf calculations for transport
Definition transport.F:746
atomic_kind_types::atomic_kind_type
Provides all information about an atomic kind.
Definition atomic_kind_types.F:45
cp2k_csr_interop_type
CP2K's C-interoperable CSR matrix This definition matches the respective type definition in the trans...
Definition libcp2k.h:277
cp2k_transport_parameters
Transport parameters read from a CP2K input file. This definition matches the respective type definit...
Definition libcp2k.h:217
cp_control_types::dft_control_type
Definition cp_control_types.F:758
cp_dbcsr_api::dbcsr_type
Definition cp_dbcsr_api.F:176
cp_log_handling::cp_logger_type
type of a logger, at the moment it contains just a print level starting at which level it should be l...
Definition cp_log_handling.F:140
input_section_types::section_vals_type
stores the values of a section
Definition input_section_types.F:127
particle_list_types::particle_list_type
represent a list of objects
Definition particle_list_types.F:46
particle_types::particle_type
Definition particle_types.F:35
pw_env_types::pw_env_type
contained for different pw related things
Definition pw_env_types.F:70
pw_pool_types::pw_pool_type
Manages a pool of grids (to be used for example as tmp objects), but can also be used to instantiate ...
Definition pw_pool_types.F:83
pw_types::pw_c1d_gs_type
Definition pw_types.F:127
pw_types::pw_r3d_rs_type
Definition pw_types.F:62
qs_environment_types::qs_environment_type
Definition qs_environment_types.F:222
qs_kind_types::qs_kind_type
Provides all information about a quickstep kind.
Definition qs_kind_types.F:190
qs_ks_types::qs_ks_env_type
calculation environment to calculate the ks matrix, holds all the needed vars. assumes that the core ...
Definition qs_ks_types.F:137
qs_linres_types::current_env_type
Definition qs_linres_types.F:94
qs_rho_types::qs_rho_type
keeps the density in various representations, keeping track of which ones are valid.
Definition qs_rho_types.F:62
qs_subsys_types::qs_subsys_type
Definition qs_subsys_types.F:56
transport_env_types::transport_env_type
Definition transport_env_types.F:95