A simple characterization of the solvability of power flow equations is of great importance in the monitoring, control, and protection of power systems. In this paper, we introduce a sufficient condition for power flow Jacobian nonsingularity. We show that this condition is second-order conic representable when load powers are fixed. Through the incorporation of the sufficient condition, we propose a voltage stability-constrained optimal power flow (VSC-OPF) formulation as a second-order cone program (SOCP). An approximate model is introduced to improve the scalability of the formulation to larger systems. Extensive computation results on MATPOWER and NESTA instances confirm the effectiveness and efficiency of the formulation.