wavefunction.h

Go to the documentation of this file.

 
#include "constants.h"
#include "lookuptable.h"
 
#ifndef WAVEFUNCTION_H 
#define WAVEFUNCTION_H
 
class Path;
 
// ========================================================================  
// WaveFunctionBase Class
// ========================================================================  
class WaveFunctionBase {
 
    public:
        WaveFunctionBase (const Path &,LookupTable &_lookup, string _name="constant");
        virtual ~WaveFunctionBase();
 
        string name;
 
        virtual double PsiTrial(const int) { return 1.0; };
        virtual double PsiTrial(const double) { return 1.0; };
        virtual double PsiTrial(const beadLocator &) { return 1.0; };
        virtual double delPsiTrial(const double ) { return 0.0; };
        virtual double delSqPsiTrial(const double ) { return 0.0; };
        virtual double gradSqPsiTrial(const int) { return 0.0; };
 
    protected:
        const Path &path;               
        LookupTable &lookup;            
 
};
 
 
// ========================================================================  
// SechWaveFunction Class
// ========================================================================  
class SechWaveFunction: public WaveFunctionBase {
 
    public:
        SechWaveFunction(const Path &,LookupTable &_lookup, string _name="SHO sech");
        ~SechWaveFunction();
 
        double PsiTrial (const int);
 
    private:
        double a;           // The parameter of the wave function 
 
};
 
// ========================================================================
// JastrowWaveFunction Class
// ========================================================================
class JastrowWaveFunction: public WaveFunctionBase {
    
public:
    JastrowWaveFunction(const Path &, LookupTable &_lookup,string _name="Jastrow");
    ~JastrowWaveFunction();
    
    double PsiTrial (const double);
    double PsiTrial (const int);
    double delPsiTrial(const double r);
    double delSqPsiTrial(const double r);
    double gradSqPsiTrial(const int);
 
    double twoBodyPsiTrial (const double);
    
private:
    double alpha;           // The parameter of the wave function
    double beta;            // The parameter of the wave function
    
};
 
 
 
// ========================================================================
// LiebLinigerWaveFunction Class
// ========================================================================
class LiebLinigerWaveFunction: public WaveFunctionBase {
    
public:
    LiebLinigerWaveFunction(const Path &, LookupTable &_lookup,string _name="LiebLiniger");
    ~LiebLinigerWaveFunction();
    
    double PsiTrial (const double);
    double PsiTrial (const beadLocator &bead1);
    double PsiTrial (const int);
    
    double delPsiTrial(const double r);
    double delSqPsiTrial(const double r);
    double gradSqPsiTrial(const int);
    
    
private:
    double R;           // The parameter length scale of the wave function
    double k;           // The wavevector of the wave function
    
};
 
// ========================================================================
// SutherlandWaveFunction Class
// ========================================================================
class SutherlandWaveFunction: public WaveFunctionBase {
    
public:
    SutherlandWaveFunction(const Path &, LookupTable &_lookup, double, string _name="Sutherland");
    ~SutherlandWaveFunction();
    
    double PsiTrial(const double r) {return pow(2.0*sin(pioL*r),lambda);}
    double PsiTrial (const int);
    
private:
    double lambda;          // Sutherland model \lambda
    double pioL;            // pi / L
};
 
 
 
#endif