Experiments on single crystal Li X FePO 4 (LFP) nanoparticles indicate rich nonequilibrium phase behavior, such as suppression of phase separation at high lithiation rates, striped patterns of coherent phase boundaries, nucleation by binarysolid surface wetting and intercalation waves. These observations have been successfully predicted (prior to the experiments) by 1D depthaveraged phase-field models, which neglect any subsurface phase separation. In this paper, using an electrochemo-mechanical phase-field model, we investigate the coherent non-equilibrium subsurface phase morphologies that develop in the abplane of platelet-like single-crystal platelet-like Li X FePO 4 nanoparticles. Finite element simulations are performed for 2D plane-stress conditions in the abplane, and validated by 3D simulations, showing similar results. We show that the anisotropy of the interfacial tension tensor, coupled with electroautocatalytic surface intercalation reactions, plays a crucial role in determining the subsurface phase morphology. With isotropic interfacial tension, subsurface phase separation is observed, independent of the reaction kinetics, but for strong anisotropy, phase separation is controlled by surface reactions, as assumed in 1D models. Moreover, the driven intercalation reaction suppresses phase separation during lithiation, while enhancing it during delithiation, by electro-autocatalysis, in quantitative agreement with in operando imaging experiments in * equal contributions † bazant@mit.edu single-crystalline nanoparticles, given measured reaction rate constants.arXiv:1802.05847v1 [cond-mat.mtrl-sci]