Stochastic Selection of Cosmological Constant in Chaotic Inflation: A String Theory Landscape Approach

n_i = np.random.randint(-10, 11, size=N_flux)  

q_i = np.random.normal(0.1, 0.02, size=N_flux)  

# Inflaton potential (KKLT-like)  

def V(phi):  

    return 1e-10 * M_Pl**4 * (1 - np.exp(-0.1 * phi/M_Pl))  

# Stochastic Langevin simulation  

def simulate_inflaton(phi0, steps=10000, dt=1e-2):  

    phi = [phi0]  

    for _ in range(steps):  

        V_prime = (V(phi[-1] + 1e-10) - V(phi[-1] - 1e-10 # Numerical derivative  

        noise = np.random.normal(0, np.sqrt(H_inf**3 * dt / (4 * np.pi**2))