Hi all,(adsbygoogle = window.adsbygoogle || []).push({});

I'm doing classical (Monte Carlo) simulations on crystal structures using so far the following tail correction where [tex]\rho(r)[/tex] is the pair correlation function and [tex]U(r)[/tex] the pair potential.

[tex]2N\pi\rho\int_{r_{c}}^{\infty}\mathrm{d}r\, r^{2}U(r)[/tex]

It is usual to assume that [tex]\rho(r)=1[/tex] for [tex]r>=r_{c}[/tex].

I'm wondering if this approximation is suitable for crystals.

For fcc and hcp Lennard Jones systems, lattice sum has been implemented (J. Chem. Phys., Vol. 115, No. 11 ; J. O. Hirschfelder, C. F. Curtiss, and R. B. Bird, Molecular Theory of Gases and Liquids, Wiley, New York, 1954, pp. 1036–1039) and I'd like to know if there is a general approach to apply it to any crystal structure.

Martin Dove talks about that as well in his book "Introduction to Lattice Dynamics" but he uses a kind of Ewald summation and I'd prefer to use the same approach as for the Lennard Jonesium.

Thanks in advance,

Eric.

**Physics Forums - The Fusion of Science and Community**

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Tail correction for crystals.

Can you offer guidance or do you also need help?

Draft saved
Draft deleted

Loading...

Similar Threads - Tail correction crystals | Date |
---|---|

A High-energy tail of H electron momentum distribution? 1/p^6? | Jul 1, 2017 |

B Is this plane correctly drawn? | May 7, 2017 |

Why does corrected Maxwell Boltzmann have decimals? | May 14, 2015 |

Hartree-Fock correction to e-e interaction | Nov 28, 2014 |

**Physics Forums - The Fusion of Science and Community**