Current research activities are strongly entangled with a “source concept” paradigm, in which all involved field quantities as well as engineering metrics are described solely in terms of electromagnetic sources. The idea of the source concept is sketched in the following figure.
CEM team collaborates with other department’s teams. Theses of various topics are supervised by members of the group.
Core Members of the Group:
Doctoral Students 
Students 
Visiting Reseachers 
Past Members 
Vit Losenicky
Michal Masek

Martin Strambach
Jakub Liska
Jonas Tucek
Vojtech Neuman

Lamyae Akrou

Projects
 Virtual Prototyping and Validation of Electromagnetic Systems (20182021, TH04010373)
 Fundamental Bounds on Electromagnetic Radiation and Scattering Phenomena and Associated Realizable Subforms (20192021, GA1906049S)
 Tools for Synthesis of Antennas and Sensors (20142017, TA04010457)
 Source Concept of Electrically Small Antenna Synthesis (20152017, GA1510280Y)
Active External Collaborations
 Brno University of Technology, Czech Republic (Petr Kadlec)
 University of Lund, Sweden (Mats Gustafson, Casimir Ehrenborg)
 University of Santa Clara, USA (Kurt Schab)
 KTH Royal Institute of Technology, Sweden (Lars Johnsson)
 Catholic University of Leuven, Belgium (Xuezhi Zhang, Guy Vandenbosch)
 Université Nice Sophia Antipolis (Fabien Ferrero)
 MECAS ESI s.r.o., Czech Republic (Jaroslav Rymus)
Series of Internal Seminars
 06012020, MATLAB on RCI (J. Tucek)
 30032020, GIT, (M. Strambach)
 20042020, AToM2TikZ (V. Neuman)
 04052020, Symmetries in AToM (M. Masek)
Journal Papers (Last Three Years)
Current and Antenna Optimization

AbstractThe optimal spectral efficiency of MIMO antennas in an ideal lineofsight 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 Qfactor a semianalytical 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 subregions 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. 

AbstractFundamental bounds on antenna gain are found via convex optimization of the current density in a prescribed region. Various constraints are considered, including selfresonance 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 Qfactor. The paper is accompanied with a series of examples practically demonstrating the relevance of the theoretical framework and entirely spanning wide range of material parameters and electrical sizes used in antenna technology. Presented results are analyzed from a perspective of effectively radiating modes. In contrast to a common approach utilizing spherical 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. 

AbstractA recently introduced technique of topology sensitivity in method of moments is extended by the possibility of adding degreesoffreedom (reconstruct) into underlying structure. The algebraic formulation is inversionfree, 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 learningrelated applications. 

AbstractConsiderable 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 radiation Qfactor (maximizing single resonance fractional bandwidth), but are far from reaching the associated physical bounds on efficiency. The relative performance of other canonical antenna designs is compared in similar ways, and the quantitative results are connected to intuitions from small antenna design, physical bounds, and matching network design. 

AbstractA 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 one hundred times faster than contemporary approaches, which allows the investigation of the sensitivity and the modification of shapes in realtime. The method is compared with known approaches and its validity and effectiveness is verified using a series of examples. The procedure can be implemented in uptodate EM simulators in a straightforward manner. It is shown that the iterative repetition of the topology sensitivity evaluation can be used for gradientbased topology synthesis. This technique can also be employed as a local step in global optimizers. 

AbstractThe tradeoff between radiation efficiency and antenna bandwidth, expressed in terms of Qfactor, 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 selfresonant current as opposed to one tuned by an external reactance. The resulting Paretooptimal sets reveal the relative cost of valuing low Qfactor over high efficiency, the cost in efficiency to require a selfresonant current, and other insights. 

AbstractExisting optimization methods are used to calculate the upperbounds on radiation efficiency with and without the constraint of selfresonance. These bounds are used for the design and assessment of small electricdipoletype antennas. We demonstrate that the assumption of lossless, lumped, 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. 

AbstractThe optimal currents on arbitrarily shaped radiators with respect to the minimum quality factor Q 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, socalled 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 nearfields. The obtained numerical results are compatible with previous results based on composition of modal currents, convex optimization, and quasistatic approximations; however, using the methodology in this paper, the class of solvable problems is significantly extended. 

AbstractAn 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 $eta$ of canonical shapes with conduction losses explicitly included. The RaoWiltonGlisson 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. 

AbstractThis work 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 Vandenboschs relations yielding stored electromagnetic energy as a function of the source current density. All calculations are based on a matrix representation of the integrodifferential 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. 
Characteristic Modes

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

AbstractA new method to improve the accuracy of characteristic modes 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 moments solution of the electric field integral equation to improve the numerical range of the real part of the impedance matrix, R, that determines the number of obtainable modes from characteristic modes decomposition. Computation speed of the R matrix and characteristic modes are improved. The method can easily be integrated in existing method of moments solvers. Several structures are investigated illustrating the improved accuracy and performance of the new method. 

AbstractCharacteristic modes of a spherical shell are found analytically as spherical harmonics normalized to radiate unitary power and to fulfill specific boundary conditions. The presented closedform formulas lead to a proposal of precise synthetic benchmarks which can be utilized to validate the method of moments matrix or performance of characteristic mode decomposition. Dependence on the mesh size, electrical size and other parameters can systematically be studied, including the performance of various mode tracking algorithms. A noticeable advantage is the independence on feeding models. Both theoretical and numerical aspects of characteristic mode decomposition are discussed and illustrated by examples. The performance of stateoftheart 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 higherorder modes. Surprisingly poor performance of the tracking routines is observed notwithstanding the recent ambitious development. 

AbstractAspects 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. 
Stored Energy and Quality Factor Q

AbstractThough commonly used to calculate Qfactor 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. 

AbstractAmbiguities in the definition of stored energy within distributed or radiating electromagnetic systems motivate the discussion of the welldefined 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 phaseshift Darlington’s synthesis procedure. 

AbstractThe commented paper [1] claims to provide a new expression for an energy stored around a general radiator. The major purpose of this comment is to show that the claim is unjustified. Alongside with this issue, it is pointed out that some of the core formulas of [1] are not completely correct, and that their correct form has in fact been derived elsewhere, though for the purpose of evaluating the quality factor QZ and not the stored energies. 

AbstractThis 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 semidefiniteness and coordinate independence, i.e. two key properties actually not met by the contemporary approaches. The proposed technique is compared with an uptodate frequency domain method that is extensively used in practice. Both concepts are discussed and compared on the basis of examples of varying complexity. 
Artificial Materials

AbstractA fullwave numerical scheme of polarizability (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 polarizability 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. 