Line data Source code
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 function that build the input sections for external [potential, density VXC]
10 : !> \par History
11 : !> 10.2005 moved out of input_cp2k [fawzi]
12 : !> 10.2020 moved out of input_cp2k_dft [JGH]
13 : !> \author fawzi
14 : ! **************************************************************************************************
15 : MODULE input_cp2k_external
16 : USE bibliography, ONLY: Tozer1996,&
17 : Zhao1994
18 : USE input_constants, ONLY: use_coulomb,&
19 : use_diff,&
20 : use_no
21 : USE input_keyword_types, ONLY: keyword_create,&
22 : keyword_release,&
23 : keyword_type
24 : USE input_section_types, ONLY: section_add_keyword,&
25 : section_add_subsection,&
26 : section_create,&
27 : section_release,&
28 : section_type
29 : USE input_val_types, ONLY: char_t,&
30 : lchar_t,&
31 : real_t
32 : USE kinds, ONLY: dp
33 : USE string_utilities, ONLY: s2a
34 : #include "./base/base_uses.f90"
35 :
36 : IMPLICIT NONE
37 : PRIVATE
38 :
39 : CHARACTER(len=*), PARAMETER, PRIVATE :: moduleN = 'input_cp2k_external'
40 :
41 : PUBLIC :: create_ext_pot_section, create_ext_den_section, create_ext_vxc_section
42 :
43 : CONTAINS
44 :
45 : ! **************************************************************************************************
46 : !> \brief Creates the section for applying an electrostatic external potential
47 : !> \param section ...
48 : !> \date 12.2009
49 : !> \author teo
50 : ! **************************************************************************************************
51 10643 : SUBROUTINE create_ext_pot_section(section)
52 : TYPE(section_type), POINTER :: section
53 :
54 : TYPE(keyword_type), POINTER :: keyword
55 : TYPE(section_type), POINTER :: subsection
56 :
57 10643 : CPASSERT(.NOT. ASSOCIATED(section))
58 : CALL section_create(section, __LOCATION__, name="EXTERNAL_POTENTIAL", &
59 : description="Section controlling the presence of an electrostatic "// &
60 : "external potential dependent on the atomic positions (X,Y,Z). "// &
61 : "As the external potential is currently applied via a grid, "// &
62 : "it only works with DFT based methods (GPW/GAPW) that already use "// &
63 : "a grid based approach to solve the Poisson equation.", &
64 10643 : n_keywords=7, n_subsections=0, repeats=.FALSE.)
65 10643 : NULLIFY (keyword, subsection)
66 :
67 : CALL keyword_create(keyword, __LOCATION__, name="FUNCTION", &
68 : description="Specifies the functional form in mathematical notation. Variables must be the atomic "// &
69 : "coordinates (X,Y,Z) of the grid.", usage="FUNCTION X^2+Y^2+Z^2+LOG(ABS(X+Y))", &
70 10643 : type_of_var=lchar_t, n_var=1)
71 10643 : CALL section_add_keyword(section, keyword)
72 10643 : CALL keyword_release(keyword)
73 :
74 : CALL keyword_create(keyword, __LOCATION__, name="PARAMETERS", &
75 : description="Defines the parameters of the functional form", &
76 : usage="PARAMETERS a b D", type_of_var=char_t, &
77 10643 : n_var=-1, repeats=.TRUE.)
78 10643 : CALL section_add_keyword(section, keyword)
79 10643 : CALL keyword_release(keyword)
80 :
81 : CALL keyword_create(keyword, __LOCATION__, name="VALUES", &
82 : description="Defines the values of parameter of the functional form", &
83 : usage="VALUES ", type_of_var=real_t, &
84 10643 : n_var=-1, repeats=.TRUE., unit_str="internal_cp2k")
85 10643 : CALL section_add_keyword(section, keyword)
86 10643 : CALL keyword_release(keyword)
87 :
88 : CALL keyword_create(keyword, __LOCATION__, name="UNITS", &
89 : description="Optionally, allows to define valid CP2K unit strings for each parameter value. "// &
90 : "It is assumed that the corresponding parameter value is specified in this unit.", &
91 : usage="UNITS angstrom eV*angstrom^-1 angstrom^1 K", type_of_var=char_t, &
92 10643 : n_var=-1, repeats=.TRUE.)
93 10643 : CALL section_add_keyword(section, keyword)
94 10643 : CALL keyword_release(keyword)
95 :
96 : CALL keyword_create(keyword, __LOCATION__, name="STATIC", &
97 : description="Specifies the external potential as STATIC or time dependent. At the moment "// &
98 : "only static potentials are implemented.", &
99 10643 : usage="STATIC T", default_l_val=.TRUE., lone_keyword_l_val=.TRUE.)
100 10643 : CALL section_add_keyword(section, keyword)
101 10643 : CALL keyword_release(keyword)
102 :
103 : CALL keyword_create(keyword, __LOCATION__, name="DX", &
104 : description="Parameter used for computing the derivative with the Ridders' method.", &
105 10643 : usage="DX <REAL>", default_r_val=0.1_dp, unit_str="bohr")
106 10643 : CALL section_add_keyword(section, keyword)
107 10643 : CALL keyword_release(keyword)
108 :
109 : CALL keyword_create(keyword, __LOCATION__, name="ERROR_LIMIT", &
110 : description="Checks that the error in computing the derivative is not larger than "// &
111 : "the value set; in case error is larger a warning message is printed.", &
112 10643 : usage="ERROR_LIMIT <REAL>", default_r_val=1.0E-12_dp)
113 10643 : CALL section_add_keyword(section, keyword)
114 10643 : CALL keyword_release(keyword)
115 :
116 : !keyword for reading the external potential from cube file
117 : CALL keyword_create(keyword, __LOCATION__, name="READ_FROM_CUBE", &
118 : description="Switch for reading the external potential from file pot.cube. The values "// &
119 : "of the potential must be on the grid points of the realspace grid. "// &
120 : "Whitespace-separated cube values are accepted. If adjacent values are written "// &
121 : "without whitespace, each value must occupy a 13-character E13.5 field.", &
122 10643 : usage="READ_FROM_CUBE T", default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
123 10643 : CALL section_add_keyword(section, keyword)
124 10643 : CALL keyword_release(keyword)
125 :
126 : !keyword for scaling the external potential that is read from file by a constant factor
127 : CALL keyword_create(keyword, __LOCATION__, name="SCALING_FACTOR", &
128 : description="A factor for scaling the the external potential that is read from file. "// &
129 : "The value of the potential at each grid point is multiplied by this factor.", &
130 10643 : usage="SCALING_FACTOR <REAL>", default_r_val=1.0_dp)
131 10643 : CALL section_add_keyword(section, keyword)
132 10643 : CALL keyword_release(keyword)
133 :
134 10643 : CALL create_maxwell_section(subsection)
135 10643 : CALL section_add_subsection(section, subsection)
136 10643 : CALL section_release(subsection)
137 :
138 10643 : END SUBROUTINE create_ext_pot_section
139 :
140 : ! **************************************************************************************************
141 : !> \brief Creates the section for applying an electrostatic external potential
142 : !> \param section ...
143 : !> \date 12.2009
144 : !> \author teo
145 : ! **************************************************************************************************
146 10643 : SUBROUTINE create_maxwell_section(section)
147 : TYPE(section_type), POINTER :: section
148 :
149 : TYPE(keyword_type), POINTER :: keyword
150 :
151 10643 : CPASSERT(.NOT. ASSOCIATED(section))
152 : CALL section_create(section, __LOCATION__, name="MAXWELL", &
153 : description="Section controlling the calculation of an electrostatic "// &
154 : "external potential calculated from Maxwell equations. ", &
155 10643 : n_keywords=1, n_subsections=0, repeats=.FALSE.)
156 10643 : NULLIFY (keyword)
157 :
158 : CALL keyword_create(keyword, __LOCATION__, name="TEST_LOGICAL", &
159 : description="Test for logical value", &
160 10643 : usage="TEST_LOGICAL T", default_l_val=.FALSE., lone_keyword_l_val=.TRUE.)
161 10643 : CALL section_add_keyword(section, keyword)
162 10643 : CALL keyword_release(keyword)
163 :
164 : CALL keyword_create(keyword, __LOCATION__, name="TEST_REAL", &
165 : description="TEST for Real", &
166 10643 : usage="TEST_REAL <REAL>", default_r_val=1.0_dp)
167 10643 : CALL section_add_keyword(section, keyword)
168 10643 : CALL keyword_release(keyword)
169 :
170 : CALL keyword_create(keyword, __LOCATION__, name="TEST_INTEGER", &
171 : description="TEST for Int", &
172 10643 : usage="TEST_INTEGER <INT>", default_i_val=0)
173 10643 : CALL section_add_keyword(section, keyword)
174 10643 : CALL keyword_release(keyword)
175 :
176 10643 : END SUBROUTINE create_maxwell_section
177 :
178 : ! **************************************************************************************************
179 : !> \brief ZMP Creates the section for reading user supplied external density
180 : !> \param section ...
181 : !> \date 03.2011
182 : !> \author D. Varsano [daniele.varsano@nano.cnr.it]
183 : ! **************************************************************************************************
184 10643 : SUBROUTINE create_ext_den_section(section)
185 : TYPE(section_type), POINTER :: section
186 :
187 : TYPE(keyword_type), POINTER :: keyword
188 :
189 10643 : CPASSERT(.NOT. ASSOCIATED(section))
190 : CALL section_create(section, __LOCATION__, name="EXTERNAL_DENSITY", &
191 : description="Section for the use of the ZMP technique on external densities.", &
192 : n_keywords=4, n_subsections=0, repeats=.FALSE., &
193 31929 : citations=[Zhao1994, Tozer1996])
194 10643 : NULLIFY (keyword)
195 :
196 : CALL keyword_create(keyword, __LOCATION__, name="FILE_DENSITY", &
197 : description="Specifies the filename containing the target density in *.cube format. "// &
198 : "In the MGRID section it must be imposed NGRID 1, as it works with only "// &
199 : "one grid. The number of points in each direction, and the spacing must "// &
200 : "be previously defined choosing the plane waves cut-off in section MGRID "// &
201 : "keyword CUTOFF, and the cube dimension in section SUBSYS / CELL / keyword ABC", &
202 : usage="FILE_DENSITY <FILENAME>", &
203 10643 : type_of_var=char_t, default_c_val="RHO_O.dat", n_var=-1)
204 10643 : CALL section_add_keyword(section, keyword)
205 10643 : CALL keyword_release(keyword)
206 :
207 : CALL keyword_create(keyword, __LOCATION__, name="LAMBDA", &
208 : description="Lagrange multiplier defined in the constraint ZMP method. When starting, use "// &
209 : "small values when starting from scratch (around 5,10). Then gradually increase "// &
210 : "the values depending, restarting from the previous calculation with the smaller "// &
211 : "value. To choose the progressive values of LAMBDA look at the convergence of the "// &
212 : "eigenvalues.", &
213 10643 : usage="LAMBDA <REAL>", default_r_val=10.0_dp)
214 10643 : CALL section_add_keyword(section, keyword)
215 10643 : CALL keyword_release(keyword)
216 :
217 : CALL keyword_create(keyword, __LOCATION__, name="ZMP_CONSTRAINT", &
218 : description="Specify which kind of constraint to solve the ZMP equation. The COULOMB default "// &
219 : "option is more stable.", &
220 : usage="ZMP_CONSTRAINT <CHAR>", &
221 : enum_c_vals=s2a("COULOMB", "DIFF", "NONE"), &
222 : enum_i_vals=[use_coulomb, use_diff, use_no], &
223 : enum_desc=s2a("Coulomb constraint, integral of [rho_0(r)-rho(r)]/|r-r'|", &
224 : "Simple constraint, [rho_0(r)-rho(r)]", &
225 : "No constrain imposed"), &
226 10643 : default_i_val=use_coulomb)
227 10643 : CALL section_add_keyword(section, keyword)
228 10643 : CALL keyword_release(keyword)
229 :
230 : CALL keyword_create(keyword, __LOCATION__, name="FERMI_AMALDI", &
231 : description="Add the Fermi-Amaldi contribution to the Hartree potential. "// &
232 : "It leads to a more stable convergence.", &
233 : usage="FERMI_AMALDI <LOGICAL>", &
234 : repeats=.FALSE., &
235 : n_var=1, &
236 10643 : default_l_val=.TRUE., lone_keyword_l_val=.TRUE.)
237 10643 : CALL section_add_keyword(section, keyword)
238 10643 : CALL keyword_release(keyword)
239 :
240 10643 : END SUBROUTINE create_ext_den_section
241 :
242 : ! **************************************************************************************************
243 : !> \brief ZMP Creates the section for creating the external v_xc
244 : !> \param section ...
245 : !> \date 03.2011
246 : !> \author D. Varsano [daniele.varsano@nano.cnr.it]
247 : ! **************************************************************************************************
248 10643 : SUBROUTINE create_ext_vxc_section(section)
249 : TYPE(section_type), POINTER :: section
250 :
251 : TYPE(keyword_type), POINTER :: keyword
252 :
253 10643 : CPASSERT(.NOT. ASSOCIATED(section))
254 : CALL section_create(section, __LOCATION__, name="EXTERNAL_VXC", &
255 : description="SCF convergence with external v_xc calculated through previous ZMP "// &
256 : "calculation", &
257 10643 : n_keywords=1, n_subsections=0, repeats=.FALSE.)
258 10643 : NULLIFY (keyword)
259 :
260 : CALL keyword_create(keyword, __LOCATION__, name="FILE_VXC", &
261 : description="The *.cube filename containing the v_xc potential. This works only "// &
262 : "with NGRID 1 imposed in the MGRID section. The number of points in each "// &
263 : "direction, and the spacing must equal to those previously used in the ZMP "// &
264 : "calculation and defined through the plane wave cut-off and the cube dimension "// &
265 : "respectively set in section MGRID / keyword CUTOFF, and in section SUBSYS / "// &
266 : "CELL / keyword ABC", &
267 : usage="FILE_VXC <FILENAME>", &
268 10643 : type_of_var=char_t, default_c_val="VXC_O.dat", n_var=-1)
269 10643 : CALL section_add_keyword(section, keyword)
270 10643 : CALL keyword_release(keyword)
271 10643 : END SUBROUTINE create_ext_vxc_section
272 :
273 : END MODULE input_cp2k_external
|
