Water-gas shift reaction on Pt(111).¶
This tutorial is intended to show how to use pyZacros from a Zacros perspective. Thus, we will literally show how to translate the Zacros input files to a pyZacros script. To do that, we use the system described in the Zacros tutorial Mapping DFT Energies to Zacros Input as an example. All physical/chemical description of the system is described in detail there. This example shows how to include gas species, transition states, as well as stable surface species, and the lateral interactions between them. All information about the energetics is obtained via density functional theory (DFT) calculations. Importantly, lateral interactions are also discussed and included in the example’s script. We invited you first to get familiar with the tutorial cited above to quickly appreciate the parallel between the Zacros input files and the pyZacros objects. This will allow you to follow line-by-line the example’s python script easily.
You can download the example’s python script from this link WaterGasShiftOnPt111.py.
If everything is working well, you should get the following information in the standard output and the figure shown at the end.
1$ amspython WaterGasShiftOnPt111.py
2PLAMS working folder: /home/user/pyzacros/examples/WaterGasShiftOnPt111/plams_workdir
3---------------------------------------------------------------------
4simulation_input.dat
5---------------------------------------------------------------------
6random_seed 123278
7temperature 500.0
8pressure 10.0
9
10snapshots on time 0.0005
11process_statistics on time 0.0005
12species_numbers on time 0.0005
13event_report off
14max_steps infinity
15max_time 0.25
16
17n_gas_species 4
18gas_specs_names CO H2 H2O CO2
19gas_energies 0.00000e+00 0.00000e+00 0.00000e+00 -6.15000e-01
20gas_molec_weights 2.79949e+01 2.01560e+00 1.80105e+01 4.39898e+01
21gas_molar_fracs 1.00000e-05 0.00000e+00 9.50000e-01 0.00000e+00
22
23n_surf_species 6
24surf_specs_names CO* H* H2O* OH* O* COOH*
25surf_specs_dent 1 1 1 1 1 1
26
27finish
28---------------------------------------------------------------------
29lattice_input.dat
30---------------------------------------------------------------------
31lattice default_choice
32 hexagonal_periodic 1.0 8 10
33end_lattice
34---------------------------------------------------------------------
35energetics_input.dat
36---------------------------------------------------------------------
37energetics
38
39cluster CO_point
40 sites 1
41 lattice_state
42 1 CO* 1
43 site_types 1
44 graph_multiplicity 1
45 cluster_eng -2.07700e+00
46end_cluster
47
48cluster H2O_point
49 sites 1
50 lattice_state
51 1 H2O* 1
52 site_types 1
53 graph_multiplicity 1
54 cluster_eng -3.62000e-01
55end_cluster
56
57cluster OH_point
58 sites 1
59 lattice_state
60 1 OH* 1
61 site_types 1
62 graph_multiplicity 1
63 cluster_eng 8.30000e-01
64end_cluster
65
66cluster O_point
67 sites 1
68 lattice_state
69 1 O* 1
70 site_types 1
71 graph_multiplicity 1
72 cluster_eng 1.29800e+00
73end_cluster
74
75cluster H_point
76 sites 1
77 lattice_state
78 1 H* 1
79 site_types 1
80 graph_multiplicity 1
81 cluster_eng -6.19000e-01
82end_cluster
83
84cluster COOH_point
85 sites 1
86 lattice_state
87 1 COOH* 1
88 site_types 1
89 graph_multiplicity 1
90 cluster_eng -1.48700e+00
91end_cluster
92
93cluster CO_pair_1NN
94 sites 2
95 neighboring 1-2
96 lattice_state
97 1 CO* 1
98 2 CO* 1
99 site_types 1 1
100 graph_multiplicity 1
101 cluster_eng 5.60000e-01
102end_cluster
103
104cluster OH_H_1NN
105 sites 2
106 neighboring 1-2
107 lattice_state
108 1 OH* 1
109 2 H* 1
110 site_types 1 1
111 graph_multiplicity 1
112 cluster_eng 2.10000e-02
113end_cluster
114
115cluster O_H_1NN
116 sites 2
117 neighboring 1-2
118 lattice_state
119 1 O* 1
120 2 H* 1
121 site_types 1 1
122 graph_multiplicity 1
123 cluster_eng 1.98000e-01
124end_cluster
125
126cluster CO_OH_1NN
127 sites 2
128 neighboring 1-2
129 lattice_state
130 1 CO* 1
131 2 OH* 1
132 site_types 1 1
133 graph_multiplicity 1
134 cluster_eng 6.60000e-02
135end_cluster
136
137cluster CO_O_1NN
138 sites 2
139 neighboring 1-2
140 lattice_state
141 1 CO* 1
142 2 O* 1
143 site_types 1 1
144 graph_multiplicity 1
145 cluster_eng 4.23000e-01
146end_cluster
147
148end_energetics
149---------------------------------------------------------------------
150mechanism_input.dat
151---------------------------------------------------------------------
152mechanism
153
154reversible_step CO_adsorption
155 gas_reacs_prods CO -1
156 sites 1
157 initial
158 1 * 1
159 final
160 1 CO* 1
161 site_types 1
162 pre_expon 2.22600e+07
163 pe_ratio 2.13700e-06
164 activ_eng 0.00000e+00
165end_reversible_step
166
167reversible_step H2_dissoc_adsorp
168 gas_reacs_prods H2 -1
169 sites 2
170 neighboring 1-2
171 initial
172 1 * 1
173 2 * 1
174 final
175 1 H* 1
176 2 H* 1
177 site_types 1 1
178 pre_expon 8.29900e+07
179 pe_ratio 7.96600e-06
180 activ_eng 0.00000e+00
181end_reversible_step
182
183reversible_step H2O_adsorption
184 gas_reacs_prods H2O -1
185 sites 1
186 initial
187 1 * 1
188 final
189 1 H2O* 1
190 site_types 1
191 pre_expon 2.77600e+02
192 pe_ratio 2.66500e-06
193 activ_eng 0.00000e+00
194end_reversible_step
195
196reversible_step H2O_dissoc_adsorp
197 sites 2
198 neighboring 1-2
199 initial
200 1 H2O* 1
201 2 * 1
202 final
203 1 OH* 1
204 2 H* 1
205 site_types 1 1
206 pre_expon 1.04200e+13
207 pe_ratio 1.00000e+00
208 activ_eng 7.77000e-01
209end_reversible_step
210
211reversible_step OH_decomposition
212 sites 2
213 neighboring 1-2
214 initial
215 1 * 1
216 2 OH* 1
217 final
218 1 O* 1
219 2 H* 1
220 site_types 1 1
221 pre_expon 1.04200e+13
222 pe_ratio 1.00000e+00
223 activ_eng 9.40000e-01
224end_reversible_step
225
226reversible_step COOH_formation
227 sites 2
228 neighboring 1-2
229 initial
230 1 CO* 1
231 2 OH* 1
232 final
233 1 * 1
234 2 COOH* 1
235 site_types 1 1
236 pre_expon 1.04200e+13
237 pe_ratio 1.00000e+00
238 activ_eng 4.05000e-01
239end_reversible_step
240
241step COOH_decomposition
242 gas_reacs_prods CO2 1
243 sites 2
244 neighboring 1-2
245 initial
246 1 COOH* 1
247 2 * 1
248 final
249 1 * 1
250 2 H* 1
251 site_types 1 1
252 pre_expon 1.04200e+13
253 activ_eng 8.52000e-01
254end_step
255
256step CO_oxidation
257 gas_reacs_prods CO2 1
258 sites 2
259 neighboring 1-2
260 initial
261 1 CO* 1
262 2 O* 1
263 final
264 1 * 1
265 2 * 1
266 site_types 1 1
267 pre_expon 1.04200e+13
268 activ_eng 9.88000e-01
269end_step
270
271end_mechanism
272[08.02|15:34:32] JOB plamsjob STARTED
273[08.02|15:34:32] JOB plamsjob RUNNING
274[08.02|15:34:57] JOB plamsjob FINISHED
275[08.02|15:34:57] JOB plamsjob SUCCESSFUL
276[08.02|15:35:10] PLAMS run finished. Goodbye
