187 |
|
GBMap%atidToGBtype(myATID) = current |
188 |
|
GBMap%GBtypes(current)%atid = myATID |
189 |
|
GBMap%GBtypes(current)%isLJ = .true. |
190 |
< |
GBMap%GBtypes(current)%d = getSigma(myATID) |
190 |
> |
GBMap%GBtypes(current)%d = getSigma(myATID) / sqrt(2.0_dp) |
191 |
|
GBMap%GBtypes(current)%l = GBMap%GBtypes(current)%d |
192 |
|
GBMap%GBtypes(current)%eps = getEpsilon(myATID) |
193 |
|
GBMap%GBtypes(current)%eps_ratio = 1.0_dp |
227 |
|
er1 = GBMap%GBtypes(i)%eps_ratio |
228 |
|
dw1 = GBMap%GBtypes(i)%dw |
229 |
|
|
230 |
< |
do j = i, nGBtypes |
230 |
> |
do j = 1, nGBtypes |
231 |
|
|
232 |
|
d2 = GBMap%GBtypes(j)%d |
233 |
|
l2 = GBMap%GBtypes(j)%l |
234 |
|
e2 = GBMap%GBtypes(j)%eps |
235 |
|
er2 = GBMap%GBtypes(j)%eps_ratio |
236 |
|
dw2 = GBMap%GBtypes(j)%dw |
237 |
< |
! write(*,*) 'd2 = ', d2,l2, d1,l2 |
238 |
< |
! GBMixingMap(i,j)%sigma0 = sqrt(d1*d1 + d2*d2) |
239 |
< |
GBMixingMap(i,j)%sigma0 = 0.5_dp*(d1 + d2) |
237 |
> |
|
238 |
> |
! Cleaver paper uses sqrt of squares to get sigma0 for |
239 |
> |
! mixed interactions. |
240 |
> |
|
241 |
> |
GBMixingMap(i,j)%sigma0 = sqrt(d1*d1 + d2*d2) |
242 |
|
GBMixingMap(i,j)%xa2 = (l1*l1 - d1*d1)/(l1*l1 + d2*d2) |
243 |
|
GBMixingMap(i,j)%xai2 = (l2*l2 - d2*d2)/(l2*l2 + d1*d1) |
244 |
|
GBMixingMap(i,j)%x2 = (l1*l1 - d1*d1) * (l2*l2 - d2*d2) / & |
257 |
|
GBMixingMap(i,j)%xp2 = xp*xp |
258 |
|
GBMixingMap(i,j)%xpap2 = xp*ap2 |
259 |
|
GBMixingMap(i,j)%xpapi2 = xp/ap2 |
258 |
– |
|
259 |
– |
if (i.ne.j) then |
260 |
– |
GBMixingMap(j,i)%sigma0 = GBMixingMap(i,j)%sigma0 |
261 |
– |
GBMixingMap(j,i)%dw = GBMixingMap(i,j)%dw |
262 |
– |
GBMixingMap(j,i)%eps0 = GBMixingMap(i,j)%eps0 |
263 |
– |
GBMixingMap(j,i)%x2 = GBMixingMap(i,j)%x2 |
264 |
– |
GBMixingMap(j,i)%xa2 = GBMixingMap(i,j)%xa2 |
265 |
– |
GBMixingMap(j,i)%xai2 = GBMixingMap(i,j)%xai2 |
266 |
– |
GBMixingMap(j,i)%xp2 = GBMixingMap(i,j)%xp2 |
267 |
– |
GBMixingMap(j,i)%xpap2 = GBMixingMap(i,j)%xpap2 |
268 |
– |
GBMixingMap(j,i)%xpapi2 = GBMixingMap(i,j)%xpapi2 |
269 |
– |
endif |
260 |
|
enddo |
261 |
|
enddo |
262 |
|
haveMixingMap = .true. |
281 |
|
gbt1 = GBMap%atidToGBtype(atomID) |
282 |
|
l = GBMap%GBtypes(gbt1)%l |
283 |
|
d = GBMap%GBtypes(gbt1)%d |
284 |
< |
cutValue = 2.5_dp*max(l,d) |
284 |
> |
|
285 |
> |
! sigma is actually sqrt(2)*l for prolate ellipsoids |
286 |
> |
|
287 |
> |
cutValue = 2.5_dp*sqrt(2.0_dp)*max(l,d) |
288 |
|
|
289 |
|
end function getGayBerneCut |
290 |
|
|