-
Notifications
You must be signed in to change notification settings - Fork 6
Expand file tree
/
Copy pathREADME.html
More file actions
140 lines (111 loc) · 8.91 KB
/
README.html
File metadata and controls
140 lines (111 loc) · 8.91 KB
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
<!DOCTYPE HTML PUBLIC "-//W3C//DTD HTML 4.01 Transitional//EN" "http://www.w3.org/TR/html4/loose.dtd">
<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=UTF-8">
<title>PyPointChargePotentialSolver</title>
<meta name='Generator' content='Zim 0.65'>
<style type='text/css'>
a { text-decoration: none }
a:hover { text-decoration: underline }
a:active { text-decoration: underline }
strike { color: grey }
u { text-decoration: none;
background-color: yellow }
tt { color: #2e3436; }
pre { color: #2e3436;
margin-left: 20px }
h1 { text-decoration: underline;
color: #4e9a06; margin-bottom: 0 }
h2 { color: #4e9a06; margin-bottom: 0 }
h3 { color: #4e9a06; margin-bottom: 0 }
h4 { color: #4e9a06; margin-bottom: 0 }
h5 { color: #4e9a06; margin-bottom: 0 }
p { margin-top: 0 }
span.zim-tag {
color: #ce5c00;
}
div.zim-object {
border-style:solid;
border-width:1px;
}
.checked-box {list-style-image: url(data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAABgAAAAYCAYAAADgdz34AAAAAXNSR0IArs4c6QAAAAZiS0dEAP8A/wD/oL2nkwAAAAlwSFlzAAANOgAADMQBiN+4gQAAAAd0SU1FB9gKGQ8sMEGsKGkAAAAZdEVYdENvbW1lbnQAQ3JlYXRlZCB3aXRoIEdJTVBXgQ4XAAAEBUlEQVRIx62V22tdRRTGf7Nn73P2ybntnNOe3NqkPTGgLTVUUZF6QatSLOKTPgqCIqLgQ0H/A1sQQbBYCBb1QfAxiC8tSO1FqHkwJVKtjdTGNraUmObsc9nXmfGh7cGYpM1D5nHWzPetteZb3wg2eB2YqYm4zSadsMtoboiNBH/3TE0awx6j+MRoxoTg/IYRvP19TQrJS0bzhdHGSyKFkLTtjSKwMjyiEz43ynhtP6bdjBCWyFobAf7eT7VhNF/q1FRbjYjmUohlCVPwnB+6FUxMTJipqSmUUhhjEGKd3bMT4ks/Y6oLBK2Yth8hHYtCJXOix7Nf7xLMzc0xOzvLzp078TyPNE3viW3QJPXzhNWbxFFKHCmMhoLn/FHodd48vGfhapdAacXQlkFK5dL6wIUm6fuTZPuvqDQhaMUYYyiVyuQr6rXDexYuAdi3tSv1ZJNs/R/CaszzT+1na88uXFnCEnJVgivBNN8uTJKmHQI/ptOOcXNZzMz9mOqFs90OHpipWcYwlo5P4ebnuOkrvr5wgrH+h3im7y36MzuwRXYZeKha/OhP0EkadFoxQSdGSotedR/+XwMc2XvKdNUFOFqZx6LKZWIiwjgkikNmLp/hm8sH+K1zjFTHXfBYdTi+eJArzXM0GxFxoBDCopLvo/fqEwi1XPkWkGqjFo2TgB1jOYZUKZTS/D1/ncmLh7jon0IbRWoiTi59ymzzJEEQE3cStNZsGxqlfPE57MBbOR8fP3hDGalOO9fq2DlBvmZw8xa2IxACGn6TydlD/O6f5OzSV/zif0cYhLQaEXGkKBbz7Ov/AOlXV1cxgBJRI3fuSTrpTawt18kWIZN1CFuaONI0w0WOXfsI43YIggh/KUSlhqxrMz74AkOZcWBm9QkH+Gw8NDLuITi+m0yzhluSyJzBLcpblUhFxywSRAEtPwQjsKVN30CNh0uvYuOubSHLtN3J0TO1j0pmBNuFbFWRK0gyPRZpktL2I5JQkclKakNlnh54g6ocvevUr/Ai2a7wineEkcJupA3S1Wg0nVZM2E6wbEF5U5G9Q++wI7sfR7h3N8HVNstykBfzH+KJEZwiWD0aIwxCgJ0R1Mu7GXOeJSuK93bZtQIle4D9pUNU5DC5jEsu55AvZakM5NicGyEj8uuz8bUCQgj67QfY671P3vEoeC69gy695U1U7NG7XV0pUwBjDJa1/JJlWWxzHuflzQe5FJ/GsgUVuZ2t8lEkTvfc0aNHb72flBhjVicQQqCUuvM3/M+WDVguWBrMVdDXEGZlBVEUrVCU9d9s5+fnaTQa2PZyPxEIhJaI1EEoZwX4ncynp6fXrmB4eJjR0VFarRbNZnP9P9rt9gohqNVq1Ov1ZbF/AZGev3hLJ2/zAAAAAElFTkSuQmCC)}
.xchecked-box {list-style-image: url(data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAABgAAAAYCAYAAADgdz34AAAAAXNSR0IArs4c6QAAAAZiS0dEAP8A/wD/oL2nkwAAAAlwSFlzAAANOgAADMQBiN+4gQAAAAd0SU1FB9gKGQ8bDYnDxEwAAAAZdEVYdENvbW1lbnQAQ3JlYXRlZCB3aXRoIEdJTVBXgQ4XAAAEK0lEQVRIx9WVS2hTWRjHf/eR3CY1nbxMH2YiZRQS6qO13YlMVxY3SnVcuNIBFezGpSADLoQqLu1sHJCqdCFSXFpw4YOCSH3BtFqttTNamabX3DS5bfO6uffMoglja3RGcDMHzuac7/z/53++//cd+L8P6VuCPQYZ8ADNgBd4J31DcDcQs+GnHByRocEDv0kfBSjAOlYCs11Q+gpwDegS8LMJ+3QIK0ATzEhV8Odnz5bzw8P4dJ25aJQ/WlvJ1df/K7hSLtOcTNI+Pk69rpMTgqIQhCDvh1/VSpw79+gRrRMTmLZNezLJJsPg+a5dmOEwQlFqg1sWG16/Jv7sGWXDwBACFQjC9HcwIMONKkGp4PGAJGEDS0IQmZlhnWnye3c3eiyGo6qr3WHbrJ+dJf7gAXI6zSIr72T7/fgzmT4FHnTBsgrQBfYvTU0km5vxz86iADnAm0rRPTWFt7cXZccOJJcLAGHb2K9ekT93jmwmwwdAAFpjI6Ntbfxw5879ag7l6o1sr5eHHR3IsRgeQK/M4sQE+YEB7JcvEY6zAj45SWFgAPPxYwzHoQxIkQjTPT0kIxE+Noj8sexFn4/xnh58iQTeSpHkHAdrbIzi0BCOrmNPTpK/eJHM3bt8sCyKQCiR4NWePWSiUZw1+ZLXJm4pFKLhzBlCsRh2RUXacVgeGaF47RrL58+zcP8+RrmMkCQinZ1EL1zAjERqmkH+tLYl1G3bCJw4QUjT0IA0MJfLMX/5MqmHD0nZNiUgtGULG/r7ccXjINWuWbnmqsuFu7sb/4EDNLlcBIEioNs2KUAFGmMxmk6dQm1tRZI+3xBqEkiShBQOox05Ql1nJ26gvuIUAWiKQnj/ftStW5Fk+YuF+NldsbBA4cYN9KdPmaso8Fc62ZJtk7l1C2t0FGdxESHE1xE4hkHh6lX0oSHSpRIewC/LrPf7CSgKNpCamkI/fZr8pUuIZBIcpyaBunahPp1mub+fDyMjGKUSChCsq6Nh717q9u2jbnSU0uAgRrFIwTThyhUCqRS+hgZKLS1fJvDm87SNjZGcnsYUAjcQ8vsJ9/Xh7u1FDgRQN20iks3iDA+zZFmYhQLqzZtsj8WY3L0baY2Sf55ICCKpFHUzM2SEQAJCHg+hY8fQDh5EDgRWDoRC1J88SePRo2geD0XAcBy8b98Sv3ePYDZbbf2rFQjLIphMsmDbaEBQVQkePox26BCSz7e6i4bDrDt+nGYhmBscpFAskheC4Js3bPR4qHphFYEnlcI7P4/jdqNpGu8TCe4oCsXr1z//F2ga3+/cSfTJE0qmSVYIsKzaOZDcbjKyzFIiwfvt21kMBLA07YsetzWNd+3tLLW0sH5igvT8PH9Go/z44kX+E4LGjg7GDYOcy4XlOEgLC//5P/5LCFzxOPLmzWyIx+m6fduu7v0NVGqyTSycKksAAAAASUVORK5CYII=)}
.unchecked-box {list-style-image: url(data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAABgAAAAYCAYAAADgdz34AAAAAXNSR0IArs4c6QAAAAZiS0dEAP8A/wD/oL2nkwAAAAlwSFlzAAANOgAADMQBiN+4gQAAAAd0SU1FB9gKGQ8qAt8h3m8AAAAZdEVYdENvbW1lbnQAQ3JlYXRlZCB3aXRoIEdJTVBXgQ4XAAAA60lEQVRIx+2VsQqDMBRF70sCLg5OLoKgjk7+lJ/hh+STXBwcnRz8ArMEkrxOFktbaC3tULzTg5e8k5vADXDq70VbobXmvu/hvQczg4heHrJfXxQFuq67blZbMc8zpmlCXddIkgTOuZcBUko45zCOI6y1Nz2xFSEEZFmGOI7fGg4A3nsQEZqmuXOu9jallACAtm3fvmutNaIoAjM/dkBECCF89KCbk4eAb+kEnIAT8EsAM0OIz3hSyrssUvss8t5fg+uIrLXPs0gIgWVZYIyBUurQyYdheO4gz3NUVQVjDNZ1PfSjpWmKsixvehfB9GBZ3NndrgAAAABJRU5ErkJggg==)}
ul {list-style-image: none}
/* ul rule needed to reset style for sub-bullets */
</style>
</head>
<body>
<!-- Wiki content -->
<h1>PyPointChargePotentialSolver</h1>
<p>
Created maanantai 19 joulukuu 2016
</p>
<h2>Description</h2>
<p>
A Python module for solving electrostatic potential of a distribution of point charges using periodic boundary conditions. Represents point charges as Gaussian charge distributions. The values of the Gaussians are gathered to points on a 3D grid and the resulting charge distribution on grid is transformed using FFT to k-space. Poisson equation is solved in k-space for the electrostatic potential and the result is inverse transformed back to real space.
</p>
<p>
Since the point charges are smeared, the electric potential near the point charges differs from the potential near actual point charges but further away these two potentials converge. <i>gaussian_width</i> parameter defines the smearing width and thus what is actually "near". The smaller the <i>gaussian_width</i> is the more dense grid of points is needed and the heavier the calculation gets.
</p>
<h2>Requirements</h2>
<p>
<ul>
<li>Cython</li>
<li>NumPy</li>
</ul>
</p>
<p>
<b>Cython</b> must be installed in order to compile the computationally heavy part of the solver. See installition instructions on<br>
<a href="http://docs.cython.org/src/quickstart/install.html" title="http://docs.cython.org/src/quickstart/install.html" class="http">http://docs.cython.org/src/quickstart/install.html</a><br>
in case you need to install it yourself.
</p>
<p>
In addition, the solver uses <b>NumPy</b>, and the examples use <b>matplotlib</b> for plotting and the 'nacl_example.py' uses <b>ASE (Atomic Simulation Environment)</b> to create the NaCl structure.
</p>
<h2>Installition</h2>
<p>
Put the directory with 'electrostatic_potential_solver.py' and 'gaussians_to_grid.pyx' to your PYTHONPATH environment variable, or run the solver in the directory where the files are.
</p>
<h2>Usage</h2>
<p>
See the folder 'examples' for Python script examples with possible use cases of the solver and read the instructions below to understand what happens in them.
</p>
<p>
The solver is a Python class called 'ElectrostaticPotentialSolver'. To initialize the solver object, use command
</p>
<p>
solver = ElectrostaticPotentialSolver(atom_coordinates, atom_charges, simulation_cell)
</p>
<p>
where 'atom_coordinates' is a NumPy array of shape (n_atoms, 3) containing the coordinates of the atoms and 'atom_charges' is an array of shape (n_atoms,) containing the charges of the atoms. 'simulation_cell' is an array of shape (3,) containing the length of the simulation cell in each dimension. The solver needs three computational parameters which determine the accuracy of the solution. These parameters are set using commands
</p>
<p>
solver.set_parameter('k_cutoff', 3.0)<br>
solver.set_parameter('gaussian_width', 0.5)<br>
solver.set_parameter('gaussian_cutoff', 6.0)
</p>
<p>
This set of parameter values has been tested to provide accurate results. 'k_cutoff' is the most important parameter, and you can try smaller values to run faster or higher ones to obtain better accuracy. To actually solve the potential, run
</p>
<p>
solver.solve_potential()
</p>
<p>
You can obtain the solution on the grid used by the solver by calling
</p>
<p>
xs, ys, zs, pot_grid = solver.get_data_grid('potential')
</p>
<p>
In addition to potential, you can get the components of the electric field by substituting 'potential' by 'efield_x', 'efield_y' or 'efield_z'. For convenience, the class provides also a method for obtaining a slice of values at a certain height along one of the axes. It can be used as
</p>
<p>
xs, ys, pot_slice = solver.get_data_slice('potential', height, normal_direction='z')
</p>
<p>
where 'height' determines the position along one of the axes at which the slice is taken and 'normal_direction' determines that axis. normal_direction='z' means slice along xy plane, for example.
</p>
<!-- End wiki content -->
</body>
</html>