Package etomica.virial

Class MCMoveClusterRingPartialRegrow

java.lang.Object

etomica.integrator.mcmove.MCMove

etomica.integrator.mcmove.MCMoveBox

etomica.virial.MCMoveClusterRingPartialRegrow

public class MCMoveClusterRingPartialRegrow
extends MCMoveBox

MCMove that partially regrows the beads of a ring polymer, accepting or rejecting the move based on the sampling weight. The move can (optionally) regrow the beads such that the beads from multiple molecules are combined to form a larger ring.

Field Summary

Fields

Modifier and Type	Field	Description
protected IAtomList	atoms
protected Vector	dcom
protected MeterPotentialEnergy	energyMeter
protected double	fac
HistogramExpanding[][]	hist
protected int	iMolecule
protected int	kStart
protected AtomIteratorLeafAtoms	leafIterator
protected int	maxNumBeads
protected AtomArrayList	myAtoms
protected int	nTrial
protected int	numBeads
protected Vector[]	oldPositions
protected double[]	pkl
protected IRandom	random
protected Vector[]	rTrial
protected Space	space
protected int[][]	tangledMolecules
protected double	uNew
protected double	uOld
protected double	weightNew
protected double	weightOld
protected double	wNew
protected double	wOld

Fields inherited from class etomica.integrator.mcmove.MCMove

moveTracker, nominalFrequency, potential

Fields inherited from class etomica.integrator.mcmove.MCMoveBox

box, perParticleFrequency

Constructor Summary

Constructors

Constructor	Description
MCMoveClusterRingPartialRegrow(PotentialMaster potentialMaster, IRandom random, Space _space)
MCMoveClusterRingPartialRegrow(PotentialMaster potentialMaster, IRandom random, Space _space, int[][] tangledMolecules)

Method Summary

Modifier and Type	Method	Description
void	acceptNotify()	Method called by IntegratorMC in the event that the most recent trial is accepted.
AtomIterator	affectedAtoms()
protected double	calcWeight()
boolean	doTrial()	Method to perform trial move.
double	energyChange()
double	getChi(double temperature)	Chi is the parameter within standard Metropolis Monte Carlo.
int	getNumBeads()
int	getNumTrial()
void	rejectNotify()	Method called by IntegratorMC in the event that the most recent trial move is rejected.
void	setBox(Box p)	Sets the box on which this move acts.
void	setEnergyFactor(double factor)
void	setNumBeads(int newNumBeads)
void	setNumTrial(int newNumTrial)

Methods inherited from class etomica.integrator.mcmove.MCMove

getNominalFrequency, getTracker, setTracker

Methods inherited from class etomica.integrator.mcmove.MCMoveBox

affectedAtoms, energyChange, getBox, isNominallyPerParticleFrequency, toString

Methods inherited from class java.lang.Object

clone, equals, finalize, getClass, hashCode, notify, notifyAll, wait, wait, wait

Field Detail

atoms

protected IAtomList atoms

space

protected final Space space

random

protected final IRandom random

oldPositions

protected Vector[] oldPositions

rTrial

protected Vector[] rTrial

nTrial

protected int nTrial

pkl

protected double[] pkl

wOld

protected double wOld

wNew

protected double wNew

weightOld

protected double weightOld

weightNew

protected double weightNew

uOld

protected double uOld

uNew

protected double uNew

dcom

protected final Vector dcom

leafIterator

protected final AtomIteratorLeafAtoms leafIterator

fac

protected double fac

hist

public HistogramExpanding[][] hist

tangledMolecules

protected final int[][] tangledMolecules

myAtoms

protected final AtomArrayList myAtoms

numBeads

protected int numBeads

maxNumBeads

protected int maxNumBeads

kStart

protected int kStart

iMolecule

protected int iMolecule

energyMeter

protected final MeterPotentialEnergy energyMeter

Constructor Detail

MCMoveClusterRingPartialRegrow

public MCMoveClusterRingPartialRegrow(PotentialMaster potentialMaster,
                                      IRandom random,
                                      Space _space)

MCMoveClusterRingPartialRegrow

public MCMoveClusterRingPartialRegrow(PotentialMaster potentialMaster,
                                      IRandom random,
                                      Space _space,
                                      int[][] tangledMolecules)

Method Detail

setNumTrial

public void setNumTrial(int newNumTrial)

getNumTrial

public int getNumTrial()

setEnergyFactor

public void setEnergyFactor(double factor)

setNumBeads

public void setNumBeads(int newNumBeads)

getNumBeads

public int getNumBeads()

setBox

public void setBox(Box p)

Description copied from class: MCMoveBox

Sets the box on which this move acts. The box itself can be changed via this method, if desired.

Overrides:

setBox in class MCMoveBox

doTrial

public boolean doTrial()

Description copied from class: MCMove

Method to perform trial move. Returns false if the trial could not be attempted, for example if there were no molecules in the box and the trial is designed to displace an atom; returns true otherwise.

Specified by:

doTrial in class MCMove

calcWeight

protected double calcWeight()

getChi

public double getChi(double temperature)

Description copied from class: MCMove

Chi is the parameter within standard Metropolis Monte Carlo. The probability of accepting the move is given by min(1,chi).

Specified by:

getChi in class MCMove

Parameters:

temperature - used to compute chi

Returns:

chi

rejectNotify

public void rejectNotify()

Description copied from class: MCMove

Method called by IntegratorMC in the event that the most recent trial move is rejected. This method should cause the system to be restored to the condition before the most recent call to doTrial.

Specified by:

rejectNotify in class MCMove

acceptNotify

public void acceptNotify()

Description copied from class: MCMove

Method called by IntegratorMC in the event that the most recent trial is accepted.

Specified by:

acceptNotify in class MCMove

energyChange

public double energyChange()

Specified by:

energyChange in class MCMoveBox

affectedAtoms

public AtomIterator affectedAtoms()

Specified by:

affectedAtoms in class MCMoveBox