dphi = (-V_prime / (3 * H_inf) * dt) + noise
phi.append(phi[-1] + dphi)
return phi
Calculate residual Lambda
phi_end = simulate_inflaton(phi0)[-1]
Lambda = Lambda0 - 0.5 * np.sum(n_i**2 * q_i**2) * M_Pl**4 + V(phi_end)
Output
print(f" = {Lambda:.2e} M_Pl^4")
```
Code Explanation:
1. Flux Parameters: The charges `q_i` and the numbers `n_i` are randomized with Gaussian and Poisson distributions.
Beri Komentar
Belum ada komentar. Jadilah yang pertama untuk memberikan komentar!