AtomCenteredSphericalFunctionBessel.h

Go to the documentation of this file.

// ---------------------------------------------------------------------
//
// Copyright (c) 2017-2025 The Regents of the University of Michigan and DFT-FE
// authors.
//
// This file is part of the DFT-FE code.
//
// The DFT-FE code is free software; you can use it, redistribute
// it, and/or modify it under the terms of the GNU Lesser General
// Public License as published by the Free Software Foundation; either
// version 2.1 of the License, or (at your option) any later version.
// The full text of the license can be found in the file LICENSE at
// the top level of the DFT-FE distribution.
//
// ---------------------------------------------------------------------
//
// @author  Vishal Subramanian, Kartick Ramakrishnan, Sambit Das
//
 
#ifndef DFTFE_ATOMCENTEREDSPHERICALFUNCTIONBESSEL_H
#define DFTFE_ATOMCENTEREDSPHERICALFUNCTIONBESSEL_H
 
#include "AtomCenteredSphericalFunctionBase.h"
 
#include <boost/math/quadrature/gauss_kronrod.hpp>
#include <dftUtils.h>
#include <vector>
#include <cmath>
 
namespace dftfe
{
  class AtomCenteredSphericalFunctionBessel
    : public AtomCenteredSphericalFunctionBase
  {
  public:
    /**
     * @brief Creates the Bessel function based on the 2 parameters
     * @param[in] RcParameter coefficent for the bessel function
     * @param[in] RmaxParameter the maximum R afterwhich the function returns 0.
     */
    AtomCenteredSphericalFunctionBessel(double      RcParameter,
                                        double      RmaxParameter,
                                        dftfe::uInt lParameter);
 
    double
    getRadialValue(double r) const override;
 
    dftfe::uInt
    getQuantumNumbern() const;
 
    double
    getrMinVal() const;
 
    std::vector<double>
    getDerivativeValue(double r) const override;
 
  private:
    double d_rMinVal;
    double d_NormalizationConstant;
    double d_Rc;
  };
 
} // end of namespace dftfe
 
#endif // DFTFE_ATOMCENTEREDSPHERICALFUNCTIONBESSEL_H