00001 /* 00002 ------------------------------------------------------------------- 00003 00004 Copyright (C) 2006, 2007, 2008, 2009, Andrew W. Steiner 00005 00006 This file is part of O2scl. 00007 00008 O2scl is free software; you can redistribute it and/or modify 00009 it under the terms of the GNU General Public License as published by 00010 the Free Software Foundation; either version 3 of the License, or 00011 (at your option) any later version. 00012 00013 O2scl is distributed in the hope that it will be useful, 00014 but WITHOUT ANY WARRANTY; without even the implied warranty of 00015 MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the 00016 GNU General Public License for more details. 00017 00018 You should have received a copy of the GNU General Public License 00019 along with O2scl. If not, see <http://www.gnu.org/licenses/>. 00020 00021 ------------------------------------------------------------------- 00022 */ 00023 00024 #ifndef O2SCL_DDC_EOS_H 00025 #define O2SCL_DDC_EOS_H 00026 00027 #include <string> 00028 #include <cmath> 00029 #include <o2scl/lib_settings.h> 00030 #include <o2scl/constants.h> 00031 #include <o2scl/collection.h> 00032 #include <o2scl/part.h> 00033 #include <o2scl/hadronic_eos.h> 00034 #include <o2scl/fermion.h> 00035 #include <o2scl/mm_funct.h> 00036 00037 #ifndef DOXYGENP 00038 namespace o2scl { 00039 #endif 00040 00041 /** 00042 \brief Relativistic mean field EOS with density dependent couplings 00043 00044 Based on \ref Typel99. 00045 00046 \future Implement the finite temperature EOS properly. 00047 */ 00048 class ddc_eos : public hadronic_eos_eden { 00049 public: 00050 00051 /// \name Masses 00052 //@{ 00053 /// nucleon mass 00054 double mnuc; 00055 /// \f$ \phi \f$ mass (in \f$ \mathrm{fm}^{-1} \f$ ) 00056 double ms; 00057 /// \f$ A_{\omega} \f$ mass (in \f$ \mathrm{fm}^{-1} \f$ ) 00058 double mw; 00059 /// \f$ A_{\rho} \f$ mass (in \f$ \mathrm{fm}^{-1} \f$ ) 00060 double mr; 00061 //@} 00062 00063 /// \name Parameters for couplings 00064 //@{ 00065 /// The coupling \f$ \Gamma_{\sigma}(\rho_{\mathrm{sat}}) \f$ 00066 double Gs; 00067 /// The coupling \f$ \Gamma_{\omega}(\rho_{\mathrm{sat}}) \f$ 00068 double Gw; 00069 /// The coupling \f$ \Gamma_{\rho}(\rho_{\mathrm{sat}}) \f$ 00070 double Gr; 00071 /// \f$ a_{\sigma} \f$ 00072 double as; 00073 /// \f$ a_{\omega} \f$ 00074 double aw; 00075 /// \f$ a_{\rho} \f$ 00076 double ar; 00077 /// \f$ b_{\sigma} \f$ 00078 double bs; 00079 /// \f$ b_{\omega} \f$ 00080 double bw; 00081 /// \f$ c_{\sigma} \f$ 00082 double cs; 00083 /// \f$ c_{\omega} \f$ 00084 double cw; 00085 /// \f$ d_{\sigma} \f$ 00086 double ds; 00087 /// \f$ d_{\omega} \f$ 00088 double dw; 00089 //@} 00090 00091 // The saturation density 00092 double rho0; 00093 00094 ddc_eos(); 00095 00096 /// Equation of state as a function of the densities 00097 virtual int calc_e(fermion &n, fermion &p, thermo &th) { 00098 return gsl_eunimpl; 00099 } 00100 00101 /** 00102 \brief Equation of state and meson field equations 00103 as a function of the density 00104 00105 This calculates the pressure and energy density as a function 00106 of \f$ \mu_n, \mu_p, \phi, A_{\omega}, A_{\rho} \f$ . When the 00107 field equations have been solved, \c f1, \c f2, and \c f3 are 00108 all zero. 00109 00110 \todo Is the thermodynamic identity is satisfied even when the 00111 field equations are not solved? Check this. 00112 */ 00113 virtual int calc_eq_e(fermion &neu, fermion &p, double sig, 00114 double ome, double rho, double &f1, 00115 double &f2, double &f3, thermo &th); 00116 00117 00118 /// Return string denoting type ("ddc_eos") 00119 virtual const char *type() { return "ddc_eos"; } 00120 00121 #ifndef DOXYGEN_INTERNAL 00122 00123 protected: 00124 00125 #endif 00126 00127 }; 00128 00129 #ifndef DOXYGENP 00130 } 00131 #endif 00132 00133 #endif
Documentation generated with Doxygen and provided under the GNU Free Documentation License. See License Information for details.
Project hosting provided by
,
O2scl Sourceforge Project Page