In this article, we focus on the development of a multiscale modeling and recurrent neural network (RNN) based optimization framework of a plasma etch process on a three-dimensional substrate with uniform thickness using the inductive coupled plasma (ICP). Specifically, the gas flow and chemical reactions of plasma are simulated by a macroscopic fluid model. In addition, the etch process on the substrate is simulated by a kinetic Monte Carlo (kMC) model. While long time horizon optimization cannot be completed due to the computational complexity of the simulation models, RNN models are applied to approximate the fluid model and kMC model. The training data of RNN models are generated by open-loop simulations of the fluid model and the kMC model. Additionally, the stochastic characteristic of the kMC model is presented by a probability function. The well-trained RNN models and the probability function are then implemented in computing an open-loop optimization problem, in which a moving optimization method is applied to overcome the error accumulation problem when using RNN models. The optimization goal is to achieve the desired average etching depth and average bottom roughness within the least amount of time. The simulation results show that our prediction model is accurate enough and the optimization objectives can be completed well.