“…W cl is the 4×4 coupled line (CL) system matrix, and E x1 , H y1 are the EM fields along line 1, while E x2 , H y2 are those along line 2. s , m and µ s , µ m are the self and mutual permittivities and permeabilities, respectively, which can be expressed in terms of CL even-and odd-mode characteristic impedances (Z 0e , Z 0o ), phase velocities (v e , v o ) and effective permittivities ( e eff , o eff ) [105][106][107]. Equations (29a)-(29b) can be solved using the M1-D CL-FDTD method [103,104], which could bypass the fine mesh for line width and spacing of coupled transmission lines. However, due to the stability constraint of such explicit method for CL, there is a time step limit that is usually more restrictive than that for the single uncoupled TL.…”