Multi-frequency experiment on determining the parameters of a dielectric layer in a waveguide Elena A. Sheinaa and Yury V. Shestopalovb a Lomonosov Moscow State University, Leninskiye Gory, 1-52, Moscow, 119991, Russia b University of Gävle, Kungsbäcksvägen 47, Gävle, SE-801 76, Sweden Abstract The inverse problem of reconstructing the real permittivity of a plane-parallel layer in a perfectly conducting rectangular waveguide from experimental data using an explicit expression for the scattering matrix is considered. This problem is ill-posed due to the presence of self-intersection points on the curves of the complex scattering coefficients. It is shown that the traditional multi-frequency measurement method used in vector network analyzers can be justified by the fact that the algorithm for processing the measurement results by the least squares method becomes stable if the number of frequencies is large enough. Keywords 1 waveguide, well-posedness of inverse problem, convergence conditions for the least squares method, multi-frequency measurements 1. Introduction Application of materials in various fields of science and technology such as material science, microwave engineering, aerospace, microelectronics, and communication industries requires the exact knowledge of material parameters such as permittivity and permeability. The most widely used methods for measuring dielectric materials parameters are cavity perturbation techniques and free- space and transmission line of waveguide methods. Each approach has its own advantages and limitations. To implement the free-space methods it is necessary to have sufficiently large sheets of the studied material. The accuracy of this method is low due to mismatch between the experimental setup and the simple mathematical model having an exact solution. The cavity perturbation techniques are more accurate, but they are applicable only over a narrow frequency band in which the resonance effect is observed. In our work we study the widely used method based on transmission line in the form of a rectangular metal (perfectly conducting) waveguide. If the material sample used in these measurements is a plane-parallel layer, then the mathematical model of the experiment is the simplest. In this case the exact explicit expression for the permittivity and permeability as a function of the S- parameters (the coefficient of transmission of the electromagnetic field through the inclusion and the reflection coefficient) is known. This procedure named Nicolson-Ross-Weir method (NRW) [1, 2]. Its disadvantage is the solution phase ambiguity for low-loss materials, except for those samples which width is less than one-half wavelength. In the general case of inclusion of an arbitrary shape numerical simulation of the experiment is necessary. The desired value of the dielectric constant of the material can be found as a solution to the inverse problem by comparing the experimental data and theoretical values of the transmission coefficient. Far Eastern Workshop on Computational Technologies and Intelligent Systems, March 2–3, 2021, Khabarovsk, Russia EMAIL: sheina_elena@mail.ru (E.A. Sheina); yury.shestopalov@hig.se (Y.V. Shestopalov) ORCID: 0000-0002-2691-2820 (Y.V. Shestopalov) © 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) We study the well-posedness condition for the inverse problem of determining the layer parameters from experimental data; namely, the existence and uniqueness of solution and its continuous dependence on the input data. Unfortunately this algorithm is improperly posed. In fact, (a) the range of the function specifying the transmission coefficient is a curve on the complex plane; therefore the probability that the experimental data belongs this curve is equal to zero; and (b) the parametric curve of the function on the complex plane has self-intersection points which means that the solution may not be unique. It is shown that the traditional multi-frequency method of measurements applied in vector network analyzers can be used to formulate a well-posed problem. For a noiseless experiment that perfectly matches the mathematical model the non-uniqueness of the solution can be eliminated if we consider the vector formulation of the problem determined by a set of frequencies in the range of single-mode waveguide. The vector function of the transmission coefficient becomes a one-to-one function of the inclusion dielectric constant if the frequency resolution is sufficiently small. For an actual physical experiment, the least squares method (LSM) can be applied for the solution of the inverse problem under study. The LSM solution converges to the desired value of the layer permittivity if the quality of the experiment (determined both by noise and defects of the measurement setup and material samples) is improved. The convergence rate is enhanced if the number of frequencies used in experiment is taken large enough. 2. The transmission coefficient of the principal mode The waveguide of rectangular cross section and perfectly conducting walls shown in Fig. 1 contains a dielectric parallel-plane diaphragm (layer). 1 = (0, a)  (0, b)  ( z−(layer ) , z+(layer ) ) with a dielectric inclusion  2 ,  2  1   . 1, r 0 \ 1 ,  The relative permittivity in the waveguide  (r ) /  0 =  teflon , r 1 \ 2 ,   cube , r 2 , where  teflon  2.04,  cube  3.8 (quartz), 10.0 (ruby), and  0 is permittivity of vacuum. Figure 1: Waveguide with a diaphragm and inclusion. We develop algorithms of reconstructing permittivity of the diaphragm (a test problem) and of a small inclusion inside the diaphragm by comparing the results of multi-frequency series of the field measurements on the waveguide flanges and the data obtained from the mathematical model corresponding to the experimental setup. The scheme of measurements using Vector Network Analyzer (VNA) are performed at the points (ports) shown in Fig. 1. Define f (1,0) = c / ( 2a ) , f (2,0) = 2 f(1,0) , f(1,0) , f (2,0) are the cutoff frequency for TE(1,0) , TE(2,0) modes, f (1,0)  6.52 (GHz). ˆ (r)e−it of It is known that in the empty parts of the waveguide the steady-state solution E(r, t ) = E Maxwell's equations is represented as a sum of harmonic and evanescent waves. If calculations are performed at the frequency f : f (1,0)  f  f (2,0) , when only one (principal) waveguide mode propagates in empty parts of the waveguide; the higher-order modes are evanescent (standing) waves decaying exponentially on both sides of the diaphragm. At the output (for z = z+( wg ) ) of the waveguide with inclusions in the diaphragm, the transmitted field has the form   Eˆ y( trans ) ( x, y, z+( wg ) ) = Eˆ y( trans ) ( wg ) ,(1,0) ( x, y , z + ) +  (1 −  ( n,1) ( m,0) ) Eˆ y( trans ) ( wg ) ,( n , m ) ( x, y , z + ), (1) m = 0 n =1 where (z) ( wg ) (z) ( wg ) Eˆ y(trans e 0,(1,0 ) + , Eˆ y(trans ik z −| k |z ) ( wg ) ,(1,0) ( x, y, z+ + ) = R(1,0) ) ( wg ) ,( n, m) ( x, y, z+ ) = R(+n , m ) X ( n ) ( x)Y( m ) ( y )e 0 ,( n ,m ) + , x  (0, a), y  (0, b) . If the waveguide is long enough then we can assume that Eˆ y(trans ) ( x, y, z+( wg ) )  Eˆ y(trans ) ( wg ) ,(1,0) ( x, y, z+ ). Here b a R(+n ,m ) =   X ( n ) ( x)Y( m ) ( y ) Eˆ y(trans ) ( x, y, z+( wg ) ) dxdy, 0 0 X ( n ) ( x) = (2 / a)1/ 2 sin(k((nx)) x) , Y( m ) ( y ) = ((2 −  ( m ,0) ) / b)1/ 2 cos(k((my )) y ) are the basis functions, k((nx)) =  n / a , k((my )) =  m / b , k0,( (z) n , m ) = ( k0 − ( k( n ) ) − ( k( m ) ) ) 2 ( x) 2 ( y ) 2 1/ 2 , k02 =  2 0 0 , n  , m  0 , and  ( n,m) = sgn | n − m | . For the sample under test situated between two antennas (ports) the measured quantities are the complex scattering S -matrix coefficients (transmission and reflection), where S S12  S =  11 ,  S 21 S11  S12 is the transmission coefficient from port 2 to port 1, calculated as the ratio of the field measured at the exit port in the presence of an inclusion to the source field measured at the input port. The input and the output (at the ports with z= z−( wg ) , z= z+( wg ) correspondingly) measurement data give the following elements of the scattering matrix describing the transmitted wave for the empty waveguide and the waveguide containing the diaphragm, respectively: Eˆ y(inc ) ( x, y, z+( wg ) ) ik0( z ) d ( wg ) S( wg ) = =e , Eˆ (inc ) ( x, y, z ( wg ) ) 12 y − (2) Eˆ y(trans ) ( x, y, z+( wg ) ) S12( wg ,layer ) = (inc ) . Eˆ ( x, y, z ( wg ) ) y − If measurements are taken at the layer boundary, then Eˆ y(trans ) ( x, y, z+(layer ) ) Eˆ y(inc ) ( x, y, z+(layer ) ) ik0( z ) d ( layer ) ( layer ) S12 = (inc ) (0,layer ) , S12 = (inc ) =e , Eˆ ( x, y, z (layer ) ) y − Eˆ ( x, y, z (layer ) ) y − d =z ( wg ) +−z ( wg ) is the waveguide length, d ( wg ) − =z −z is the layer width. ( layer ) (layer ) + (layer ) − Define the transmission coefficient of the principal waveguide mode as Eˆ y(trans ) ,(1,0) ( x, y, z+ ( wg ) ) F(1,0) = + = R(1,0) / A. (3) Eˆ ( inc ) ( ( x, y, z )wg ) y ,(1,0) + The waveguide transmission coefficient F1,0 defined according to (1), (2), and (3) can be represented as Eˆ y(trans ) ( x, y, z+( wg ) ) S12( wg ,layer ) F(1,0)  (inc ) = . Eˆ ( x, y, z ( wg ) ) S12( wg ) y + Thus the transmission coefficient of the principal mode through a layer in a waveguide is calculated as the ratio of the field measured at the output port far enough from inclusion to the field measured at the output port in the empty waveguide. Problem for the Maxwell's equations in the unbounded waveguide containing a dielectric layer with permittivity  admits the closed-form solution for transmission coefficient of the principal mode TE(1,0) : ( z ) ( layer ) 4t ( f )e − ik0 d F ( , f ) = ( z ) ( layer ) ( z ) ( layer ) , (1 + t ( f )) 2 e − ik d − (1 − t ( f )) 2 eik d or in another form [3] 1 1 F ( , f ) = , (4) S (0,layer ) 12 ( f ) g ( , f ) where g ( , f ) = c ( f ) − i H ( t ( f ) ) s ( f ) , (5) s ( f ) = sin ( k( z ) ( f )d (layer ) ) , c ( f ) = cos ( k( z ) ( f )d (layer ) ) , (6) k( z ) ( f ) ( ) ( ) 1/2 1/2 t ( f ) = (z) , k( z ) ( f ) = k( ,(1,0) z) ( f ) = k2 ( f ) − (k ( x ) )2 , k0( z ) ( f ) = k0,(1,0) ( z) ( f ) = k02 ( f ) − (k ( x ) )2 , k0 ( f ) k ( f ) =  1/2 k0 ( f ), k0 ( f ) = ( 0 0 )1/2 ,  ( f ) = 2 / k ( f ) ,  =  0 ,  1 H ( x) = 0.5  x +  ,  x with x  0. 3. Algorithms of experimental data processing Let us assume that the relative dielectric constant of the sample is real and does not depend on the frequency of oscillations of the electromagnetic field in a certain range within the single-mode range of the waveguide. If the material properties are a priori unknown, we will look for the effective dielectric constant of the inclusion for the selected bandwidth. Introduce the vectors (exp) N (exp) f = ( f1 ,..., f N (exp) )  ℝ N , F (exp) = ( F1(exp) ,..., FN(exp) (exp) )  ℂ of the frequency and complex-valued measurement data of N (exp) experiments. Consider the equation ( g  (layer ) , f = g (exp) , ) (7) for the (unknown) dielectric constant of the layer  (layer )  1 , where ( g( , f ) = g ( , f1 ),..., g ( , f N (exp) )  ℂ N ) (exp) with g defined in (5), (6), (exp) g(exp) = ( g1(exp) ,..., g N(exp) exp ) ℂ N , 1 1 g n(exp) = ( layer ) (exp) , Z 0 ( fn ) F n n = 1,..., N (exp) . We formulate inverse problems that constitute different permittivity reconstruction scenarios of the layer in the waveguide. To this end, let ( ) = { :   1}, ( ) = { :1    },   1 , and by G ( f , ( ) ) =  g( , f )  ℂ N ,  ( )  (exp) denote the set of values of function g( , f ) for the fixed frequency vector f (it is a curve in N (exp) - dimensional complex space). Formulate two inverse problems considered in the study. Problem 1 ( ) Find a real  (layer ) ( ) satisfying relation (7) for a given complex vector g(exp)  G f , ( ) with the fixed frequency vector f . Problem 2 Find a real  (layer ) ( ) satisfying relation (7) for a given complex vector g(exp) ℂ N (exp) with the fixed frequency vector f . Check the fulfillment of the well-posedness condition for these problems; namely, the existence and uniqueness of solution and its continuous dependence on the input data. Problem 1 describes a perfect experiment exactly corresponding to the mathematical model, it is ( ) solvable by the definition of the set G f , ( ) . However, its uniqueness may be violated. In fact, if N (exp) = 1 for any chosen frequency the solution is not unique due to the existence of a countable set { m }m=1,... , satisfying the relation sin(k m ( f ) d (layer ) ) = 0 that specifies self-intersections points of curve G( f , ( ) ) (Fig. 2). Using a priori information about  (layer ) we can achieve the uniqueness by adjusting domain ( ) and a frequency range [ f1 , f N (exp) ]. However, the formally well-posed problem may be ill-conditioned in the vicinity of the intersection points mentioned above where the parameter values are such that the quantity sin(k ( f )d (layer ) ) in the denominator virtually vanishes. Proposition 1 below demonstrates that for N (exp)  1 the solution to Problem 1 is unique if the frequency resolution is sufficiently small. In fact, g( , f ) becomes a one-to-one vector function of real variable  for a fixed set of frequency values f . Problem 2 simulates the processing of noisy experimental data. This problem is also ill-posed since it may be unsolvable: in actual experiments, it is typical that g(exp)  G(f , ( ) ) because the set (a curve) has the zero measure on the complex plane. We will replace Problem 2 with an LSM problem such that its solution approximates the sought solution of perfect Problem 1 when the defects of the experimental setup and measurement error decrease. 4. One-to-one correspondence Let us consider the problem of determining the value of the dielectric constant of a plane-parallel non-magnetic layer without losses. One can show that the transition from a single-frequency experiment to a multi-frequency experiment improves the properties of the inverse problem, ensuring its unique solvability. We present without proof the following statement [4]. Proposition 1 Set N (exp) = 1 . For any   1 there is one-to-one correspondence between ( ) and G( f , ( ) ) for the fixed frequency f if d (layer )  1. (8) 0.5E ( f ) Assume that N (exp)  2 . For any   1 there is one-to-one correspondence between ( ) and G(f , ( ) ) for the fixed frequency vector f if the following condition is satisfied in at least one of the two equivalent forms: d (layer ) d (layer ) − 1, (9) 0.5E ( f n +1 ) 0.5E ( f n ) n = 1,..., N (exp) − 1 , or c 1 h( f )  hE( f ) = ( layer ) . (10) 2d E1/2 Figure 2: The branches of the curve G f , ( ) , ( ) f = 9.25 GHz, ( ) = { :1.0    9.0} , 1( ) = (2.05,2.13) , (2 ) = (3.06,3.12) , their intersection point g (1 , f ) = g ( 2 , f ) , 1 = 2.09 1( ) ,  2 = 3.12 (2 ) , red point is experimental value g (exp) . 5. Multi-frequency least squares method Let us show that the solution of ill-posed Problem 2 can be found approximately with an accuracy determined by the quality of the experiment. Here along with the problem for determining the real dielectric constant of a low-loss sample we formulate a problem for determining the complex parameter of a loss sample without discussing its properties; the results of the numerical solution of both problems are shown in Fig. 3. Problem 3 (4) (LSM) Find  ( LS ,1)  ℝ N (  ( LS ,2)  ℂ N ) satisfying the condition (exp) (exp) ( g( ( LS ) , f ) − g(exp) = min g( , f ) − g(exp) ,  (E ) ) (exp) for a given vector g(exp) ℂ N with the fixed frequency vector f . We present without proof the following statements [5]. Proposition 2 The Problems 3 is solvable. Proposition 3 If the conditions of Proposition 1 are satisfied and (exp) g ( ( layer ) , f ) − g (exp) →0, (11) then  ( LS ,1) →  (layer ) . (12) Figure 3: Application of least squares to the determination of permittivity  min     max from transmission coefficient experimental data {Fn(exp) }n=1,... Nexp for frequency set { f n(exp) }n=1,... Nexp ; F ( min , f n(exp) )n=1,... Nexp , F ( max , f n(exp) )n=1,... Nexp are test points; F ( ( LS ,1 ) , f n(exp)  n =1,... Nexp , F ( ( LS ,2 ) , f n(exp)  n =1,... Nexp are theoretical transmission coefficient values for least squares solutions. 6. Conclusion Our findings have demonstrated that conducting an experiment in a multi-frequency mode makes it possible to turn a technical possibility into a mathematical achievement. The inverse problem of determining the sample parameter for one perfect experiment fully corresponding to the mathematical model may have a non-unique solution. Considering the results of several experiments as a vector, we arrive at an inverse problem which is well-posed if the frequency step is small enough. For a non- perfect multi-frequency experiment it is shown that the solution can be found by the least squares method. It converges to the solution of the perfect problem if the experimental data approach the theoretically predicted values for the unknown dielectric constant of the inclusion. 7. References [1] Nicolson, A. M., Ross, G. F. (1970), Measurement of the intrinsic properties of materials by time-domain techniques, IEEE Trans. Instrum. Meas., 19(7), 377 – 382. [2] Weir, W. B. (1974), Automatic Measurement of Complex Dielectric Constant and Permeability at Microwave Frequencies, Proceedings of the IEEE, 62(1), 33 – 36. [3] Shestopalov, Y., Smirnov, Y., & Derevyanchuk E. (2015), Permittivity Reconstruction of a Diaphragm in a Rectangular Waveguide: Unique Solvability of Benchmark Inverse Problems, PIERS Proc., Prague, 1528 – 1532. [4] Sheina, E., Shestopalov, Y., & Smirnov, A. (2019), Stability of least squares for the solution of an ill-posed inverse problem of reconstructing real value of the permittivity of a dielectric layer in a rectangular waveguide, PIERS Proc., Roma, 1732 – 1738. [5] Sheina, E., Shestopalov, Y., & Smirnov, A. (2021), Advantages of a Multi-Frequency Experiment for Determining the Dielectric Constant of a Layer in a Rectangular Waveguide and Free Space, Radio Science, DOI: 10.1029/2020RS007115