Merge pull request #104 from gdsfactory/add_heater_model

add heater models

Author: joamatab

CHANGELOG.md

@@ -2,6 +2,16 @@
 
 <!-- towncrier release notes start -->
 
+## 0.14.0
+
+- fix release flow [#102](https://github.com/gdsfactory/cspdk/pull/102)
+- update gdsfactory9 [#100](https://github.com/gdsfactory/cspdk/pull/100)
+
+## 0.13.0
+
+- add radius_loopback [#99](https://github.com/gdsfactory/cspdk/pull/99)
+- update gdsfactory9 [#100](https://github.com/gdsfactory/cspdk/pull/100)
+
 ## 0.12.4
 
 - update to gdsfactory 8.31.0

cspdk/si220/models.py

@@ -373,6 +373,42 @@ def heater() -> sax.SDict:
     raise NotImplementedError("No model for 'heater'")
 
 
+def straight_heater_metal_sc(
+    wl: float = 1.55,
+    neff: float = 2.34,
+    voltage: float = 0,
+    vpi: float = 1.0,  # Voltage required for π-phase shift
+    length: float = 10,
+    loss: float = 0.0,
+) -> sax.SDict:
+    """Returns simple phase shifter model.
+
+    Args:
+        wl: wavelength.
+        neff: effective index.
+        voltage: applied voltage.
+        vpi: voltage required for a π-phase shift.
+        length: length.
+        loss: loss.
+    """
+    # Calculate additional phase shift due to applied voltage.
+    deltaphi = (voltage / vpi) * jnp.pi
+    phase = 2 * jnp.pi * neff * length / wl + deltaphi
+    amplitude = jnp.asarray(10 ** (-loss * length / 20), dtype=complex)
+    transmission = amplitude * jnp.exp(1j * phase)
+    return sax.reciprocal(
+        {
+            ("o1", "o2"): transmission,
+            ("l_e1", "r_e1"): 0.0,
+            ("l_e2", "r_e2"): 0.0,
+            ("l_e3", "r_e3"): 0.0,
+            ("l_e4", "r_e4"): 0.0,
+        }
+    )
+
+
+straight_heater_metal_so = straight_heater_metal_sc
+
 crossing_so = sm.crossing
 crossing_rc = sm.crossing
 crossing_sc = sm.crossing