HardRodPotential Class Reference

Computes the effective potential from the exact two-body density matrix for hard rods in 1D. More...

#include <potential.h>

Inheritance diagram for HardRodPotential:

Collaboration diagram for HardRodPotential:

Public Member Functions
	HardRodPotential (double)
	Constructor. More...
	~HardRodPotential ()
	Destructor.
virtual double	V (const dVec &r)
	The classical potential.
double	V (const dVec &, const dVec &)
	The effective potential. More...
double	dVdlambda (const dVec &, const dVec &)
	The derivative of the effective potential with respect to lambda. More...
double	dVdtau (const dVec &, const dVec &)
	The derivative of the effective potential with respect to tau. More...
Public Member Functions inherited from PotentialBase
	PotentialBase ()
	Constructor.
virtual	~PotentialBase ()
	Destructor.
virtual dVec	gradV (const dVec &)
	The gradient of the potential.
virtual double	grad2V (const dVec &)
	Grad^2 of the potential.
virtual Array< dVec, 1 >	initialConfig (const Container *, MTRand &, const int)
	Default Initial configuration of particles. More...
void	output (const double)
	A debug method that output's the potential to a supplied separation. More...
virtual Array< double, 1 >	getExcLen ()
	Array to hold data elements. More...

Additional Inherited Members
Data Fields inherited from PotentialBase
double	tailV
	Tail correction factor.
Protected Member Functions inherited from PotentialBase
double	deltaSeparation (double sep1, double sep2) const
	Return the minimum image difference for 1D separations.

Detailed Description

Computes the effective potential from the exact two-body density matrix for hard rods in 1D.

Definition at line 1102 of file potential.h.

Constructor & Destructor Documentation

HardRodPotential()

HardRodPotential::HardRodPotential

(

double

_a

)

Constructor.

Parameters

_a	The radius of the hard rod (also the scattering length)

Definition at line 2367 of file potential.cpp.

                                            : 
    PotentialBase(),
    a(_a) {
 
}

Member Function Documentation

dVdlambda()

double HardRodPotential::dVdlambda ( const dVec &  sep1, const dVec &  sep2  )

virtual

The derivative of the effective potential with respect to lambda.

Computes the non-local two-body effective pair potential.

Tested and working with Mathematica on 2013-06-17.

Parameters

sep1	the first separation
sep2	the second separation

Returns

the derivative of the effective potential with respect to lambda

Reimplemented from PotentialBase.

Definition at line 2426 of file potential.cpp.

{
 
    double r1 = sqrt(dot(sep1,sep1));
    double r2 = sqrt(dot(sep2,sep2));
    double d1 = deltaSeparation(r1,a);
    double d2 = deltaSeparation(r2,a);
 
    double t1 = d1*d2;
    double t2 = t1/(2.0*constants()->lambda()*constants()->lambda()*constants()->tau());
 
    return ((0.5*t1/(exp(t2)-1.0))/(constants()->lambda()*constants()->lambda()*constants()->tau()));
}

Here is the call graph for this function:

dVdtau()

double HardRodPotential::dVdtau ( const dVec &  sep1, const dVec &  sep2  )

virtual

The derivative of the effective potential with respect to tau.

Computes the non-local two-body effective pair potential.

Tested and working with Mathematica on 2013-06-17.

Parameters

sep1	the first separation
sep2	the second separation

Returns

the derivative of the effective potential with respect to tau

Reimplemented from PotentialBase.

Definition at line 2451 of file potential.cpp.

{
 
    double r1 = sqrt(dot(sep1,sep1));
    double r2 = sqrt(dot(sep2,sep2));
    double d1 = deltaSeparation(r1,a);
    double d2 = deltaSeparation(r2,a);
 
    double t1 = d1*d2;
    double t2 = t1/(2.0*constants()->lambda()*constants()->tau());
 
    return ((0.5*t1/(exp(t2)-1.0))/(constants()->lambda()*constants()->tau()*constants()->tau()));
}

Here is the call graph for this function:

V()

double HardRodPotential::V ( const dVec &  sep1, const dVec &  sep2  )

virtual

The effective potential.

Computes the non-local two-body effective pair potential.

Tested and working with Mathematica on 2013-06-17.

Parameters

sep1	The first separation
sep2	The second separation

Returns

the two-body effective pair potential

Reimplemented from PotentialBase.

Definition at line 2391 of file potential.cpp.

{
    double r1 = sqrt(dot(sep1,sep1));
    double r2 = sqrt(dot(sep2,sep2));
 
    /* We need to enforce the distinguishable particle constraint at short
     * imaginary times */
    /* if ( (sep1[0]*sep2[0] < 0.0) || (r1 <= a ) || (r2 <= a) ) */ 
    /*     return LBIG; */
 
    if ( (r1 <= a ) || (r2 <= a) ) 
        return LBIG;
 
    double d1 = deltaSeparation(r1,a);
    double d2 = deltaSeparation(r2,a);
 
    double t1 = -d1*d2/(2.0*constants()->lambda()*constants()->tau());
 
    /* communicate()->file("debug")->stream() << sep1[0] << "\t" << sep2[0] << "\t" */ 
    /*     << -log(1.0-exp(t1)) << endl; */
 
    return (-log(1.0 - exp(t1)));
}

Here is the call graph for this function:

---

The documentation for this class was generated from the following files:

• include/potential.h
• src/potential.cpp