Cómo trazar la capacidad calorífica específica en función de la temperatura, en una simulación MC

Hay dos métodos para calcular la capacidad calorífica. El método uno es diferenciar numéricamente la energía, es decir

C_V (T) = [E (T + (\ delta T) / 2) – E (T – (\ delta T) / 2)] / (\ delta T).

Para obtener el calor específico, que es la capacidad calorífica por partícula, divida por el número de partículas, N.

El método dos es utilizar la fórmula del teorema de fluctuación-disipación, como sugirió en su pregunta.

Para un algoritmo correcto, estos dos métodos deben estar de acuerdo, dentro de los errores estadísticos. Si no están de acuerdo, significa que una de las dos posibilidades es verdadera. Ya sea:

A) Hay un error en su código de computadora. Por ejemplo, puede estar usando un generador de números pseudoaleatorio incorrecto, o puede estar usando un número pseudoaleatorio bueno de una manera incorrecta.

o

B) No ha ejecutado su simulación el tiempo suficiente.

En la práctica, el Método Uno parece converger numéricamente más rápido que el Método Dos. Por lo tanto, siempre debe calcular el calor específico utilizando ambos métodos y asegurarse de que estén de acuerdo. Si no están de acuerdo, sus resultados no son válidos.

NOTA: Hay un error en su fórmula FDT, como se indica en la pregunta. Debería decir (E – (E1)) ^ 2, no (E – (E1) ^ 2).