df = pd.DataFrame(np.random.randint(0,10,(10,2)), columns=["random", "a"]) df[["sq_a","cube_a"]] = df.apply(lambda x: [x.a**2, x.a**3], axis=1, result_type="expand")