search

Dougal-s / Baryogenesis-Stability

1 stars 1 forks source link

Correct thermal spectrum #1

Open j824h opened 5 years ago

j824h commented 5 years ago

In initialise.cu, the initial condition of a field, that is, Phi and PhiDot are being set with PreFFTInitPhi and PreFFTInitPhiDot. We already have the option thermal_spectrum but it is insufficient, since

  1. In the BE distribution, the mass of the field and the temperature are naively set to 1,
  2. Normalization is naively set to get the dimension in units right,
  3. PhiDot gets the same scale as Phi.

To improve,

  1. Mass and temperature should be obtained as arguments, to be suitably used to calculate the distribution,
  2. Normalization needs to be better considered.
  3. Use the correct dispersion relation to get omega for given kx (assuming free field for now?). Then, instead of setting Phi and PhiDot separately, set the two modes of the oscillation with BE distribution. (Maybe this is just same as giving PhiDot additional factor of omega?)