CEM group concentrates on theoretical research in classical electromagnetic theory with an emphasis on the computational aspects. Currently developed topics include radiation properties of electrically small radiators and scatterers, fundamental bounds, and inverse design. Members of the group have for many years been also involved in classical antenna theory and field propagation in artificial materials.

Current research activities are strongly entangled with a “source concept” paradigm, in which field observables operate over electromagnetic sources, which are the system variables often represented in finite-dimensional bases. The emphasis is placed on the development of original computational codes based on this paradigm.

CEM group members offer variety of topics for final projects and maintain extensive library of books related to classical and quantum electromagnetism, antenna theory, and numerical modeling. The team periodically participates in workshops. We like MATLAB, LaTeX, Beamer, and Ti*k*Z.

### Miloslav Čapek

Associate professor

- Antenna theory
- Numerical methods
- Optimization
- MATLAB, LaTeX, Ti
*k*Z

### Lukáš Jelínek

Associate professor

- Electromagnetic field theory
- Antenna theory
- Artificial electromagnetic materials

### Enrique Moreno Perez

Postdoctoral researcher

- Numerical methods
- Theory of electromagnetic field
- Multiphysics

### Martin Štrambach

Student

- Machine learning
- GPU computing

Kurt Schab

- Visit March-June 2022
- Electromagnetic theory, Numerical methods, Optimization.

Lamye Akrou

- Visit 2019
- Precise characteristic mode decomposition.

- Fundamental Bounds in EM
- Optimal Inverse Design
- Small Antennas
- Modal Decompositions
- Method of Moments

Fundamental bounds determine the best attainable values of physical metrics and are typically evaluated using tools of convex optimization and matrix operators. The fundamental bounds delimit performance of a hypothetical device. In line with the “source concept”, fundamental bounds result in optimal current densities defined in a prescribed region. Formulation of fundamental bounds commonly include physically motivated constraints such as enforcement of self-resonance or complex power balance. Among others, bounds are for example known on minimum antenna Q-factor, maximum antenna gain, minimum dissipation factor, minimum or maximum scattering cross sections.

Quite a few constraints can be used, *e.g.*, a constraint on self-resonance, on complex power balance, or on only partly controllable region. All bounds can be transformed into multi-criteria form describing the mutual trade-offs.

*Figure: Optimal current densities for an L-shape plate (electric size ka = 1/2). The sequence corresponds to a specific trade-off point between dissipation factor and Q-factor, maximal radiation efficiency (externally tuned and self-resonant), and minimal Q-factor, respectively.*

In the current state of the art, a direct shape synthesis is not possible, the main obstacle being the combinatorial nature of an associated optimization problem which is non-polynomial in time. On the other hand, there is a strong evidence that a skilled designer can provide designs with performance close to fundamental bounds...

How close can we go with clever, albeit brute-force techniques? Can we go closer? Is it possible to improve the design performance with the knowledge of the first differences or utilization of machine learning? To address these questions, a novel memetic framework combining local and global optimization routines is developed, combining the advantages of an adjoint formulation of topology optimization and of an evolutional algorithm. Various geometry- and topology-based metrics like shape regularity are being incorporated as well.

*Figure: One run of topology optimization based on the exact re-analysis.*

Electrically small antennas are ubiquitous, yet their true importance is still about to be revealed with the upcoming Internet of Things (IoT). Other technologies relying on small radiators such as implantable antennas, directional nanoscatterers, or field concentrators and absorbers are developed too. For all of them, the performance nontrivially depends on shape of the structure, material distribution, and excitation.

Various metrics relevant mainly in the electrically small regime are studied and formulated in terms of method of moments operators. A few examples are stored energy, ohmic losses, electric and magnetic moments, total active reflection coefficient and others.

*Figure: A meanderline situated in a radiansphere of radius a as a typical representant of planar electrically small antenna.*

The most salient features of many complicated EM phenomena can be revealed by proper modal decomposition, which also reduces the complexity of problem and often offers additional physical insight. Most commonly, the modal decomposition is achieved via eigenvalue problem the modes of which are used as a new basis for the engineering problem at hand.

For example, characteristic modes diagonalize the impedance matrix, result in orthogonal far-fields, and are thus excellent for a design of electrically small MIMO antennas. Other bases, like radiation modes, are perfect for reducing the numerical complexity of fundamental bounds evaluation via convex optimization routines.

*Figure: First two dominant modes on a rectangular plate.*

Method of moments is used to convert linear operator equations into linear algebraic equation systems by representing the solution in a suitable basis of expansion functions. The crucial step of identifying associated Green’s function makes this method ideal for open (radiating) problems since only the region containing sources is discretized.

Application of method of moments to integral equations leads to dense and relatively small (as compared the application to differential equations also known as finite element method) system matrix. Direct solvers are used to invert these matrices up to dimensions of several thousands. For bigger matrices, the efficient indirect solvers are available.

Conventionally, the inversion of a system matrix yields the contrast current flowing in the solution domain. This is, however, only one of many possibilities how to employ this method. Other techniques utilize the system matrix and its derivatives directly.

*Figure: Surface current density (3rd characteristic mode) on a triangularized obstacle.*

- Virtual Prototyping of EM Systems
- Fundamental Bounds
- AToM
- Source Concept

Fast and precise virtual prototyping of EM systems is a strong prerequisite for massive expansion of IoT devices and is one of the final goals of Industry 4.0. This project aims at the development of numerical tools capable to design and simulate wireless systems, to expand existing optimization routines, and to characterize the connectors and calibration kits.

This project has been supported by Technology Agency of the Czech Republic within the frame of project TH04010373, 2018-2021 ("Virtual Prototyping and Validation of Electromagnetic Systems").

*Figure: Parametric sweep of a double cylindrical helix.*

The optimality of passive electromagnetic structures is studied within this project. The main goal is to investigate the existence of and build up an understanding of the physical bounds on the primary physical quantities as well as system metrics related to radiation and scattering phenomena arising in wireless communication and power transfer. This project strives to fill a gap in the knowledge of classical electromagnetism as the investigation of the optimal source current distributions is, in many cases, absent.

This project has been supported by Czech Science Foundation within the frame of project GACR 19-06049S, 2019-2021 ("Fundamental Bounds on Electromagnetic Radiation and Scattering Phenomena and Associated Realizable Subforms").

*Figure: Pareto-optimal solutions with respect to three antenna parameters for three arrangements of dipole arrays.*

AToM (Antenna Toolbox for MATLAB) originated from project TA04010457 (2014-2017) of Technology Agency of the Czech Republic and is written in MATLAB as a semi-open code offering user-friendly operation through GUI or direct access to low level functions. The simulation core of AToM is based on Method of moments solution to field integral equations on wire, planar, and volumetric structures. Together with modal decomposition (characteristic modes), the source concept, feeding synthesis and powerful optimization package (FOPS optimizer), AToM presents unique tool for analysis and synthesis of antennas.

AToM also contains several add-ons, notable examples being the determination of fundamental bounds on antenna metrics, the design of optimal shapes via topology (shape) optimization and parametric strip generator. An extensive list of integrodifferential operators is available both in RWG MoM (2D surface code) and piecewise constant MoM (3D volumetric code).

Visit AToM webpage for more details.

Evaluating antenna and scattering characteristics solely by means of contrast currents flowing in a region occupied by the radiator opens a new paradigm called “source concept”, which also involves diverse methods of calculation, decomposition, and electric and magnetic currents modification. Within the method of moments, the source concept typically leads to observables, the current-quadratic forms of which provide the performance metrics and allows for their efficient analysis or for setting up their fundamental bounds. Automated search for sub-optimal designs via topology optimization is the latest application of the source concept.

A precise and efficient numerical implementation of the source concept is an inevitable part of its practically-oriented development. With respect to this, the MATLAB toolbox AToM is currently being developed at CTU in Prague.

- Posters
- Ti
*k*Z-ed Science - Other materials
- Season greeting cards
- Internal seminars

Author: J. Tucek

Presented on: January 06, 2020

Author: M. Strambach

Presented on: March 03, 2020

Author: V. Neuman

Presented on: April 20, 2020

AToM2TikZ is an extension for AToM package, which allows conversion of AToM plots into LaTeX/Ti*k*Z macros. Package offers additional functionality such as animations, shading, opacity settings and possibility to easily integrate another graphical objects. The presentation contains basic description of plot conversion, examples, how to use package and restrictions of current version of AToM2TikZ. The direct visibility of triangles algorithm is presented at the end of seminar.

Author: M. Masek

Presented on: May 04, 2020

Symmetries are an integral part of our everyday lives. The presentation describes the utilization of symmetries and application of point group theory into MoM framework. Their detailed implementation into Antenna Toolbox for Matlab (AToM) package is explained and thoroughly illustrated. Later on, application examples are presented. Namely: the solution to the modal tracking, block-diagonalizing of any arbitrary MoM operator, or design of orthogonal feeding schemes.

Authors: M. Capek, V. Neuman, and J. Tucek

Presented on: February 22, 2022

Authors: L. Jelinek and M. Capek

Presented on: February 22, 2022

Authors: M. Masek and M. Capek

Presented on: February 22, 2022

Authors: V. Losenicky, L. Jelinek, and M. Capek

Presented on: February 25, 2020

Authors: L. Jelinek and M. Capek

Presented on: February 25, 2020

Authors: M. Capek, *et al.*

Presented on: February 26, 2019

Authors: M. Capek, L. Jelinek, and V. Losenicky

Presented on: February 26, 2019

Authors: M. Capek, *et al.*

Presented on: February 26, 2019

Authors: M. Capek and L. Jelinek

Presented on: February 27, 2018

## Templates & LaTeX packages

Authors: M. Masek, M. Capek, and P. Chvojka

A departmental beamer template to create a presentations in LaTeX.

The template is available in the free use at Overleaf.

Authors: M. Masek and M. Capek

A departmental template to create a poster in LaTeX.

The template is available in the free use at Overleaf.

Author: M. Masek

The mtbpars package provides a parsing algorithm which is able to print MATLAB commands in LaTeX with the same style as they are written in the original software. Direct inline commands and loading the code from the .m files are included.

The package is still under development for some corner-cases, however, it works correctly for standard usage. Is in the free use at Overleaf.

## Games

Author: M. Masek

Can you find a pair of objects having the same symmetries, *i.e.*, belonging into the same point group?

Try it online or download print&play version.

May all your ports be perfectly matched in the entire year 2022.

May your joy and success be unbounded in 2021.

May the stars shine upon your road in 2020.

Let the left and right hand sides of your equations match for the entire year 2019.

Dominant characteristic mode of IFS pre-fractal of the 3rd iteration.

Author: M. Capek

Dimensions: 77.5 cm × 83.7 cm

Optimal excitation of an array of 10 thin-strip dipoles for maximum gain, *l=λ/2*, *d = λ/4*, *f = 1* GHz, copper (radiation pattern is in log scale).

Author: M. Capek

Dimensions: 94.5 cm × 284.4 cm

Bistatic RCS of an RPG-7 missile. Induced current density is depicted together with *φ = 0* cut of scattering radiation characteristic and values of the RCS on a Lebedev sphere of 110 equidistantly spaced samples. Radiation characteristic and values of the RCS are shown in logarithmic scale.

Author: M. Capek

Dimensions: 86.4 cm × 72.0 cm

Optimal current densities for various trade-offs between radiation Q-factor, *Q/η*, and normalized dissipation factor, *δ/R*_{s}. L-shape plate, *ka = 1/2*.

Author: M. Capek

Dimensions: 94.1 cm × 54.1 cm

Poynting streamlines and magnitude of Poynting vector for a plane wave impinging from left to a receiving short dipole (*L ≈ 0.03 λ*) made of lossy shell and with conjugate matched load placed in the middle. Two relevant perpendicular cut planes are shown.

Author: M. Capek

Dimensions: 92.0 cm × 67.0 cm

RF spectrum in terms of spectral density as measured on the roof of the CTU building in Dejvice Campus. The measured peaks are associated with various services.

Authors: M. Capek and J. Spacil

Dimensions: 258.0 cm × 78.0 cm

The arts hang on walls at the Department.

Casimir Ehrenborg

Guy Vandenbosch

Raúl Rodríguez Berral

Ville Viikari

- Fundamental Bounds
- Optimal Design
- Modal Decomposition
- Others

## Submitted:

Knowledge of the fundamental limitations on a magnetic trap for neutral particles is of paramount interest to designers as it allows for the rapid assessment of the feasibility of specific trap requirements or the quality of a given design. In this paper, performance limitations are defined for convexity of magnetic trapping potential and bias field using a local approximation in the trapping center. As an example, the fundamental bounds are computed for current supporting regions in the form of a spherical shell, a cylindrical region, and a box. A Pareto-optimal set considering both objectives is found and compared with known designs of the Baseball trap and Ioffe-Pritchard trap. The comparison reveals a significant gap in the performance of classical trap designs from the fundamental limitations. This indicates a possibility of improved trap designs and modern techniques of shape synthesis are applied in order to prove their existence. The topologically optimized traps perform almost two times better as compared to conventional designs. Last, but not least, the developed framework might serve as a prototype for the formulation of fundamental limitations on plasma confinement in a wider sense.

Available online at arXiv.

xx

## Published:

*Optics Express, vol. 29, pp. 24068 - 24082, Jul. 2021,*

**Editor's pick**Fundamental bounds on the performance of monochromatic scattering-cancellation and fieldzeroing cloaks made of prescribed linear passive materials occupying a predefined design region are formulated by projecting field quantities onto a sub-sectional basis and applying quadratically constrained quadratic programming. Formulations are numerically tested revealing key physical trends as well as advantages and disadvantages between the two classes of cloaks. Results show that the use of low-loss materials with high dielectric contrast affords the highest potential for effective cloaking.

*IEEE Transactions on Antennas and Propagation*, Apr. 2021

A problem of the erroneous duality gap caused by the presence of symmetries is solved in this paper utilizing point group theory. The optimization problems are first divided into two classes based on their predisposition to suffer from this deficiency. Then, the classical problem of Q-factor minimization is shown in an example where the erroneous duality gap is eliminated by combining solutions from orthogonal subspaces. Validity of this treatment is demonstrated in a series of subsequent examples of increasing complexity spanning the wide variety of optimization problems, namely minimum Qfactor, maximum antenna gain, minimum total active reflection coefficient, or maximum radiation efficiency with self-resonant constraint. They involve problems with algebraic and geometric multiplicities of the eigenmodes, and are completed by an example introducing the selective modification of modal currents falling into one of the symmetry conformal sub-spaces. The entire treatment is accompanied with a discussion of finite numerical precision, and mesh grid imperfections and their influence on the results. Finally, the robust and unified algorithm is proposed and discussed, including advanced topics such as the uniqueness of the optimal solutions, dependence on the number of constraints, or an interpretation of the qualitative difference between the two classes of the optimization problems.

*IEEE Transactions on Antennas and Propagation*, Feb. 2021

The optimal spectral efficiency of MIMO antennas in an ideal line-of-sight channel is investigated when bandwidth requirements are placed on the antenna. By posing the problem as a convex optimization problem restricted by the input port Q-factor a semi-analytical expression is formed for its solution. It is shown that this solution is solely dependent on energy modes of the antenna. These modes are compared to the characteristic modes and the ability to induce them through sub-regions of a plate is investigated. The position of these regions is also investigated when they are raised above the ground plane. Their performance is illustrated through spectral efficiency over Q, a quantity that is connected to the true capacity. It is demonstrated that the spatial diversity of the controlled regions correlates with the number of significant energy modes.

*Optics Express, vol. 28, pp. 36584 - 36599, Nov. 2020*

Trade-offs between feasible absorption and scattering cross sections of obstacles confined to an arbitrarily shaped volume are formulated as a multi-objective optimization problem solvable by Lagrangian-dual methods. Solutions to this optimization problem yield a Pareto-optimal set, the shape of which reveals the feasibility of achieving simultaneously extremal absorption and scattering. Two forms of the trade-off problems are considered involving both loss and reactive material parameters. Numerical comparisons between the derived multi-objective bounds and several classes of realized structures are made. Additionally, low-frequency (electrically small, long wavelength) limits are examined for certain special cases.

*IEEE Transactions on Antennas and Propagation*, Oct. 2020

A numerically effective description of the total active reflection coefficient and realized gain are studied for multi-port antennas. Material losses are fully considered. The description is based on operators represented in an entire-domain port-mode basis, *i.e.*, on matrices with favorably small dimensions. Optimal performance is investigated and conditions on optimal excitation and matching are derived. The solution to the combinatorial problem of optimal ports' placement and optimal feeding synthesis is also accomplished. Four examples of various complexity are numerically studied, demonstrating the advantages of the proposed method. The final formulas can easily be implemented in existing electromagnetic simulators using integral equation solver.

*New Journal of Physics*, vol. 22, 073013, Sept. 2020.

A general framework for determining fundamental bounds in nanophotonics is introduced in this paper. The theory is based on convex optimization of dual problems constructed from operators generated by electromagnetic integral equations. The optimized variable is a contrast current defined within a prescribed region of a given material constitutive relations. Two power conservation constraints analogous to optical theorem are utilized to tighten the bounds and to prescribe either losses or material properties. Thanks to the utilization of matrix rank-1 updates, modal decompositions, and model order reduction techniques, the optimization procedure is computationally efficient even for complicated scenarios. No dual gaps are observed. The method is well-suited to accommodate material anisotropy and inhomogeneity. To demonstrate the validity of the method, bounds on scattering, absorption, and extinction cross sections are derived first and evaluated for several canonical regions. The tightness of the bounds is verified by comparison to optimized spherical nanoparticles and shells. The next metric investigated is bi-directional scattering studied closely on a particular example of an electrically thin slab. Finally, the bounds are established for Purcell's factor and local field enhancement where a dimer is used as a practical example.

*IEEE Transactions on Antennas and Propagation*, vol. 67, no. 8, pp. 5282 – 5293, Aug. 2019.

Fundamental bounds on antenna gain are found via convex optimization of the current density in a prescribed region. Various constraints are considered, including self-resonance and only partial control of the current distribution. Derived formulas are valid for arbitrarily shaped radiators of a given conductivity. All the optimization tasks are reduced to eigenvalue problems, which are solved efficiently. The second part of the paper deals with superdirectivity and its associated minimal costs in efficiency and Q-factor. The paper is accompanied by a series of examples practically demonstrating the relevance of the theoretical framework and entirely spanning a wide range of material parameters and electrical sizes used in antenna technology. Presented results are analyzed from the perspective of effectively radiating modes. In contrast to a common approach utilizing modes, the radiating modes of a given body are directly evaluated and analyzed here. All crucial mathematical steps are reviewed in the appendices, including a series of important subroutines to be considered making it possible to reduce the computational burden associated with the evaluation of electrically large structures and structures of high conductivity.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 67, no. 4, pp. 2482-2493, Apr. 2019.

The tradeoff between radiation efficiency and antenna bandwidth, expressed in terms of Q-factor for small antennas, is formulated as a multiobjective optimization problem in current distributions of predefined support. Variants on the problem are constructed to demonstrate the consequences of requiring a self-resonant current as opposed to the one tuned by an external reactance. The tradeoffs are evaluated for sample problems and the resulting Pareto-optimal sets reveal the relative cost of valuing low radiation Q-factor over high efficiency, the cost in efficiency to require a self-resonant current, the effects of lossy parasitic loading, and other insights. Observations are drawn from the sample problems selected, though the tradeoff evaluation method is valid for studying arbitrary antenna geometries. In the examples considered here, we observe that small increases in Q-factor away from its lower bound allow for dramatic increases in efficiency toward its upper bound.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 66, no. 12, pp. 6716-6723, Dec. 2018.

Existing optimization methods are used to calculate the upper bounds on radiation efficiency with and without the constraint on self-resonance. These bounds are used for the design and assessment of small electric-dipole-type antennas. We demonstrate that the assumption of lossless, lumped, and external tuning skews the true nature of radiation efficiency bounds when practical material characteristics are used in the tuning network. A major result is that, when realistic (*e.g.*, finite conductivity) materials are used, small antenna systems exhibit dissipation factors which scale as (ka)^{-4} , rather than (ka)^{-2} as previously predicted under the assumption of lossless external tuning.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 66, no. 9, pp. 4948-4952, Sept. 2018.

Radiation efficiencies of modal current densities distributed on a spherical shell are evaluated in terms of dissipation factor. The presented approach is rigorous, yet simple and straightforward, leading to closed-form expressions. The same approach is utilized for a two-layered shell and the results are compared with other models existing in the literature. Discrepancies in this comparison are reported and reasons are analyzed. Finally, it is demonstrated that radiation efficiency potentially benefits from the use of internal volume, which contrasts with the case of the radiation Q-factor.

Published version available online at IEEEXplore.

*IEEE Transactions on Antennas and Propagation*, vol. 65, no. 8, pp. 4115-4123, Aug. 2017.

Optimal currents on arbitrarily shaped radiators with respect to the minimum quality factor are found using a simple and efficient procedure. The solution starts with a reformulation of the problem of minimizing quality factor Q as an alternative, so-called dual, problem. Taking advantage of modal decomposition and group theory, it is shown that the dual problem can easily be solved and always results in minimal quality factor Q. Moreover, the optimization procedure is generalized to minimize quality factor Q for embedded antennas, with respect to the arbitrarily weighted radiation patterns, or with prescribed magnitude of the electric and magnetic near fields. The obtained numerical results are compatible with previous results based on the composition of modal currents, convex optimization, and quasi-static approximations; however, using the methodology in this paper, the class of solvable problems is significantly extended.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 65, no. 1, pp. 329-341, Jan. 2017.

An optimization problem has been formulated to find a resonant current extremizing various antenna parameters. The method is presented on, but not limited to, particular cases of gain G, quality factor Q, gain to quality factor ratio G/Q, and radiation efficiency η of canonical shapes with conduction losses explicitly included. The Rao-Wilton-Glisson basis representation is used to simplify the underlying algebra while still allowing surface current regions of arbitrary shape to be treated. By switching to another basis generated by a specific eigenvalue problem, it is finally shown that the optimal current can, in principle, be found as a combination of a few eigenmodes. The presented method constitutes a general framework in which the antenna parameters, expressed as bilinear forms, can automatically be extremized.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

## Submitted:

Memetics for shape synthesis, introduced in Part 1, is examined on antenna design examples. It combines local and global techniques to accelerate convergence and to maintain robustness. Method-of-moments matrices are used to evaluate objective functions. By applying the Shermann-Morrison-Woodbury identity, the repetitively performed structural update is inversion-free yet full-wave in nature. The technique can easily be combined with additional features often required in practice, e.g., only a part of the structure is controllable or evaluation of an objective function is required in a subdomain only. The framework supports multi-frequency and multi-port optimization and offers many other advantages, such as an actual shape being known at every moment of the optimization. The performance of the method is assessed, including its convergence and computational cost.

Available online at arXiv.

A novel inverse design framework over a fixed discretization grid using method of moments for the shape optimization of radiation devices is proposed. The formulation is discrete and solves the original combinatorial problem. Fixed discretization provides an inversion-free, exact-reanalysis-based evaluation of the smallest topology perturbations. The fixed grid also allows to precalculate all expensive evaluations prior to the optimization. Local greedy-based structural updates are combined with global, heuristic steps, resulting in a rapidly converging and robust memetic algorithm applicable to various problems of electromagnetic inverse design. The method is robust and versatile, capable of evaluating a vast amount of possible shapes in a full-wave manner. It shows superb performance over widely used heuristic approaches, such as genetic algorithms or swarm optimization. Thanks to the discrete formulation, no filtering or rounding is needed as compared to the adjoint formulation of topology optimization.

Available online at arXiv.

## Published:

*International Journal of RF and Microwave Computer-Aided Engineering*, pp. 1-8, May 2020.

A single‐ and multi‐objective optimization package is presented and described in detail. It contains an ensemble of local and global optimization routines. Procedures controlling variable number of dimensions are implemented as well, which is a rare feature among optimization oriented packages. The package is provided as a MATLAB toolbox. It excels in versatility and extensibility, which is demonstrated on a series of examples covering classical electromagnetism and antenna design. It is taken for granted that defining parameters of the optimization method can be set prior to the simulation run. However, its effective performance can be changed during the optimization run thanks to the full control feature. Moreover, it opens new possibilities in merging various algorithms into hybrids, performing complex dynamic programming tasks, or exploiting third party software. These advantages render the package as a perfect tool to deal with nowadays challenging engineering tasks.

*IEEE Transactions on Antennas and Propagation*, vol. 18, no. 11, pp. 2311-2315, Nov. 2019.

A recently introduced technique of topology sensitivity in method of moments is extended by the possibility of adding degrees of freedom (reconstruct) into the underlying structure. The algebraic formulation is inversion-free, suitable for parallelization, and scales favorably with the number of unknowns. The reconstruction completes the nearest neighbors procedure for an evaluation of the smallest shape perturbation. The performance of the method is studied with a greedy search over a Hamming graph representing the structure in which initial positions are chosen from a random set. The method is shown to be an effective data mining tool for machine learning-related applications.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 61, no. 4, pp. 19-29, Aug. 2019.

Considerable time is often spent optimizing antennas to meet specific design metrics. Rarely, however, are the resulting antenna designs compared to rigorous physical bounds on those metrics. Here, we study the performance of optimized planar meander line antennas with respect to such bounds. Results show that these simple structures meet the lower bound on the radiation quality factor (Q-factor) (maximizing single-resonance fractional bandwidth) but are far from reaching the associated physical bounds for efficiency. The relative performance of other canonical antenna designs is comparable in similar ways, and the quantitative results are connected to intuitions from small antenna design, physical bounds, and matching network design.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 67, no. 6, pp. 3889-3901, Jun. 2019.

A method evaluating the sensitivity of a given parameter to topological changes is proposed within the method of moments paradigm. The basis functions are used as degrees of freedom which, when compared to the classical pixeling technique, provide important advantages, one of them being impedance matrix inversion free evaluation of the sensitivity. The devised procedure utilizes port modes and their superposition which, together with only a single evaluation of all matrix operators, leads to a computationally effective procedure. The proposed method is approximately 100 times faster than the contemporary approaches, which allows the investigation of the sensitivity and the modification of shapes in real time. The method is compared with the known approaches and its validity and effectiveness are verified using a series of examples. The procedure can be implemented in up-to-date electromagnetic (EM) simulators in a straightforward manner. It is shown that the iterative repetition of the topology sensitivity evaluation can be used for gradient-based topology synthesis. This technique can also be employed as a local step in global optimizers.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

## Submitted:

An iterative algorithm is adopted to construct approximate representations of matrices describing the scattering properties of arbitrary objects. The method is based on the implicit evaluation of scattering responses from iteratively generated excitations. The method does not require explicit knowledge of any system matrices (e.g., stiffness or impedance matrices) and is well-suited for use with matrix-free and iterative full-wave solvers, such as FDTD, FEM, and MLFMA. The proposed method allows for significant speed-up compared to the direct construction of a full transition matrix or scattering dyadic. The method is applied to the characteristic mode decomposition of arbitrarily shaped obstacles of arbitrary material distribution. Examples demonstrating the speed-up and complexity of the algorithm are studied with several commercial software packages.

Available online at arXiv.

## Published:

*IEEE Transactions on Antennas and Propagation*, 2022

This is the second component of a two-part paper dealing with a unification of characteristic mode decomposition. This second part addresses modal tracking and losses and presents several numerical examples for both surface- and volume-based method-of-moment formulations. A new tracking algorithm based on algebraic properties of the transition matrix is developed, achieving excellent precision and requiring a very low number of frequency samples as compared to procedures previously reported in the literature. The transition matrix is further utilized to show that characteristic mode decomposition of lossy objects fails to deliver orthogonal far fields and to demonstrate how characteristic modes can be evaluated using the finite element method.

*IEEE Transactions on Antennas and Propagation*, 2022

A unification of characteristic mode decomposition for all method-of-moment formulations of field integral equations describing free-space scattering is derived. The work is based on an algebraic link between impedance and transition matrices, the latter of which was used in early definitions of characteristic modes and is uniquely defined for all scattering scenarios. This also makes it possible to extend the known application domain of characteristic mode decomposition to any other frequency-domain solver capable of generating transition matrices, such as finite difference or finite element methods. The formulation of characteristic modes using a transition matrix allows for the decomposition of induced currents and scattered fields from arbitrarily shaped objects, providing high numerical dynamics and increased stability, removing the issue of spurious modes, and offering good control of convergence. This first part of a two-part paper introduces the entire theory, extensively discusses its properties and offers its basic numerical validation.

*IEEE Transactions on Antennas and Propagation*, 2022

Characteristic modes are formulated using the scattering dyadic, which maps incident plane waves to scattered far fields generated by an object of arbitrary material composition. Numerical construction of the scattering dyadic using arbitrary full-wave electromagnetic solvers is demonstrated in examples involving a variety of dielectric and magnetic materials. Wrapper functions for computing characteristic modes in method-ofmoments, finite-difference time domain, and finite element solvers are provided as supplementary material.

*IEEE Transactions on Antennas and Propagation*, Feb. 2021

A technique of designing antenna excitation realizing orthogonal free-space channels is presented. It is shown that a symmetric antenna geometry is required in order to achieve orthogonality with respect to all physical quantities. A maximal number of reachable uncorrelated channels and a minimal number of feeders required to excite them are rigorously determined from the knowledge of an antenna's symmetries. The number of channels and number of feeders are summarized for the commonly used point groups (a rectangle, a square, etc.). The theory is applied to an example of a rectangular rim where the positions of feeders providing the best total active reflection coefficient, an important metric in multi-port systems, are determined. The described technique can easily be implemented in existing integral solvers.

*IEEE Transactions on Antennas and Propagation*, vol. 64, no. 2, pp. 23-31, Apr. 2022.

Nearly all practical applications of the theory of characteristic modes (CMs) involve the use of computational tools. Here in Paper 2 of this Series on CMs, we review the general transformations that move CMs from a continuous theoretical framework to a discrete representation compatible with numerical methods. We also review several key topics related to computational CMs, including modal tracking, dynamic range, code validation, electrically large problems, and non-PEC techniques.

*IEEE Transactions on Antennas and Propagation*, vol. 64, no. 2, pp. 14-22, Apr. 2022.

Over the past decade, characteristic mode analysis (CMA) research has grown from a niche topic to a mainstream topic, warranting a tutorial-style special issue to survey the signiﬁcant progress that has been made in this ﬁeld. In this introductory article (PAPER 1), the focus is on providing the big picture. We start with a simple description of characteristic modes. Next, we examine the trends in this ﬁeld, followed by providing further insights into CMA’s historical development. We will also address common myths surrounding the subject. Then, leaving the detailed coverage of major topics to the following papers, we summarize recent applications of CMA in scattering and other emerging topics. Finally, we conclude with some future perspectives on this ﬁeld.

*IEEE Transactions on Antennas and Propagation*, vol. 68, no. 2, pp. 927-937, Feb. 2020.

Issues in modal tracking in the presence of crossings and crossing avoidances between eigenvalue traces are solved via the theory of point groups. The von Neumann-Wigner theorem is used as a key factor in predictively determining mode behavior over arbitrary frequency ranges. The implementation and capabilities of the proposed procedure are demonstrated using characteristic mode decomposition as a motivating example. The procedure is, nevertheless, general and can be applied to an arbitrarily parametrized eigenvalue problem. A treatment of modal degeneracies is included and several examples are presented to illustrate modal tracking improvements and the immediate consequences of improper modal tracking. An approach leveraging a symmetry-adapted basis to accelerate computation is also discussed. A relationship between geometrical and physical symmetries is demonstrated on a practical example.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 68, no. 2, pp. 927-937, Dec. 2018.

A new method to improve the accuracy and efficiency of characteristic mode (CM) decomposition for perfectly conducting bodies is presented. The method uses the expansion of the Green dyadic in spherical vector waves. This expansion is utilized in the method-of-moment (MoM) solution of the electric-field integral equation to factorize the real part of the impedance matrix. The factorization is then employed in the computation of CMs, which improves the accuracy as well as the computational speed. An additional benefit is a rapid computation of far fields. The method can easily be integrated into existing MoM solvers. Several structures are investigated, illustrating the improved accuracy and performance of the new method.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 65, no. 8, pp. 4134-4145, Aug. 2017.

Characteristic modes (CMs) of a spherical shell are found analytically as spherical harmonics normalized to radiate unitary power and to fulfill specific boundary conditions. The presented closed-form formulas lead to a proposal of precise synthetic benchmarks that can be utilized to validate the method-of-moments matrix or performance of CM decomposition. Dependence on the mesh size, electrical size, and other parameters can systematically be studied, including the performance of various mode tracking algorithms. A notable advantage is the independence on feeding models. Both theoretical and numerical aspects of CM decomposition are discussed and illustrated by examples. The performance of state-of-the-art commercial simulators and academic packages having been investigated, we can conclude that all contemporary implementations are capable of identifying the first dominant modes while having severe difficulties with higher order modes. Surprisingly poor performance of the tracking routines is observed notwithstanding the recent ambitious development.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 65, no. 2, pp. 713-720, Feb. 2017.

Aspects of the theory of characteristic modes (CMs), based on their variational formulation, are presented and an explicit form of a related functional, involving only currents in a spatial domain, is derived. The new formulation leads to deeper insight into the modal behavior of radiating structures as demonstrated by a detailed analysis of three canonical structures: a dipole, an array of two dipoles and a loop, cylinder, and a sphere. It is demonstrated that knowledge of the analytical functional can be utilized to solve important problems related to the theory of CM decomposition such as the resonance of inductive modes or the benchmarking of method of moments code.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 64, no. 12, pp. 5230-5242, Dec. 2016.

This paper describes a powerful, yet simple, procedure how to acquire a current approaching the lower bound of quality factor Q. This optimal current can be determined for an arbitrarily shaped electrically small radiator made of a perfect conductor. Quality factor Q is evaluated by Vandenbosch's relations yielding stored electromagnetic energy as a function of the source current density. All calculations are based on a matrix representation of the integro-differential operators. This approach simplifies the entire development and results in a straightforward numerical evaluation. The optimal current is represented in a basis of modal currents suitable for solving the optimization problem so that the minimum is approached by either one mode tuned to the resonance, or, by two properly combined modes. An overview of which modes should be selected and how they should be combined is provided and results concerning rectangular plate, spherical shell, capped dipole antenna, and fractal shapes of varying geometrical complexity are presented. The reduction of quality factor Q and the G/ Q ratio are studied and, thanks to the modal decomposition, the physical interpretation of the results is discussed in conjunction with the limitations of the proposed procedure.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

## Pre-prints:

Three contradictory but state-of-the-art concepts for defining and evaluating stored electromagnetic energy are treated in this communication, and are collated with the widely accepted definition of stored energy, which is the total energy minus the radiated energy. All three concepts are compared, and the results are discussed on an example of a dominant spherical mode, which is known to yield dissimilar results for the concepts dealt with here. It is shown that various definitions of stored energy density immanently imply diverse meanings of the term "radiation".

Available online at arXiv.

Ambiguities in the definition of stored energy within distributed or radiating electromagnetic systems motivate the discussion of the well-defined concept of recoverable energy. This concept is commonly overlooked by the community and the purpose of this communication is to recall its existence and to discuss its relationship to fractional bandwidth. Using a rational function approximation of a system's input impedance, the recoverable energy of lumped and radiating systems is calculated in closed form and is related to stored energy and fractional bandwidth. Lumped circuits are also used to demonstrate the relationship between recoverable energy and the energy stored within equivalent circuits produced by the minimum phase-shift Darlington's synthesis procedure.

Available online at arXiv.

## Published:

*IEEE Transactions on Antennas and Propagation*, Dec. 2021

Hybrid computational schemes combining the advantages of a method of moments formulation of a field integral equation and T-matrix method are developed in this paper. The hybrid methods are particularly efficient when describing the interaction of electrically small complex objects and electrically large objects of canonical shapes such as, spherical multi-layered bodies where the T-matrix method is reduced to the Mie series making the method an interesting alternative in the design of implantable antennas or exposure evaluations. Method performance is tested on a spherical multi-layer model of the human head. Along with the hybrid method, an evaluation of the transition matrix of an arbitrarily shaped object is presented and the characteristic mode decomposition is performed, exhibiting fourfold numerical precision as compared to conventional approaches.

*IEEE Access*, vol. 6, pp. 10553-10568, Feb. 2018.

Though commonly used to calculate Q-factor and fractional bandwidth, the energy stored by radiating systems (antennas) is a subtle and challenging concept that has perplexed researchers for over half a century. Here, the obstacles in defining and calculating stored energy in general electromagnetic systems are presented from first principles as well as using demonstrative examples from electrostatics, circuits, and radiating systems. Along the way, the concept of unobservable energy is introduced to formalize such challenges. Existing methods of defining stored energy in radiating systems are then reviewed in a framework based on technical commonalities rather than chronological order. Equivalences between some methods under common assumptions are highlighted, along with the strengths, weaknesses, and unique applications of certain techniques. Numerical examples are provided to compare the relative margin between methods on several radiating structures.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IET Microwaves, Antennas and Propagation*, vol. 11, no. 6, pp. 852-858, May 2017.

A full-wave numerical scheme of polarisability (polarisability) tensors evaluation is presented. The method accepts highly conducting bodies of arbitrary shape and explicitly accounts for the radiation as well as ohmic losses. The method is verified on canonical bodies with known polarisability tensors, such as a sphere and a cube, as well as on realistic scatterers. The theoretical developments are followed by a freely available code whose sole user input is the triangular mesh covering the surface of the body under consideration.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*IEEE Transactions on Antennas and Propagation*, vol. 64, no. 10, pp. 4575-4576, Oct. 2016.

This paper discusses the methods for evaluating the stored electromagnetic energies and the radiation Q for an arbitrary lossless antenna. New expressions for the stored electromagnetic energies are derived by using the Poynting theorem in the complex frequency domain, and they are compared with previous theory and are validated by numerical examples. The minimization of radiation Q for small antenna is also investigated. There exists an optimal current distribution that minimizes the radiation Q for specified small antenna geometry. The optimized Q and the optimal current distribution for small antenna may be determined by solving a generalized eigenvalue equation obtained from the Rayleigh quotient for the radiation Q.

Published version available online at IEEEXplore.

Final manuscript available online at arXiv.

*Proceedings of Royal Society A*, vol. 472, no. 2188, pp. 1-16, Apr. 2016.

This paper deals with the old yet unsolved problem of defining and evaluating the stored electromagnetic energy—a quantity essential for calculating the quality factor, which reflects the intrinsic bandwidth of the considered electromagnetic system. A novel paradigm is proposed to determine the stored energy in the time domain leading to the method, which exhibits positive semi-definiteness and coordinate independence, i.e. two key properties actually not met by the contemporary approaches. The proposed technique is compared with an up-to-date frequency domain method that is extensively used in practice. Both concepts are discussed and compared on the basis of examples of varying complexity.

Published version available online at Proceedings of Royal Society A.

Final manuscript available online at arXiv.

*Radioengineering*(Feature Article), vol. 2016, no. 1, pp. 12-17, Apr. 2016.

Antenna parameters particularly relevant to electrically small antenna design are reviewed in this paper. Source current definitions are accentuated leading to the introduction of the source concept which advantageously utilize only spatially bounded quantities. The framework of the source concept incorporates powerful techniques such as structural and modal decomposition, operator’s inversion and current optimization, thus opening new, challenging possibilities for antenna design, analysis and synthesis.

Published version available online at Radioengineering.