1
vote
0answers
20 views

How do you measure numerically the central charge of a system?

Let's say that you are doing some Monte-Carlo simulations of a statistical system on a lattice and you observe scale invariance, meaning that you are at a conformal point. Can you get a numerical ...
1
vote
1answer
71 views

I'm getting weird autocorrelations when simulating an Ising model below the critical temperature

So I'm simulating an Ising model using Monte Carlo and the Metropolis algorithm. After letting it reach equilibrium, I try to calculate the autocorrelation of the magnetization. As long as the system ...
0
votes
1answer
68 views

Failure of a simple stat mech simulation

so I did a simple simulation that I thought would yield a Boltzmann distribution, but it failed to, and I was wondering if anyone has insight into why it failed. Ok, so I had a simple discrete system ...
1
vote
0answers
61 views

How does the Lennard Jones Potential changes for interaction between particles of different sizes?

I am interested in incorporating a Lennard-Jones potential in a simulation. When the interaction only involves the same type of particle, with same characteristics, we can use reduced units, scaling ...
1
vote
0answers
47 views

Find out ground sates for large 2D classical spin model

Reaching the ground state of a large 2D classical spin model (e.g. classical Heisenberg model) might be a relatively difficult task while using conventional "flip/reject" Monte Carlo method. The ...
0
votes
1answer
83 views

Justifications for different Monte Carlo trial moves

Why do we perform different trial moves in Monte Carlo simulations in Statistical Mechanics? For example, in NVT ensemble simulations, Why only atom displacement moves? In NPT ensemble, why do we need ...
0
votes
2answers
108 views

Grand canonical ensemble with interaction, simulation doubts

First of all, this is a follow-up of my first question. The idea is the same, every state would consist of some particles in a line (with an energy associated) and the particles can not be nearer than ...
1
vote
1answer
110 views

Why NPT ensemble is used for solid state phase transitions?

In Monte Carlo simulations of solid state phase transitions, why often Isobaric Isothermal ensemble (NPT) is used ? Why not NVT ? Here, N is number of atoms, P is pressure, T is temperature and V is ...
1
vote
1answer
113 views

Deriving the change in the Helmholtz free energy in the context of the free energy perturbation method

I am reading Free Energy Calculations: Theory and Applications in Chemistry and Biology by Chipot and Pohorille. At the beginning of the text (page 19, for example), the authors define the Helmholtz ...
2
votes
1answer
178 views

Local minima in Ising model in a Monte Carlo simulation

Is there any way to check whether in a Monte Carlo simulation using Ising model is stuck in any (false) local minima of energy or not, particularly in 3D system ?
3
votes
2answers
284 views

Erogodicity in a Monte Carlo simulation

Q1: What is the ergodicity and ergodicity breaking in a Monte Carlo simulation of a statistical physics problem? Q2: How does one ensure that the ergodicity is maintained ?
1
vote
1answer
242 views

Numerical algorithms to generate a random wavefunction from a thermal ensemble

I am seeking an algorithm to generate a random wavefunction = $\sum {c_i |\varphi _i\rangle }$ from a thermal ensemble, whose density matrix $\rho \sim e^{-\beta H}$, without the need to diagonalize ...
14
votes
2answers
54 views

Sampling typical clusters between distant points in subcritical percolation

I have on several occasions wondered how one might proceed in order to sample large subcritical Bernoulli bond-percolation clusters, say on the square lattice. More precisely, let's consider the ...