pyedb 0.54.0__py3-none-any.whl → 0.55.0__py3-none-any.whl

pyedb/libraries/rf_libraries/base_functions.py

@@ -0,0 +1,1358 @@
+# Copyright (C) 2023 - 2024 ANSYS, Inc. and/or its affiliates.
+# SPDX-License-Identifier: MIT
+#
+#
+# Permission is hereby granted, free of charge, to any person obtaining a copy
+# of this software and associated documentation files (the "Software"), to deal
+# in the Software without restriction, including without limitation the rights
+# to use, copy, modify, merge, publish, distribute, sublicense, and/or sell
+# copies of the Software, and to permit persons to whom the Software is
+# furnished to do so, subject to the following conditions:
+#
+# The above copyright notice and this permission notice shall be included in all
+# copies or substantial portions of the Software.
+#
+# THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR
+# IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
+# FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE
+# AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
+# LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING FROM,
+# OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER DEALINGS IN THE
+# SOFTWARE.
+
+from __future__ import annotations
+
+import math
+from typing import List, Optional, Tuple, Union
+
+from pyedb import Edb
+from pyedb.libraries.common import Substrate
+
+
+class HatchGround:
+    """
+    Create a square demo board whose ground layer is made of an
+    orthogonal hatched copper pattern.  Any requested copper fill
+    ratio between 10 % and 90 % can be realised.
+
+    Parameters
+    ----------
+    pitch : float, default 17.07 mm
+        Centre-to-centre distance of the hatch bars [m].
+    width : float, default 5 mm
+        Width of each copper bar [m].
+    fill_target : float, default 50 %
+        Desired copper area in percent (10 … 90).
+    board_size : float, default 100 mm
+        Edge length of the square board [m].
+    layer_gnd : str, default "GND"
+        Name of the layer that receives the hatch pattern.
+
+    Examples
+    --------
+    >>> hatch = HatchGround(pitch=0.5e-3, width=0.2e-3,
+    ...                     fill_target=70, board_size=5e-3)
+    >>> edb = Edb("demo.aedb")
+    >>> hatch._edb = edb
+    >>> hatch.create()
+    >>> round(hatch.copper_fill_ratio, 1)
+    70.0
+    """
+
+    def __init__(
+        self,
+        edb_cell: Optional[Edb] = None,
+        pitch: Union[str, float] = 17.07e-3,
+        width: Union[str, float] = 5.0e-3,
+        fill_target: Union[str, float] = 50.0,
+        board_size: Union[str, float] = 100e-3,
+        layer_gnd: str = "GND",
+    ):
+        """Initialize the hatch ground object."""
+        self._edb = edb_cell
+        self.pitch = self._edb.value(pitch)
+        self.width = self._edb.value(width)
+        self.fill_target = self._edb.value(fill_target)
+        self.board_size = self._edb.value(board_size)
+        self.layer_gnd = layer_gnd
+        self._outline = None
+
+    @property
+    def copper_fill_ratio(self) -> float:
+        """
+        Return the **actual** copper fill ratio in percent.
+
+        Returns
+        -------
+        float
+            Percentage of the board area that is copper after the hatch
+            has been generated.
+        """
+        cu_area = self._edb.modeler.polygons[0].area()
+        return 100.0 * cu_area / (self.board_size**2)
+
+    def _generate_hatch(self) -> None:
+        """Draw orthogonal stripes, then punch gaps for the requested fill."""
+        # ---------- horizontal bars ----------
+        y = 0.0
+        while y < self.board_size:
+            self._add_stripe(0.0, y, self.board_size, y + self.width)
+            y += self.pitch
+
+        # ---------- vertical bars ------------
+        x = 0.0
+        while x < self.board_size:
+            self._add_stripe(x, 0.0, x + self.width, self.board_size)
+            x += self.pitch
+
+        # ---------- punch square gaps --------
+        gaps: List[List[Tuple[float, float]]] = []
+        x = 0.0
+        while x < self.board_size:
+            y = 0.0
+            while y < self.board_size:
+                gaps.append(
+                    [
+                        (x, y),
+                        (x + self.width, y),
+                        (x + self.width, y + self.width),
+                        (x, y + self.width),
+                        (x, y),
+                    ]
+                )
+                y += self.pitch
+            x += self.pitch
+        polygons = self._edb.modeler.polygons
+        polygons[0].unite(polygons[1:])
+
+    def _add_stripe(self, x0: float, y0: float, x1: float, y1: float) -> None:
+        """Create one rectangular copper bar on the GND layer."""
+        points = [(x0, y0), (x1, y0), (x1, y1), (x0, y1), (x0, y0)]
+        self._edb.modeler.create_polygon(points, layer_name=self.layer_gnd, net_name="GND")
+
+    def create(self) -> bool:
+        """
+        Generate the stack-up, board outline and hatch pattern.
+
+        Returns
+        -------
+        bool
+            True when geometry has been created successfully.
+        """
+        self._generate_hatch()
+        return True
+
+
+class Meander:
+    """
+    Fully-parametric micro-strip meander line.
+
+    Parameters
+    ----------
+    pitch : float, default 1 mm
+        Vertical spacing between successive meander rows [m].
+    trace_width : float, default 0.3 mm
+        Width of the micro-strip [m].
+    amplitude : float, default 5 mm
+        Horizontal excursion of each U-turn [m].
+    num_turns : int, default 8
+        Number of 180° bends.
+    layer : str, default "TOP"
+        EDB metal layer.
+    net : str, default "SIG"
+        Net name assigned to the trace.
+
+    Examples
+    --------
+    >>> m = Meander(pitch=0.2e-3, trace_width=0.15e-3,
+    ...             amplitude=2e-3, num_turns=4)
+    >>> edb = Edb("meander.aedb")
+    >>> m._pedb = edb
+    >>> m.create()
+    >>> f"{m.analytical_z0:.1f} Ω"
+    '50.1 Ω'
+    >>> m.electrical_length_deg(1e9)
+    59.8
+    """
+
+    def __init__(
+        self,
+        edb_cell: Edb,
+        pitch: Union[str, float] = 1e-3,
+        trace_width: Union[str, float] = 0.3e-3,
+        amplitude: Union[str, float] = 5e-3,
+        num_turns: int = 8,
+        layer: str = "TOP",
+        net: str = "SIG",
+    ):
+        self._edb = edb_cell
+        self.pitch = self._edb.value(pitch)
+        self.trace_width = self._edb.value(trace_width)
+        self.amplitude = self._edb.value(amplitude)
+        self.num_turns = num_turns
+        self.layer = layer
+        self.net = net
+        self.substrate = Substrate()
+        self.length = 0.0
+
+    # ------------------------------------------------------------------ #
+    # Analytical models
+    # ------------------------------------------------------------------ #
+    @property
+    def analytical_z0(self) -> float:
+        """
+        Micro-strip characteristic impedance using the Hammerstad & Jensen
+        closed-form expression.
+
+        Returns
+        -------
+        float
+            Z0 in Ohm.
+        """
+        return 60 / math.sqrt(self.substrate.er) * math.log(5.98 * 1.6e-3 / (0.8 * self.trace_width + self.trace_width))
+
+    def electrical_length_deg(self, freq: float) -> float:
+        """
+        Electrical length of the meander at the specified frequency.
+
+        Parameters
+        ----------
+        freq : float
+            Frequency in Hz.
+
+        Returns
+        -------
+        float
+            Phase shift in degrees.
+        """
+        c = 299_792_458
+        v = c / math.sqrt(self.substrate.er)
+        beta = 2 * math.pi * freq / v
+        return math.degrees(beta * self.length)
+
+    # ------------------------------------------------------------------ #
+    # EDB creation
+    # ------------------------------------------------------------------ #
+    def create(self) -> bool:
+        """
+        Draw the meander in the attached EDB cell and calculate its
+        physical length.
+
+        Returns
+        -------
+        bool
+            True on success.
+        """
+
+        # Parameters
+        self._edb.add_design_variable("w", self.trace_width)  # trace width
+        self._edb.add_design_variable("p", self.pitch)  # pitch (centre-to-centre)
+        self._edb.add_design_variable("a", self.amplitude)  # meander amplitude
+        self._edb.add_design_variable("n_turns", self.num_turns)  # number of U-turns (integer)
+
+        # ----------------------------------------------------------
+        # 3.  Build the point list
+        # ----------------------------------------------------------
+        pts = [(0, 0)]  # start on the axis
+
+        for i in range(self.num_turns):
+            y = f"{i + 1}*p"  # next row
+            if i % 2 == 0:  # even → left
+                pts.extend([("-a/2", f"{i}*p"), ("-a/2", y), (0, y)])  # step left  # vertical  # back to axis
+            else:  # odd → right
+                pts.extend([(" a/2", f"{i}*p"), (" a/2", y), (0, y)])  # step right  # vertical  # back to axis
+
+        self._edb.modeler.create_trace(
+            path_list=pts, layer_name=self.layer, width=self.trace_width, net_name=self.net, corner_style="Round"
+        )
+        self.length = self._edb.modeler.paths[0].length
+
+
+class MIMCapacitor:
+    """
+    Metal–Insulator–Metal parallel-plate capacitor.
+
+    Parameters
+    ----------
+    area : float, default 0.1 mm²
+        Plate area [m²].
+    gap : float, default 1 µm
+        Dielectric thickness between plates [m].
+    er : float, default 7
+        Relative permittivity of the dielectric.
+    layer_top : str, default "M1"
+        Top plate layer.
+    layer_bottom : str, default "M2"
+        Bottom plate layer.
+    net : str, default "RF"
+        Net name for both plates.
+
+    Examples
+    --------
+    >>> cap = MIMCapacitor(area=200e-12, gap=0.5e-6, er=4.1)
+    >>> edb = Edb("mim.aedb")
+    >>> cap._pedb = edb
+    >>> cap.create()
+    >>> f"{cap.capacitance_f*1e12:.2f} pF"
+    '1.45 pF'
+    """
+
+    def __init__(
+        self,
+        edb_cell: Edb,
+        area: Union[str, float] = 0.1e-6,
+        gap: Union[str, float] = 1e-6,
+        layer_top: str = "M1",
+        layer_bottom: str = "M2",
+        net: str = "RF",
+    ):
+        self._edb = edb_cell
+        self.area = self._edb.value(area)
+        self.gap = self._edb.value(gap)
+        self.layer_top = layer_top
+        self.layer_bottom = layer_bottom
+        self.net = net
+        self.substrate = Substrate()
+
+    @property
+    def capacitance_f(self) -> float:
+        """
+        Analytical parallel-plate capacitance.
+
+        Returns
+        -------
+        float
+            Capacitance in Farads.
+        """
+        eps0 = 8.854e-12
+        return eps0 * self.substrate.er * self.area / self.gap
+
+    def create(self) -> bool:
+        """
+        Create the top plate, bottom plate and assign variables.
+
+        Returns
+        -------
+        bool
+            True on success.
+        """
+        self._edb["area"] = self.area
+        self._edb["gap"] = self.gap
+        side = math.sqrt(self.area)
+        self._edb["side"] = side
+
+        self._edb.modeler.create_rectangle(self.layer_top, self.net, [0, -side / 2, side, side])
+        self._edb.modeler.create_rectangle(self.layer_bottom, self.net, [0, -side / 2, side, side])
+        return True
+
+
+class SpiralInductor:
+    """
+    Square spiral inductor with an optional under-pass bridge.
+
+    Parameters
+    ----------
+    turns : float, default 4.5
+        Number of half-turns (4.5 = 4 full turns + 1 half turn).
+    trace_width : float, default 20 µm
+        Width of the spiral trace.
+    spacing : float, default 12 µm
+        Gap between successive turns.
+    inner_diameter : float, default 60 µm
+        Side length of the innermost square.
+    layer : str, default "M1"
+        Layer on which the spiral is drawn.
+    bridge_layer : str, default "M2"
+        Layer used for the under-pass.
+    via_layer : str, default "M3"
+        Via layer connecting spiral end to the under-pass.
+    net : str, default "IN"
+        Net name.
+    inductor_center : tuple[float, float], default (0, 0)
+        Absolute centre coordinates of the structure.
+    via_size : float, default 25 µm
+        Side length of the square via pad.
+    bridge_width : float, default 12 µm
+        Width of the under-pass trace.
+    bridge_clearance : float, default 6 µm
+        Dielectric clearance under the bridge.
+    bridge_length : float, default 200 µm
+        Length of the under-pass beyond the via.
+    ground_layer : str, default "GND"
+        Layer on which the ground plane is drawn.
+
+    Examples
+    --------
+    >>> sp = SpiralInductor(turns=3.5, trace_width=25e-6,
+    ...                     inner_diameter=80e-6)
+    >>> edb = Edb("spiral.aedb")
+    >>> sp._pedb = edb
+    >>> sp.create()
+    >>> f"{sp.inductance_nh:.1f} nH"
+    '3.4 nH'
+    """
+
+    def __init__(
+        self,
+        edb_cell: Optional[Edb] = None,
+        turns: Union[int, float] = 4.5,
+        trace_width: Union[str, float] = 20e-6,
+        spacing: Union[str, float] = 12e-6,
+        inner_diameter: Union[str, float] = 60e-6,
+        layer: str = "M1",
+        bridge_layer: str = "M2",
+        via_layer: str = "M3",
+        net: str = "IN",
+        inductor_center: Tuple[Union[str, float], Union[str, float]] = (0, 0),  # centre of spiral
+        via_size: Union[str, float] = 25e-6,  # via metal pad
+        bridge_width: Union[str, float] = 12e-6,  # under-pass trace width
+        bridge_clearance: Union[str, float] = 6e-6,  # dielectric gap under bridge
+        bridge_length: Union[str, float] = 200e-6,  # how far the bridge extends
+        ground_layer: str = "GND",
+    ):
+        self._edb = edb_cell
+        self.turns = turns  # half-turns → 4.5 = 4 full + 1 half
+        self.trace_width = self._edb.value(trace_width)
+        self.spacing = self._edb.value(spacing)
+        self.inner_diameter = self._edb.value(inner_diameter)  # first inner square side
+        self.via_size = self._edb.value(via_size)  # centre via finished hole
+        self.inductor_center = [self._edb.value(v) for v in inductor_center]  # via centre position
+        self.bridge_width = self._edb.value(bridge_width)  # under-pass trace width
+        self.bridge_clearance = self._edb.value(bridge_clearance)  # dielectric gap under bridge
+        self.bridge_length = self._edb.value(bridge_length)  # how far the bridge extends
+        self.layer = layer
+        self.bridge_layer = bridge_layer
+        self.via_layer = via_layer  # layer for the centre via
+        self.substrate = Substrate()  # default substrate
+        self.net_name = net
+        self.ground_layer = ground_layer
+
+    @property
+    def inductance_nh(self) -> float:
+        """
+        Accurate inductance calculated with the improved Wheeler formula
+        for square spirals.
+
+        Returns
+        -------
+        float
+            Inductance in nano-Henries.
+        """
+        w = self.trace_width
+        s = self.spacing
+        N = self.turns
+        d_in = self.inner_diameter
+        d_out = d_in + 4 * N * (w + s)  # outer side length
+
+        # 2. Parameters for the improved Wheeler formula (square)
+        d_avg = (d_out + d_in) / 2.0
+        rho = (d_out - d_in) / (d_out + d_in)
+
+        C1, C2, C3, C4 = 1.27, 2.07, 0.18, 0.13  # square coefficients
+        mu0 = 4e-7 * math.pi  # H·m⁻¹
+
+        # 3. Wheeler formula in SI
+        L_h = (mu0 * N**2 * d_avg * C1 / 2.0) * (math.log(C2 / rho) + C3 * rho**2 + C4 * rho**3)
+
+        return L_h * 1e9  # → nH
+
+    def create(self) -> bool:
+        cx = self.inductor_center[0]
+        cy = self.inductor_center[1]
+        pts = []
+        x, y = cx, cy
+        side = self.inner_diameter
+        direction = 0  # 0=+y, 1=-x, 2=-y, 3=+x (vertical first)
+
+        for i in range(int(math.ceil(2 * self.turns))):
+            # choose direction and length
+            if direction == 0:
+                dx, dy = 0, side / 2
+            elif direction == 1:
+                dx, dy = -side / 2, 0
+            elif direction == 2:
+                dx, dy = 0, -side / 2
+            else:
+                dx, dy = side / 2, 0
+
+            x += dx
+            y += dy
+            pts.append((x, y))
+
+            direction = (direction + 1) % 4
+            if i % 2 == 1:  # every two sides increase pitch
+                side += 2 * (self.trace_width + self.spacing)
+
+        via_pos = pts[0]
+        bridge_dir = -1  # -1 = down, +1 = up
+        bridge_end = (via_pos[0] + bridge_dir * self.bridge_length, via_pos[1])
+
+        self._edb.modeler.create_trace(
+            path_list=pts,
+            layer_name=self.layer,
+            width=self.trace_width,
+            net_name=self.net_name,
+            start_cap_style="Flat",
+            end_cap_style="Flat",
+        )
+
+        # Centre via (spiral inner end → bottom metal)
+        self._edb.modeler.create_rectangle(
+            layer_name=self.via_layer,
+            center_point=pts[0],
+            width=self.via_size,
+            height=self.via_size,
+            representation_type="CenterWidthHeight",
+            net_name=self.net_name,
+        )
+
+        # bridge trace on bottom metal
+        self._edb.modeler.create_trace(
+            path_list=[via_pos, bridge_end],
+            layer_name=self.bridge_layer,
+            width=self.bridge_width,
+            net_name=self.net_name,
+            start_cap_style="Flat",
+            end_cap_style="Flat",
+        )
+
+        # ground plane
+        self._edb.modeler.create_rectangle(
+            layer_name=self.ground_layer,
+            center_point=self.inductor_center,
+            width=self.substrate.size[0],
+            height=self.substrate.size[1],
+            representation_type="CenterWidthHeight",
+            net_name="ref",
+        )
+        return True
+
+
+class CPW:
+    """
+    Coplanar waveguide with side ground planes.
+
+    Parameters
+    ----------
+    length : float, default 1 mm
+        Physical length of the line.
+    width : float, default 0.3 mm
+        Width of the centre conductor.
+    gap : float, default 0.1 mm
+        Gap between centre conductor and ground planes.
+    layer : str, default "TOP"
+        Layer on which the CPW is drawn.
+    ground_net : str, default "GND"
+        Net name for the ground planes.
+    ground_width : float, default 0.1 mm
+        Width of the side ground strips.
+    ground_layer : str, default "GND"
+        Layer for the underlying ground plane.
+    net : str, default "SIG"
+        Net name for the centre conductor.
+    substrate : Substrate, default 100 µm FR4
+        Substrate definition.
+
+    Examples
+    --------
+    >>> cpw = CPW(length=5e-3, width=0.4e-3, gap=0.2e-3)
+    >>> edb = Edb("cpw.aedb")
+    >>> cpw._pedb = edb
+    >>> cpw.create()
+    >>> f"{cpw.analytical_z0:.1f} Ω"
+    '46.5 Ω'
+    """
+
+    def __init__(
+        self,
+        edb_cell: Optional[Edb] = None,
+        length: Union[str, float] = 1e-3,
+        width: Union[str, float] = 0.3e-3,
+        gap: Union[str, float] = 0.1e-3,
+        layer: str = "TOP",
+        ground_net: str = "GND",
+        ground_width: float = 0.1e-3,
+        ground_layer: str = "GND",
+        net: str = "SIG",
+        substrate: Substrate = Substrate(100e-6, 4.4, 0.02, "SUB", (0.001, 0.001)),
+    ):
+        self._edb = edb_cell
+        self.length = self._edb.value(length)
+        self.width = self._edb.value(width)
+        self.gap = self._edb.value(gap)
+        self.layer = layer
+        self.ground_net = ground_net
+        self.ground_width = ground_width  # width of ground plane
+        self.ground_layer = ground_layer  # layer for the ground plane
+        self.net = net
+        self.substrate = substrate
+
+    @property
+    def analytical_z0(self) -> float:
+        """
+        Characteristic impedance obtained with the conformal-mapping
+        formula for CPW.
+
+        Returns
+        -------
+        float
+            Z0 in Ohm.
+        """
+        a = self.width / 2
+        b = a + self.gap
+        k = a / b
+        kpr = math.sqrt(1 - k**2)
+        # Complete elliptic integral ratio approx.
+        k_ratio = math.pi / (2 * math.log(2 * (1 + math.sqrt(kpr)) / (1 - math.sqrt(kpr))))
+        return 30 * math.pi / math.sqrt(self.substrate.er) * k_ratio
+
+    def create(self) -> bool:
+        """
+        Draw the centre strip, side grounds and bottom ground plane.
+
+        Returns
+        -------
+        bool
+            True on success.
+        """
+        self._edb["l"] = self.length
+        self._edb["w"] = self.width
+        self._edb["g"] = self.gap
+
+        # signal
+        self._edb.modeler.create_rectangle(
+            self.layer,
+            net_name=self.net,
+            lower_left_point=[-self.width / 2, 0],
+            upper_right_point=[self.width / 2, self.length],
+        )
+        # grounds
+        self._edb.modeler.create_rectangle(
+            self.layer,
+            net_name=self.ground_net,
+            lower_left_point=[-self.width / 2 - self.gap - self.ground_width, 0],
+            upper_right_point=[-self.width / 2 - self.gap, self.length],
+        )
+        self._edb.modeler.create_rectangle(
+            self.layer,
+            net_name=self.ground_net,
+            lower_left_point=[self.width / 2 + self.gap, 0],
+            upper_right_point=[self.width / 2 + self.gap + self.ground_width, self.length],
+        )
+        self._edb.modeler.create_rectangle(
+            self.ground_layer,
+            lower_left_point=[-self.width / 2 - self.gap - self.ground_width, 0],
+            upper_right_point=[self.width / 2 + self.gap + self.ground_width, self.length],
+        )
+        return True
+
+
+class RadialStub:
+    """
+    Radial (fan) open stub for RF matching.
+
+    Parameters
+    ----------
+    radius : float, default 500 µm
+        Radius of the fan [m].
+    angle_deg : float, default 60°
+        Opening angle of the sector.
+    width : float, default 0.2 mm
+        Width of the feeding micro-strip line [m].
+    layer : str, default "TOP"
+        Metal layer.
+    net : str, default "RF"
+        Net name.
+
+    Examples
+    --------
+    >>> stub = RadialStub(radius=1e-3, angle_deg=90)
+    >>> edb = Edb("radial.aedb")
+    >>> stub._pedb = edb
+    >>> stub.create()
+    >>> f"{stub.electrical_length_deg(2e9):.1f}°"
+    '108.0°'
+    """
+
+    def __init__(
+        self,
+        edb_cell,
+        radius: Union[str, float] = 500e-6,
+        angle_deg: Union[str, float] = 60,
+        width: Union[str, float] = 0.2e-3,
+        layer: str = "TOP",
+        net: str = "RF",
+    ):
+        self._edb = edb_cell
+        self.radius = self._edb.value(radius)
+        self.angle_deg = self._edb.value(angle_deg)
+        self.width = self._edb.value(width)
+        self.layer = layer
+        self.net = net
+        self.substrate: Substrate = Substrate(100e-6, 4.4, 0.02, "SUB", (0.001, 0.001))
+
+    @property
+    def electrical_length_deg(self, freq: float = 2e9) -> float:
+        """
+        Electrical length of the radial stub at a given frequency.
+
+        Parameters
+        ----------
+        freq : float, default 2 GHz
+            Frequency in Hz.
+
+        Returns
+        -------
+        float
+            Phase shift in degrees contributed by the stub.
+        """
+        c = 299_792_458
+        v = c / math.sqrt(self.substrate.er)
+        beta = 2 * math.pi * freq / v
+        return math.degrees(beta * self.radius)
+
+    def create(self) -> bool:
+        """
+        Draw the fan-shaped polygon and the feeding line.
+
+        Returns
+        -------
+        bool
+            True on success.
+        """
+        self._edb["r"] = self.radius
+        self._edb["ang"] = self.angle_deg
+        self._edb["w"] = self.width
+
+        # Create wedge polygon
+        theta = math.radians(self.angle_deg)
+        pts = [
+            [0, 0],
+            [self.radius * math.cos(-theta / 2), self.radius * math.sin(-theta / 2)],
+            [self.radius * math.cos(theta / 2), self.radius * math.sin(theta / 2)],
+        ]
+        self._edb.modeler.create_polygon(main_shape=pts, layer_name=self.layer, net_name=self.net)
+        # feed
+        self._edb.modeler.create_rectangle(
+            layer_name=self.layer,
+            net_name=self.net,
+            lower_left_point=[0, -self.width / 2],
+            upper_right_point=[self.width, self.width / 2],
+        )
+        return True
+
+
+class RatRace:
+    """
+    180° rat-race (ring) hybrid coupler.
+
+    Parameters
+    ----------
+    z0 : float, default 50 Ω
+        Characteristic impedance of the ring.
+    freq : float, default 10 GHz
+        Centre frequency.
+    layer : str, default "TOP"
+        Layer on which the ring is drawn.
+    bottom_layer : str | None
+        Layer for the ground plane (if None, no ground is drawn).
+    net : str, default "RR"
+        Net name.
+    width : float, default 0.2 mm
+        Micro-strip width for the ring and port stubs.
+    nr_segments : int, default 32
+        Number of straight segments per 90° arc.
+
+    Examples
+    --------
+    >>> rr = RatRace(freq=5e9)
+    >>> edb = Edb("ratrace.aedb")
+    >>> rr._pedb = edb
+    >>> rr.create()
+    >>> f"{rr.circumference*1e3:.2f} mm"
+    '45.00 mm'
+    """
+
+    def __init__(
+        self,
+        edb_cell: Optional[Edb] = None,
+        z0: Union[float, str] = 50,
+        freq: Union[float, str] = 10e9,
+        layer: str = "TOP",
+        bottom_layer: Optional[str] = None,
+        net: str = "RR",
+        width: Union[float, str] = 0.2e-3,
+        nr_segments: int = 32,
+    ):
+        self._edb = edb_cell
+        self.z0 = self._edb.value(z0)
+        self.freq = self._edb.value(freq)
+        self.layer = layer
+        self.bottom_layer = bottom_layer
+        self.net = net
+        self.width = self._edb.value(width)
+        self.nr_segments = nr_segments
+        self.substrate = Substrate()
+
+    @property
+    def circumference(self) -> float:
+        """
+        Physical circumference of the ring.
+
+        Returns
+        -------
+        float
+            Circumference in metres (1.5 guided wavelengths).
+        """
+        c = 299_792_458
+        v = c / math.sqrt(self.substrate.er)
+        return 1.5 * v / self.freq
+
+    @property
+    def radius(self) -> float:
+        """
+        Mean radius of the ring.
+
+        Returns
+        -------
+        float
+            Radius in metres.
+        """
+        return self.circumference / (2 * math.pi)
+
+    # ------------------------------------------------------------------
+    # Geometry builders
+    # ------------------------------------------------------------------
+    def _arc_points(
+        self,
+        centre: Tuple[float, float],
+        radius: float,
+        start_angle: float,  # rad
+        delta_angle: float,  # rad
+    ) -> List[Tuple[float, float]]:
+        """Return a list of (x,y) for a discretised arc."""
+        points = []
+        step = delta_angle / self.nr_segments
+        for i in range(self.nr_segments + 1):
+            ang = start_angle + i * step
+            x = centre[0] + radius * math.cos(ang)
+            y = centre[1] + radius * math.sin(ang)
+            points.append((x, y))
+        return points
+
+    def _port_stub(self, start: Tuple[float, float], length: float, angle: float):
+        """Return a two-point list for a straight stub."""
+        dx = length * math.cos(angle)
+        dy = length * math.sin(angle)
+        return [start, (start[0] + dx, start[1] + dy)]
+
+    # ------------------------------------------------------------------
+    # Main creation routine
+    # ------------------------------------------------------------------
+    def create(self) -> bool:
+        """
+        Draw the discretised ring and four 50 Ω port stubs.
+
+        Returns
+        -------
+        bool
+            True on success.
+        """
+        self._edb["c"] = self.circumference
+        self._edb["r"] = self.radius
+        self._edb["w"] = self.width
+
+        r = self.radius
+        w = self.width
+        stub_len = 1e-3  # 1 mm straight sections at the four ports
+
+        # Build the ring as four connected arcs
+        # ----------------------------------------------------------------------
+        # We go counter-clockwise starting at the right-hand port (0°, port 1).
+        #
+        # port 1 (0°) -------- arc A (90°) ------- port 2 (90°)
+        #      |                                      |
+        #      |                                      |
+        #      |             270° arc D               |
+        #      |                                      |
+        # port 4 (270°) ----- arc C (90°) ------ port 3 (180°)
+
+        # Arc A: 0° -> 90°
+        pts_a = self._arc_points((0, 0), r, math.radians(0), math.radians(90))
+
+        # Arc B: 90° -> 180°
+        pts_b = self._arc_points((0, 0), r, math.radians(90), math.radians(90))
+
+        # Arc C: 180° -> 270°
+        pts_c = self._arc_points((0, 0), r, math.radians(180), math.radians(90))
+
+        # Arc D: 270° -> 360° (=0°)  (the long 270° section)
+        pts_d = self._arc_points((0, 0), r, math.radians(270), math.radians(90))
+
+        # Stitch them together into a single closed path
+        ring_points = pts_a + pts_b[1:] + pts_c[1:] + pts_d[1:]
+        self._edb.modeler.create_trace(
+            path_list=ring_points,
+            layer_name=self.layer,
+            net_name=self.net,
+            width=w,
+        )
+
+        # Add the four 50 Ω port stubs
+        # ----------------------------------------------------------------------
+        # Angles of the ports on the ring
+        port_angles = [0, 90, 180, 270]  # degrees
+
+        for idx, ang_deg in enumerate(port_angles, start=1):
+            ang = math.radians(ang_deg)
+            x_ring = r * math.cos(ang)
+            y_ring = r * math.sin(ang)
+
+            # Direction vector pointing outwards
+            stub_pts = self._port_stub((x_ring, y_ring), stub_len, ang)
+            self._edb.modeler.create_trace(
+                path_list=stub_pts,
+                layer_name=self.layer,
+                net_name=f"{self.net}_P{idx}",
+                width=w,
+                start_cap_style="flat",
+                end_cap_style="flat",
+            )
+
+        self._edb.modeler.create_rectangle(
+            lower_left_point=[-self.circumference / 4, -self.circumference / 4],
+            upper_right_point=[self.circumference / 4, self.circumference / 4],
+            layer_name=self.bottom_layer,
+            net_name=self.net,
+        )
+        return True
+
+
+class InterdigitalCapacitor:
+    """
+    Inter-digitated comb capacitor with fully parametric fingers.
+
+    All dimensions are stored as native EDB variables so they remain
+    editable inside AEDT after the library cell is imported.
+
+    Parameters
+    ----------
+    fingers : int, default 8
+        Number of finger pairs.
+    finger_length : str, default "0.9 mm"
+        Length of each finger (string with units).
+    finger_width : str, default "0.08 mm"
+        Width of each finger.
+    gap : str, default "0.04 mm"
+        Gap between adjacent fingers.
+    comb_gap : str, default "0.06 mm"
+        Gap between the two combs at the open end.
+    bus_width : str, default "0.25 mm"
+        Width of the top/bottom bus bars.
+    layer : str, default "TOP"
+        Layer on which the capacitor is drawn.
+    net_a : str, default "PORT1"
+        Net for the bottom comb.
+    net_b : str, default "PORT2"
+        Net for the top comb.
+
+    Examples
+    --------
+    >>> idc = InterdigitalCapacitor(fingers=10,
+    ...                             finger_length="0.5mm",
+    ...                             gap="0.03mm")
+    >>> edb = Edb("idc.aedb")
+    >>> idc._pedb = edb
+    >>> idc.create()
+    >>> f"{idc.capacitance_pf:.2f} pF"
+    '0.74 pF'
+    """
+
+    VAR_PREFIX = "IDC"  # prefix for every EDB variable
+
+    def __init__(
+        self,
+        edb_cell: Optional[Edb] = None,
+        fingers: int = 8,
+        finger_length: Union[float, str] = "0.9mm",
+        finger_width: Union[float, str] = "0.08mm",
+        gap: Union[float, str] = "0.04mm",
+        comb_gap: Union[float, str] = "0.06mm",
+        bus_width: Union[float, str] = "0.25mm",
+        layer: str = "TOP",
+        net_a: str = "PORT1",
+        net_b: str = "PORT2",
+    ):
+        self._edb = edb_cell
+        self.layer = layer
+        self.net_a = net_a
+        self.net_b = net_b
+        self.substrate = Substrate()
+
+        pfx = self.VAR_PREFIX
+        self._vars = {
+            "fingers": int(fingers),
+            "finger_length": self._edb.value(finger_length),
+            "finger_width": self._edb.value(finger_width),
+            "gap": self._edb.value(gap),
+            "comb_gap": self._edb.value(comb_gap),
+            "bus_width": self._edb.value(bus_width),
+        }
+        for k, v in self._vars.items():
+            var_name = f"{pfx}_{k}"
+            # Use design variable (no $ prefix) -> stored in cell
+            self._edb[var_name] = v
+
+    @property
+    def capacitance_pf(self) -> float:
+        """
+        Quick parallel-plate estimate of the total capacitance.
+
+        Returns
+        -------
+        float
+            Capacitance in pico-Farads.
+        """
+        pfx = self.VAR_PREFIX
+        eps0 = 8.854e-12
+        er = self.substrate.er
+        N = self._edb[f"{pfx}_fingers"]
+        L = self._edb[f"{pfx}_finger_length"]
+        W = self._edb[f"{pfx}_finger_width"]
+        g = self._edb[f"{pfx}_gap"]
+        return (eps0 * er * N * L * W / g) * 1e12
+
+    def create(self) -> bool:
+        """
+        Draw the two bus bars and interleaved fingers.
+
+        Returns
+        -------
+        bool
+            True on success.
+        """
+        pfx = self.VAR_PREFIX
+        pitch = f"({pfx}_finger_width + {pfx}_gap)"
+        total_width = f"(2*{pfx}_bus_width + (2*{pfx}_fingers-1)*{pitch} + {pfx}_finger_width)"
+        total_height = f"(2*{pfx}_bus_width + 2*{pfx}_finger_length + {pfx}_comb_gap)"
+
+        # 1. Bottom bus bar (NET_A)
+        self._edb.modeler.create_rectangle(
+            layer_name=self.layer,
+            net_name=self.net_a,
+            lower_left_point=["0", "0"],
+            upper_right_point=[total_width, f"{pfx}_bus_width"],
+        )
+
+        # 2. Top bus bar (NET_B)
+        right_bar_y = f"({total_height} - {pfx}_bus_width)"
+        self._edb.modeler.create_rectangle(
+            layer_name=self.layer,
+            net_name=self.net_b,
+            lower_left_point=["0", right_bar_y],
+            upper_right_point=[total_width, total_height],
+        )
+
+        # 3. Fingers (interleaved)
+        y_a = f"{pfx}_bus_width"  # base of up fingers
+        y_b = f"({total_height} - {pfx}_bus_width - {pfx}_finger_length)"  # base of down fingers
+
+        for i in range(self._vars["fingers"] * 2):
+            x = f"({pfx}_bus_width + {i}*{pitch})"
+            if i % 2 == 0:  # NET_A finger (points up)
+                self._edb.modeler.create_rectangle(
+                    layer_name=self.layer,
+                    net_name=self.net_a,
+                    lower_left_point=[x, y_a],
+                    upper_right_point=[f"{x} + {pfx}_finger_width", f"{y_a} + {pfx}_finger_length"],
+                )
+            else:  # NET_B finger (points down)
+                self._edb.modeler.create_rectangle(
+                    layer_name=self.layer,
+                    net_name=self.net_b,
+                    lower_left_point=[x, y_b],
+                    upper_right_point=[f"{x} + {pfx}_finger_width", f"{y_b} + {pfx}_finger_length"],
+                )
+
+        return True
+
+
+class DifferentialTLine:
+    """
+    Edge-coupled differential pair with fully parametric geometry.
+
+    Parameters
+    ----------
+    length : float, default 10 mm
+        Total length of the pair.
+    width : float, default 0.2 mm
+        Width of each individual trace.
+    spacing : float, default 0.2 mm
+        Edge-to-edge separation between the traces.
+    x0 : float, default 0
+        Start x-coordinate (metres).
+    y0 : float, default 0
+        Start y-coordinate (metres).
+    angle : float, default 0 rad
+        Rotation angle of the pair (radians, CCW).
+    layer : str, default "TOP"
+        Layer on which the traces are drawn.
+    net_p : str, default "P"
+        Net name for the positive trace.
+    net_n : str, default "N"
+        Net name for the negative trace.
+
+    Examples
+    --------
+    >>> diff = DifferentialTLine(Edb("diff.aedb"),
+    ...                          length=5e-3,
+    ...                          width=0.15e-3,
+    ...                          spacing=0.1e-3,
+    ...                          angle=math.pi/4)
+    >>> traces = diff.create()
+    >>> f"{diff.diff_impedance:.1f} Ω"
+    '95.6 Ω'
+    """
+
+    def __init__(
+        self,
+        edb: Edb,
+        length: Union[float, str] = 10e-3,
+        width: Union[float, str] = 0.20e-3,
+        spacing: Union[float, str] = 0.20e-3,
+        x0: Union[float, str] = 0.0,
+        y0: Union[float, str] = 0.0,
+        angle: Union[float, str] = 0.0,
+        layer: str = "TOP",
+        net_p: str = "P",
+        net_n: str = "N",
+    ):
+        self._edb = edb
+        self.length = self._edb.value(length)
+        self.width = self._edb.value(width)
+        self.spacing = self._edb.value(spacing)
+        self.x0 = self._edb.value(x0)
+        self.y0 = self._edb.value(y0)
+        self.layer = layer
+        self.net_p = net_p
+        self.net_n = net_n
+
+        self._edb["diff_len"] = self.length  # total length
+        self._edb["diff_w"] = self.width  # single trace width
+        self._edb["diff_s"] = self.spacing  # edge-to-edge spacing
+        self._edb["diff_x0"] = self.x0  # start-x
+        self._edb["diff_y0"] = self.y0  # start-y
+        self._edb["diff_angle"] = math.degrees(angle)  # rotation in degrees
+
+    @property
+    def diff_impedance(self) -> float:
+        """
+        Rough odd-mode impedance estimate for the differential pair.
+
+        Returns
+        -------
+        float
+            Differential impedance in Ohms.
+        """
+        w = self._edb["diff_w"]
+        s = self._edb["diff_s"]
+        z0_single = 60.0
+        return 2 * z0_single * (1 - 0.48 * math.exp(-0.96 * s / w))
+
+    def create(self) -> List[float]:
+        """
+        Create the two traces using only parameter strings so the
+        geometry stays fully editable in AEDT.
+
+        Returns
+        -------
+        list[float]
+            EDB object IDs of the positive and negative traces.
+        """
+        pos_trace = self._edb.modeler.create_trace(
+            path_list=[
+                ["diff_x0", "diff_y0"],
+                ["diff_x0 + diff_len*cos(diff_angle*1deg)", "diff_y0 + diff_len*sin(diff_angle*1deg)"],
+            ],
+            layer_name=self.layer,
+            width="diff_w",
+            net_name=self.net_p,
+            start_cap_style="Flat",
+            end_cap_style="Flat",
+        )
+
+        neg_y_expr = "diff_y0 - diff_w - diff_s"
+        neg_trace = self._edb.modeler.create_trace(
+            path_list=[
+                ["diff_x0", neg_y_expr],
+                ["diff_x0 + diff_len*cos(diff_angle*1deg)", f"{neg_y_expr} + diff_len*sin(diff_angle*1deg)"],
+            ],
+            layer_name=self.layer,
+            width="diff_w",
+            net_name=self.net_n,
+            start_cap_style="Flat",
+            end_cap_style="Flat",
+        )
+        return [pos_trace, neg_trace]
+
+    import math
+    from typing import List, Optional
+
+
+class MicroStripLine:
+    """
+    Single-ended micro-strip line with fixed geometry and editable stack-up.
+
+    Parameters
+    ----------
+    edb : object
+        EDB application / API handle.
+    layer : str
+        Conductor layer on which the trace is drawn.
+    net : str
+        Net name assigned to the trace.
+    x0, y0 : float
+        Global start coordinates (mm, m, or whatever length unit is active in AEDT).
+    length : float
+        Physical length of the line (same unit as x0/y0).
+    angle : float
+        Rotation angle in degrees (0° = along +X).
+    freq : float, optional
+        Reference frequency (Hz) used to compute the electrical length.
+        If None, electrical_length will be None.
+    """
+
+    def __init__(
+        self,
+        edb_cell,
+        layer: str,
+        net: str,
+        x0: Union[float, str] = 0.0,
+        y0: Union[float, str] = 0.0,
+        length: Union[float, str] = "1mm",
+        width: Union[float, str] = "0.2mm",
+        angle: Union[float, str] = 0.0,
+        freq: Optional[Union[float, str]] = None,
+    ):
+        self._edb = edb_cell
+        self.layer = layer
+        self.net = net
+        self.x0 = self._edb.value(x0)
+        self.y0 = self._edb.value(y0)
+        self._length = self._edb.value(length)
+        self._width = self._edb.value(width)  # width of the trace
+        self.angle = self._edb.value(angle)  # degrees
+        self.freq = freq
+        self.substrate = Substrate()
+        self.id = None
+
+        # Edb variable
+        self._edb["w"] = self.width
+        self._edb["l"] = self.length
+
+    @property
+    def width(self):
+        return self._width
+
+    @width.setter
+    def width(self, value: Union[float, str]):
+        """Set the trace width and update the EDB variable."""
+        self._width = self._edb.value(value)
+        self._edb["w"] = self._width
+
+    @property
+    def length(self) -> float:
+        """Physical length of the micro-strip line."""
+        return self._length
+
+    @length.setter
+    def length(self, value: Union[float, str]):
+        """Set the trace length and update the EDB variable."""
+        self._length = self._edb.value(value)
+        self._edb["l"] = self._length
+
+    @property
+    def _ereff(self) -> float:
+        w = self._edb["w"]
+        er = self.substrate.er
+        h = self.substrate.h
+
+        u = w / h
+        a = (
+            1
+            + (1 / 49) * math.log((u**4 + (u / 52) ** 2) / (u**4 + 0.432))
+            + (1 / 18.7) * math.log(1 + (u / 18.1) ** 3)
+        )
+        b = 0.564 * ((er - 0.9) / (er + 3)) ** 0.053
+        return (er + 1) / 2 + (er - 1) / 2 * (1 + 10 / u) ** (-a * b)
+
+    # ------------------------------------------------------------------
+    # Electrical properties
+    # ------------------------------------------------------------------
+    @property
+    def impedance(self) -> float:
+        """
+        Accurate characteristic impedance (Ω) for a micro-strip line.
+        Valid for 0.05 ≤ w/h ≤ 1000 and 1 ≤ εr ≤ 128.
+        Source: G. Wheeler, E. Hammerstad & Jensen, IEEE-MTT 1980.
+        """
+        w = self._edb["w"]
+        h = self.substrate.h
+        er = self.substrate.er
+
+        u = w / h
+
+        # Effective permittivity (εeff)
+        eta0 = 376.73  # Free-space impedance, Ω
+        a0 = (
+            1
+            + (1 / 49) * math.log((u**4 + (u / 52) ** 2) / (u**4 + 0.432))
+            + (1 / 18.7) * math.log(1 + (u / 18.1) ** 3)
+        )
+        b0 = 0.564 * ((er - 0.9) / (er + 3)) ** 0.053
+        er_eff = (er + 1) / 2 + (er - 1) / 2 * (1 + 10 / u) ** (-a0 * b0)
+
+        # Impedance (Z0)
+        # Two regions: u ≤ 1 and u > 1
+        if u <= 1:
+            z0 = eta0 / (2 * math.pi * math.sqrt(er_eff + 1.41)) * math.log(8 / u + u / 4)
+        else:
+            f = 6 + (2 * math.pi - 6) * math.exp(-((30.666 / u) ** 0.7528))
+            z0 = eta0 / (math.sqrt(er_eff) * (u + 1.98 * u**0.172))
+
+        return round(z0, 2)
+
+    @property
+    def electrical_length(self) -> Optional[float]:
+        """
+        Electrical length in degrees at self.freq.
+        Returns None if freq is None.
+        """
+        if self.freq is None:
+            return None
+        c0 = 299_792_458  # m/s
+        beta_l = 2 * math.pi * self.freq * math.sqrt(self._ereff) / c0 * self.length
+        return round(math.degrees(beta_l), 3)
+
+    # ------------------------------------------------------------------
+    # Geometry creation
+    # ------------------------------------------------------------------
+    def create(self) -> List[float]:
+        """Create the trace and return its EDB object."""
+        angle_rad = math.radians(self.angle)
+        trace = self._edb.modeler.create_trace(
+            path_list=[
+                [self.x0, self.y0],
+                [
+                    self.x0 + self.length * math.cos(angle_rad),
+                    self.y0 + self.length * math.sin(angle_rad),
+                ],
+            ],
+            layer_name=self.layer,
+            width="w",
+            net_name=self.net,
+            start_cap_style="Flat",
+            end_cap_style="Flat",
+        )
+        self.id = trace.id
+        return trace