Design of modern antennas relies-for reliability reasons-on full-wave electromagnetic simulation tools. In addition, increasingly stringent specifications pertaining to electrical and field performance, growing complexity of antenna topologies, along with the necessity for handling multiple objectives, make numerical optimization of antenna geometry parameters a highly recommended design procedure. Conventional algorithms, particularly global ones, entail often-unmanageable computational costs, so alternative approaches are needed. This work proposes a novel method for cost-efficient globalized design optimization of multi-band antennas incorporating the response feature technology into the trust-region framework. It allows for unequivocal allocation of the antenna resonances even for poor initial designs, where conventional local algorithms fail. Furthermore, the algorithm is accelerated by means of Jacobian variability tracking, which reduces the number of expensive finite-differentiation updates. Two real-world antenna design cases are used for demonstration purposes. The optimization cost is comparable to that of local routines while ensuring nearly global search capabilities. INDEX TERMS Antenna design, input characteristics, EM-driven design, trust-region methods, response features.