gw.h

// Copyright (C) 2024 Lisa Biermann, Margarete Mühlleitner, Rui Santos, João
// Viana
//
// SPDX-FileCopyrightText: 2024 Lisa Biermann, Margarete Mühlleitner, Rui
// Santos, João Viana
//
// SPDX-License-Identifier: GPL-3.0-or-later
 
#pragma once
 
#include <BSMPT/bounce_solution/bounce_solution.h> // BounceSolution
#include <BSMPT/models/SMparam.h>
#include <BSMPT/utility/Logger.h>
#include <gsl/gsl_math.h>
 
namespace BSMPT
{
 
struct GravitationalWaveData
{
  bool swON             = true;  // enable sound wave contribution by default
  bool turbON           = true;  // enable turbulence contribution by default
  double transitionTemp = false; // transition temperature
  double PTStrength     = false; // strength of EW phase transition
  double InvTimeScale   = false; // inverse time scale beta/H
  double vb             = false; // bubble wall velocity
 
  double kappa_sw = false;   // efficiency factor for sound waves
  double K_sw     = false;   // kinetic energy fraction in fluid of total bubble
                             // energy for sound waves
  double HR         = false; // time scale times max. velocity for sound waves
  double kappa_turb = false; // efficiency factor for turbulence
  double K_turb     = false; // kinetic energy fraction in fluid of total bubble
                             // energy for turbulence
  double gstar = false;      // number of eff. d.o.f.
  double Hstar = false;      // Hubble rate at percolation temperature
 
  double fPeakSoundWave        = false; // peak frequency for sound wave
  double h2OmegaPeakSoundWave  = false; // peak amplitude for sound wave
  double fPeakTurbulence       = false; // peak frequency for turbulence
  double h2OmegaPeakTurbulence = false; // peak amplitude for turbulence
 
  StatusGW status = StatusGW::NotSet; // gw calculation status
};
struct GravitationalWaveData {…};
 
class GravitationalWave
{
private:
public:
  GravitationalWaveData data;
  GravitationalWave(BounceSolution &BACalc,
                    const int &which_transition_temp = 3);
  ~GravitationalWave();
 
  double CalcEpsTurb(double epsturb_in);
 
  const double AbsErr = 0;
 
  const double RelErr = 1e-6;
 
  double h = 0.674;
 
  void CalcPeakFrequencySoundWave();
 
  void CalcPeakAmplitudeSoundWave();
 
  void CalcPeakFrequencyTurbulence();
 
  void CalcPeakAmplitudeTurbulence();
 
  double CalcGWAmplitude(double f, bool swON, bool turbON);
 
  double GetSNR(const double fmin, const double fmax, const double T = 3);
 
  friend double snr_integrand(double freq, void *params);
};
class GravitationalWave {…};
 
double SIfunc(const double f);
 
double Rfunc(const double f);
 
double powspec_density(const double f);
 
double h2OmSens(const double f);
 
struct resultErrorPair
Nintegrate_SNR(GravitationalWave &obj, const double fmin, const double fmax);
 
double
Getkappa_sw(const double &alpha, const double &vwall, const double &Csound);
 
double GetK_sw(const double &alpha, const double &vwall, const double &Csound);
 
double
GetHR(const double &invTimeScale, const double &vwall, const double &Csound);
 
double GetK_turb(const double &alpha, const double &kappa);
 
bool IsFluidTurnoverApproxOne(const double &HR, const double &K);
 
} // namespace BSMPT