35 double CS_VS80(
long int ipISO,
long int nelem,
long int ipHi,
long int ipLo,
double Aul,
double temp,
long int Collider )
83 double Apn, bp, Bpn, delta;
88 long ipISO, nelem, ipHi, ipLo, Collider;
120 Apn = 2.*ryd/Epn*abs_osc_str;
122 bp = 1.4*log(p)/p - .7/p - .51/p/p + 1.16/p/p/p - 0.55/p/p/p/p;
124 Bpn = 4.*ryd*ryd/n/n/n*(1./Epn/Epn + 4./3.*Epi/
POW3(Epn) + bp*Epi*Epi/
powi(Epn,4));
126 delta = exp(-Bpn/Apn) - 0.4*Epn/ryd;
128 gamma = ryd*(8. + 23.*
POW2(s/p))/
129 ( 8. + 1.1*n*s + 0.8/s/s + 0.4*(pow(n,1.5))/(pow(s,0.5))*fabs(s-1.0) );
135 if( energy/2./ryd+delta <= 0 )
141 cross_section = 2.*ryd/(energy + gamma)*(Apn*log(energy/2./ryd+delta) + Bpn);
142 cross_section =
MAX2( cross_section, 0. );
146 coll_str = cross_section * gLo * (energy/
EVRYD);
181 denom = pow(epi,2.33) + 4.38*pow(epi,1.72) + 1.32*epi;
184 coef = 9.56e-6 * pow(t_eV,-1.5) *
StatesElem[ipISO][nelem][level].ConBoltz / denom;
221 static const double arrH[4] = {1.48,3.64,5.93,8.32};
222 static const double arrRnp[8] = {2.20,1.90,1.73,1.65,1.60,1.56,1.54,1.52};
223 static const double arrg[10] = {0.8675,0.932,0.952,0.960,0.965,0.969,0.972,0.975,0.978,0.981};
225 static double small = 0.;
234 charge = nelem - ipISO;
278 etwo =
StatesElem[ipISO][nelem][level].ConBoltz - y*eone;
279 ethree = (
StatesElem[ipISO][nelem][level].ConBoltz - y*etwo)/2.;
282 t2 = 1./3./xn*(
StatesElem[ipISO][nelem][level].ConBoltz - y*ethree);
283 t3 = 3.*H/xn/(3. - Rnp)*(y*etwo - 2.*y*eone +
StatesElem[ipISO][nelem][level].ConBoltz);
284 t4 = 3.36*y*(eone - etwo);
285 rate = 7.69415e-9*(double)
phycon.
sqrte*9.28278e-3*
powi(xn/(charge+1),4)*g*y;
286 rate *= t1 - t2 + t3 + t4;
288 rate =
MAX2(rate,small);
293 rate2 =
MAX2(rate2,small);
298 if( rate==0. || rate2==0. )
299 HydColIon_v =
MAX2(rate,rate2);
301 HydColIon_v =
MIN2(rate,rate2);
303 ASSERT( HydColIon_v >= 0. );
304 return( HydColIon_v );
311 double Anp, bn, Bnp, delta_np;
313 double Gamma_np, p, n, g_p, g_n;
314 double ryd, s, kT_eV, rate, col_str, abs_osc_str;
345 Anp = 2.*ryd/Enp*abs_osc_str;
347 bn = 1.4*log(n)/n - .7/n - .51/n/n + 1.16/n/n/n - 0.55/n/n/n/n;
349 Bnp = 4.*ryd*ryd/p/p/p*(1./Enp/Enp + 4./3.*Eni/
POW3(Enp) + bn*Eni*Eni/
powi(Enp,4));
351 delta_np = exp(-Bnp/Anp) + 0.1*Enp/ryd;
353 Gamma_np = ryd*log(1. + n*n*n*kT_eV/ryd) * (3. + 11.*s*s/n/n) /
354 ( 6. + 1.6*p*s + 0.3/s/s + 0.8*(pow(p,1.5))/(pow(s,0.5))*fabs(s-0.6) );
356 if( 0.3*kT_eV/ryd+delta_np <= 0 )
362 rate = 1.6E-7 * pow(kT_eV, 0.5) * g_n / g_p / ( kT_eV + Gamma_np ) *
363 ( Anp * log(0.3*kT_eV/ryd + delta_np) + Bnp );