file models/Models/src/models/Axions.cpp

[No description available] More…

Defines

Name
MODEL
FRIEND
MODEL
MODEL
MODEL
MODEL
MODEL
MODEL
MODEL
MODEL

Detailed Description

Author: Sebastian Hoof (hoof@uni-goettingen.de)

Date:

  • 2019 Feb
  • 2020 June-July, Oct

Models for QCD axions and axion-like particles.

Authors (add name and date if you modify):

Macros Documentation

define MODEL 

#define MODEL GeneralCosmoALP

define FRIEND 

#define FRIEND DecayingDM_photon

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

define MODEL 

#define MODEL GeneralCosmoALP

Source code 

//   GAMBIT: Global and Modular BSM Inference Tool
//   *********************************************
///  \file
///
///  Models for QCD axions and axion-like particles.
///
///  *********************************************
///
///  Authors (add name and date if you modify):
///
///  \author Sebastian Hoof
///          (hoof@uni-goettingen.de)
///  \date 2019 Feb
///  \date 2020 June-July, Oct
///
///  *********************************************

#include <cmath>

#include "gambit/Models/model_macros.hpp"
#include "gambit/Models/model_helpers.hpp"
#include "gambit/Logs/logger.hpp"
#include "gambit/Utils/numerical_constants.hpp"

#include "gambit/Models/models/Axions.hpp"

#define MODEL GeneralCosmoALP
#define FRIEND DecayingDM_photon
void MODEL_NAMESPACE::GeneralCosmoALP_to_DecayingDM_photon (const ModelParameters &myparams, ModelParameters &friendparams)
{
    USE_MODEL_PIPE(FRIEND) // get pipe for "interpret as friend" function
    logger()<<"Running interpret_as_friend calculations for GeneralCosmoALP -> DecayingDM_photon ..."<<EOM;

    friendparams.setValue("lifetime", *Dep::lifetime);
    friendparams.setValue("mass", 1e-9*myparams["ma0"]); // Convert units from eV (GeneralCosmoALP) to GeV (DecayingDM_photon)
    friendparams.setValue("fraction", *Dep::DM_fraction);
}
#undef FRIEND
#undef MODEL

#define MODEL CosmoALP
void MODEL_NAMESPACE::CosmoALP_to_GeneralCosmoALP (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for CosmoALP -> GeneralCosmoALP ..."<<EOM;

    const double alpha_red = alpha_EM/(2.0*pi);
    double fa  = myparams["fa"];

    parentparams.setValue("gagg", alpha_red*myparams["Cagg"]/fa);
    parentparams.setValue("gaee", 0);
    parentparams.setValue("gaN", 0);
    parentparams.setValue("fa", fa);
    parentparams.setValue("ma0", myparams["ma0"]);
    parentparams.setValue("Tchi", 1E99);
    parentparams.setValue("beta", 0);
    parentparams.setValue("thetai", myparams["thetai"]);
    parentparams.setValue("f0_thermal", myparams["f0_thermal"]);
    parentparams.setValue("T_R", myparams["T_R"]);
}
#undef MODEL

#define MODEL GeneralALP
void MODEL_NAMESPACE::GeneralALP_to_GeneralCosmoALP (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for GeneralALP -> GeneralCosmoALP ..."<<EOM;

    parentparams.setValue("gagg", myparams["gagg"]);
    parentparams.setValue("gaee", myparams["gaee"]);
    parentparams.setValue("gaN", myparams["gaN"]);
    parentparams.setValue("fa", myparams["fa"]);
    parentparams.setValue("ma0", myparams["ma0"]);
    parentparams.setValue("Tchi", myparams["Tchi"]);
    parentparams.setValue("beta", myparams["beta"]);
    parentparams.setValue("thetai", myparams["thetai"]);
    // Set f0_thermal = 0, i.e. no thermal component.
    parentparams.setValue("f0_thermal", 0);
    // Use extremely high reheating temperature to guarantee that Tosc < T_R
    parentparams.setValue("T_R", 1.0E99);
}
#undef MODEL

#define MODEL CosmoALP_gg_tau
void MODEL_NAMESPACE::CosmoALP_gg_tau_to_GeneralCosmoALP (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for CosmoALP_gg_tau -> GeneralCosmoALP ..."<<EOM;

    // Compute coupling out of the lifetime and mass
    parentparams.setValue("gagg", sqrt(64e27 * pi * hbar/myparams["tau"] / pow(myparams["ma0"],3)) );
    // Set matter couplings to 0
    parentparams.setValue("gaee", 0.0);
    parentparams.setValue("gaN", 0.0);
    parentparams.setValue("fa", myparams["fa"]);
    parentparams.setValue("ma0", myparams["ma0"]);
    parentparams.setValue("Tchi", myparams["Tchi"]);
    parentparams.setValue("beta", myparams["beta"]);
    parentparams.setValue("thetai", myparams["thetai"]);
    parentparams.setValue("f0_thermal", myparams["f0_thermal"]);
    parentparams.setValue("T_R", myparams["T_R"]);
}
#undef MODEL

#define MODEL QCDAxion
void MODEL_NAMESPACE::QCDAxion_to_GeneralALP (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for QCDAxion -> GeneralALP ..."<<EOM;

    const double alpha_red = alpha_EM/(2.0*pi);
    double fa  = myparams["fa"];
    double L2  = myparams["LambdaChi"]*myparams["LambdaChi"];
    double EoN = myparams["EoverN"];
    double CG  = myparams["CaggQCD"];

    parentparams.setValue("gagg", alpha_red*std::fabs(EoN-CG)/fa);
    parentparams.setValue("gaee", m_electron*myparams["Caee"]/fa);
    parentparams.setValue("gaN", myparams["CaN"]/fa); // N.B. the energy scale of the 'nucleon mass' is ~ 1 GeV
    parentparams.setValue("fa", fa);
    parentparams.setValue("ma0", 1E+3*L2/fa);
    parentparams.setValue("Tchi", myparams["Tchi"]);
    parentparams.setValue("beta", myparams["beta"]);
    parentparams.setValue("thetai", myparams["thetai"]);
}
#undef MODEL

#define MODEL KSVZAxion
void MODEL_NAMESPACE::KSVZAxion_to_QCDAxion (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for KSVZAxion -> QCDAxion ..."<<EOM;

    const double prefactor = 3.0*alpha_EM*alpha_EM/(2.0*pi);
    const double scale     = 1.0;

    double EoN     = myparams["EoverN"];
    double CaggQCD = myparams["CaggQCD"];
    double fa      = myparams["fa"];

    parentparams.setValue("EoverN", EoN);
    parentparams.setValue("CaggQCD", CaggQCD);
    parentparams.setValue("Caee", prefactor*(EoN*std::log(fa/m_electron) - CaggQCD*std::log(scale/m_electron)));
    parentparams.setValue("CaN", myparams["CaN"]);
    parentparams.setValue("fa", fa);
    parentparams.setValue("LambdaChi", myparams["LambdaChi"]);
    parentparams.setValue("Tchi", myparams["Tchi"]);
    parentparams.setValue("beta", myparams["beta"]);
    parentparams.setValue("thetai", myparams["thetai"]);
}
#undef MODEL

#define MODEL DFSZAxion_I
void MODEL_NAMESPACE::DFSZAxion_I_to_QCDAxion (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for DFSZAxion I -> QCDAxion ..."<<EOM;

    double angle = std::atan(myparams["tanbeta"]);
    double s2    = std::sin(angle);
           s2    = s2*s2;

    parentparams.setValue("EoverN", myparams["EoverN"]);
    parentparams.setValue("CaggQCD", myparams["CaggQCD"]);
    parentparams.setValue("Caee", s2/3.0);
    parentparams.setValue("CaN", myparams["CaN"]);
    parentparams.setValue("fa", myparams["fa"]);
    parentparams.setValue("LambdaChi", myparams["LambdaChi"]);
    parentparams.setValue("Tchi", myparams["Tchi"]);
    parentparams.setValue("beta", myparams["beta"]);
    parentparams.setValue("thetai", myparams["thetai"]);
}
#undef MODEL

#define MODEL DFSZAxion_II
void MODEL_NAMESPACE::DFSZAxion_II_to_QCDAxion (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for DFSZAxion II -> QCDAxion ..."<<EOM;

    double angle = std::atan(myparams["tanbeta"]);
    double s2    = std::sin(angle);
           s2    = s2*s2;

    parentparams.setValue("EoverN", myparams["EoverN"]);
    parentparams.setValue("CaggQCD", myparams["CaggQCD"]);
    parentparams.setValue("Caee", (1.0-s2)/3.0);
    parentparams.setValue("CaN", myparams["CaN"]);
    parentparams.setValue("fa", myparams["fa"]);
    parentparams.setValue("LambdaChi", myparams["LambdaChi"]);
    parentparams.setValue("Tchi", myparams["Tchi"]);
    parentparams.setValue("beta", myparams["beta"]);
    parentparams.setValue("thetai", myparams["thetai"]);
}
#undef MODEL

#define MODEL ConstantMassALP
void MODEL_NAMESPACE::ConstantMassALP_to_GeneralALP (const ModelParameters &myparams, ModelParameters &parentparams)
{
    logger()<<"Running interpret_as_parent calculations for ConstantMassALP -> GeneralALP ..."<<EOM;

    const double alpha_red = alpha_EM/(2.0*pi);

    double L2 = myparams["Lambda"]*myparams["Lambda"];
    double fa = myparams["fa"];

    parentparams.setValue("gagg", alpha_red*myparams["Cagg"]/fa);
    parentparams.setValue("gaee", m_electron*myparams["Caee"]/fa);
    parentparams.setValue("gaN", myparams["CaN"]/fa); // N.B. the energy scale of the 'nucleon mass' is ~ 1 GeV
    parentparams.setValue("fa", fa);
    parentparams.setValue("ma0", 1E+3*L2/fa);
    parentparams.setValue("Tchi", 1.0E99);
    parentparams.setValue("beta", 0);
    parentparams.setValue("thetai", myparams["thetai"]);
}
#undef MODEL

Updated on 2023-06-26 at 21:36:54 +0000