Abstract: A novel implicit frequency-domain finite volume method for directly solving the scattered-field formulation of frequency-domain Maxwell’s equations is proposed to obtain the precise spatial distribution of electromagnetic fields and scattering characteristics of targets under harmonic wave incidence. This method is an extension of the time-domain finite volume method from temporal to spectral domains, transforming solution variables from real-valued 4D spacetime coordinates into complex-valued 3D spatial representations while enabling applications of various steady-state acceleration techniques. It is a semi-discrete scheme, whose solution process involves two steps: steady-state virtual time integration and implicit LU-SGS (lower upper symmetric Gauss Siedel) solving of the spatial flux residual. Calculations with Large Courant number show that this implicit algorithm has the advantage of unconditional stability. Numerical examples demonstrate that the frequency-domain finite volume method exhibits comparable computational accuracy to methods of the moments and Mie series solutions. It indicates that this method has extensive applicability across various scenarios.