Abstract-We propose eP lace-M S, an electrostatics based placement algorithm for large-scale mixed-size circuits. ePlace-MS is generalized, flat, analytic and nonlinear. The density modeling method eDensity is extended to handle the mixedsize placement. We conduct detailed analysis on the correctness of the gradient formulation and the numerical solution, as well as the rationale of direct-current removal and the advantages over prior density functions. Nesterov's method is used as the nonlinear solver, which shows high yet stable performance over mixed-size circuits. The steplength as the inverse of Lipschitz constant of the gradient function, while we develop a backtracking method to prevent overestimation. An approximated nonlinear preconditioner is developed to minimize the topological and physical differences between large macros and standard cells. Besides, we devise a simulated annealer to legalize the layout of macros and use a second-phase global placement to reoptimize the standard cell layout. All the above innovations are integrated into our mixed-size placement prototype ePlace-MS, which outperforms all the related works in literature with better quality and efficiency. Compared to the leading-edge mixed-size placer NTUplace3 [13], ePlace-MS produces up to 22.98% and on average 8.22% shorter wirelength over all the sixteen modern mixed-size (MMS) benchmark circuits with the same runtime.