Damghan University Press Analytical and Numerical Solutions for Nonlinear Equations 3060-785X 10 2 2025 12 30 Nonlinear Eigenvalue Methods for Quantifying Quantum Entanglement 300 313 2037 10.22128/ansne.2025.3141.1175 EN Abrar Ahmed Naqash Department of Physics, National Institute of Technology Srinagar, Jammu and Kashmir, 190006, India 0000-0003-3891-4740 Fardeen Ahmad Sofi Department of Physics, University of Kashmir, Srinagar 190006, India 0009-0004-0792-5353 Mohammad Haris Khan Department of Physics, University of Kashmir, Srinagar 190006, India 0009-0007-8927-484X Sundus Abdi Department of Physics, University of Toronto, Toronto, Ontario M5S 1A7, Canada 0000-0002-7810-653X Journal Article 2025 11 11 We present a hybrid analytical numerical method to evaluate the geometric measure of entanglement for pure multipartite states by formulating the closest separable state problem as a coupled nonlinear eigenvalue condition. We develop a hybrid analytical numerical framework in which a formal perturbative linearization around a reference product state is combined with an iterative fixed-point scheme. The approach combines a Gauss-Seidel block fixed-point iteration with a controlled first order linearization about a stationary reference product state. The perturbative analysis provides local structural insight and initialization guidance, while the iterative method yields accurate numerical estimates of the geometric measure of entanglement. We make explicit and prove an equal multiplier stationarity identity showing that, at an optimum, all block Lagrange multipliers coincide and are fixed by the optimal fidelity to the target state. A normalization preserving linearization is obtained by projecting onto local tangent spaces, which yields an explicit first order correction and a corresponding scalar shift in the effective eigenvalue. We further establish a monotonic block ascent property: the overlap with the target state increases at every update, remains bounded, and converges to a stationary value. For standard three qubit benchmarks, the hybrid solver converges smoothly and reproduces the known exact optima for the GHZ and W states. https://ansne.du.ac.ir/article_2037_26abd9c4509320a2089c51fab82388d6.pdf