Inheritance diagram for EwaldHandler2D:

Collaboration diagram for EwaldHandler2D:

Public Member Functions
	EwaldHandler2D (ParticleSet &ref, mRealType kc_in=-1.0)

LRHandlerBase *	makeClone (ParticleSet &ref) const override
	make clone More...

mRealType	evaluate (mRealType r, mRealType rinv) const override

mRealType	evaluateLR_r0 () const override
	evaluate for the self-interaction term More...

mRealType	evaluateLR (mRealType r) const override

mRealType	evaluateSR_k0 () const override
	evaluate $v_{s}(k=0) = \frac{4\pi}{V}\int_0^{r_c} r^2 v_s(r) dr$ More...

void	fillFk (const KContainer &KList)

mRealType	srDf (mRealType r, mRealType rinv) const override

mRealType	lrDf (mRealType r) const override

virtual mRealType	evaluate_vlr_k (mRealType k) const override

void	initBreakup (ParticleSet &ref) override

void	Breakup (ParticleSet &ref, mRealType rs_in) override

void	resetTargetParticleSet (ParticleSet &ref) override

Public Member Functions inherited from LRHandlerBase
	LRHandlerBase (mRealType kc)

virtual	~LRHandlerBase ()=default

mRealType	get_rc () const

mRealType	get_kc () const

mRealType	evaluate_w_sk (const std::vector< int > &kshell, const pRealType *restrict sk) const

mRealType	evaluate (const std::vector< int > &kshell, const pRealType restrict rk1_r, const pRealType restrict rk1_i, const pRealType restrict rk2_r, const pRealType restrict rk2_i) const
	evaluate $\sum_k F_{k} \rho^1_{-{\bf k}} \rho^2_{\bf k}$ More...

virtual mRealType	evaluate_slab (pRealType z, const std::vector< int > &kshell, const pRealType restrict rk1_r, const pRealType restrict rk1_i, const pRealType restrict rk2_r, const pRealType restrict rk2_i) const
	Evaluate the long-range potential with the open BC for the D-1 direction. More...

void	evaluateGrad (const ParticleSet &A, const ParticleSet &B, int specB, std::vector< pRealType > &Zat, std::vector< TinyVector< pRealType, OHMMS_DIM >> &grad1) const
	evaluate $\sum_k F_{k} \rho^1_{-{\bf k}} \rho^2_{\bf k}$ and $\sum_k F_{k} \rho^1_{-{\bf k}} \rho^2_{\bf k}$ More...

SymTensor< pRealType, OHMMS_DIM >	evaluateStress (const std::vector< int > &kshell, const pRealType rhokA_r, const pRealType rhokA_i, const pRealType rhokB_r, const pRealType rhokB_i) const
	FIX_PRECISION. More...

virtual SymTensor< mRealType, OHMMS_DIM >	evaluateLR_r0_dstrain () const
	These functions return the strain derivatives of all corresponding quantities in total energy. More...

virtual SymTensor< mRealType, OHMMS_DIM >	evaluateSR_k0_dstrain () const

virtual SymTensor< mRealType, OHMMS_DIM >	evaluateLR_dstrain (TinyVector< pRealType, OHMMS_DIM > k, pRealType kmag) const

virtual SymTensor< mRealType, OHMMS_DIM >	evaluateSR_dstrain (TinyVector< pRealType, OHMMS_DIM > r, pRealType rmag) const

Private Attributes
mRealType	alpha

mRealType	area

Additional Inherited Members
Public Attributes inherited from LRHandlerBase
DECLARE_COULOMB_TYPES int	MaxKshell
	Maxkimum Kshell for the given Kc. More...

mRealType	LR_kc
	Maximum k cutoff. More...

mRealType	LR_rc
	Maximum r cutoff. More...

Vector< mRealType >	Fk
	Fourier component for all the k-point. More...

Vector< mRealType >	Fkg
	Fourier component of the LR part, fit to optimize the gradients. More...

std::vector< SymTensor< mRealType, OHMMS_DIM > >	dFk_dstrain
	Fourier component of the LR part of strain tensor, by optimized breakup. More...

Vector< mRealType >	Fkgstrain
	Vector of df_k/dk, fit as to optimize strains. More...

Vector< mRealType >	Fk_symm
	Fourier component for each k-shell. More...

std::vector< mRealType >	coefs
	Fourier component for each k-shell Coefficient. More...

std::vector< mRealType >	gcoefs
	Coefficient for gradient fit. More...

std::vector< mRealType >	gstraincoefs
	Coefficient for strain fit. More...

Protected Attributes inherited from LRHandlerBase
std::string	ClassName

Detailed Description

Definition at line 24 of file EwaldHandler2D.h.

Constructor & Destructor Documentation

◆ EwaldHandler2D()

EwaldHandler2D	(	ParticleSet &	ref,
		mRealType	kc_in = `-1.0`
	)

Definition at line 17 of file EwaldHandler2D.cpp.

References EwaldHandler2D::alpha, qmcplusplus::app_log(), EwaldHandler2D::area, EwaldHandler2D::fillFk(), ParticleSet::getLattice(), ParticleSet::getSimulationCell(), LRHandlerBase::LR_kc, LRHandlerBase::LR_rc, and qmcplusplus::sqrt().

Referenced by EwaldHandler2D::makeClone().

   : LRHandlerBase(kc_in)
 {
   if (ref.getLattice().ndim != 2)
     throw std::runtime_error("2D Ewald requires 2D Lattice");
   LR_rc = ref.getLattice().LR_rc; // CoulombPBC needs get_rc() to createSpline4RbyVs
   LR_kc = ref.getLattice().LR_kc; // get_kc() is used in QMCFiniteSize
   alpha = std::sqrt(LR_kc/2.0/LR_rc);
   area = ref.getLattice().Volume/ref.getLattice().R(2,2);
   // report
   app_log() << "    alpha = " << alpha << " area = " << area << std::endl;
   fillFk(ref.getSimulationCell().getKLists());
 }

Member Function Documentation

◆ Breakup()

void Breakup	(	ParticleSet &	ref,
		mRealType	rs_in
	)

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 55 of file EwaldHandler2D.h.

References EwaldHandler2D::initBreakup().

55 { initBreakup(ref); }

qmcplusplus::EwaldHandler2D::initBreakup

void initBreakup(ParticleSet &ref) override

Definition: EwaldHandler2D.h:54

◆ evaluate()

mRealType evaluate	(	mRealType	r,
		mRealType	rinv
	)		const

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 33 of file EwaldHandler2D.h.

References EwaldHandler2D::alpha.

33 { return erfc(alpha*r) * rinv; }

qmcplusplus::EwaldHandler2D::alpha

mRealType alpha

Definition: EwaldHandler2D.h:59

◆ evaluate_vlr_k()

virtual mRealType evaluate_vlr_k ( mRealType k ) const

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 50 of file EwaldHandler2D.h.

   {
     throw std::runtime_error("2D Ewald vlr_k not implemented");
   }

◆ evaluateLR()

mRealType evaluateLR ( mRealType r ) const

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 37 of file EwaldHandler2D.h.

References EwaldHandler2D::alpha.

37 { return erf(alpha*r) / r; }

qmcplusplus::EwaldHandler2D::alpha

mRealType alpha

Definition: EwaldHandler2D.h:59

◆ evaluateLR_r0()

mRealType evaluateLR_r0 ( ) const

inlineoverridevirtual

evaluate $ v_s(r=0) $ for the self-interaction term

Reimplemented from LRHandlerBase.

Definition at line 34 of file EwaldHandler2D.h.

References EwaldHandler2D::alpha, and qmcplusplus::sqrt().

34 { return 2.0 * alpha / std::sqrt(M_PI); }

qmcplusplus::EwaldHandler2D::alpha

mRealType alpha

Definition: EwaldHandler2D.h:59

qmcplusplus::sqrt

MakeReturn< UnaryNode< FnSqrt, typename CreateLeaf< Vector< T1, C1 > >::Leaf_t > >::Expression_t sqrt(const Vector< T1, C1 > &l)

Definition: OhmmsVectorOperators.h:136

◆ evaluateSR_k0()

mRealType evaluateSR_k0 ( ) const

inlineoverridevirtual

evaluate $v_{s}(k=0) = \frac{4\pi}{V}\int_0^{r_c} r^2 v_s(r) dr$

Reimplemented from LRHandlerBase.

Definition at line 38 of file EwaldHandler2D.h.

References EwaldHandler2D::alpha, EwaldHandler2D::area, and qmcplusplus::sqrt().

38 { return 2.0 * std::sqrt(M_PI) / (alpha*area); }

qmcplusplus::EwaldHandler2D::alpha

mRealType alpha

Definition: EwaldHandler2D.h:59

qmcplusplus::sqrt

MakeReturn< UnaryNode< FnSqrt, typename CreateLeaf< Vector< T1, C1 > >::Leaf_t > >::Expression_t sqrt(const Vector< T1, C1 > &l)

Definition: OhmmsVectorOperators.h:136

qmcplusplus::EwaldHandler2D::area

mRealType area

Definition: EwaldHandler2D.h:60

◆ fillFk()

void fillFk ( const KContainer & KList )

Definition at line 31 of file EwaldHandler2D.cpp.

References EwaldHandler2D::alpha, EwaldHandler2D::area, LRHandlerBase::Fk, LRHandlerBase::Fk_symm, KContainer::kpts_cart, KContainer::kshell, KContainer::ksq, LRHandlerBase::MaxKshell, Vector< T, Alloc >::resize(), Vector< T, Alloc >::size(), and qmcplusplus::sqrt().

Referenced by EwaldHandler2D::EwaldHandler2D().

 {
   const mRealType knorm = 2.0*M_PI / area;
   const mRealType kalpha = 1.0 / (2.0*alpha);
   mRealType kmag, uk;
 
   Fk.resize(KList.kpts_cart.size());
   MaxKshell = KList.kshell.size() - 1;
   Fk_symm.resize(MaxKshell);
 
   for (int ks = 0, ki = 0; ks < Fk_symm.size(); ks++)
   {
     kmag = std::sqrt(KList.ksq[ki]);
     uk = knorm * erfc(kalpha*kmag)/kmag;
     Fk_symm[ks] = uk;
     while (ki < KList.kshell[ks + 1] && ki < Fk.size())
       Fk[ki++] = uk;
   }
 }

◆ initBreakup()

void initBreakup ( ParticleSet & ref )

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 54 of file EwaldHandler2D.h.

Referenced by EwaldHandler2D::Breakup().

54 {}

◆ lrDf()

mRealType lrDf ( mRealType r ) const

inlineoverridevirtual

Reimplemented from LRHandlerBase.

Definition at line 46 of file EwaldHandler2D.h.

   {
     throw std::runtime_error("2D Ewald lr derivative not implemented");
   }

◆ makeClone()

LRHandlerBase* makeClone ( ParticleSet & ref ) const

inlineoverridevirtual

make clone

Implements LRHandlerBase.

Definition at line 30 of file EwaldHandler2D.h.

References EwaldHandler2D::EwaldHandler2D().

30 { return new EwaldHandler2D(*this); }

qmcplusplus::EwaldHandler2D::EwaldHandler2D

EwaldHandler2D(ParticleSet &ref, mRealType kc_in=-1.0)

Definition: EwaldHandler2D.cpp:17

◆ resetTargetParticleSet()

void resetTargetParticleSet ( ParticleSet & ref )

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 56 of file EwaldHandler2D.h.

56 {}

◆ srDf()

mRealType srDf	(	mRealType	r,
		mRealType	rinv
	)		const

inlineoverridevirtual

Implements LRHandlerBase.

Definition at line 42 of file EwaldHandler2D.h.

   {
     throw std::runtime_error("2D Ewald sr derivative not implemented");
   }

Member Data Documentation

◆ alpha

mRealType alpha

private

Definition at line 59 of file EwaldHandler2D.h.

Referenced by EwaldHandler2D::evaluate(), EwaldHandler2D::evaluateLR(), EwaldHandler2D::evaluateLR_r0(), EwaldHandler2D::evaluateSR_k0(), EwaldHandler2D::EwaldHandler2D(), and EwaldHandler2D::fillFk().

◆ area

mRealType area

private

Definition at line 60 of file EwaldHandler2D.h.

Referenced by EwaldHandler2D::evaluateSR_k0(), EwaldHandler2D::EwaldHandler2D(), and EwaldHandler2D::fillFk().

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

/home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/Particle/LongRange/EwaldHandler2D.h
/home/pk7/projects/qmc/for_cron_doxygen/qmcpack/src/Particle/LongRange/EwaldHandler2D.cpp

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ EwaldHandler2D()

Member Function Documentation

◆ Breakup()

◆ evaluate()

◆ evaluate_vlr_k()

◆ evaluateLR()

◆ evaluateLR_r0()

◆ evaluateSR_k0()

◆ fillFk()

◆ initBreakup()

◆ lrDf()

◆ makeClone()

◆ resetTargetParticleSet()

◆ srDf()

Member Data Documentation

◆ alpha

◆ area