Skip to content

RingModel

RingModel(
*,
kappa1: complex
| collections.abc.Sequence[complex]
| collections.abc.Sequence[collections.abc.Sequence[complex]]
| photonforge.Interpolator,
kappa2: complex
| collections.abc.Sequence[complex]
| collections.abc.Sequence[collections.abc.Sequence[complex]]
| photonforge.Interpolator
| None = None,
n_eff: complex
| collections.abc.Sequence[complex]
| photonforge.Interpolator,
length: Annotated[float, units='μm'],
propagation_loss: Annotated[float, minimum=0, units='dB/μm']
| collections.abc.Sequence[
Annotated[float, minimum=0, units='dB/μm']
] = 0.0,
n_group: float
| collections.abc.Sequence[float]
| photonforge.Interpolator
| None = None,
dispersion: Annotated[float, units='s/μm²']
| collections.abc.Sequence[Annotated[float, units='s/μm²']] = 0.0,
dispersion_slope: Annotated[float, units='s/μm³']
| collections.abc.Sequence[Annotated[float, units='s/μm³']] = 0.0,
reference_frequency: Annotated[float, minimum=0, units='Hz']
| None = None,
dn_dT: Annotated[
complex | collections.abc.Sequence[complex], units='1/K'
] = 0.0,
dL_dT: Annotated[
float | collections.abc.Sequence[float], units='dB/μm/K'
] = 0.0,
temperature: Annotated[float, minimum=0, units='K'] = 293.0,
reference_temperature: Annotated[float, minimum=0, units='K'] = 293.0,
voltage: Annotated[float, units='V'] = 0.0,
dn_dv: Annotated[complex, units='1/V'] = 0.0,
dn_dv2: Annotated[complex, units='1/V²'] = 0.0,
dloss_dv: Annotated[float, units='dB/μm/V'] = 0.0,
dloss_dv2: Annotated[float, units='dB/μm/V²'] = 0.0,
ports: Annotated[collections.abc.Sequence[str], maxItems=2, minItems=2]
| Annotated[collections.abc.Sequence[str], maxItems=4, minItems=4]
| None = None,
) -> None

Analytic ring resonator/modulator model with single or double bus.

The S matrix for each mode is given by:

S=[0t1d0t100dd00t20dt20]S = \begin{bmatrix} 0 & t_1 & d & 0 \\ t_1 & 0 & 0 & d \\ d & 0 & 0 & t_2 \\ 0 & d & t_2 & 0 \\ \end{bmatrix}

in which t1t_1 and t2t_2 are the thru-port transmissions, and dd is the drop-port transmission. Temperature- and voltage-dependent effective index and ring loss are introduced analogously to the AnalyticWaveguideModel.

For single-bus models, τ2=1\tau_2 = 1 and the S matrix becomes:

S=[0t1t10]S = \begin{bmatrix} 0 & t_1 \\ t_1 & 0 \\ \end{bmatrix}

Args

  • kappa1: Coupling coefficient for the first bus waveguide. It must include magnitude and phase of the field coupling coefficient.
  • kappa2: Coupling coefficient for the second bus waveguide. It must include magnitude and phase for the field coupling coeffcient. If None, models a single-bus ring.
  • n_eff: Effective refractive index (loss can be included here by using complex values).
  • length: Round-trip length of the ring.
  • propagation_loss: Propagation loss.
  • n_group: Group index. If None, the value of n_eff is used.
  • dispersion: Chromatic dispersion coefficient.
  • dispersion_slope: Chromatic dispersion slope.
  • reference_frequency: Reference frequency for dispersion coefficients. If None, the central frequency is used.
  • dn_dT: Temperature sensitivity for n_eff.
  • dL_dT: Temperature sensitivity for propagation_loss.
  • temperature: Operating temperature.
  • reference_temperature: Reference temperature.
  • voltage: Operating voltage.
  • dn_dv: Linear voltage-dependent effective index coefficient.
  • dn_dv2: Quadratic voltage-dependent effective index coefficient.
  • dloss_dv: Linear voltage-dependent propagation loss coefficient.
  • dloss_dv2: Quadratic voltage-dependent propagation loss coefficient.
  • ports: List of port names. If not set, the sorted list of port names from the component is used.

Note

For multimode ports, mixed-mode coefficients are 0. Parameters can be specified per mode by using sequences of values. Dispersion for kappa1 and kappa2 can also be manually included by setting the coefficients to an Interpolator (with multiple values for multimode ports), or a 2D array with shape (M, N), in which M is the number of modes, and N the length of the frequency sequence used in the S matrix computation.

The internal transmission coefficients for directional couplers assume lossless coupling and add a 90° phase shift with respect to the coupling coefficients.

MethodDescription
autograd_smatrixCompute an autograd-compatible S matrix for traced parameters.
black_box_componentCreate a black-box component using this model.
estimate_costEstimate the cost for S matrix computation.
parameters_from_qRing loss and coupling from the quality factors.
s_matrixCompute the S matrix for a component using this model.
setup_time_stepperObtain a time stepper for a component using this model.
startStart computing the S matrix response from a component.
transformApply a transformation to this model.
updateUpdate this model.

autograd_smatrix(
self,
*,
component,
traced_params,
frequencies,
component_for_params,
param_names,
show_progress=True,
model_kwargs=None,
)

Compute an autograd-compatible S matrix for traced parameters.

Important

This is an abstract method that must be implemented by derived classes that support autograd-based S-matrix computation.

Args

  • component (Component): Nominal component associated with this autograd evaluation.
  • traced_params (ndarray): Traced parameter vector used by autograd.
  • frequencies (Sequence[float]): Frequency values at which to calculate the scattering parameters (in Hz).
  • component_for_params (Callable[[ndarray], Component]): Callback that returns the component for a given traced parameter vector.
  • param_names (Sequence[str]): Flattened traced parameter names aligned with traced_params.
  • show_progress (bool): Progress display hint from Component.s_matrix.
  • model_kwargs (dict[str, Any] | None): Model-specific keyword arguments forwarded from Component.s_matrix.

Returns

AutogradSMatrix: Autograd-compatible S matrix.

black_box_component(
self,
port_spec: str | photonforge.PortSpec | None = None,
technology: photonforge.Technology | None = None,
name: str | None = None,
) -> photonforge.Component

Create a black-box component using this model.

The layout is the model’s schematic symbol in a box, sized by the port width, with the port profiles drawn as stubs.

Args

  • port_spec: Port specification used in the component. If None, look for "port_spec" in config.default_kwargs.
  • technology: Component technology. If None, the default technology is used.
  • name: Component name. If None a default is used.

Returns

Component with ports and model.

estimate_cost(self, *args, **kwargs)

Estimate the cost for S matrix computation.

All arguments are passed to the model’s start method with the addition of the keyword argument cost_estimation=True.

Returns

float: Cost estimation.

parameters_from_q(
resonance_frequency: Annotated[float, minimum=0, units='Hz'],
*,
q_intrinsic: Annotated[float, exclusiveMinimum=0],
q_coupling: Annotated[float, exclusiveMinimum=0]
| collections.abc.Sequence[float],
length: Annotated[float, units='μm'],
n_group: Annotated[float, exclusiveMinimum=0],
) -> dict

Ring loss and coupling from the quality factors.

Coupled-mode papers describe a ring by its intrinsic and coupling quality factors; RingTimeStepper takes the geometry, a propagation loss and the field coupling of each bus. The two are related through the round-trip time T=ngL/cT = n_g L / c:

γ=ω02Q∣κ∣2=2γeTα=2γivg\begin{gathered} \gamma = \frac{\omega_0}{2Q} \\ |\kappa|^2 = 2\gamma_e T \\ \alpha = \frac{2\gamma_i}{v_g} \end{gathered}

Args

  • resonance_frequency: Resonance the quality factors describe.
  • q_intrinsic: Intrinsic (unloaded) quality factor.
  • q_coupling: Coupling quality factor, one value per bus.
  • length: Round-trip length.
  • n_group: Group index of the ring waveguide.

Returns

Dict with propagation_loss and kappa1 (kappa2 as well for two buses).

s_matrix(self, component, frequencies, show_progress=True, model_kwargs={})

Compute the S matrix for a component using this model.

Args

  • component (Component): Component to perform the calculation.
  • frequencies (Sequence[float]): Frequency values at which to calculate the scattering parameters (in Hz).
  • show_progress (bool): If True, show computation progress.
  • model_kwargs (dict[str, object]): Keyword arguments passed to Model.start. Accepted arguments depend on the model type.

Returns

SMatrix: Computed S matrix.

Example

Compute only the responses from mode 0 of P0 with Tidy3D to reduce the number and cost of simulations:

>>> model = pf.Tidy3DModel()
>>> s_matrix = model.s_matrix( # doctest: +SKIP
... component,
... frequencies,
... model_kwargs={"inputs": ("P0@0",)},
... )

A port name such as ("P0",) selects all modes of that port. All outputs for the selected inputs are returned; responses from unselected inputs are not computed. Use Model.estimate_cost with the same inputs passed directly to estimate the reduced run.

setup_time_stepper(
self,
component,
time_step,
carrier_frequency=0,
show_progress=True,
time_stepper_kwargs={},
)

Obtain a time stepper for a component using this model.

Args

  • component (Component): Component for the time stepper.
  • time_step (float): The interval between time steps (in seconds).
  • carrier_frequency (float): The carrier frequency used to construct the time stepper. The carrier should be omitted from the input signals, as it is handled automatically by the time stepper.
  • show_progress (bool): If True, show computation progress.
  • time_stepper_kwargs (dict[str, object]): Keyword arguments passed to TimeStepper.setup. Each time stepper type can accept a different set of arguments. Check the documentation for the each time stepper’s setup_state method for more information.

Returns

TimeStepper: Initialized time stepper.

start(
self,
component: photonforge.Component,
frequencies: collections.abc.Sequence[
typing.Annotated[float, minimum=0, units='Hz']
],
temperature: Annotated[float, minimum=0, units='K'] | None = None,
voltage: Annotated[float, units='V'] | None = None,
**kwargs: object,
) -> photonforge.SMatrix

Start computing the S matrix response from a component.

Args

  • component: Component from which to compute the S matrix.
  • frequencies: Sequence of frequencies at which to perform the computation.
  • temperature: Operating temperature override.
  • voltage: Operating voltage override.
  • **kwargs: Unused.

Returns

Result object with attributes status and s_matrix.

transform(
self, translation=(0, 0), rotation=0, scaling=1, x_reflection=False
)

Apply a transformation to this model.

The base implementation is a no-op and returns self. Derived classes can override this method to transform model-specific geometry settings.

Args

  • translation (Sequence[float] | complex): Translation vector.
  • rotation (float): Rotation angle in degrees.
  • scaling (float): Magnification factor.
  • x_reflection (bool): Mirror across the horizontal axis.

Returns

Model: This object with transformations applied.

update(self, *args, **kwargs)

Update this model.

All contents of the model are updated with the contents from the updated version. object contents or modifications introduced after the model creation (outside of the original parametric function) are not propagated.

The model’s time stepper does not participate in the update, unless the __init__ method of the derived model class replaces it via set_time_stepper.

Args

  • *args (object): Positional arguments to the parametric model.
  • **kwargs (object): Keyword arguments to the parametric model.

Returns

Model: This model.

Note

The original keyword arguments used to generate the model are stored and updated with the passed kwargs. This updated version is used in the function call and will be stored for future updates. Positional arguments are not stored.

AttributeDescription
parametric_functionstr: Function used to update the model.
parametric_kwargsdict[str, object]: Keyword arguments used to update the model.
propertiesProperties: Object properties.
random_variableslist[RandomVariable]: Random variables associated to the model’s parameters.
time_stepperTimeStepper | None: Time stepper associated with this model.

str: Function used to update the model.

dict[str, object]: Keyword arguments used to update the model.

Properties: Object properties.

list[RandomVariable]: Random variables associated to the model’s parameters.

TimeStepper | None: Time stepper associated with this model.