Révision a581d3c1 modules/screening.py
b/modules/screening.py | ||
---|---|---|
311 | 311 |
|
312 | 312 |
def correct_coll(molec, slab, ctr_coord, site_coord, num_pts, |
313 | 313 |
min_coll_height, norm_vect, slab_nghbs, molec_nghbs, |
314 |
coll_coeff): |
|
314 |
coll_coeff, height=2.5):
|
|
315 | 315 |
# TODO Rethink this function |
316 | 316 |
"""Tries to adsorb a molecule on a slab trying to avoid collisions by doing |
317 | 317 |
small rotations. |
... | ... | |
331 | 331 |
@param coll_coeff: The coefficient that multiplies the covalent radius of |
332 | 332 |
atoms resulting in a distance that two atoms being closer to that is |
333 | 333 |
considered as atomic collision. |
334 |
@param height: Height on which to try adsorption |
|
334 | 335 |
@return collision: bool, whether the structure generated has collisions |
335 | 336 |
between slab and adsorbate. |
336 | 337 |
""" |
... | ... | |
345 | 346 |
slab_molec = deepcopy(slab) |
346 | 347 |
molec.euler_rotate(k * d_angle, k * d_angle / 2, k * d_angle, |
347 | 348 |
center=ctr_coord) |
348 |
add_adsorbate(slab_molec, molec, site_coord, ctr_coord, 2.5,
|
|
349 |
add_adsorbate(slab_molec, molec, site_coord, ctr_coord, height,
|
|
349 | 350 |
norm_vect=norm_vect) |
350 | 351 |
collision = check_collision(slab_molec, slab_num_atoms, min_coll_height, |
351 | 352 |
norm_vect, slab_nghbs, molec_nghbs, |
Formats disponibles : Unified diff