EPRPABreakup< T > Struct Template Reference

Collaboration diagram for EPRPABreakup< T >:

Public Member Functions
	EPRPABreakup ()
void	reset (ParticleSet &ref)
void	reset (ParticleSet &ref, T rs)
T	operator() (T r, T rinv) const
T	df (T r) const
T	Fk (T k, T rc) const
T	Xk (T k, T rc) const
T	integrate_r2 (T rc) const
T	Uk (T kk) const
	return RPA value at |k| More...
T	derivUk (T kk) const
	return d u(k)/d rs More...

Public Attributes
T	Rs
T	Kf
T	Density
T	NormFactor

Detailed Description

template<class T = double>
struct qmcplusplus::EPRPABreakup< T >

Definition at line 334 of file LRBreakupUtilities.h.

Constructor & Destructor Documentation

EPRPABreakup()

EPRPABreakup ( )

inline

Definition at line 340 of file LRBreakupUtilities.h.

{}

Member Function Documentation

derivUk()

T derivUk ( T  kk ) const

inline

return d u(k)/d rs

Implement a correct one

Definition at line 396 of file LRBreakupUtilities.h.

{ return 0.0; }

df()

T df ( T  r ) const

inline

Definition at line 365 of file LRBreakupUtilities.h.

{ return 0.0; }

Fk()

T Fk ( T  k, T  rc  ) const

inline

Definition at line 367 of file LRBreakupUtilities.h.

References EPRPABreakup< T >::Xk().

{ return -Xk(k, rc); }

integrate_r2()

T integrate_r2 ( T  rc ) const

inline

Definition at line 385 of file LRBreakupUtilities.h.

{ return 0.0; }

operator()()

T operator() ( T  r, T  rinv  ) const

inline

Definition at line 363 of file LRBreakupUtilities.h.

{ return 0.0; }

reset() [1/2]

void reset ( ParticleSet &  ref )

inline

Definition at line 342 of file LRBreakupUtilities.h.

References EPRPABreakup< T >::Density, ParticleSet::getLattice(), ParticleSet::getTotalNum(), EPRPABreakup< T >::Kf, EPRPABreakup< T >::NormFactor, qmcplusplus::pow(), and EPRPABreakup< T >::Rs.

  {
    NormFactor = 1.0 / ref.getTotalNum();
    Density    = ref.getTotalNum() / ref.getLattice().Volume;
    Rs         = std::pow(3.0 / (4.0 * M_PI * Density), 1.0 / 3.0);
    //unpolarized K_f
    Kf = std::pow(2.25 * M_PI, 1.0 / 3.0) / Rs;
  }

reset() [2/2]

void reset ( ParticleSet &  ref, T  rs  )

inline

Definition at line 351 of file LRBreakupUtilities.h.

References EPRPABreakup< T >::Density, ParticleSet::getLattice(), ParticleSet::getTotalNum(), EPRPABreakup< T >::Kf, EPRPABreakup< T >::NormFactor, qmcplusplus::pow(), and EPRPABreakup< T >::Rs.

  {
    //       NormFactor=4.0*M_PI/ref.getLattice().Volume;
    NormFactor = 1.0 / ref.getTotalNum();
    //       NormFactor=4.0*M_PI/ref.getTotalNum();
    Density = ref.getTotalNum() / ref.getLattice().Volume;
    Rs      = rs;
    //unpolarized
    Kf = std::pow(2.25 * M_PI, 1.0 / 3.0) / Rs;
  }

Uk()

T Uk ( T  kk ) const

inline

return RPA value at |k|

Parameters

kk

|k|^2

Definition at line 390 of file LRBreakupUtilities.h.

References EPRPABreakup< T >::NormFactor, and EPRPABreakup< T >::Rs.

{ return NormFactor * Rs / kk; }

Xk()

T Xk ( T  k, T  rc  ) const

inline

Definition at line 369 of file LRBreakupUtilities.h.

References EPRPABreakup< T >::Kf, EPRPABreakup< T >::NormFactor, qmcplusplus::pow(), and EPRPABreakup< T >::Rs.

Referenced by EPRPABreakup< T >::Fk().

  {
    T y = 0.5 * k / Kf;
    T Sy;
    if (y >= 1.0)
    {
      Sy = 1.0;
    }
    else
    {
      Sy = 1.5 * y - 0.5 * y * y * y;
    };
    T val = 12.0 / (k * k * k * k * Rs * Rs * Rs);
    return -0.5 * NormFactor * val * std::pow(1.0 / (Sy * Sy) + val, -0.5);
  }

Member Data Documentation

Density

T Density

Definition at line 338 of file LRBreakupUtilities.h.

Referenced by EPRPABreakup< T >::reset().

Kf

T Kf

Definition at line 337 of file LRBreakupUtilities.h.

Referenced by EPRPABreakup< T >::reset(), and EPRPABreakup< T >::Xk().

NormFactor

T NormFactor

Definition at line 339 of file LRBreakupUtilities.h.

Referenced by EPRPABreakup< T >::reset(), EPRPABreakup< T >::Uk(), and EPRPABreakup< T >::Xk().

Rs

T Rs

Definition at line 336 of file LRBreakupUtilities.h.

Referenced by EPRPABreakup< T >::reset(), EPRPABreakup< T >::Uk(), and EPRPABreakup< T >::Xk().

---

The documentation for this struct was generated from the following file:

• /home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/QMCWaveFunctions/Jastrow/LRBreakupUtilities.h