Wave Impedance Determination for PCB SIWs Considering Metal Roughness Loss Effects

This paper presents an expression to determine the complex wave impedance of a substrate-integrated waveguide for the dominant TE10 propagation mode, notably enhancing the accuracy in modeling the corresponding imaginary part. This was accomplished by systematically identifying the need to consider additional conductor losses caused by the interaction of the propagating fields with the conductor material. In fact, by using the proposed expression, the complex impedance can straightforwardly be determined by combining propagation constant data, and the resistance that represents the loss associated with longitudinal currents occurring at the top and bottom walls, which are influenced by the conductor surface roughness. This allows for completely describing the characteristics of the waveguide when assuming uniform propagation along its length. Furthermore, since the voltage–current, the power–current, and the power–voltage impedances can easily be obtained from the wave impedance, the proposal enables representing the matching characteristics of the waveguide for circuit design purposes. An agreement between simulated and experimental insertion loss is achieved for substrate-integrated waveguides of two different widths when reconstructing the corresponding S-parameters using the determined wave impedance.


Introduction
Substrate-integrated waveguides (SIWs) are widely used in high-speed electronics due to their compatibility with printed circuit boards (PCBs) [1][2][3][4].In fact, the corresponding structure and performance resemble those of a rectangular waveguide (RWG), enabling their application from interconnects [5,6] to resonators used for determining material properties [7][8][9] in microwave circuitry.In this regard, the full characterization of a uniform section of SIW propagating signals in the single mode can be achieved once the complex propagation constant (γ) and the wave impedance (Z wave ) are known.Furthermore, from these parameters, the properties of the constituting dielectric and conductor materials can be inferred [10], the effect of the SIW when used as an access component can be deembedded [11], and even equivalent circuit models can be implemented to represent these structures in practical circuits [12,13].
For the purpose of carrying out the abovementioned characterization, it is a common practice to obtain the attenuation (α) and phase delay (β) versus frequency (f ) curves by applying the eigensolution of the thru-reflect-line (TRL) algorithm to data measured to SIWs varying only in length [10], which is simple and accurate.Conversely, obtaining Z wave is problematic due to the significant effect of the signal launchers on the experimental return loss, even after de-embedding.Alternatively, treating the SIW as a RWG with effective dimensions [14] and assuming propagation in the dominant TE 10 mode, Z wave can be indirectly obtained from γ considering that the propagation occurs within a nonmagnetic medium characterized by the permeability of vacuum [15].Notwithstanding, since the conductor walls of the SIW are made of copper, the corresponding losses at microwave frequencies are noticeable and accentuated by the skin effect; thus, it is inferred that an effective complex permeability needs to be considered, rather than µ 0 , when calculating Z wave .
For the case of RWGs, a modeling approach to account for the losses occurring within the conductor material was proposed in [16] for obtaining Z wave , which provides good results for hollow waveguides.Evidently, the lossless dielectric assumption required to simplify the corresponding analysis yields substantial errors for SIWs constructed, even with low-loss PCB dielectric materials.Therefore, to take into account both dielectric and conductor losses, a transmission line representation can be employed to indirectly determine Z wave [12].Bear in mind, however, that in current technologies, the copper layers forming the top and bottom walls are intentionally roughened to promote adherence with the dielectric laminate [6].Therefore, it is essential to account for the increase in loss introduced by the rough conductor surface.In this regard, knowledge of the surface topography eases the model's implementation [17].Nevertheless, we demonstrate that using the nominal rms roughness (R q ) specified by the PCB manufacturer provides acceptable accuracy for characterizing SIWs operating at tens of gigahertz, and avoids the additional inspection of the conductor surfaces, either before or after fabrication.
Here, Z wave is obtained by combining experimental γ data and a parameter extraction for the conductor losses, including the effect of the roughness of the metal foils used to form the waveguide.In fact, we demonstrate that, even though the real part of Z wave can be computed from γ using the traditional approach [15], its imaginary part is strongly sensitive to losses occurring within the conductor material due to the time-varying fields.This sensitivity underscores the need for improvement in determining Im(Z wave ) to ensure accuracy in characterizing SIWs operating at microwave frequencies.For this reason, a new expression explicitly incorporating the resistance of the top and bottom waveguide walls for calculating Z wave is presented in this paper, along with the corresponding implementation for experimentally analyzing prototypes in practical PCB technologies.Moreover, since the increase in the resistance with frequency is substantial when these walls are rough, this effect is also considered during the development of the proposal.
Once Z wave is determined, it can be used to analyze the intrinsic properties of the waveguide propagating in the single TE 10 mode in a similar manner, as the characteristic impedance concept is associated with transmission lines operating in the quasi-transverse electromagnetic mode.Moreover, the power-current, voltage-current, and power-voltage impedances (i.e., Z vi , Z pi , and Z pv , respectively) can directly be obtained depending on the desired description for the SIW [18].Additionally, although the resistance associated with the top and bottom walls is demonstrated to significantly impact the response of PCB SIWs, it is also observed that the effect of the corresponding internal inductance is negligible.This enables the determination of the complex relative dielectric permittivity (ε r ) and, consequently, the loss tangent (tanδ) of the laminate used as the PCB substrate within the measurement frequency range.

Relationship between Z wave and γ in PCB SIWs
Consider the sketch shown in Figure 1, where a simplified representation of a PCB SIW is presented.This structure can be assumed to behave as a RWG, provided that the adjacent metallic posts forming the side walls exhibit a separation much smaller than the wavelength of the propagating signals [10], which can easily be achieved using commercial PCB fabrication processes and for circuit operation at microwave frequencies.Hence, this section is dedicated to establishing a relationship that allows for the determination of Z wave from γ, assuming that the single TE 10 propagation mode as in a RWG is taking place.In this case, Z wave is usually calculated from the following simplified expression [15]: where j 2 = −1, ω = 2πf is the angular frequency and µ 0 is the permeability of vacuum.To distinguish the real and imaginary parts of Z wave , Equation ( 1) is expanded to yield: where γ = α + jβ was used, and α 2 ≪ β 2 was assumed, given that the phase delay significantly outweighs the attenuation in practical scenarios.Notice in these expressions that the magnetic field is assumed to be propagating through a vacuum, hence, µ 0 is used.To circumvent this inaccuracy when the interaction of the magnetic field with the conductor material causes loss, the well-known RLGC transmission line model in Figure 2a is useful for analyzing the waveguide [16], which can be modified into the equivalent representation shown in Figure 2b [12].Notice in these figures that, since the per-unit-length series elements R and L exhibit no variation between the two models, Z wave can be expressed as for a generic transmission line by [19]: Micromachines 2024, 15, x FOR PEER REVIEW 3 of 13 where j 2 = −1, ω = 2πf is the angular frequency and μ0 is the permeability of vacuum.To distinguish the real and imaginary parts of Zwave, Equation ( 1) is expanded to yield: and where γ = α + jβ was used, and α 2 ≪ β 2 was assumed, given that the phase delay significantly outweighs the attenuation in practical scenarios.Notice in these expressions that the magnetic field is assumed to be propagating through a vacuum, hence,  0 is used.To circumvent this inaccuracy when the interaction of the magnetic field with the conductor material causes loss, the well-known RLGC transmission line model in Figure 2a is useful for analyzing the waveguide [16], which can be modified into the equivalent representation shown in Figure 2b [12].Notice in these figures that, since the per-unit-length series elements R and L exhibit no variation between the two models, Zwave can be expressed as for a generic transmission line by [19]: It is important to notice that Equation (4) becomes Equation (1) in case R = 0 and L = μ0 are assumed, which is equivalent to neglecting the loss and delay associated with longitudinal currents occurring on the top and bottom walls of the SIW.To assess the impact of neglecting these effects, here, an expression for the complex Zwave is developed, considering R > 0 and L the total series inductance of the waveguide.This inductance is provided by the sum of the inductances accounting for the interaction of the magnetic field with the non-magnetic media external to the conductor material (i.e., represented by μ0) and internal to it (i.e., Lint).For this purpose, several assumptions are made.For instance, although  > 0 should be considered for SIWs, in PCB technology, these waveguides are made of low-loss dielectric materials, and copper.Hence,  ≪  at microwave frequencies.In this case, after substituting L = μ0 + Lint and γ = α + jβ into Equation (4), it is possible to obtain: where j 2 = −1, ω = 2πf is the angular frequency and μ0 is the permeability of vacuum.To distinguish the real and imaginary parts of Zwave, Equation ( 1) is expanded to yield: and where γ = α + jβ was used, and α 2 ≪ β 2 was assumed, given that the phase delay significantly outweighs the attenuation in practical scenarios.Notice in these expressions that the magnetic field is assumed to be propagating through a vacuum, hence,  0 is used.To circumvent this inaccuracy when the interaction of the magnetic field with the conductor material causes loss, the well-known RLGC transmission line model in Figure 2a is useful for analyzing the waveguide [16], which can be modified into the equivalent representation shown in Figure 2b [12].Notice in these figures that, since the per-unit-length series elements R and L exhibit no variation between the two models, Zwave can be expressed as for a generic transmission line by [19]: It is important to notice that Equation (4) becomes Equation (1) in case R = 0 and L = μ0 are assumed, which is equivalent to neglecting the loss and delay associated with longitudinal currents occurring on the top and bottom walls of the SIW.To assess the impact of neglecting these effects, here, an expression for the complex Zwave is developed, considering R > 0 and L the total series inductance of the waveguide.This inductance is provided by the sum of the inductances accounting for the interaction of the magnetic field with the non-magnetic media external to the conductor material (i.e., represented by μ0) and internal to it (i.e., Lint).For this purpose, several assumptions are made.For instance, although  > 0 should be considered for SIWs, in PCB technology, these waveguides are made of low-loss dielectric materials, and copper.Hence,  ≪  at microwave frequencies.In this case, after substituting L = μ0 + Lint and γ = α + jβ into Equation (4), it is possible to obtain:  It is important to notice that Equation (4) becomes Equation (1) in case R = 0 and L = µ 0 are assumed, which is equivalent to neglecting the loss and delay associated with longitudinal currents occurring on the top and bottom walls of the SIW.To assess the impact of neglecting these effects, here, an expression for the complex Z wave is developed, considering R > 0 and L the total series inductance of the waveguide.This inductance is provided by the sum of the inductances accounting for the interaction of the magnetic field with the non-magnetic media external to the conductor material (i.e., represented by µ 0 ) and internal to it (i.e., L int ).For this purpose, several assumptions are made.For instance, although α > 0 should be considered for SIWs, in PCB technology, these waveguides are made of low-loss dielectric materials, and copper.Hence, α ≪ β at microwave frequencies.In this case, after substituting L = µ 0 + L int and γ = α + jβ into Equation (4), it is possible to obtain: and Bear in mind that R and L are series elements in the transmission line model and thus are associated with longitudinal currents, which only take place in the top and bottom walls of an SIW operating in the TE 10 mode.Therefore, to quantify the impacts of R and L in Z wave through Equations ( 5) and ( 6), firstly, the surface impedance (Z S ) for a smooth conductor is expressed as [20]: In this equation, the resistance under a smooth-conductor assumption (R smooth ) can be calculated from the skin depth (δ), the conductor conductivity (σ), and the height (h) of the waveguide in the direction parallel to the electric field by means of [12]: where the numerator is multiplied by 2 to consider the resistance of the top and bottom walls, and δ is given by: Now, consider that the top and bottom walls of the SIW exhibit roughness, as illustrated in Figure 3.In this case, the imperfect surface introduces additional wave scattering that increases the loss experienced by a signal propagating along the waveguide.For incorporating this effect into the surface impedance model, the complex and frequencydependent roughness factor ( krough = k 1 + jk 2 ) is combined with Equation ( 7) to write [21]: Micromachines 2024, 15, x FOR PEER REVIEW 4 of 13 and Bear in mind that R and L are series elements in the transmission line model and thus are associated with longitudinal currents, which only take place in the top and bottom walls of an SIW operating in the TE10 mode.Therefore, to quantify the impacts of R and L in Zwave through Equations ( 5) and ( 6), firstly, the surface impedance (ZS) for a smooth conductor is expressed as [20]: In this equation, the resistance under a smooth-conductor assumption (Rsmooth) can be calculated from the skin depth (δ), the conductor conductivity (σ), and the height (h) of the waveguide in the direction parallel to the electric field by means of [12]: where the numerator is multiplied by 2 to consider the resistance of the top and bottom walls, and δ is given by: Now, consider that the top and bottom walls of the SIW exhibit roughness, as illustrated in Figure 3.In this case, the imperfect surface introduces additional wave scattering that increases the loss experienced by a signal propagating along the waveguide.For incorporating this effect into the surface impedance model, the complex and frequency-dependent roughness factor ( ̂rough =  1 +  2 ) is combined with Equation ( 7) to write [21]: From Equation (10), considering that the surfaces of the top and bottom walls are rough, R and  int can respectively be expressed as: and where  R =  1 −  2 and  L =  1 +  2 .
To continue with the analysis corresponding with the wave impedance, substituting Equations ( 11) and ( 12) into ( 5) and ( 6) and making some simplifications yield: From Equation (10), considering that the surfaces of the top and bottom walls are rough, R and L int can respectively be expressed as: and where To continue with the analysis corresponding with the wave impedance, substituting Equations ( 11) and ( 12) into ( 5) and ( 6) and making some simplifications yield: and In these equations, (αK R + βK L ) ≈ βK L was assumed, since not only α ≪ β, but also K L is a few times larger than K R for typical copper foils on PCB technology [22].
In order to verify the validity of the presented analysis, a RWG uniform along its length and exhibiting the cross-section depicted in Figure 4a was simulated using the commercial full-wave electromagnetic (EM) solver Ansys HFSS, employing the modal solution.In this regard, when establishing the simulation conditions for the RWG, waveguide ports were defined at the terminations of the structure.The material properties were specified for the filling dielectric, as depicted in Figure 4a, while copper was assumed to constitute the walls.The simulation covers the frequency range from 25 GHz to 80 GHz with a discrete sweep (i.e., without interpolations).Additionally, the dimensions of the waveguide are set to achieve a cutoff frequency of the TE 10 mode, f TE10 ≈ 27 GHz.Two cases were simulated, one considering that all the inner surfaces of the walls are smooth, and another where the top and bottom walls exhibit a standard profile encountered in PCB technology.For this latter case, Huray's causal model was defined to represent the surface impedance during the EM simulation by establishing a nodule radius n r = 0.5 µm and the Hall-Huray surface ratio SR = 3 [23].In these simulations, the complex Z wave was directly obtained from one of the excitation waveguide ports.
In these equations, ( R +  L ) ≈  L was assumed, since not only  ≪ , but also  L is a few times larger than  R for typical copper foils on PCB technology [22].
In order to verify the validity of the presented analysis, a RWG uniform along its length and exhibiting the cross-section depicted in Figure 4a was simulated using the commercial full-wave electromagnetic (EM) solver Ansys HFSS, employing the modal solution.In this regard, when establishing the simulation conditions for the RWG, waveguide ports were defined at the terminations of the structure.The material properties were specified for the filling dielectric, as depicted in Figure 4a, while copper was assumed to constitute the walls.The simulation covers the frequency range from 25 GHz to 80 GHz with a discrete sweep (i.e., without interpolations).Additionally, the dimensions of the waveguide are set to achieve a cutoff frequency of the TE10 mode, fTE10 ≈ 27 GHz.Two cases were simulated, one considering that all the inner surfaces of the walls are smooth, and another where the top and bottom walls exhibit a standard profile encountered in PCB technology.For this latter case, Huray's causal model was defined to represent the surface impedance during the EM simulation by establishing a nodule radius nr = 0.5 μm and the Hall-Huray surface ratio SR = 3 [23].In these simulations, the complex Zwave was directly obtained from one of the excitation waveguide ports.As illustrated in Figure 4b, no noticeable change was observed in the simulated Re(Zwave) curves with and without considering the conductor surface roughness, which suggests that the effect of the internal inductance for this structure is negligible for PCB materials and within this frequency range.This is expected since Lint significantly decreases with frequency, as expressed in Equation (12).Hence, including the internal inductance term in Equation ( 13) is necessary only at frequencies where the skin effect allows the current to flow within a significant portion of the cross-section of the top and bottom walls.Considering copper, this occurs well below 10 GHz.In this case, SIWs would need to be excessively wide to sufficiently lower their cutoff frequency, rendering them impractical for PCB applications.Hence,  0 ≫  L  smooth can be assumed in Equation (13), making it possible to apply Equation (2) for obtaining Re(Zwave) with accuracy in PCB SIWs operating at tens of gigahertz, as also shown in Figure 4b.Conversely, even when omitting the term  L  smooth in Equation ( 14), the resulting expression for the imaginary part of the wave impedance differs from Equation (3); this is: As illustrated in Figure 4b, no noticeable change was observed in the simulated Re(Z wave ) curves with and without considering the conductor surface roughness, which suggests that the effect of the internal inductance for this structure is negligible for PCB materials and within this frequency range.This is expected since L int significantly decreases with frequency, as expressed in Equation (12).Hence, including the internal inductance term in Equation ( 13) is necessary only at frequencies where the skin effect allows the current to flow within a significant portion of the cross-section of the top and bottom walls.Considering copper, this occurs well below 10 GHz.In this case, SIWs would need to be excessively wide to sufficiently lower their cutoff frequency, rendering them impractical for PCB applications.Hence, ωµ 0 ≫ K L R smooth can be assumed in Equation (13), making it possible to apply Equation (2) for obtaining Re(Z wave ) with accuracy in PCB SIWs operating at tens of gigahertz, as also shown in Figure 4b.Conversely, even when omitting the term K L R smooth in Equation ( 14), the resulting expression for the imaginary part of the wave impedance differs from Equation (3); this is: where R is defined in Equation (11).This is the novel definition of Im(Z wave ) reported in this paper.In fact, Figure 4c illustrates that using the simulated α, β, and R and applying Equation ( 15) reproduces, with accuracy, Im(Z wave ) for the smooth and rough cases.Conversely, when comparing the imaginary part of the simulated Z wave with the calculation using Equation (1), significant discrepancy is observed in Figure 4c.Also note that, in this figure, determining Z wave from the simulated γ and involving the effective dimensions of the waveguide, as in [16], is not appropriate as the dielectric losses are assumed to be negligible for the corresponding formulation to remain valid.In fact, the approach in [16] is more suitable for hollow waveguides.Therefore, as explained in this section, accounting for the conductor effects in the determination of Im(Z wave ) requires considering the βR term as proposed in this work.As β = Im(γ) can be easily extracted from experiments, the only unknown in Equation ( 15) is R. A simple strategy to obtain this parameter is explained hereafter.

Experimental Details
Several SIW structures were built on the PCB prototype shown in Figure 5.The dielectric laminates used in this board exhibited a thickness, h = 1.4 mm, and nominal relative permittivity ε ′ r = 2.2 and loss tangent tanδ = 0.0009 at 10 GHz [24].On the other hand, the foils forming the top and bottom walls of the SIWs are made of low-profile copper with a root-mean-square surface roughness, R q = 1.198 µm ± 0.053 µm, as specified by the PCB manufacturer.The sidewalls are implemented using arrays of copper via holes, as described in Figure 5. Now, considering these dimensions and equivalences with a RWG [25], from the transverse center-to-center width w = 4.1 mm, a cutoff frequency f TE10 ≈ 27 GHz was obtained, whereas for w = 2.3 mm, the result is f TE10 ≈ 47 GHz.Since these SIWs are constructed on the same board layer, the propagating waves can be assumed to experience identical properties of the dielectric and conductor media.Thus, including SIWs of different widths within the prototype provides a way to experimentally verify consistency in the parameter extraction methodology for structures exhibiting different electrical responses.The microstrip lines used to apply the signals to the SIWs are shown in Figure 5.These lines were designed to excite odd electric modes and feature coplanar pad terminations for landing probes with a pitch of 250 µm.This allows for measuring S-parameters up to 80 GHz using a vector network analyzer (VNA), which was calibrated to shift the measurement plane up to the end of the probes through a line-reflect-reflect-match (LRRM) algorithm.where R is defined in Equation (11).This is the novel definition of Im(Zwave) reported in this paper.In fact, Figure 4c illustrates that using the simulated α, β, and R and applying Equation ( 15) reproduces, with accuracy, Im(Zwave) for the smooth and rough cases.Conversely, when comparing the imaginary part of the simulated Zwave with the calculation using Equation (1), significant discrepancy is observed in Figure 4c.Also note that, in this figure, determining Zwave from the simulated γ and involving the effective dimensions of the waveguide, as in [16], is not appropriate as the dielectric losses are assumed to be negligible for the corresponding formulation to remain valid.In fact, the approach in [16] is more suitable for hollow waveguides.Therefore, as explained in this section, accounting for the conductor effects in the determination of Im(Zwave) requires considering the  term as proposed in this work.As β = Im(γ) can be easily extracted from experiments, the only unknown in Equation ( 15) is R. A simple strategy to obtain this parameter is explained hereafter.

Experimental Details
Several SIW structures were built on the PCB prototype shown in Figure 5.The dielectric laminates used in this board exhibited a thickness, h = 1.4 mm, and nominal relative permittivity   ′ = 2.2 and loss tangent tanδ = 0.0009 at 10 GHz [24].On the other hand, the foils forming the top and bottom walls of the SIWs are made of low-profile copper with a root-mean-square surface roughness, Rq = 1.198 μm ± 0.053 μm, as specified by the PCB manufacturer.The sidewalls are implemented using arrays of copper via holes, as described in Figure 5. Now, considering these dimensions and equivalences with a RWG [25], from the transverse center-to-center width w = 4.1 mm, a cutoff frequency fTE10 ≈ 27 GHz was obtained, whereas for w = 2.3 mm, the result is fTE10 ≈ 47 GHz.Since these SIWs are constructed on the same board layer, the propagating waves can be assumed to experience identical properties of the dielectric and conductor media.Thus, including SIWs of different widths within the prototype provides a way to experimentally verify consistency in the parameter extraction methodology for structures exhibiting different electrical responses.The microstrip lines used to apply the signals to the SIWs are shown in Figure 5.These lines were designed to excite odd electric modes and feature coplanar pad terminations for landing probes with a pitch of 250 μm.This allows for measuring S-parameters up to 80 GHz using a vector network analyzer (VNA), which was calibrated to shift the measurement plane up to the end of the probes through a line-reflect-reflect-match (LRRM) algorithm.Since the core of the methodology is using experimental γ to characterize the uniform section of an SIW, these data are obtained free of the effect of the signal launchers (i.e., including the microstrip sections) using a line-line method based on the TRL formulation [26].For this purpose, SIWs of two lengths, l1 = 76.2 mm and l2 = 254 mm, were implemented in the prototype for each width.Since the core of the methodology is using experimental γ to characterize the uniform section of an SIW, these data are obtained free of the effect of the signal launchers (i.e., including the microstrip sections) using a line-line method based on the TRL formulation [26].For this purpose, SIWs of two lengths, l 1 = 76.2 mm and l 2 = 254 mm, were implemented in the prototype for each width.

Parameter Determination
Electromagnetic waves propagate through dielectric and conductor media within SIWs and RWGs; thus, the overall energy loss and delay involve the interaction of waves with these two types of materials.Consequently, a simultaneous determination of model parameters for the effects associated with these interactions is necessary to assess their impact on the waveguide performance.In this regard, although obtaining the dielectric complex permittivity is not necessary to implement the model for Z wave proposed in this paper, it is nonetheless included as the first step of the parameter-extraction methodology.This allows for presenting a complete modeling approach for a uniform section of SIW and facilitates comparing the dielectric and conductor-related losses occurring at microwave frequencies.
As demonstrated in Section 2, the effect of the internal inductance on the performance of the SIW can be neglected for PCB SIWs operating at tens of gigahertz; thus, since β = Im(γ) is already known from experiments, the relative permittivity ε ′ r = Re(ε r ) can be obtained from [27]: where w e f f = w − d 2 /0.95s is the SIW effective width [10], ε 0 is the vacuum permittivity, and d and s are defined in Figure 5.
For PCB laminates, Djordjevic's model can be used to consider the causal relationship between ε ′ r and tanδ [28].In the corresponding mathematical representation, ε ′ r is assumed to exhibit a variation ∆ε′ between a lower and an upper frequency limit, ω 1 and ω 2 , respectively.Furthermore, a useful simplification can be made without significantly losing accuracy at frequencies where ω ≫ ω 1 and ω ≪ ω 2 .For the real part of the permittivity, the approximation is as follows [28]: where ε ′ ∞ is the asymptotic value toward which ε ′ r approaches at high frequencies, whereas 1 = log10(ω 1 ) and m 2 = log10(ω 2 ).For the case studied here, the intention is represent- ing ε ′ r well above the cutoff frequency of the TE 10 mode for the SIW with w 1 = 4.1 mm, and below the cutoff frequency of the TE 30 mode; this allows for assuming single-mode operation.In round numbers, this range is defined from 30 GHz to 80 GHz.Hence, for Equation (17) to remain valid, ω 1 = 10 4 rad/s and ω 2 = 10 12 rad/s are arbitrarily defined.
6 illustrates the excellent model-experiment correlation for the relative permittivity for the SIW with w 1 = 4.1 mm, where ε ∞ = 2.204 and ∆ε ′ = 0.026 were obtained through curve fitting.Furthermore, Figure 6 also shows that, even though the model was implemented from data corresponding to one of the SIW widths, the ε ′ r curves for the two structures with different widths barely differed from each other well above the cutoff frequency of their fundamental modes.This verifies the usefulness of using more than one SIW for obtaining experimental broadband relative permittivity data.Now, employing Equation ( 17) under the same assumptions explained in [28], the loss tangent is approximately given by: which allows for determining tanδ ≈ 0.001 as shown in Figure 6.
The next step is to calculate R, which is associated with the longitudinal currents taking place along the top and bottom walls of the SIW and thus suffers from the additional loss introduced by the surface roughness effect on these walls [29].In this regard, R, can be obtained using Equation (11), which requires the knowledge of the roughness correction factor [30]. Owing to its simplicity and ability to maintain accuracy for low-profile copper foils, the mathematical form presented in reference [31] for K R is utilized here.This model assumes that a pyramidal array of spheres of fixed radius r allows for defining a unit cell for representing the surface roughness, and is expressed as: In this equation, c x = 2πµ 0 σr 2 , and σ is assumed here to be that of copper, σ ≈ 5.8 × 10 7 S/m.Furthermore, considering an array of 14 spheres, ∆K = 8.33 and r ≈ R q /4.8 [32]; thus, r ≈ 0.25 µm is obtained by rounding R q = 1.2 µm from the PCB manufacturer's specification detailed in Section 3.  17) under the same assumptions explained in [28], the loss tangent is approximately given by: tan ≈ π 2 • ∆′ ln (10)( 2 − 1 )  ′ (18) which allows for determining tanδ ≈ 0.001 as shown in Figure 6.
The next step is to calculate R, which is associated with the longitudinal currents taking place along the top and bottom walls of the SIW and thus suffers from the additional loss introduced by the surface roughness effect on these walls [29].In this regard, R, can be obtained using Equation (11), which requires the knowledge of the roughness correction factor [30]. Owing to its simplicity and ability to maintain accuracy for low-profile copper foils, the mathematical form presented in reference [31] for KR is utilized here.This model assumes that a pyramidal array of spheres of fixed radius r allows for defining a unit cell for representing the surface roughness, and is expressed as: In this equation,   = 2 0  2 , and σ is assumed here to be that of copper, σ ≈ 5.8×10 7 S/m.Furthermore, considering an array of 14 spheres, ΔK = 8.33 and r ≈ Rq/4.8 [32]; thus, r ≈ 0.25 μm is obtained by rounding Rq = 1.2 μm from the PCB manufacturer's specification detailed in Section 3.
To verify the suitability of using Equation ( 19), a 3D model of a uniform section of the SIW with w1 = 4.1 mm was simulated in Ansys HFSS.The simulation setup was defined as described in Section 2. However, in this case, the sidewalls were modeled using metallic posts, consistent with the fabricated structure shown in Figure 5.Moreover, the frequency-dependent permittivity and loss tangent of the dielectric material were defined based on the results shown in Figure 6.Afterward, a model-experiment correlation was performed for the attenuation curve by only varying the surface impedance on the top and bottom walls in the 3D model.The result of this correlation is shown in Figure 7a, demonstrating excellent agreement up to approximately 60 GHz.Beyond this frequency, the trend is well-predicted by the EM model, but some data dispersion is observed.This is attributed to the increased difficulty in accurately determining α at higher frequencies, even after rigorous de-embedding of the signal launchers, which involves multiple transitions.From this result, it is possible to illustrate in Figure 7b the concordance between KR obtained from the simulation and calculated using Equation (19).Hence, using the nominal value of Rq provided by the manufacturer enables a very approximate calculation To verify the suitability of using Equation ( 19), a 3D model of a uniform section of the SIW with w 1 = 4.1 mm was simulated in Ansys HFSS.The simulation setup was defined as described in Section 2. However, in this case, the sidewalls were modeled using metallic posts, consistent with the fabricated structure shown in Figure 5.Moreover, the frequencydependent permittivity and loss tangent of the dielectric material were defined based on the results shown in Figure 6.Afterward, a model-experiment correlation was performed for the attenuation curve by only varying the surface impedance on the top and bottom walls in the 3D model.The result of this correlation is shown in Figure 7a, demonstrating excellent agreement up to approximately 60 GHz.Beyond this frequency, the trend is well-predicted by the EM model, but some data dispersion is observed.This is attributed to the increased difficulty in accurately determining α at higher frequencies, even after rigorous de-embedding of the signal launchers, which involves multiple transitions.From this result, it is possible to illustrate in Figure 7b the concordance between K R obtained from the simulation and calculated using Equation (19).Hence, using the nominal value of R q provided by the manufacturer enables a very approximate calculation of K R , useful for the purposes of this paper.This fact is further corroborated in the Results Section.

Results
Figure 8a shows the Re(Zwave) versus frequency curve obtained from an EM simulation including the surface roughness effect for the prototyped SIW with w = 4.1 mm.This simulation is accurately reproduced by the calculation performed using Equation ( 2) and the experimental β.As previously mentioned, at frequencies of gigahertz, the inductance

Results
Figure 8a shows the Re(Z wave ) versus frequency curve obtained from an EM simulation including the surface roughness effect for the prototyped SIW with w = 4.1 mm.This simulation is accurately reproduced by the calculation performed using Equation ( 2) and the experimental β.As previously mentioned, at frequencies of gigahertz, the inductance increase introduced by the surface roughness effect is negligible for practical PCB SIWs, and thus no correction in this equation is necessary.Nonetheless, in Figure 8b, it is observed that Equation (3) fails representing Im(Z wave ), since the losses associated with the longitudinal currents occurring in the top and bottom conductor walls are neglected, underestimating the SIW's overall loss.Likewise, as expected, the curve corresponding to neglecting the dielectric losses to determine Im(Z wave ) from the experimental γ using the approach in [16] shows a significant discrepancy with the EM simulation.Conversely, the proposed Equation (15) and the parameter extraction strategy that led to its implementation allow to obtain a curve that correlates the EM simulation with accuracy.Complementary, Figure 9a shows the Re(Z wave ) curves corresponding to the SIW that exhibits a higher cutoff frequency for the TE 10 mode.For this structure, bear in mind that its reduced cross-section introduces higher insertion loss, which contributes to increasing the measurement uncertainty for loss-related effects.In this case, despite the data dispersion introduced by involving the experimental α for the calculation of Im(Z wave ), there is a noticeable downward shift in the corresponding curve when considering the impact of R on the performance of this narrower SIW; this is illustrated in Figure 9b.

Results
Figure 8a shows the Re(Zwave) versus frequency curve obtained from an EM simulation including the surface roughness effect for the prototyped SIW with w = 4.1 mm.This simulation is accurately reproduced by the calculation performed using Equation ( 2) and the experimental β.As previously mentioned, at frequencies of gigahertz, the inductance increase introduced by the surface roughness effect is negligible for practical PCB SIWs, and thus no correction in this equation is necessary.Nonetheless, in Figure 8b, it is observed that Equation (3) fails representing Im(Zwave), since the losses associated with the longitudinal currents occurring in the top and bottom conductor walls are neglected, underestimating the SIW's overall loss.Likewise, as expected, the curve corresponding to neglecting the dielectric losses to determine Im(Zwave) from the experimental γ using the approach in [16] shows a significant discrepancy with the EM simulation.Conversely, the proposed Equation (15) and the parameter extraction strategy that led to its implementation allow to obtain a curve that correlates the EM simulation with accuracy.Complementary, Figure 9a shows the Re(Zwave) curves corresponding to the SIW that exhibits a higher cutoff frequency for the TE10 mode.For this structure, bear in mind that its reduced crosssection introduces higher insertion loss, which contributes to increasing the measurement uncertainty for loss-related effects.In this case, despite the data dispersion introduced by involving the experimental α for the calculation of Im(Zwave), there is a noticeable downward shift in the corresponding curve when considering the impact of R on the performance of this narrower SIW; this is illustrated in Figure 9b.To demonstrate the application of the proposed model for Zwave, the curves for the equivalent circuit elements in Figure 2b are obtained versus frequency using the following equations: To demonstrate the application of the proposed model for Z wave , the curves for the equivalent circuit elements in Figure 2b are obtained versus frequency using the following equations: R = Re(γZ wave ) where Y RL is the shunt admittance associated with the series connection of R′ and L′.Thus, when using the experimentally determined data for γ, ε ′ r , and tanδ in the previous equations, the curves shown in Figure 10 are obtained for the R and G′ elements, which represent the more significant losses occurring in the SIW.A correlation is observed between the experimentally determined data and EM simulations performed using the 3D model for the SIW.In contrast, unexpected curve trends for R and G′ are obtained when neglecting the dielectric losses in the calculation of Z wave : R is clearly overestimated while G′ drops with frequency.This can be seen in Figure 10a,b, respectively.Likewise, ignoring the longitudinal currents in the top and bottom walls when obtaining this impedance as in Equation ( 1) yields R = 0, whereas R′, which is predicted to be negligible by the EM simulation, is overestimated.Finally, the S-parameters of an SIW were reconstructed using the experimental γ and Zwave obtained from Equations (3) and (15).To do so, the following transmission linebased equation was employed to firstly obtain the ABCD parameters of a uniform section of SIW of length (l): Afterward, an ABCD-to-S-parameter transformation is performed by considering Zwave as the reference impedance.This implies that perfect matching is assumed at the SIW terminations, and the signal transmission is maximized.Hence, the model can be directly assessed by analyzing the insertion loss (i.e., |S21|).Figure 11a,b show that the reconstructed magnitude and phase for S21 corresponding to the uniform section of the SIW agree with the EM simulations for the two considered widths.Furthermore, considering the experimental γ and Zwave obtained through applying the proposal, it is possible to implement the equivalent circuit for the SIW using Equations ( 20)-( 26).This enables the representation of the insertion loss of the SIW in a circuit simulator.The corresponding curves are also shown in Figure 11a,b.This verifies the usefulness of the proposal for obtaining Finally, the S-parameters of an SIW were reconstructed using the experimental γ and Z wave obtained from Equations (3) and (15).To do so, the following transmission linebased equation was employed to firstly obtain the ABCD parameters of a uniform section of SIW of length (l): T = cosh(γl) Z wave sinh(γl) Z −1 wave sinh(γl) cosh(γl) (27) Afterward, an ABCD-to-S-parameter transformation is performed by considering Z wave as the reference impedance.This implies that perfect matching is assumed at the SIW terminations, and the signal transmission is maximized.Hence, the model can be directly assessed by analyzing the insertion loss (i.e., |S 21 |). Figure 11a,b show that the reconstructed magnitude and phase for S 21 corresponding to the uniform section of the SIW agree with the EM simulations for the two considered widths.Furthermore, considering the experimental γ and Z wave obtained through applying the proposal, it is possible to implement the equivalent circuit for the SIW using Equations ( 20)-( 26).This enables the representation of the insertion loss of the SIW in a circuit simulator.The corresponding curves are also shown in Figure 11a,b.This verifies the usefulness of the proposal for obtaining Z wave , which can also be converted to other impedance definitions (i.e., Z pi , Z vi , and Z pv ) [18,33].
wave −1 sinh () cosh () Afterward, an ABCD-to-S-parameter transformation is performed by considering Zwave as the reference impedance.This implies that perfect matching is assumed at the SIW terminations, and the signal transmission is maximized.Hence, the model can be directly assessed by analyzing the insertion loss (i.e., |S21|).Figure 11a,b show that the reconstructed magnitude and phase for S21 corresponding to the uniform section of the SIW agree with the EM simulations for the two considered widths.Furthermore, considering the experimental γ and Zwave obtained through applying the proposal, it is possible to implement the equivalent circuit for the SIW using Equations ( 20)- (26).This enables the representation of the insertion loss of the SIW in a circuit simulator.The corresponding curves are also shown in Figure 11a,b.This verifies the usefulness of the proposal for obtaining Zwave, which can also be converted to other impedance definitions (i.e., Zpi, Zvi, and Zpv) [18,33].

Conclusions
An expression for obtaining the complex wave impedance of an SIW operating in the TE10 propagation mode was proposed.It was shown that the effect of the conductor surface roughness exhibited by the top and bottom walls of the SIW substantially impacts the imaginary part of this parameter, which was taken into consideration in the proposed expression and was successfully used to represent Zwave for SIWs of different widths up to 80 GHz.Furthermore, the model implementation based on the proposed expression can be applied to other SIW dimensions, provided that single-mode operation is maintained and the equivalence with a RWG remains valid.However, it is important to bear in mind that extending the application to higher frequencies necessitates verifying the appropriateness of the representation used for the series resistance, which must consider the specific metal surface profile required to keep metal losses within practical limits.

Figure 1 .Figure 2 .
Figure 1.Sketch showing the structure of a PCB substrate-integrated waveguide.

Figure 1 . 13 𝑍
Figure 1.Sketch showing the structure of a PCB substrate-integrated waveguide.

Figure 1 .Figure 2 .
Figure 1.Sketch showing the structure of a PCB substrate-integrated waveguide.

Figure 2 .
Figure 2. RLGC model for a RWG: (a) generic transmission line model and (b) modified topology reported in [12].

Figure 3 .
Figure 3. Simplified sketch of the cross-section of a PCB SIW illustrating the rough surfaces exhibited by the top and bottom walls.

Figure 3 .
Figure 3. Simplified sketch of the cross-section of a PCB SIW illustrating the rough surfaces exhibited by the top and bottom walls.

Figure 4 .
Figure 4.  wave obtained from EM simulated data for a RWG assuming PCB materials considering smooth and rough conductor conditions: (a) cross-section of the simulated structure, (b) real part, and (c) imaginary part.

Figure 4 .
Figure 4. Z wave obtained from EM simulated data for a RWG assuming PCB materials considering smooth and rough conductor conditions: (a) cross-section of the simulated structure, (b) real part, and (c) imaginary part.

Figure 5 .
Figure 5. Photographs showing a perspective view of the prototype, the details of the microstrip-to-SIW transition, and a zoomed-in image of the probing pads.

Figure 5 .
Figure 5. Photographs showing a perspective view of the prototype, the details of the microstrip-to-SIW transition, and a zoomed-in image of the probing pads.

Micromachines 2024 , 13 Figure 6 .
Figure 6.Permittivity and loss tangent determination using a causal model.Now, employing Equation (17) under the same assumptions explained in[28], the loss tangent is approximately given by:

Figure 6 .
Figure 6.Permittivity and loss tangent determination using a causal model.

Figure 7 .
Figure 7. Verification of suitability of using Equation (20): (a) α from an EM simulation correlated with experimental data.(b) Comparison between KR = R/Rsmooth obtained from the correlated simulation and from the causal model expressed by Equation (19).

Figure 7 .
Figure 7. Verification of suitability of using Equation (20): (a) α from an EM simulation correlated with experimental data.(b) Comparison between K R = R/R smooth obtained from the correlated simulation and from the causal model expressed by Equation (19).

Figure 7 .
Figure 7. Verification of suitability of using Equation (20): (a) α from an EM simulation correlated with experimental data.(b) Comparison between KR = R/Rsmooth obtained from the correlated simulation and from the causal model expressed by Equation (19).

Figure 8 .Figure 8 .Figure 9 .
Figure 8. Zwave curves calculated from experimental α and β data using different approaches.EM simulations accounting for the surface roughness effect on the SIW top and bottom walls are also depicted for comparison: (a) real part, and (b) imaginary part.The curves correspond to the SIW with w = 4.1 mm.

Figure 9 .
Figure 9. Z wave curves calculated from experimental α and β data using different approaches.EM simulations accounting for the surface roughness effect on the SIW top and bottom walls are also depicted for comparison: (a) real part, and (b) imaginary part.The curves correspond to the SIW with w = 2.3 mm.

Figure 10 .
Figure 10.Curves for the circuit elements associated with the higher losses in the model shown in Figure 2b: (a) R, and (b) G'.The cases where Zwave was calculated with simplified approaches are also shown.

Figure 10 .
Figure 10.Curves for the circuit elements associated with the higher losses in the model shown in Figure 2b: (a) R, and (b) G'.The cases where Z wave was calculated with simplified approaches are also shown.

Figure 11 .
Figure 11.Model-experiment correlation for the insertion loss obtained through a reconstruction of S-parameters using γ and Zwave for the SIW with: (a) w = 4.1 mm and (b) w = 2.3 mm; l = 15 mm was assumed for both cases.For the calculation of S21, the S-parameters are referenced to Zwave for minimizing the return loss.