(git:3f0d10b)
Loading...
Searching...
No Matches
input_cp2k_negf.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 Input section for NEGF based quantum transport calculations.
10! **************************************************************************************************
11
12MODULE input_cp2k_negf
13 USE bibliography, ONLY: bailey2006,&
14 papior2017
15 USE cp_output_handling, ONLY: add_last_numeric,&
16 cp_print_key_section_create,&
17 debug_print_level,&
18 high_print_level,&
19 low_print_level,&
20 medium_print_level,&
21 silent_print_level
22 USE input_constants, ONLY: negfint_method_cc,&
23 negfint_method_simpson
24 USE input_keyword_types, ONLY: keyword_create,&
25 keyword_release,&
26 keyword_type
27 USE input_section_types, ONLY: section_add_keyword,&
28 section_add_subsection,&
29 section_create,&
30 section_release,&
31 section_type
32 USE input_val_types, ONLY: char_t,&
33 enum_t,&
34 integer_t,&
35 real_t
36 USE kinds, ONLY: dp
37 USE physcon, ONLY: kelvin
38 USE qs_density_mixing_types, ONLY: create_mixing_section
39 USE string_utilities, ONLY: s2a
40#include "./base/base_uses.f90"
41
42 IMPLICIT NONE
43 PRIVATE
44
45 CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'input_cp2k_negf'
46
47 PUBLIC :: create_negf_section
48
49CONTAINS
50
51! **************************************************************************************************
52!> \brief Create NEGF input section.
53!> \param section input section
54!> \par History
55!> * 02.2017 created [Sergey Chulkov]
56! **************************************************************************************************
57 SUBROUTINE create_negf_section(section)
58 TYPE(section_type), POINTER :: section
59
60 TYPE(keyword_type), POINTER :: keyword
61 TYPE(section_type), POINTER :: print_key, subsection
62
63 cpassert(.NOT. ASSOCIATED(section))
64 CALL section_create(section, __location__, name="NEGF", &
65 description="Parameters which control quantum transport calculation"// &
66 " based on Non-Equilibrium Green Function method.", &
67 citations=[bailey2006, papior2017], &
68 n_keywords=18, n_subsections=6, repeats=.false.)
69
70 NULLIFY (keyword, print_key, subsection)
71
72 CALL create_contact_section(subsection)
73 CALL section_add_subsection(section, subsection)
74 CALL section_release(subsection)
75
76 CALL create_scattering_region_section(subsection)
77 CALL section_add_subsection(section, subsection)
78 CALL section_release(subsection)
79
80 ! mixing section
81 CALL create_mixing_section(subsection, ls_scf=.false.)
82 CALL section_add_subsection(section, subsection)
83 CALL section_release(subsection)
84
85 CALL keyword_create(keyword, __location__, name="DISABLE_CACHE", &
86 description="Do not keep contact self-energy matrices for future reuse", &
87 default_l_val=.false., lone_keyword_l_val=.true.)
88 CALL section_add_keyword(section, keyword)
89 CALL keyword_release(keyword)
90
91 ! convergence thresholds
92 CALL keyword_create(keyword, __location__, name="EPS_DENSITY", &
93 description="Target accuracy for electronic density.", &
94 n_var=1, type_of_var=real_t, default_r_val=1.0e-5_dp)
95 CALL section_add_keyword(section, keyword)
96 CALL keyword_release(keyword)
97
98 CALL keyword_create(keyword, __location__, name="EPS_GREEN", &
99 description="Target accuracy for surface Green's functions.", &
100 n_var=1, type_of_var=real_t, default_r_val=1.0e-5_dp)
101 CALL section_add_keyword(section, keyword)
102 CALL keyword_release(keyword)
103
104 CALL keyword_create(keyword, __location__, name="EPS_SCF", &
105 description="Target accuracy for SCF convergence.", &
106 n_var=1, type_of_var=real_t, default_r_val=1.0e-5_dp)
107 CALL section_add_keyword(section, keyword)
108 CALL keyword_release(keyword)
109
110 CALL keyword_create(keyword, __location__, name="EPS_GEO", &
111 description="Accuracy in mapping atoms between different force environments.", &
112 n_var=1, type_of_var=real_t, unit_str="angstrom", &
113 default_r_val=1.0e-6_dp)
114 CALL section_add_keyword(section, keyword)
115 CALL keyword_release(keyword)
116
117 CALL keyword_create(keyword, __location__, name="ENERGY_LBOUND", &
118 description="Lower bound energy of the conductance band.", &
119 n_var=1, type_of_var=real_t, unit_str="hartree", &
120 default_r_val=-5.0_dp)
121 CALL section_add_keyword(section, keyword)
122 CALL keyword_release(keyword)
123
124 CALL keyword_create(keyword, __location__, name="ETA", &
125 description="Infinitesimal offset from the real axis.", &
126 n_var=1, type_of_var=real_t, unit_str="hartree", &
127 default_r_val=1.0e-5_dp)
128 CALL section_add_keyword(section, keyword)
129 CALL keyword_release(keyword)
130
131 CALL keyword_create(keyword, __location__, name="HOMO_LUMO_GAP", &
132 description="The gap between the HOMO and some fictitious LUMO. This option is used as"// &
133 " an initial offset to determine the actual Fermi level of bulk contacts."// &
134 " It does not need to be exact HOMO-LUMO gap, just some value to start with.", &
135 n_var=1, type_of_var=real_t, unit_str="hartree", &
136 default_r_val=0.2_dp)
137 CALL section_add_keyword(section, keyword)
138 CALL keyword_release(keyword)
139
140 CALL keyword_create(keyword, __location__, name="DELTA_NPOLES", &
141 description="Number of poles of Fermi function to consider.", &
142 n_var=1, type_of_var=integer_t, &
143 default_i_val=4)
144 CALL section_add_keyword(section, keyword)
145 CALL keyword_release(keyword)
146
147 CALL keyword_create(keyword, __location__, name="GAMMA_KT", &
148 description="Offset from the axis (in terms of k*T)"// &
149 " where poles of the Fermi function reside.", &
150 n_var=1, type_of_var=integer_t, &
151 default_i_val=20)
152 CALL section_add_keyword(section, keyword)
153 CALL keyword_release(keyword)
154
155 CALL keyword_create(keyword, __location__, name="INTEGRATION_METHOD", &
156 description="Method to integrate Green's functions along a closed-circuit contour.", &
157 default_i_val=negfint_method_cc, &
158 enum_c_vals=s2a("CLENSHAW-CURTIS", "SIMPSON"), &
159 enum_desc=s2a( &
160 "Adaptive Clenshaw-Curtis quadrature method. Requires FFTW3 library.", &
161 "Adaptive Simpson method. Works without FFTW3."), &
162 enum_i_vals=[negfint_method_cc, negfint_method_simpson])
163 CALL section_add_keyword(section, keyword)
164 CALL keyword_release(keyword)
165
166 CALL keyword_create(keyword, __location__, name="INTEGRATION_MIN_POINTS", &
167 description="Initial (minimal) number of grid point for adaptive numerical integration.", &
168 n_var=1, type_of_var=integer_t, &
169 default_i_val=16)
170 CALL section_add_keyword(section, keyword)
171 CALL keyword_release(keyword)
172
173 CALL keyword_create(keyword, __location__, name="INTEGRATION_MAX_POINTS", &
174 description="Maximal number of grid point for adaptive numerical integration.", &
175 n_var=1, type_of_var=integer_t, &
176 default_i_val=512)
177 CALL section_add_keyword(section, keyword)
178 CALL keyword_release(keyword)
179
180 CALL keyword_create(keyword, __location__, name="MAX_SCF", &
181 description="Maximum number of SCF iterations to be performed.", &
182 n_var=1, type_of_var=integer_t, &
183 default_i_val=30)
184 CALL section_add_keyword(section, keyword)
185 CALL keyword_release(keyword)
186
187 CALL keyword_create(keyword, __location__, name="NPROC_POINT", &
188 description="Number of MPI processes to be used per energy point."// &
189 " Default is to use all processors (0).", &
190 n_var=1, type_of_var=integer_t, &
191 default_i_val=0)
192 CALL section_add_keyword(section, keyword)
193 CALL keyword_release(keyword)
194
195 CALL keyword_create(keyword, __location__, name="V_SHIFT", &
196 description="Initial value of the Hartree potential shift", &
197 n_var=1, type_of_var=real_t, unit_str="hartree", &
198 default_r_val=0.0_dp)
199 CALL section_add_keyword(section, keyword)
200 CALL keyword_release(keyword)
201
202 CALL keyword_create(keyword, __location__, name="V_SHIFT_OFFSET", &
203 description="Initial offset to determine the optimal shift in Hartree potential.", &
204 n_var=1, type_of_var=real_t, default_r_val=0.10_dp)
205 CALL section_add_keyword(section, keyword)
206 CALL keyword_release(keyword)
207
208 CALL keyword_create(keyword, __location__, name="V_SHIFT_MAX_ITERS", &
209 description="Maximal number of iteration to determine the optimal shift in Hartree potential.", &
210 n_var=1, type_of_var=integer_t, default_i_val=30)
211 CALL section_add_keyword(section, keyword)
212 CALL keyword_release(keyword)
213
214 ! PRINT subsection
215 CALL section_create(subsection, __location__, "PRINT", "Printing of information during the NEGF.", &
216 repeats=.false.)
217
218 CALL create_print_program_run_info_section(print_key)
219 CALL section_add_subsection(subsection, print_key)
220 CALL section_release(print_key)
221
222 CALL create_print_dos_section(print_key, "DOS", "the Density of States (DOS) in the scattering region")
223 CALL section_add_subsection(subsection, print_key)
224 CALL section_release(print_key)
225
226 CALL create_print_dos_section(print_key, "TRANSMISSION", "the transmission function")
227 CALL section_add_subsection(subsection, print_key)
228 CALL section_release(print_key)
229
230 CALL keyword_create(keyword, __location__, name="ENERGY_UNIT", &
231 description="Energy unit used for the printed Transmission/DOS energy axis. "// &
232 "Density of states (DOS) values are converted consistently to the selected energy unit.", &
233 usage="ENERGY_UNIT HARTREE", type_of_var=enum_t, &
234 enum_c_vals=s2a("HARTREE", "EV"), &
235 enum_i_vals=[1, 2], &
236 enum_desc=s2a("Print energies in Hartree (a.u.).", &
237 "Print energies in electronvolt (eV)."), &
238 default_i_val=2)
239 CALL section_add_keyword(subsection, keyword)
240 CALL keyword_release(keyword)
241
242 ! RESTART subsection
243 CALL cp_print_key_section_create(print_key, __location__, "RESTART", &
244 description="Controls the dumping of the restart files during NEGF SCF.", &
245 print_level=low_print_level, common_iter_levels=3, &
246 each_iter_names=s2a("NEGF_SCF"), each_iter_values=[1], &
247 add_last=add_last_numeric, filename="RESTART", do_backup=.true.)
248 CALL section_add_subsection(subsection, print_key)
249 CALL section_release(print_key)
250
251 CALL section_add_subsection(section, subsection)
252 CALL section_release(subsection)
253
254 ! SCF subsection
255 CALL section_create(subsection, __location__, "SCF", "Additional parameters for an NEGF SCF run.", &
256 repeats=.false.)
257
258 CALL keyword_create(keyword, __location__, name="UPDATE_HS", &
259 description="Change the atomic Hamiltonian during NEGF self-consistent cycle."// &
260 " If .FALSE., only the density matrix is updated.", &
261 default_l_val=.true., lone_keyword_l_val=.true.)
262 CALL section_add_keyword(subsection, keyword)
263 CALL keyword_release(keyword)
264
265 CALL keyword_create(keyword, __location__, name="RESTART_SCF", &
266 description="Restart SCF from files saved as defined in NEGF%PRINT%RESTART", &
267 default_l_val=.true., lone_keyword_l_val=.true.)
268 CALL section_add_keyword(subsection, keyword)
269 CALL keyword_release(keyword)
270
271 CALL section_add_subsection(section, subsection)
272 CALL section_release(subsection)
273
274 END SUBROUTINE create_negf_section
275
276! **************************************************************************************************
277!> \brief Create NEGF%CONTACT input section.
278!> \param section input section
279!> \par History
280!> * 09.2017 split from create_negf_section() [Sergey Chulkov]
281! **************************************************************************************************
282 SUBROUTINE create_contact_section(section)
283 TYPE(section_type), POINTER :: section
284
285 TYPE(keyword_type), POINTER :: keyword
286 TYPE(section_type), POINTER :: print_key, subsection, subsection2
287
288 cpassert(.NOT. ASSOCIATED(section))
289
290 CALL section_create(section, __location__, name="CONTACT", &
291 description="Section defining the contact region of NEGF setup.", &
292 n_keywords=5, n_subsections=3, repeats=.true.)
293
294 NULLIFY (keyword, print_key, subsection, subsection2)
295
296 CALL create_atomlist_section(subsection, "BULK_REGION", &
297 "the bulk contact adjacent to the screening region.", .false.)
298 CALL section_add_subsection(section, subsection)
299 CALL create_atomlist_section(subsection2, "CELL", &
300 "a single bulk contact unit cell. Bulk Hamiltonian will be contstructed "// &
301 "using two such unit cells instead of performing k-point bulk calculation. "// &
302 "FORCE_EVAL_SECTION must be 0.", .true.)
303 CALL section_add_subsection(subsection, subsection2)
304 CALL section_release(subsection2)
305 CALL section_release(subsection)
306
307 CALL create_atomlist_section(subsection, "SCREENING_REGION", &
308 "the given contact adjacent to the scattering region.", .false.)
309 CALL section_add_subsection(section, subsection)
310 CALL section_release(subsection)
311
312 CALL keyword_create(keyword, __location__, name="FORCE_EVAL_SECTION", &
313 description=" Index of the FORCE_EVAL section which will be used for bulk calculation.", &
314 n_var=1, type_of_var=integer_t, default_i_val=0)
315 CALL section_add_keyword(section, keyword)
316 CALL keyword_release(keyword)
317
318 CALL keyword_create(keyword, __location__, name="ELECTRIC_POTENTIAL", &
319 description="External electrostatic potential applied to the given contact.", &
320 n_var=1, type_of_var=real_t, unit_str="hartree", &
321 default_r_val=0.0_dp)
322 CALL section_add_keyword(section, keyword)
323 CALL keyword_release(keyword)
324
325 CALL keyword_create(keyword, __location__, name="FERMI_LEVEL", &
326 description="Contact Fermi level at the given temperature."// &
327 " If this keyword is not given explicitly, the Fermi level"// &
328 " will be automatically computed prior the actual NEGF calculation.", &
329 n_var=1, type_of_var=real_t, unit_str="hartree", &
330 default_r_val=0.0_dp)
331 CALL section_add_keyword(section, keyword)
332 CALL keyword_release(keyword)
333
334 CALL keyword_create(keyword, __location__, name="REFINE_FERMI_LEVEL", &
335 description="Compute the Fermi level using the value from the FERMI_LEVEL keyword"// &
336 " as a starting point. By default the Fermi level is computed only"// &
337 " when the keyword FERMI_LEVEL is not given explicitly.", &
338 default_l_val=.false., lone_keyword_l_val=.true.)
339 CALL section_add_keyword(section, keyword)
340 CALL keyword_release(keyword)
341
342 CALL keyword_create(keyword, __location__, name="FERMI_LEVEL_SHIFTED", &
343 description="Used to shift the zero-energy level of an electrode to the common zero-energy level."// &
344 " If this keyword is specified, the Fermi level,"// &
345 " calculated by standard DFT or NEGF (using the REFINE_FERMI_LEVEL keyword),"// &
346 " or previously specified using the FERMI_LEVEL keyword,"// &
347 " is changed to this value. All diagonal elements of the Hamiltonian are shifted accordingly.", &
348 n_var=1, type_of_var=real_t, unit_str="hartree", &
349 default_r_val=0.0_dp)
350 CALL section_add_keyword(section, keyword)
351 CALL keyword_release(keyword)
352
353 CALL keyword_create(keyword, __location__, name="TEMPERATURE", &
354 description="Electronic temperature.", &
355 n_var=1, type_of_var=real_t, unit_str="K", &
356 default_r_val=300.0_dp/kelvin)
357 CALL section_add_keyword(section, keyword)
358 CALL keyword_release(keyword)
359
360 ! PRINT subsection
361 CALL section_create(subsection, __location__, "PRINT", "Print properties for the given contact.", &
362 repeats=.false.)
363
364 CALL create_print_dos_section(print_key, "DOS", "the Density of States (DOS)")
365 CALL section_add_subsection(subsection, print_key)
366 CALL section_release(print_key)
367
368 CALL section_add_subsection(section, subsection)
369 CALL section_release(subsection)
370
371 ! RESTART subsection
372 CALL section_create(subsection, __location__, "RESTART", &
373 "Controls reading and writing of the restart files for the specified contact.", repeats=.false.)
374
375 CALL keyword_create(keyword, __location__, name="FILENAME", &
376 description=' Controls part of the filename for output. '// &
377 ' Use filename to obtain projectname-filename. '// &
378 ' Use ./filename to get filename.'// &
379 ' A middle name (contact number, type of matrix, spin if 2 spins) and extension '// &
380 ' are always added to the filename.', &
381 usage="FILENAME ./filename ", &
382 default_lc_val="")
383 CALL section_add_keyword(subsection, keyword)
384 CALL keyword_release(keyword)
385
386 CALL keyword_create(keyword, __location__, name="READ_WRITE_HS", &
387 description="Requests reading of the electrode Hamiltonian and overlap matrices from a file."// &
388 " If at least one of these files doesn't exist, all Hamiltonian and overlap matrices"// &
389 " are calculated and saved. If no name is specified by FILENAME, the default file names are"// &
390 " projectname-Nn-H00, projectname-Nn-H01, projectname-Nn-S00, projectname-Nn-S01"// &
391 " for restricted calculations with identical spin population, or"// &
392 " projectname-Nn-H00-Ss, projectname-Nn-H01-Ss, projectname-Nn-S00-Ss, projectname-Nn-S01-Ss"// &
393 " for unrestricted calculations with two spin components,"// &
394 " where n is the number of the contact and s is the spin index."// &
395 " Otherwise, projectname is modified. Note that the code does not distinguish"// &
396 " between the files created by a separate electrode or the entire system.", &
397 default_l_val=.false., lone_keyword_l_val=.true.)
398 CALL section_add_keyword(subsection, keyword)
399 CALL keyword_release(keyword)
400
401 CALL section_add_subsection(section, subsection)
402 CALL section_release(subsection)
403
404 END SUBROUTINE create_contact_section
405
406! **************************************************************************************************
407!> \brief Create a scattering region section.
408!> \param section NEGF section
409!> \par History
410!> * 02.2026 created [Dmitry Ryndyk]
411! **************************************************************************************************
412 SUBROUTINE create_scattering_region_section(section)
413 TYPE(section_type), POINTER :: section
414
415 TYPE(keyword_type), POINTER :: keyword
416 TYPE(section_type), POINTER :: subsection
417
418 cpassert(.NOT. ASSOCIATED(section))
419
420 CALL section_create(section, __location__, name="SCATTERING_REGION", &
421 description="Section defining the scattering region of NEGF setup.", &
422 n_keywords=2, n_subsections=1, repeats=.false.)
423
424 NULLIFY (keyword, subsection)
425
426 CALL keyword_create(keyword, __location__, name="LIST", &
427 description="Specifies a list of atoms.", &
428 usage="LIST {integer} {integer} .. {integer}", repeats=.true., &
429 n_var=-1, type_of_var=integer_t)
430 CALL section_add_keyword(section, keyword)
431 CALL keyword_release(keyword)
432
433 CALL keyword_create(keyword, __location__, name="MOLNAME", &
434 description="Specifies a list of named molecular fragments.", &
435 usage="MOLNAME WAT MEOH", repeats=.true., &
436 n_var=-1, type_of_var=char_t)
437 CALL section_add_keyword(section, keyword)
438 CALL keyword_release(keyword)
439
440 ! RESTART subsection
441 CALL section_create(subsection, __location__, "RESTART", &
442 "Controls reading and writing of the initial restart files for the scattering region,"// &
443 " including the coupling to the electrodes.", repeats=.false.)
444 CALL keyword_create(keyword, __location__, name="FILENAME", &
445 description=' Controls part of the filename. '// &
446 ' Use filename to obtain projectname-filename. '// &
447 ' Use ./filename to get filename.'// &
448 ' A middle name (type of matrix, contact number for coupling matrices, spin if 2 spins)'// &
449 ' and extension are always added to the filename.', &
450 usage="FILENAME ./filename ", &
451 default_lc_val="")
452 CALL section_add_keyword(subsection, keyword)
453 CALL keyword_release(keyword)
454 CALL keyword_create(keyword, __location__, name="READ_WRITE_HS", &
455 description="Requests reading of the Hamiltonian and overlap matrices from a file."// &
456 " For the scattering region, including the coupling to the electrodes."// &
457 " If at least one of these files doesn't exist, all Hamiltonian and overlap matrices"// &
458 " are calculated and saved. If no name is specified by FILENAME, the default file names are"// &
459 " projectname-H, projectname-H-Nn, projectname-S, projectname-S-Nn"// &
460 " for restricted calculations with identical spin population, or"// &
461 " projectname-H-Ss, projectname-H-Nn-Ss, projectname-S-Ss, projectname-S-Nn-Ss"// &
462 " for unrestricted calculations with two spin components,"// &
463 " where n is the number of the contact and s is the spin index."// &
464 " Otherwise, projectname is modified.", &
465 default_l_val=.false., lone_keyword_l_val=.true.)
466 CALL section_add_keyword(subsection, keyword)
467 CALL keyword_release(keyword)
468 CALL section_add_subsection(section, subsection)
469 CALL section_release(subsection)
470
471 END SUBROUTINE create_scattering_region_section
472
473! **************************************************************************************************
474!> \brief Create an atomic list section.
475!> \param section NEGF section
476!> \param name name of the new section
477!> \param description section description
478!> \param repeats whether the section can be repeated
479!> \par History
480!> * 02.2017 created [Sergey Chulkov]
481! **************************************************************************************************
482 SUBROUTINE create_atomlist_section(section, name, description, repeats)
483 TYPE(section_type), POINTER :: section
484 CHARACTER(len=*), INTENT(in) :: name, description
485 LOGICAL, INTENT(in) :: repeats
486
487 TYPE(keyword_type), POINTER :: keyword
488
489 cpassert(.NOT. ASSOCIATED(section))
490
491 CALL section_create(section, __location__, name=trim(adjustl(name)), &
492 description="Atoms belonging to "//trim(adjustl(description)), &
493 n_keywords=2, n_subsections=0, repeats=repeats)
494
495 NULLIFY (keyword)
496
497 CALL keyword_create(keyword, __location__, name="LIST", &
498 description="Specifies a list of atoms.", &
499 usage="LIST {integer} {integer} .. {integer}", repeats=.true., &
500 n_var=-1, type_of_var=integer_t)
501 CALL section_add_keyword(section, keyword)
502 CALL keyword_release(keyword)
503
504 CALL keyword_create(keyword, __location__, name="MOLNAME", &
505 description="Specifies a list of named molecular fragments.", &
506 usage="MOLNAME WAT MEOH", repeats=.true., &
507 n_var=-1, type_of_var=char_t)
508 CALL section_add_keyword(section, keyword)
509 CALL keyword_release(keyword)
510 END SUBROUTINE create_atomlist_section
511
512! **************************************************************************************************
513!> \brief Create the PROGRAM_RUN_INFO print section.
514!> \param section section to create
515!> \par History
516!> * 11.2020 created [Dmitry Ryndyk]
517! **************************************************************************************************
518 SUBROUTINE create_print_program_run_info_section(section)
519
520 TYPE(section_type), POINTER :: section
521
522 TYPE(keyword_type), POINTER :: keyword
523
524 CALL cp_print_key_section_create(section, __location__, "PROGRAM_RUN_INFO", &
525 description="Controls the printing of basic information during the NEGF.", &
526 print_level=low_print_level, filename="__STD_OUT__")
527 NULLIFY (keyword)
528
529 CALL keyword_create(keyword, __location__, name="_SECTION_PARAMETERS_", &
530 description="Level starting at which this property is printed", &
531 usage="_SECTION_PARAMETERS_", &
532 default_i_val=low_print_level, lone_keyword_i_val=low_print_level, &
533 enum_c_vals=s2a("on", "off", "silent", "low", "medium", "high", "debug"), &
534 enum_i_vals=[silent_print_level - 1, debug_print_level + 1, &
535 silent_print_level, low_print_level, &
536 medium_print_level, high_print_level, debug_print_level])
537 CALL section_add_keyword(section, keyword)
538 CALL keyword_release(keyword)
539
540 CALL keyword_create(keyword, __location__, name="PRINT_LEVEL", &
541 variants=["IOLEVEL"], &
542 description="Determines the verbose level for this section "// &
543 "additionally to GLOBAL%PRINT_LEVEL and SECTION_PARAMETERS, "// &
544 "which switch on printing.", &
545 usage="PRINT_LEVEL HIGH", &
546 default_i_val=low_print_level, enum_c_vals= &
547 s2a("SILENT", "LOW", "MEDIUM", "HIGH", "DEBUG"), &
548 enum_desc=s2a("No output", &
549 "Little output", "Quite some output", "Lots of output", &
550 "Everything is written out, useful for debugging purposes only"), &
551 enum_i_vals=[silent_print_level, low_print_level, medium_print_level, &
552 high_print_level, debug_print_level])
553 CALL section_add_keyword(section, keyword)
554 CALL keyword_release(keyword)
555
556 END SUBROUTINE create_print_program_run_info_section
557
558! **************************************************************************************************
559!> \brief Create the DOS/TRANSMISSION print section.
560!> \param section section to create
561!> \param name name of the new section
562!> \param description section description
563!> \par History
564!> * 11.2017 created [Sergey Chulkov]
565! **************************************************************************************************
566 SUBROUTINE create_print_dos_section(section, name, description)
567 TYPE(section_type), POINTER :: section
568 CHARACTER(len=*), INTENT(in) :: name, description
569
570 TYPE(keyword_type), POINTER :: keyword
571
572 CALL cp_print_key_section_create(section, __location__, trim(adjustl(name)), &
573 description="Controls the printing of "//trim(adjustl(description))//".", &
574 print_level=high_print_level, filename="__STD_OUT__")
575 NULLIFY (keyword)
576
577 CALL keyword_create(keyword, __location__, name="FROM_ENERGY", &
578 description="Energy point to start with.", &
579 n_var=1, type_of_var=real_t, unit_str="hartree", &
580 default_r_val=-1.0_dp)
581 CALL section_add_keyword(section, keyword)
582 CALL keyword_release(keyword)
583
584 CALL keyword_create(keyword, __location__, name="TILL_ENERGY", &
585 description="Energy point to end with.", &
586 n_var=1, type_of_var=real_t, unit_str="hartree", &
587 default_r_val=1.0_dp)
588 CALL section_add_keyword(section, keyword)
589 CALL keyword_release(keyword)
590
591 CALL keyword_create(keyword, __location__, name="N_GRIDPOINTS", &
592 description="Number of points to compute.", &
593 n_var=1, type_of_var=integer_t, default_i_val=201)
594 CALL section_add_keyword(section, keyword)
595 CALL keyword_release(keyword)
596 END SUBROUTINE create_print_dos_section
597END MODULE input_cp2k_negf
string_utilities::s2a
Definition string_utilities.F:49
bibliography
collects all references to literature in CP2K as new algorithms / method are included from literature...
Definition bibliography.F:21
bibliography::papior2017
integer, save, public papior2017
Definition bibliography.F:36
bibliography::bailey2006
integer, save, public bailey2006
Definition bibliography.F:36
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_section_create
subroutine, public cp_print_key_section_create(print_key_section, location, name, description, print_level, each_iter_names, each_iter_values, add_last, filename, common_iter_levels, citations, unit_str, do_backup)
creates a print_key section
Definition cp_output_handling.F:146
cp_output_handling::debug_print_level
integer, parameter, public debug_print_level
Definition cp_output_handling.F:77
cp_output_handling::low_print_level
integer, parameter, public low_print_level
Definition cp_output_handling.F:77
cp_output_handling::medium_print_level
integer, parameter, public medium_print_level
Definition cp_output_handling.F:77
cp_output_handling::high_print_level
integer, parameter, public high_print_level
Definition cp_output_handling.F:77
cp_output_handling::add_last_numeric
integer, parameter, public add_last_numeric
Definition cp_output_handling.F:74
cp_output_handling::silent_print_level
integer, parameter, public silent_print_level
Definition cp_output_handling.F:77
input_constants
collects all constants needed in input so that they can be used without circular dependencies
Definition input_constants.F:18
input_constants::negfint_method_simpson
integer, parameter, public negfint_method_simpson
Definition input_constants.F:1360
input_constants::negfint_method_cc
integer, parameter, public negfint_method_cc
Definition input_constants.F:1360
input_cp2k_negf
Input section for NEGF based quantum transport calculations.
Definition input_cp2k_negf.F:12
input_cp2k_negf::create_negf_section
subroutine, public create_negf_section(section)
Create NEGF input section.
Definition input_cp2k_negf.F:58
input_keyword_types
represents keywords in an input
Definition input_keyword_types.F:14
input_keyword_types::keyword_release
subroutine, public keyword_release(keyword)
releases the given keyword (see doc/ReferenceCounting.html)
Definition input_keyword_types.F:448
input_keyword_types::keyword_create
subroutine, public keyword_create(keyword, location, name, description, usage, type_of_var, n_var, repeats, variants, default_val, default_l_val, default_r_val, default_lc_val, default_c_val, default_i_val, default_l_vals, default_r_vals, default_c_vals, default_i_vals, lone_keyword_val, lone_keyword_l_val, lone_keyword_r_val, lone_keyword_c_val, lone_keyword_i_val, lone_keyword_l_vals, lone_keyword_r_vals, lone_keyword_c_vals, lone_keyword_i_vals, enum_c_vals, enum_i_vals, enum, enum_strict, enum_desc, unit_str, citations, deprecation_notice, removed)
creates a keyword object
Definition input_keyword_types.F:156
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_create
subroutine, public section_create(section, location, name, description, n_keywords, n_subsections, repeats, citations, deprecation_notice)
creates a list of keywords
Definition input_section_types.F:158
input_section_types::section_add_keyword
subroutine, public section_add_keyword(section, keyword)
adds a keyword to the given section
Definition input_section_types.F:514
input_section_types::section_add_subsection
subroutine, public section_add_subsection(section, subsection)
adds a subsection to the given section
Definition input_section_types.F:570
input_section_types::section_release
recursive subroutine, public section_release(section)
releases the given keyword list (see doc/ReferenceCounting.html)
Definition input_section_types.F:238
input_val_types
a wrapper for basic fortran types.
Definition input_val_types.F:14
input_val_types::real_t
integer, parameter, public real_t
Definition input_val_types.F:41
input_val_types::char_t
integer, parameter, public char_t
Definition input_val_types.F:41
input_val_types::integer_t
integer, parameter, public integer_t
Definition input_val_types.F:41
input_val_types::enum_t
integer, parameter, public enum_t
Definition input_val_types.F:41
kinds
Defines the basic variable types.
Definition kinds.F:23
kinds::dp
integer, parameter, public dp
Definition kinds.F:34
physcon
Definition of physical constants:
Definition physcon.F:68
physcon::kelvin
real(kind=dp), parameter, public kelvin
Definition physcon.F:165
qs_density_mixing_types
module that contains the definitions of the scf types
Definition qs_density_mixing_types.F:14
qs_density_mixing_types::create_mixing_section
subroutine, public create_mixing_section(section, ls_scf, active_space)
Create CP2K input section for SCF or active-space density mixing.
Definition qs_density_mixing_types.F:492
string_utilities
Utilities for string manipulations.
Definition string_utilities.F:16
input_keyword_types::keyword_type
represent a keyword in the input
Definition input_keyword_types.F:89
input_section_types::section_type
represent a section of the input file
Definition input_section_types.F:102