The Choice of Line Lengths in Multiline Thru-Reflect-Line Calibration
Pith reviewed 2026-05-16 21:06 UTC · model grok-4.3
The pith
Optimal line lengths for multiline TRL calibration are selected by nonlinear constrained optimization of the eigenvalue problem to spread uncertainty evenly across frequencies.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The optimal lengths of the line standards are obtained through nonlinear constrained optimization of the eigenvalue problem that underlies multiline TRL calibration; a practical shortcut uses sparse rulers to select near-optimal sets. When these lengths are implemented, measurement-based Monte Carlo analysis that employs error boxes extracted from impedance-standard-substrate data shows that calibration uncertainty is distributed more evenly across the lines than occurs with a commercial calibration kit.
What carries the argument
Nonlinear constrained optimization of the multiline TRL eigenvalue problem, with sparse rulers as a fallback length-selection rule.
If this is right
- The required number of lines for a target bandwidth follows directly from the optimization constraints.
- The same length-selection procedure applies without change to lossy and dispersive transmission lines.
- Banded length solutions exist for waveguide calibrations.
- The method has been demonstrated on printed-circuit-board kits of varied materials and stackups up to 150 GHz.
Where Pith is reading between the lines
- The even distribution of uncertainty could allow engineers to cover the same bandwidth with fewer physical line standards than current practice requires.
- The sparse-ruler shortcut offers a quick manual design rule that avoids the need for numerical optimization software in routine work.
- The approach may be adapted to other multi-standard VNA calibration techniques that rely on eigenvalue solutions.
Load-bearing premise
The error boxes obtained from impedance standard substrate measurements capture the dominant uncertainty sources present in the multiline TRL measurements.
What would settle it
If a Monte Carlo simulation rerun with the same error boxes but using the commercial-kit lengths instead of the optimized lengths yields lower or equal peak uncertainty values across the band, the claim of more even distribution would be falsified.
Figures
read the original abstract
This paper presents an analysis and rigorous procedure for determining the optimal lengths of line standards in multiline thru-reflect-line (TRL) calibration of vector network analyzers (VNAs). The solution is obtained through nonlinear constrained optimization of the eigenvalue problem in multiline TRL calibration. Additionally, we propose a simplified approach for near-optimal length selection based on predefined sparse rulers. Alongside the length calculation, we discuss the required number of lines to meet bandwidth requirements. The proposed methods are validated through measurements of multiple multiline TRL calibration kits on printed circuit boards of different materials and stackups, covering frequencies up to 150 GHz. A measurement-based Monte Carlo uncertainty analysis, using error boxes derived from impedance standard substrate measurements, demonstrates that the proposed line lengths distribute calibration uncertainty more evenly across lines compared to a commercial calibration kit. Practical examples are provided for various applications, including lossy and dispersive lines, as well as banded solutions for waveguides.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper claims to provide a nonlinear constrained optimization method for choosing optimal line lengths in multiline TRL calibration by optimizing the eigenvalue problem, along with a sparse ruler based heuristic for near-optimal selection. It addresses the number of lines required for sufficient bandwidth and validates the approach with measurements on various PCB kits up to 150 GHz. A Monte Carlo uncertainty analysis using error boxes from impedance standard substrate measurements is used to show that the proposed lengths distribute calibration uncertainty more evenly than a commercial calibration kit. Examples for lossy, dispersive lines and waveguides are included.
Significance. If the uncertainty model holds, the result would be significant for improving VNA calibration accuracy in high-frequency PCB applications. The measurement validation across multiple materials and the Monte Carlo analysis provide empirical grounding, while the heuristic offers a practical alternative to full optimization; these elements could help standardize better line selection practices.
major comments (1)
- [Monte Carlo uncertainty analysis] Monte Carlo uncertainty analysis: the demonstration that proposed lengths distribute uncertainty more evenly than the commercial kit rests on error boxes derived from ISS measurements as the uncertainty model. The manuscript does not show that these boxes capture all dominant sources (probe contact repeatability, fabrication tolerances, substrate variations) up to 150 GHz across the tested materials, which directly affects the strength of the empirical superiority claim.
minor comments (2)
- [Abstract] The abstract refers to a 'rigorous procedure' without naming the objective function or constraints; adding one sentence on the optimization formulation would improve clarity.
- [Results figures] Figures comparing uncertainty distributions should explicitly label the proposed lengths versus the commercial kit and include error bars from the Monte Carlo runs for direct visual assessment.
Simulated Author's Rebuttal
We thank the referee for the detailed review and constructive comments on our manuscript. We provide a point-by-point response to the major comment below.
read point-by-point responses
-
Referee: [Monte Carlo uncertainty analysis] Monte Carlo uncertainty analysis: the demonstration that proposed lengths distribute uncertainty more evenly than the commercial kit rests on error boxes derived from ISS measurements as the uncertainty model. The manuscript does not show that these boxes capture all dominant sources (probe contact repeatability, fabrication tolerances, substrate variations) up to 150 GHz across the tested materials, which directly affects the strength of the empirical superiority claim.
Authors: We appreciate the referee's observation regarding the scope of the uncertainty model. The Monte Carlo analysis employs error boxes extracted from repeated ISS measurements to model the calibration uncertainty, following established practices in the literature for VNA calibration studies. This approach captures the combined effects observed in the measurement setup, including contributions from probe positioning and contact variations during the calibration process itself. Nevertheless, we acknowledge that explicitly quantifying every possible source, such as long-term fabrication tolerances or material variations across different batches, would require additional extensive experiments beyond the scope of this work. In the revised manuscript, we will expand the discussion in the uncertainty analysis section to explicitly state the assumptions underlying the error model and qualify the comparison as being valid under the measured variability from the ISS standards. This will ensure the empirical results are presented with appropriate caveats, while preserving the demonstration that the optimized lengths provide a more uniform uncertainty distribution within this framework. We believe this addresses the concern without necessitating changes to the optimization procedure or the measurement data. revision: partial
Circularity Check
Optimization on standard multiline TRL eigenvalue problem; Monte Carlo uses independent ISS error boxes
full rationale
The derivation applies nonlinear constrained optimization directly to the established multiline TRL eigenvalue problem to select line lengths, then validates via Monte Carlo sampling of error boxes obtained from separate ISS measurements. Neither the optimization criterion nor the uncertainty propagation reduces to a fitted parameter defined by the target distribution itself, nor does any load-bearing step collapse to a self-citation chain or ansatz smuggled from prior author work. The result remains self-contained against the standard TRL formulation and external measurement data.
Axiom & Free-Parameter Ledger
axioms (1)
- domain assumption The multiline TRL calibration can be expressed as a well-posed eigenvalue problem whose solution depends on line lengths.
Reference graph
Works this paper leans on
-
[1]
G. Engen and C. Hoer, “Thru-reflect-line: An improved technique for calibrating the dual six-port automatic network analyzer,”IEEE Transactions on Microwave Theory and Techniques, vol. 27, no. 12, pp. 987–993, 1979, doi: 10.1109/TMTT.1979.1129778
-
[2]
Uncertainty of vna s-parameter measurement due to nonideal trl calibration items,
U. Stumper, “Uncertainty of vna s-parameter measurement due to nonideal trl calibration items,”IEEE Transactions on Instrumentation and Measurement, vol. 54, no. 2, pp. 676–679, 2005
work page 2005
-
[3]
A multiline method of network analyzer calibration,
R. Marks, “A multiline method of network analyzer calibration,”IEEE Transactions on Microwave Theory and Techniques, vol. 39, no. 7, pp. 1205–1215, 1991, doi: 10.1109/22.85388
-
[4]
Improving the reliability of the multiline trl calibration algorithm,
Z. Hatab, M. Gadringer, and W. B ¨osch, “Improving the reliability of the multiline trl calibration algorithm,” in2022 98th ARFTG Mi- crowave Measurement Conference (ARFTG). IEEE, Jan. 2022, doi: 10.1109/ARFTG52954.2022.9844064. pp. 1–5
-
[5]
TRL–calibration — microwaves101.com,
Microwaves101, “TRL–calibration — microwaves101.com,” https:// www.microwaves101.com/encyclopedias/trl-calibration, [Accessed 15- 01-2024]
work page 2024
-
[6]
N. M. Ridler, “Choosing line lengths for calibrating waveguide vector network analysers at millimetre and sub-millimetre wavelengths.” NPL Report, March 2009. [Online]. Available: http://eprintspublications.npl. co.uk/4346/
work page 2009
-
[7]
Strategies for traceable submillimeter-wave vector network analyzer,
N. M. Ridler, R. G. Clarke, C. Li, and M. J. Salter, “Strategies for traceable submillimeter-wave vector network analyzer,”IEEE Transac- tions on Terahertz Science and Technology, vol. 9, no. 4, pp. 392–398, Jul. 2019, doi: 10.1109/TTHZ.2019.2911870
-
[8]
N. M. Ridler, S. Johny, M. J. Salter, X. Shang, W. Sun, and A. Wilson, “Establishing waveguide lines as primary standards for scattering param- eter measurements at submillimetre wavelengths,”Metrologia, vol. 58, no. 1, p. 015015, jan 2021, doi: 10.1088/1681-7575/abd371
-
[9]
H.-J. Eul and B. Schiek, “Thru-match-reflect: One result of a rig- orous theory for de-embedding and network analyzer calibration,” in 1988 18th European Microwave Conference. IEEE, Sep. 1988, doi: 10.1109/euma.1988.333924. pp. 909–914
-
[10]
D. DeGroot, J. Jargon, and R. Marks, “Multiline trl revealed,” in60th ARFTG Conference Digest, Fall 2002., ser. ARFTGF-02. IEEE, 2002, doi: 10.1109/ARFTGF.2002.1218696. pp. 131–155
-
[11]
An optimal multiline trl cal- ibration algorithm,
D. Williams, C. Wang, and U. Arz, “An optimal multiline trl cal- ibration algorithm,” inIEEE MTT-S International Microwave Sym- posium Digest, 2003, ser. MWSYM-03, vol. 3. IEEE, 2003, doi: 10.1109/MWSYM.2003.1210494. pp. 1819–1822 vol.3
-
[12]
Choosing line lengths for calibrating network analyzers,
C. Hoer, “Choosing line lengths for calibrating network analyzers,”IEEE Transactions on Microwave Theory and Techniques, vol. 31, no. 1, pp. 76–78, 1983, doi: 10.1109/TMTT.1983.1131433
-
[13]
Some questions and answers concerning air lines as impedance standards,
C. A. Hoer, “Some questions and answers concerning air lines as impedance standards,” in29th ARFTG Conference Digest, vol. 11, 1987, doi: 10.1109/ARFTG.1987.323862. pp. 161–173
-
[14]
Reliable method for computing the phase shift of multiline lrl calibration technique,
J. Reynoso-Hernandez, “Reliable method for computing the phase shift of multiline lrl calibration technique,”IEEE Microwave and Wireless Components Letters, vol. 12, no. 10, pp. 395–397, Oct. 2002, doi: 10.1109/lmwc.2002.804572
-
[15]
Broadband coaxial s-parameter measurements for cryogenic quantum technologies,
S.-H. Shin, M. Stanley, J. Skinner, S. E. de Graaf, and N. M. Ridler, “Broadband coaxial s-parameter measurements for cryogenic quantum technologies,”IEEE Transactions on Microwave Theory and Techniques, vol. 72, no. 4, pp. 2193–2201, Apr. 2024, doi: 10.1109/tmtt.2023.3322909
-
[16]
On-line accuracy assessment for the dual six-port ana: Extension to nonmating connectors,
C. A. Hoer and G. F. Engen, “On-line accuracy assessment for the dual six-port ana: Extension to nonmating connectors,”IEEE Transactions on Instrumentation and Measurement, vol. IM-36, no. 2, pp. 524–529, 1987, doi: 10.1109/TIM.1987.6312732
-
[17]
A general waveguide circuit theory,
R. Marks and D. Williams, “A general waveguide circuit theory,” Journal of Research (NIST JRES), National Institute of Standards and Technology, Gaithersburg, MD, no. 97, 1992, doi: 10.6028/jres.097.024
-
[18]
Steer,Microwave and RF Design: Transmission Lines, Volume 2
M. Steer,Microwave and RF Design: Transmission Lines, Volume 2. NC State University, Jul. 2019. ISBN 9781469656939
work page 2019
-
[19]
Kronecker products and matrix calculus in system theory,
J. Brewer, “Kronecker products and matrix calculus in system theory,” IEEE Transactions on Circuits and Systems, vol. 25, no. 9, pp. 772–781, Sep. 1978, doi: 10.1109/tcs.1978.1084534
-
[20]
Propagation of linear uncertainties through multiline thru-reflect-line calibration,
Z. Hatab, M. E. Gadringer, and W. B ¨osch, “Propagation of linear uncertainties through multiline thru-reflect-line calibration,”IEEE Trans- actions on Instrumentation and Measurement, vol. 72, pp. 1–9, 2023, doi: 10.1109/TIM.2023.3296123
-
[21]
J. H. Wilkinson,The algebraic eigenvalue problem, reprinted from corr. sheets of the 1. ed. ed., ser. Monographs on numerical analysis. Oxford: Clarendon Press, 1978. ISBN 0198534035
work page 1978
-
[22]
Hogben, Ed.,Handbook of linear algebra, second edition
L. Hogben, Ed.,Handbook of linear algebra, second edition. ed., ser. Discrete mathematics and its applications. Boca Raton: CRC Press, 2014. ISBN 9781466507296 Includes bibliographical references. Description based on print version record. ”A Chapman & Hall book.”
work page 2014
-
[23]
Tell me the method, i’ll give you the mean,
A. Lann and R. Falk, “Tell me the method, i’ll give you the mean,” The American Statistician, vol. 60, no. 4, pp. 322–327, Nov. 2006, doi: 10.1198/000313006x151460
-
[24]
On differentiating eigenvalues and eigenvectors,
J. R. Magnus, “On differentiating eigenvalues and eigenvectors,” Econometric Theory, vol. 1, no. 2, pp. 179–191, 1985. [Online]. Available: http://www.jstor.org/stable/3532409
-
[25]
R. Storn and K. Price, “Differential evolution – a simple and efficient heuristic for global optimization over continuous spaces,”Journal of Global Optimization, vol. 11, no. 4, pp. 341–359, Dec. 1997, doi: 10.1023/a:1008202821328
-
[26]
A constraint handling approach for the differential evolu- tion algorithm,
J. Lampinen, “A constraint handling approach for the differential evolu- tion algorithm,” inProceedings of the 2002 Congress on Evolutionary Computation. CEC’02 (Cat. No.02TH8600), ser. CEC-02, vol. 2. IEEE, 2002, doi: 10.1109/cec.2002.1004459. pp. 1468–1473
-
[27]
Difference bases and sparse sensor arrays,
D. Linebarger, I. Sudborough, and I. Tollis, “Difference bases and sparse sensor arrays,”IEEE Transactions on Information Theory, vol. 39, no. 2, pp. 716–721, Mar. 1993, doi: 10.1109/18.212309
-
[28]
H. Xiao, B. Zhou, Y . Zhang, and G. Xiao, “On the foundation of spar- sity constrained sensing—part ii: Diophantine sampling with arbitrary temporal and spatial sparsity,”IEEE Transactions on Signal Processing, vol. 71, pp. 1277–1292, 2023, doi: 10.1109/tsp.2023.3248493
-
[29]
Compressive co- variance sensing: Structure-based compressive sensing beyond sparsity,
D. Romero, D. D. Ariananda, Z. Tian, and G. Leus, “Compressive co- variance sensing: Structure-based compressive sensing beyond sparsity,” IEEE Signal Processing Magazine, vol. 33, no. 1, pp. 78–93, Jan. 2016, doi: 10.1109/msp.2015.2486805
-
[30]
On the representation of 1, 2, . . . , n by differences,
J. Leech, “On the representation of 1, 2, . . . , n by differences,”Journal of the London Mathematical Society, vol. s1-31, no. 2, pp. 160–169, Apr. 1956, doi: 10.1112/jlms/s1-31.2.160
-
[31]
A note on restricted difference bases,
B. Wichmann, “A note on restricted difference bases,”Journal of the London Mathematical Society, vol. s1-38, no. 1, pp. 465–466, 1963, doi: 10.1112/jlms/s1-38.1.465
-
[32]
Ein satz ueber trigonometrische polynome und seine anwen- dung in der theorie der fourier-reihen,
S. Sidon, “Ein satz ueber trigonometrische polynome und seine anwen- dung in der theorie der fourier-reihen,”Mathematische Annalen, vol. 106, no. 1, pp. 536–539, Dec. 1932, doi: 10.1007/bf01455900. 16
-
[33]
On the complexity of construct- ing golomb rulers,
C. Meyer and P. A. Papakonstantinou, “On the complexity of construct- ing golomb rulers,”Discrete Applied Mathematics, vol. 157, no. 4, pp. 738–748, Feb. 2009, doi: 10.1016/j.dam.2008.07.006
-
[34]
A thru-free multiline calibration,
Z. Hatab, M. E. Gadringer, and W. B ¨osch, “A thru-free multiline calibration,”IEEE Transactions on Instrumentation and Measurement, vol. 72, pp. 1–9, 2023, doi: 10.1109/TIM.2023.3308226
-
[35]
A quasi-tem approach for designing microvias for pcb layer transition with minimal return loss,
Z. Hatab, H. Takahashi, A. B. A. Alterkawi, M. E. Gadringer, and W. B¨osch, “A quasi-tem approach for designing microvias for pcb layer transition with minimal return loss,” in2023 53rd European Microwave Conference (EuMC), 2023, doi: 10.23919/EuMC58039.2023.10290259. pp. 62–65
-
[36]
Miniaturized ka-band metasur- face filter with wide out-of-band rejection up to the 5th harmonic,
A. Arsanjani, A. Abdi, Z. Hatab, A. B. Alothman Alterkawi, M. Ernst Gadringer, and W. B ¨osch, “Miniaturized ka-band metasur- face filter with wide out-of-band rejection up to the 5th harmonic,” in2024 54th European Microwave Conference (EuMC), 2024, doi: 10.23919/EuMC61614.2024.10732666. pp. 188–191
-
[37]
A. Arsanjani, Z. Hatab, A. B. A. Alterkawi, M. E. Gadringer, and W. B¨osch, “Modified semi-additive manufacturing of pcbs for enabling accurate device measurements at millimeter-wave and sub-terahertz fre- quencies,” in2024 103rd ARFTG Microwave Measurement Conference (ARFTG), 2024, doi: 10.1109/ARFTG61196.2024.10660794. pp. 1–4
-
[38]
Impedance standard substrate (iss) 104-783: 75-150 um pitch, up to 145 ghz,
FormFactor Inc., “Impedance standard substrate (iss) 104-783: 75-150 um pitch, up to 145 ghz,” 2023, [Accessed 15-01- 2024]. [Online]. Available: https://www.formfactor.com/product/probes/ calibration-tools-probes/impedance-standard-substrates/
work page 2023
-
[39]
Ac2-2 calibration substrate: 100-250 um pitch, up to 110 ghz,
MPI Corporation, “Ac2-2 calibration substrate: 100-250 um pitch, up to 110 ghz,” 2023, [Accessed 15-01-2024]. [Online]. Avail- able: https://www.mpi-corporation.com/ast/mpi-rf-probes-accessories/ mpi-rf-calibration-substrates/
work page 2023
-
[40]
On the influence of thru- and line-length- related effects in cpw- based multiline trl calibrations,
G. N. Phung and U. Arz, “On the influence of thru- and line-length- related effects in cpw- based multiline trl calibrations,” in2021 97th ARFTG Microwave Measurement Conference (ARFTG), 2021, doi: 10.1109/ARFTG52261.2021.9639909. pp. 1–4
-
[41]
Modeling dispersion and radiation characteristics of conductor-backed cpw with finite ground width,
F. Schnieder, T. Tischler, and W. Heinrich, “Modeling dispersion and radiation characteristics of conductor-backed cpw with finite ground width,”IEEE Transactions on Microwave Theory and Techniques, vol. 51, no. 1, pp. 137–143, Jan. 2003, doi: 10.1109/tmtt.2002.806926
-
[42]
Quasi-tem description of mmic coplanar lines including conductor-loss effects,
W. Heinrich, “Quasi-tem description of mmic coplanar lines including conductor-loss effects,”IEEE Transactions on Microwave Theory and Techniques, vol. 41, no. 1, pp. 45–52, 1993, doi: 10.1109/22.210228
-
[43]
P. Virtanenet al., “SciPy 1.0: Fundamental Algorithms for Scientific Computing in Python,”Nature Methods, vol. 17, pp. 261–272, 2020, doi: 10.1038/s41592-019-0686-2
-
[44]
Project ogr - distributed.net,
distributed.net, “Project ogr - distributed.net,” https://www.distributed. net/OGR, 2024, accessed: 2025-01-15. [Online]. Available: https: //www.distributed.net/OGR
work page 2024
-
[45]
Model of thin-film microstrip line for cir- cuit design,
F. Schnieder and W. Heinrich, “Model of thin-film microstrip line for cir- cuit design,”IEEE Transactions on Microwave Theory and Techniques, vol. 49, no. 1, pp. 104–110, 2001, doi: 10.1109/22.899967
-
[46]
Investigating the impact of final finishes on the insertion loss in as received and after aging,
B. Schafsteller, M. Schwaemmlein, M. Rosin, G. Ramos, Z. Hatab, M. E. Gadringer, and E. Schlaffer, “Investigating the impact of final finishes on the insertion loss in as received and after aging,”IMAPSource Proceed- ings, vol. 2023, no. Symposium, Feb. 2024, doi: 10.4071/001c.94519
-
[47]
Terahertz micromachined on-wafer probes: Re- peatability and reliability,
L. Chenet al., “Terahertz micromachined on-wafer probes: Re- peatability and reliability,”IEEE Transactions on Microwave Theory and Techniques, vol. 60, no. 9, pp. 2894–2902, Sep. 2012, doi: 10.1109/tmtt.2012.2205016
-
[48]
A 1.1 thz micromachined on-wafer probe,
M. F. Bauwens, N. Alijabbari, A. W. Lichtenberger, N. S. Barker, and R. M. Weikle, “A 1.1 thz micromachined on-wafer probe,” in2014 IEEE MTT-S International Microwave Symposium (IMS2014). IEEE, Jun. 2014, doi: 10.1109/mwsym.2014.6848607
-
[49]
Thz characterization and modeling of sige hbts: Review (invited),
S. Fregonese, M. Deng, M. Cabbia, C. Yadav, M. De Matos, and T. Zimmer, “Thz characterization and modeling of sige hbts: Review (invited),”IEEE Journal of the Electron Devices Society, vol. 8, pp. 1363–1372, 2020, doi: 10.1109/jeds.2020.3036135
-
[50]
M. Spirito, U. Arz, G. N. Phung, F. J. Schm ¨uckle, W. Heinrich, and R. Lozar, “Guidelines for the design of calibration substrates, including the suppression of parasitic modes for frequencies up to and including 325 ghz : Empir 14ind02,” 2019, doi: 10.7795/530.20190424A
-
[51]
3-d printed thz waveguide components,
L. Zhu, S.-H. Shin, R. Payapulli, I. W. Rossuck, N. Klein, N. M. Ridler, and S. Lucyszyn, “3-d printed thz waveguide components,”IEEE Access, vol. 11, pp. 79 073–79 086, 2023, doi: 10.1109/access.2023.3297271
-
[52]
“Ieee standard for rectangular metallic waveguides and their interfaces for frequencies of 110 ghz and above—part 1: Frequency bands and waveguide dimensions,”IEEE Std 1785.1-2012, pp. 1–22, 2013, doi: 10.1109/IEEESTD.2013.6471987
-
[53]
Transmission line model for rectangular waveguides accurately incorporating loss effects,
K. Lomakin, G. Gold, and K. Helmreich, “Transmission line model for rectangular waveguides accurately incorporating loss effects,” in2017 IEEE 21st Workshop on Signal and Power Integrity (SPI). IEEE, May 2017, doi: 10.1109/sapiw.2017.7944024. pp. 1–4
-
[54]
Analytical waveguide model precisely predicting loss and delay including surface roughness,
——, “Analytical waveguide model precisely predicting loss and delay including surface roughness,”IEEE Transactions on Microwave The- ory and Techniques, vol. 66, no. 6, pp. 2649–2662, Jun. 2018, doi: 10.1109/tmtt.2018.2827383
-
[55]
T. Needham and R. Penrose,Visual Complex Analysis: 25th Anniver- sary Edition. Oxford University PressOxford, Feb. 2023. ISBN 9780191964947
work page 2023
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.