class CustomDrude[source]#

Bases: CustomDispersiveMedium, Drude

A spatially varying dispersive medium described by the Drude model.

  • name (Optional[str] = None) – Optional unique name for medium.

  • frequency_range (Optional[Tuple[float, float]] = None) – [units = (Hz, Hz)]. Optional range of validity for the medium.

  • allow_gain (bool = False) – Allow the medium to be active. Caution: simulations with a gain medium are unstable, and are likely to diverge.Simulations where ‘allow_gain’ is set to ‘True’ will still be charged even if diverged. Monitor data up to the divergence point will still be returned and can be useful in some cases.

  • nonlinear_spec (Union[NonlinearSpec, NonlinearSusceptibility] = None) – Nonlinear spec applied on top of the base medium properties.

  • modulation_spec (Optional[ModulationSpec] = None) – Modulation spec applied on top of the base medium properties.

  • heat_spec (Union[FluidSpec, SolidSpec, NoneType] = None) – Specification of the medium heat properties. They are used for solving the heat equation via the HeatSimulation interface. Such simulations can be used for investigating the influence of heat propagation on the properties of optical systems. Once the temperature distribution in the system is found using HeatSimulation object, Simulation.perturbed_mediums_copy() can be used to convert mediums with perturbation models defined into spatially dependent custom mediums. Otherwise, the heat_spec does not directly affect the running of an optical Simulation.

  • eps_inf (SpatialDataArray) – [units = None (relative permittivity)]. Relative permittivity at infinite frequency (\(\epsilon_\infty\)).

  • coeffs (Tuple[Tuple[tidy3d.components.data.data_array.SpatialDataArray, tidy3d.components.data.data_array.SpatialDataArray], ...]) – [units = (Hz, Hz)]. List of (\(f_i, \delta_i\)) values for model.

  • interp_method (Literal['nearest', 'linear'] = nearest) – Interpolation method to obtain permittivity values that are not supplied at the Yee grids; For grids outside the range of the supplied data, extrapolation will be applied. When the extrapolated value is smaller (greater) than the minimal (maximal) of the supplied data, the extrapolated value will take the minimal (maximal) of the supplied data.

  • subpixel (bool = False) – If True and simulation’s subpixel is also True, applies subpixel averaging of the permittivity on the interface of the structure, including exterior boundary and intersection interfaces with other structures.


The frequency-dependence of the complex-valued permittivity is described by:

\[\epsilon(f) = \epsilon_\infty - \sum_i \frac{ f_i^2}{f^2 + jf\delta_i}\]


>>> x = np.linspace(-1, 1, 5)
>>> y = np.linspace(-1, 1, 6)
>>> z = np.linspace(-1, 1, 7)
>>> coords = dict(x=x, y=y, z=z)
>>> eps_inf = SpatialDataArray(np.ones((5, 6, 7)), coords=coords)
>>> f1 = SpatialDataArray(np.random.random((5, 6, 7)), coords=coords)
>>> delta1 = SpatialDataArray(np.random.random((5, 6, 7)), coords=coords)
>>> drude_medium = CustomDrude(eps_inf=eps_inf, coeffs=[(f1,delta1),])
>>> eps = drude_medium.eps_model(200e12)

See also


A dispersive medium described by the Drude model.





Permittivity array at frequency.

Inherited Common Usage


Permittivity array at frequency.


frequency (float) – Frequency to evaluate permittivity at (Hz).


The permittivity evaluated at frequency.

Return type:

Tuple[SpatialDataArray, SpatialDataArray, SpatialDataArray]


Hash method.