tidy3d.CustomMedium#

class CustomMedium[source]#

Bases: AbstractCustomMedium

Medium with user-supplied permittivity distribution.

Parameters:

attrs (dict = {}) – Dictionary storing arbitrary metadata for a Tidy3D object. This dictionary can be freely used by the user for storing data without affecting the operation of Tidy3D as it is not used internally. Note that, unlike regular Tidy3D fields, attrs are mutable. For example, the following is allowed for setting an attr obj.attrs['foo'] = bar. Also note that Tidy3D` will raise a TypeError if attrs contain objects that can not be serialized. One can check if attrs are serializable by calling obj.json().
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.
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, apply the subpixel averaging method specified by Simulation’s field subpixel for this type of material on the interface of the structure, including exterior boundary and intersection interfaces with other structures.
eps_dataset (Optional[PermittivityDataset] = None) – [To be deprecated] User-supplied dataset containing complex-valued permittivity as a function of space. Permittivity distribution over the Yee-grid will be interpolated based on interp_method.
permittivity (Union[SpatialDataArray, Annotated[Union[tidy3d.components.data.dataset.TriangularGridDataset, tidy3d.components.data.dataset.TetrahedralGridDataset], FieldInfo(default=PydanticUndefined, discriminator='type', extra={})], NoneType] = None) – [units = None (relative permittivity)]. Spatial profile of relative permittivity.
conductivity (Union[SpatialDataArray, Annotated[Union[tidy3d.components.data.dataset.TriangularGridDataset, tidy3d.components.data.dataset.TetrahedralGridDataset], FieldInfo(default=PydanticUndefined, discriminator='type', extra={})], NoneType] = None) – [units = S/um]. Spatial profile Electric conductivity. Defined such that the imaginary part of the complex permittivity at angular frequency omega is given by conductivity/omega.

Example

>>> Nx, Ny, Nz = 10, 9, 8
>>> X = np.linspace(-1, 1, Nx)
>>> Y = np.linspace(-1, 1, Ny)
>>> Z = np.linspace(-1, 1, Nz)
>>> coords = dict(x=X, y=Y, z=Z)
>>> permittivity= SpatialDataArray(np.ones((Nx, Ny, Nz)), coords=coords)
>>> conductivity= SpatialDataArray(np.ones((Nx, Ny, Nz)), coords=coords)
>>> dielectric = CustomMedium(permittivity=permittivity, conductivity=conductivity)
>>> eps = dielectric.eps_model(200e12)

Attributes

`freqs`	float array of frequencies.
`is_isotropic`	Check if the medium is isotropic or anisotropic.
`is_spatially_uniform`	Whether the medium is spatially uniform.
`n_cfl`	This property computes the index of refraction related to CFL condition, so that the FDTD with this medium is stable when the time step size that doesn't take material factor into account is multiplied by n_cfl`.
`attrs`

Methods

`compute_derivatives`(derivative_info)	Compute the adjoint derivatives for this object.
`eps_dataarray_freq`(frequency)	Permittivity array at `frequency`.
`eps_diagonal`(frequency)	Main diagonal of the complex-valued permittivity tensor at `frequency`.
`eps_diagonal_on_grid`(frequency, coords)	Spatial profile of main diagonal of the complex-valued permittivity at `frequency` interpolated at the supplied coordinates.
`eps_model`(frequency)	Spatial and polarizaiton average of complex-valued permittivity as a function of frequency.
`from_eps_raw`(eps[, freq, interp_method])	Construct a `CustomMedium` from datasets containing raw permittivity values.
`from_nk`(n[, k, freq, interp_method])	Construct a `CustomMedium` from datasets containing n and k values.
`grids`(bounds)	Make a `Grid` corresponding to the data in each `eps_ii` component.

Inherited Common Usage

eps_dataset#

permittivity#

conductivity#

property is_isotropic#: Check if the medium is isotropic or anisotropic.

property is_spatially_uniform#: Whether the medium is spatially uniform.

property freqs#: float array of frequencies. This field is to be deprecated in v3.0.

property n_cfl#

This property computes the index of refraction related to CFL condition, so that the FDTD with this medium is stable when the time step size that doesn’t take material factor into account is multiplied by n_cfl`.

For dispersiveless custom medium, it equals min[sqrt(eps_inf)], where min is performed over all components and spatial points.

eps_dataarray_freq(frequency)[source]#

Permittivity array at frequency. ()

Parameters:

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

Returns:

Tuple[ –

Union[
SpatialDataArray, TriangularGridDataset, TetrahedralGridDataset,

], Union[

SpatialDataArray, TriangularGridDataset, TetrahedralGridDataset,

], Union[

SpatialDataArray, TriangularGridDataset, TetrahedralGridDataset,

],
] – The permittivity evaluated at frequency.

eps_diagonal_on_grid(frequency, coords)[source]#

Spatial profile of main diagonal of the complex-valued permittivity at frequency interpolated at the supplied coordinates.

Parameters:

frequency (float) – Frequency to evaluate permittivity at (Hz).
coords (Coords) – The grid point coordinates over which interpolation is performed.

Returns:

The complex-valued permittivity tensor at frequency interpolated at the supplied coordinate.

Return type:

Tuple[ArrayComplex3D, ArrayComplex3D, ArrayComplex3D]

eps_diagonal(frequency)[source]#: Main diagonal of the complex-valued permittivity tensor at frequency. Spatially, we take max{|eps|}, so that autoMesh generation works appropriately.

eps_model(frequency)[source]#: Spatial and polarizaiton average of complex-valued permittivity as a function of frequency.

classmethod from_eps_raw(eps, freq=None, interp_method='nearest', **kwargs)[source]#

Construct a CustomMedium from datasets containing raw permittivity values.

Parameters:

eps (Union[) –

SpatialDataArray,
ScalarFieldDataArray, TriangularGridDataset, TetrahedralGridDataset,

] Dataset containing complex-valued permittivity as a function of space.
freq (float, optional) – Frequency at which eps are defined.
interp_method (InterpMethod, optional) – Interpolation method to obtain permittivity values that are not supplied at the Yee grids.

Notes

For lossy medium that has a complex-valued eps, if eps is supplied through SpatialDataArray, which doesn’t contain frequency information, the freq kwarg will be used to evaluate the permittivity and conductivity. Alternatively, eps can be supplied through ScalarFieldDataArray, which contains a frequency coordinate. In this case, leave freq kwarg as the default of None.

Returns:: Medium containing the spatially varying permittivity data.
Return type:: CustomMedium

classmethod from_nk(n, k=None, freq=None, interp_method='nearest', **kwargs)[source]#

Construct a CustomMedium from datasets containing n and k values.

Parameters:

n (Union[) –

SpatialDataArray,
ScalarFieldDataArray, TriangularGridDataset, TetrahedralGridDataset,

] Real part of refractive index.
k (Union[) –

SpatialDataArray,
ScalarFieldDataArray, TriangularGridDataset, TetrahedralGridDataset,

], optional Imaginary part of refrative index for lossy medium.
freq (float, optional) – Frequency at which n and k are defined.
interp_method (InterpMethod, optional) – Interpolation method to obtain permittivity values that are not supplied at the Yee grids.

Note

For lossy medium, if both n and k are supplied through SpatialDataArray, which doesn’t contain frequency information, the freq kwarg will be used to evaluate the permittivity and conductivity. Alternatively, n and k can be supplied through ScalarFieldDataArray, which contains a frequency coordinate. In this case, leave freq kwarg as the default of None.

Returns:: Medium containing the spatially varying permittivity data.
Return type:: CustomMedium

grids(bounds)[source]#: Make a Grid corresponding to the data in each eps_ii component. The min and max coordinates along each dimension are bounded by bounds.

compute_derivatives(derivative_info)[source]#: Compute the adjoint derivatives for this object.

__hash__()#: Hash method.

tidy3d.CustomMedium

Contents

tidy3d.CustomMedium#