Photonic band gaps (PBG), photonic analogues of electronic semiconductor band gaps, have attracted much attention recently because of numerous potential applications in communications and computing. Aközbek and John (Phys. Rev. E, 57, 2287 (1998).) developed a variational model of such band gaps, using action functionals, where solitary waves are expanded in terms of a finite orthonormal basis. These expansions to finite order N converged to solitary waves. The nonlinear polynomial equations for the coefficients in the expansions, have nonunique solutions. Our paper, makes a study of the multiplicity of the solutions for one-dimensional photonic band-gap structures. It is found that the nonuniqueness grows dramatically with the order of the expansion N. We use homotopy, which continuously deforms the solutions of exactly solvable systems, into the solutions of the systems to be solved with new results in numeric algebraic geometry, such that all solutions are determined. We used Maple 7 to obtain the polynomial equations for the variational coefficients, extending Aközbek and John's approach. A homotopy-based package PHCpack was used to solve the systems for N ≤ 4 and a linearization-extrapolation method was developed to find real solutions for N ≥ 5. The results are compared with the exact soliton solutions and their convergence behavior is discussed. The interplay of geometrical, topological and variational methods is seen in these interesting physical band-gap structures. PACS Nos.: 42.65.Tg, 42.70.Qs, 02.30.Xx, 02.70.Wz