Adds "FieldFinal" option to InflationEvolution and related functions. This allows one to stop inflation reliably if the position of the minimum of the potential is known a priori.
Tests
First, evaluate a natural inflation model with default options: InflationEvolution[1/2 a'[t]^2 - (1 - Cos[a[t]/5]), {a[t], 15., 0}, t], observe the inflation-end time is 209.
Plot the field evolution, Plot[%["Field"][t], {t, 0, 209}], observe that the field goes all the way to zero.
Now, evaluate the same model with "FieldFinal" set to 10: InflationEvolution[1/2 a'[t]^2 - (1 - Cos[a[t]/5]), {a[t], 15., 0}, t, "FieldFinal" -> 10], note inflation ends at t = 133 this time.
Plot the evolution: Plot[%["Field"][t], {t, 0, 133}], note the field stops at 10.
Changes
"FieldFinal"
option toInflationEvolution
and related functions. This allows one to stop inflation reliably if the position of the minimum of the potential is known a priori.Tests
InflationEvolution[1/2 a'[t]^2 - (1 - Cos[a[t]/5]), {a[t], 15., 0}, t]
, observe the inflation-end time is 209.Plot[%["Field"][t], {t, 0, 209}]
, observe that the field goes all the way to zero.InflationEvolution[1/2 a'[t]^2 - (1 - Cos[a[t]/5]), {a[t], 15., 0}, t, "FieldFinal" -> 10]
, note inflation ends at t = 133 this time.Plot[%["Field"][t], {t, 0, 133}]
, note the field stops at 10.