Many computational methods have been developed and used for modeling, understanding, and tailoring extreme optical effects at the nanoscale. Among them, this review focuses on the integral equation-based methods: within the local response limit, a potential-based boundary integral equation (BIE) formalism and a field-based volume integral equation (VIE) formalism; within the nonlocal hydrodynamic model, a potential-based BIE formalism. These formalisms are derived from macroscopic electrodynamics (together with appropriate constitutive relations). The derivations are based on three pillars: the Green function, the field relation(s) (for the VIE formalism, the incident-scattered-total field relation; for the BIE formalism, the interface conditions connecting the fields at two sides of the interface), and the field equivalence principle (for the VIE formalism, the volume equivalence principle; for the BIE formalism, the Huygens principle). By applying the method of moments (MoM) algorithm, the derived integral equations are converted into matrix equations, with possible problems in the implementation being discussed. Levels of solutions, including the eigenmode and natural mode solutions, and group representation theory are introduced as powerful post-processing steps. Many examples are shown to demonstrate the effectiveness of the reviewed algorithms.developed for microwave frequency applications, for example, antennas, radars, etc., and is crucial for the design of modern telecommunication systems, for example, smartphones, WIFI networks, etc. Hence, it is undoubtedly playing an essential role in forging modern society.Prompted by the grand leap in nanotechnology over the last decades, the classical techniques in CEM have been applied to smaller and smaller scales and have significantly assisted in the understanding of the EM wave-matter interaction in the micrometric and nanometric worlds. One of the main applications is the research domain of plasmonics where the interaction of light with collective oscillations of free electrons in subwavelength nanostructured metallic objects is studied. [1] The research on plasmonics has discovered many extreme optical effects at the nanoscale. [2][3][4] These effects are of fundamental importance for future technologies, for example, cancer therapy, single molecule sensing, assisting chemical reactions, efficient light sources, etc., which span all the domains in physical sciences, from biology over chemistry to physics. Theoretically, the interaction of light with metal-based nanostructures can be abstracted as a scattering problem (see Figure 1). In the problem, there is an external source, for example, externally imposed current and charge sources ρ(r) and j(r), occupying a volume in the space. This source radiates and generates an incident field propagating freely in the space until it meets an inhomogeneity, that is, the scatterer. The scatterer reflects and deflects the incident field and induces the scattered field which together with the incident field forms the tot...