Global modeling of atmospheric chemistry is a grand computational challenge due to the large number of coupled chemical species, the nonlinearity and numerical stiffness of chemical mechanisms, and the interactions with transport on all scales. The U.S. National Research Counci's National Strategy for Advancing Climate Modeling identifies atmospheric chemistry as a priority frontier for Earth System Model (ESM) development (National Research Council, 2012). Current atmospheric chemistry models integrate the coupled chemical kinetic equations for the mechanism species over model time steps by using high-order implicit numerical solvers, but these solvers are expensive (Sandu et al., 1997) and often dominate the cost of an atmospheric simulation (Eastham et al., 2018). Here, we explore the potential of machine learning (ML) neural network algorithms to dramatically reduce the computational intensity of atmospheric chemistry in global simulations.