Using Machine Learning to Compute Constrained Optimal Carbon Tax Rules

We develop a computational framework for deriving Pareto-improving and constrained optimal carbon tax rules in a stochastic overlapping generations (OLG) model with climate change. By integrating Deep Equilibrium Networks for fast policy evaluation and Gaussian process surrogate modeling with Bayesian active learning, the framework systematically locates optimal carbon tax schedules for heterogeneous agents exposed to climate risk. We apply our method to a 12-period OLG model in which exogenous shocks affect the carbon intensity of energy production, as well as the damage function. Constrained optimal carbon taxes consist of tax rates that are simple functions of observables and revenue-sharing rules that guarantee that the introduction of the taxes is Pareto improving. This reveals that a straightforward policy is highly effective: a Pareto-improving linear tax on cumulative emissions alone yields a 0.42% aggregate welfare gain in consumption-equivalent terms while adding further complexity to the tax provides only a marginal increase to 0.45%. The application demonstrates that the proposed approach produces scalable tools for macro-policy design in complex stochastic settings. Beyond climate economics, the framework offers a template for systematically analyzing welfare-improving policies in various heterogeneous-agent problems.