PyOptik 1.9.8__0-py3-none-any.whl

PyOptik/data/tabulated/zinc.yml

@@ -0,0 +1,162 @@
+# this file is part of refractiveindex.info database
+# refractiveindex.info database is in the public domain
+# copyright and related rights waived via CC0 1.0
+
+REFERENCES: "W. S. M. Werner, K. Glantschnig, C. Ambrosch-Draxl. Optical constants and inelastic electron-scattering data for 17 elemental metals, <a href=\"https://doi.org/10.1063/1.3243762\"><i>J. Phys Chem Ref. Data</i> <b>38</b>, 1013-1092 (2009)</a>"
+COMMENTS: "Experimental data: Derived from reflection electron energy-loss spectroscopy (REELS) spectra."
+DATA:
+  - type: tabulated nk
+    data: |
+        0.017586 0.9522 0.0688
+        0.017712 0.9523 0.0693
+        0.017839 0.9518 0.0704
+        0.017969 0.9519 0.0709
+        0.018100 0.9514 0.0720
+        0.018233 0.9515 0.0725
+        0.018368 0.9510 0.0736
+        0.018505 0.9510 0.0741
+        0.018644 0.9506 0.0752
+        0.018785 0.9507 0.0763
+        0.018929 0.9502 0.0768
+        0.019074 0.9498 0.0779
+        0.019222 0.9499 0.0790
+        0.019373 0.9494 0.0800
+        0.019525 0.9495 0.0811
+        0.019680 0.9490 0.0817
+        0.019837 0.9491 0.0827
+        0.019997 0.9487 0.0838
+        0.020160 0.9488 0.0848
+        0.020325 0.9484 0.0859
+        0.020493 0.9479 0.0870
+        0.020664 0.9480 0.0881
+        0.020838 0.9476 0.0892
+        0.021014 0.9477 0.0902
+        0.021194 0.9473 0.0913
+        0.021377 0.9469 0.0924
+        0.021562 0.9470 0.0940
+        0.021752 0.9466 0.0951
+        0.021944 0.9467 0.0961
+        0.022140 0.9463 0.0972
+        0.022339 0.9459 0.0988
+        0.022543 0.9460 0.0999
+        0.022749 0.9457 0.1015
+        0.022960 0.9453 0.1026
+        0.023175 0.9454 0.1037
+        0.023393 0.9450 0.1053
+        0.023616 0.9452 0.1069
+        0.023843 0.9448 0.1080
+        0.024075 0.9445 0.1096
+        0.024311 0.9441 0.1112
+        0.024551 0.9442 0.1123
+        0.024797 0.9439 0.1139
+        0.025047 0.9436 0.1155
+        0.025303 0.9438 0.1171
+        0.025564 0.9434 0.1187
+        0.025830 0.9431 0.1203
+        0.026102 0.9428 0.1220
+        0.026380 0.9430 0.1235
+        0.026663 0.9427 0.1252
+        0.026953 0.9424 0.1268
+        0.027249 0.9421 0.1290
+        0.027552 0.9418 0.1306
+        0.027862 0.9421 0.1322
+        0.028178 0.9418 0.1343
+        0.028502 0.9416 0.1365
+        0.028834 0.9413 0.1381
+        0.029173 0.9411 0.1403
+        0.029520 0.9409 0.1424
+        0.029876 0.9406 0.1441
+        0.030240 0.9404 0.1462
+        0.030613 0.9402 0.1484
+        0.030996 0.9400 0.1505
+        0.031388 0.9399 0.1532
+        0.031791 0.9398 0.1554
+        0.032204 0.9396 0.1575
+        0.032627 0.9395 0.1602
+        0.033062 0.9393 0.1623
+        0.033509 0.9388 0.1651
+        0.033968 0.9386 0.1673
+        0.034440 0.9386 0.1699
+        0.034925 0.9385 0.1726
+        0.035424 0.9380 0.1754
+        0.035937 0.9379 0.1781
+        0.036466 0.9374 0.1808
+        0.037010 0.9375 0.1840
+        0.037571 0.9370 0.1868
+        0.038149 0.9365 0.1901
+        0.038745 0.9367 0.1932
+        0.039360 0.9362 0.1960
+        0.039995 0.9358 0.1993
+        0.040651 0.9354 0.2026
+        0.041328 0.9352 0.2064
+        0.042029 0.9349 0.2097
+        0.042753 0.9340 0.2131
+        0.043503 0.9338 0.2169
+        0.044280 0.9331 0.2208
+        0.045085 0.9328 0.2240
+        0.045920 0.9322 0.2280
+        0.046786 0.9315 0.2319
+        0.047686 0.9305 0.2364
+        0.048621 0.9294 0.2405
+        0.049594 0.9288 0.2444
+        0.050606 0.9274 0.2491
+        0.051660 0.9263 0.2531
+        0.052759 0.9244 0.2580
+        0.053906 0.9230 0.2627
+        0.055104 0.9212 0.2676
+        0.056356 0.9189 0.2726
+        0.057667 0.9161 0.2778
+        0.059040 0.9128 0.2832
+        0.060480 0.9097 0.2891
+        0.061992 0.9055 0.2949
+        0.063582 0.9010 0.3013
+        0.065255 0.8960 0.3080
+        0.067018 0.8901 0.3151
+        0.068880 0.8834 0.3232
+        0.070848 0.8751 0.3314
+        0.072932 0.8662 0.3411
+        0.075142 0.8560 0.3517
+        0.077490 0.8440 0.3637
+        0.079990 0.8304 0.3775
+        0.082656 0.8150 0.3939
+        0.085506 0.7981 0.4135
+        0.088560 0.7791 0.4370
+        0.091840 0.7589 0.4658
+        0.095372 0.7388 0.5008
+        0.099187 0.7209 0.5438
+        0.103320 0.7093 0.5942
+        0.107812 0.7080 0.6483
+        0.112713 0.7192 0.6966
+        0.118080 0.7342 0.7273
+        0.123984 0.7318 0.7400
+        0.130510 0.6955 0.7541
+        0.137760 0.6304 0.7971
+        0.145864 0.5561 0.8811
+        0.154980 0.4922 1.0066
+        0.165312 0.4516 1.1670
+        0.177120 0.4444 1.3570
+        0.190745 0.4937 1.5648
+        0.206640 0.6184 1.7124
+        0.225426 0.6196 1.7464
+        0.247968 0.4994 1.9795
+        0.261019 0.4873 2.1720
+        0.275520 0.5524 2.3940
+        0.291728 0.7446 2.4832
+        0.309960 0.6411 2.4613
+        0.330625 0.5362 2.7034
+        0.354241 0.5240 2.9907
+        0.381490 0.5508 3.2986
+        0.413281 0.6029 3.6376
+        0.450852 0.6794 4.0209
+        0.495937 0.7858 4.4636
+        0.551041 0.9322 4.9850
+        0.619921 1.1356 5.6116
+        0.708481 1.4247 6.3811
+        0.826561 1.8489 7.3486
+        0.991874 2.4976 8.5981
+        1.239842 3.5417 10.2617
+        1.653123 5.3318 12.5605
+        2.479684 8.6709 15.9470
+SPECS:
+    n_is_absolute: true
+    wavelength_is_vacuum: true

PyOptik/directories.py

@@ -0,0 +1,44 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+from pathlib import Path
+import PyOptik
+
+
+__all__ = [
+    'root_path',
+    'project_path',
+    'doc_path',
+    'doc_css_path',
+    'logo_path'
+]
+
+root_path = Path(PyOptik.__path__[0])
+
+project_path = root_path.parents[0]
+
+example_directory = root_path.joinpath('examples')
+
+doc_path = project_path.joinpath('docs')
+
+doc_css_path = doc_path.joinpath('source/_static/default.css')
+
+logo_path = doc_path.joinpath('images/logo.png')
+
+examples_path = root_path.joinpath('examples')
+
+sellmeier_data_path = root_path.joinpath('data/sellmeier')
+tabulated_data_path = root_path.joinpath('data/tabulated')
+
+data_path = root_path.joinpath('data')
+
+libraries_path = root_path.joinpath('libraries')
+
+
+if __name__ == '__main__':
+    for path_name in __all__:
+        path = locals()[path_name]
+        print(path)
+        assert path.exists(), f"Path {path_name} do not exists"
+
+# -

PyOptik/libraries/classics.yml

@@ -0,0 +1,28 @@
+sellmeier:
+  ZBLAN: https://refractiveindex.info/database/data-nk/glass/misc/ZBLAN/Gan.yml
+  BK7: https://refractiveindex.info/database/data-nk/glass/schott/N-BK7.yml
+  fused_silica: https://refractiveindex.info/database/data-nk/main/SiO2/Malitson.yml
+  BAK1: https://refractiveindex.info/database/data-nk/glass/schott/N-BAK1.yml
+  SF5: https://refractiveindex.info/database/data-nk/glass/schott/N-SF5.yml
+  water: https://refractiveindex.info/database/data-nk/main/H2O/Daimon-19.0C.yml
+  polystyren: https://refractiveindex.info/database/data-nk/organic/(C8H8)n%20-%20polystyrene/Sultanova.yml
+  lithium_niobate: https://refractiveindex.info/database/data-nk/main/LiNbO3/Zelmon-o.yml
+  germanium: https://refractiveindex.info/database/data-nk/main/Ge/Icenogle.yml
+  silicon: https://refractiveindex.info/database/data-nk/main/Si/Salzberg.yml
+  soda_lime_glass: https://refractiveindex.info/database/data-nk/glass/misc/soda-lime/Rubin-clear.yml
+  crown: https://refractiveindex.info/database/data-nk/glass/schott/K7.yml
+  flint: https://refractiveindex.info/database/data-nk/glass/schott/F2.yml
+  argon: https://refractiveindex.info/database/data-nk/main/Ar/Bideau-Mehu.yml
+  air: https://refractiveindex.info/database/data-nk/other/mixed%20gases/air/Ciddor.yml
+  water: https://refractiveindex.info/database/data-nk/main/H2O/Daimon-19.0C.yml
+
+tabulated:
+  silver: https://refractiveindex.info/database/data-nk/main/Ag/Johnson.yml
+  aluminium: https://refractiveindex.info/database/data-nk/main/Al/Rakic.yml
+  gold: https://refractiveindex.info/database/data-nk/main/Au/Johnson.yml
+  calcium: https://refractiveindex.info/database/data-nk/main/Ca/Mathewson.yml
+  copper: https://refractiveindex.info/database/data-nk/main/Cu/Johnson.yml
+  iron: https://refractiveindex.info/database/data-nk/main/Fe/Johnson.yml
+  sodium: https://refractiveindex.info/database/data-nk/main/Na/Smith.yml
+  nickel: https://refractiveindex.info/database/data-nk/main/Ni/Johnson.yml
+  zinc: https://refractiveindex.info/database/data-nk/main/Zn/Werner.yml

PyOptik/libraries/dielectrics.yml

@@ -0,0 +1,20 @@
+sellmeier:
+  ZBLAN: https://refractiveindex.info/database/data-nk/glass/misc/ZBLAN/Gan.yml
+  BK7: https://refractiveindex.info/database/data-nk/glass/schott/N-BK7.yml
+  fused_silica: https://refractiveindex.info/database/data-nk/main/SiO2/Malitson.yml
+  BAK1: https://refractiveindex.info/database/data-nk/glass/schott/N-BAK1.yml
+  BAF10: https://refractiveindex.info/database/data-nk/glass/schott/N-BAF10.yml
+  SF5: https://refractiveindex.info/database/data-nk/glass/schott/N-SF5.yml
+  SF10: https://refractiveindex.info/database/data-nk/glass/schott/N-SF10.yml
+  SF11: https://refractiveindex.info/database/data-nk/glass/schott/N-SF11.yml
+  LASF9: https://refractiveindex.info/database/data-nk/glass/schott/N-LASF9.yml
+  FK51A: https://refractiveindex.info/database/data-nk/glass/schott/N-FK51A.yml
+  water: https://refractiveindex.info/database/data-nk/main/H2O/Daimon-19.0C.yml
+  polystyren: https://refractiveindex.info/database/data-nk/organic/(C8H8)n%20-%20polystyrene/Sultanova.yml
+  lithium_niobate: https://refractiveindex.info/database/data-nk/main/LiNbO3/Zelmon-o.yml
+  germanium: https://refractiveindex.info/database/data-nk/main/Ge/Icenogle.yml
+  silicon: https://refractiveindex.info/database/data-nk/main/Si/Salzberg.yml
+  soda_lime_glass: https://refractiveindex.info/database/data-nk/glass/misc/soda-lime/Rubin-clear.yml
+  crown: https://refractiveindex.info/database/data-nk/glass/schott/K7.yml
+  flint: https://refractiveindex.info/database/data-nk/glass/schott/F2.yml
+

PyOptik/libraries/metals.yml

@@ -0,0 +1,20 @@
+tabulated:
+  silver: https://refractiveindex.info/database/data-nk/main/Ag/Johnson.yml
+  aluminium: https://refractiveindex.info/database/data-nk/main/Al/Rakic.yml
+  gold: https://refractiveindex.info/database/data-nk/main/Au/Johnson.yml
+  calcium: https://refractiveindex.info/database/data-nk/main/Ca/Mathewson.yml
+  copper: https://refractiveindex.info/database/data-nk/main/Cu/Johnson.yml
+  iron: https://refractiveindex.info/database/data-nk/main/Fe/Johnson.yml
+  sodium: https://refractiveindex.info/database/data-nk/main/Na/Smith.yml
+  nickel: https://refractiveindex.info/database/data-nk/main/Ni/Johnson.yml
+  zinc: https://refractiveindex.info/database/data-nk/main/Zn/Werner.yml
+  chromium: https://refractiveindex.info/database/data-nk/main/Cr/Johnson.yml
+  titanium: https://refractiveindex.info/database/data-nk/main/Ti/Johnson.yml
+  tungsten: https://refractiveindex.info/database/data-nk/main/W/Werner.yml
+  lead: https://refractiveindex.info/database/data-nk/main/Pb/Werner.yml
+  bismuth: https://refractiveindex.info/database/data-nk/main/Bi/Werner.yml
+  magnesium: https://refractiveindex.info/database/data-nk/main/Mg/Hagemann.yml
+  paladium: https://refractiveindex.info/database/data-nk/main/Pd/Johnson.yml
+  rhodium: https://refractiveindex.info/database/data-nk/main/Rh/Arndt.yml
+  vanadium: https://refractiveindex.info/database/data-nk/main/V/Johnson.yml
+  tantalum: https://refractiveindex.info/database/data-nk/main/Ta/Werner.yml

PyOptik/libraries/minimal.yml

@@ -0,0 +1,9 @@
+sellmeier:
+  argon: https://refractiveindex.info/database/data-nk/main/Ar/Bideau-Mehu.yml
+  air: https://refractiveindex.info/database/data-nk/other/mixed%20gases/air/Ciddor.yml
+  water: https://refractiveindex.info/database/data-nk/main/H2O/Daimon-19.0C.yml
+  soda_lime_glass: https://refractiveindex.info/database/data-nk/glass/misc/soda-lime/Rubin-clear.yml
+  silicon: https://refractiveindex.info/database/data-nk/main/Si/Salzberg.yml
+
+tabulated:
+  zinc: https://refractiveindex.info/database/data-nk/main/Zn/Werner.yml

PyOptik/libraries/organics.yml

@@ -0,0 +1,11 @@
+sellmeier:
+  ethanol: https://refractiveindex.info/database/data-nk/organic/C2H6O%20-%20ethanol/Rheims.yml
+  methanol: https://refractiveindex.info/database/data-nk/organic/CH4O%20-%20methanol/El-Kashef.yml
+  propanol: https://refractiveindex.info/database/data-nk/organic/C3H8O%20-%20propanol/Kozma.yml
+  ethylene: https://refractiveindex.info/database/data-nk/organic/C2H4%20-%20ethylene/Loria.yml
+  acetylene: https://refractiveindex.info/database/data-nk/organic/C2H2%20-%20acetylene/Loria.yml
+  methane: https://refractiveindex.info/database/data-nk/organic/CH4%20-%20methane/Loria.yml
+  ethane: https://refractiveindex.info/database/data-nk/organic/C2H6%20-%20ethane/Loria.yml
+  acetone: https://refractiveindex.info/database/data-nk/organic/C3H6O%20-%20acetone/Rheims.yml
+  glycerol: https://refractiveindex.info/database/data-nk/organic/C3H8O3%20-%20glycerol/Rheims.yml
+

PyOptik/libraries/others.yml

@@ -0,0 +1,4 @@
+sellmeier:
+  air: https://refractiveindex.info/database/data-nk/other/mixed%20gases/air/Ciddor.yml
+  water: https://refractiveindex.info/database/data-nk/main/H2O/Daimon-19.0C.yml
+  argon: https://refractiveindex.info/database/data-nk/main/Ar/Bideau-Mehu.yml

PyOptik/libraries/polymers.yml

@@ -0,0 +1,8 @@
+sellmeier:
+  polyetylene: https://refractiveindex.info/database/data-nk/organic/(C2H4)n%20-%20polyethylene/Smith.yml
+  pva: https://refractiveindex.info/database/data-nk/organic/(C2H3Cl)n%20-%20polyvinyl%20chloride/Zhang.yml
+  pvc: https://refractiveindex.info/database/data-nk/organic/(C2H4O)n%20-%20polyvinyl%20alcohol/Schnepf.yml
+  pdms: https://refractiveindex.info/database/data-nk/organic/(C2H6OSi)n%20-%20polydimethylsiloxane/Schneider-RTV615.yml
+  polystyrene: https://refractiveindex.info/database/data-nk/organic/(C8H8)n%20-%20polystyrene/Sultanova.yml
+  polycarbonate: https://refractiveindex.info/database/data-nk/organic/(C16H14O3)n%20-%20polycarbonate/Sultanova.yml
+  cellulose: https://refractiveindex.info/database/data-nk/organic/(C6H10O5)n%20-%20cellulose/Sultanova.yml

PyOptik/material/__init__.py

@@ -0,0 +1,3 @@
+from .base_class import BaseMaterial
+from .sellmeier_class import SellmeierMaterial
+from .tabulated_class import TabulatedMaterial

PyOptik/material/base_class.py

@@ -0,0 +1,90 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+from typing import Callable
+from PyOptik.units import Quantity, meter
+import numpy
+import warnings
+
+
+class BaseMaterial:
+    def __eq__(self, other) -> bool:
+        if not isinstance(other, self.__class__):
+            return False
+
+        if self.filename != other.filename:
+            return False
+
+        return True
+
+    def __str__(self) -> str:
+        """
+        Provides an informal string representation of the Material object.
+
+        Returns
+        -------
+        str
+            Informal representation of the Material object.
+        """
+        return f"Material: {self.filename}"
+
+    def __repr__(self) -> str:
+        """
+        Provides a formal string representation of the Material object, including key attributes.
+
+        Returns
+        -------
+        str
+            Formal representation of the Material object.
+        """
+        return self.__str__()
+
+    def _check_wavelength(self, wavelength: Quantity) -> None:
+        """
+        Checks if a wavelength is within the material's allowable range and raises a warning if it is not.
+
+        Parameters
+        ----------
+        wavelength : Quantity
+            The wavelength to check, in micrometers.
+
+        Raises
+        ------
+        UserWarning
+            If the wavelength is outside the allowable range.
+        """
+        if self.wavelength_bound is not None:
+            min_value, max_value = self.wavelength_bound
+
+            if numpy.any((wavelength < min_value) | (wavelength > max_value)):
+                warnings.warn(
+                    f"Wavelength range goes from {wavelength.min().to_compact()} to {wavelength.max().to_compact()} "
+                    f"which is outside the allowable range of {min_value.to_compact()} to {max_value.to_compact()} µm. "
+                    f"[Material: {self.filename}]"
+                )
+
+    def ensure_units(func) -> Callable:
+        """
+        Decorator to ensure that the wavelength parameter has the correct units.
+
+        Parameters
+        ----------
+        func : Callable
+            The function to wrap.
+
+        Returns
+        -------
+        Callable
+            The wrapped function that ensures wavelength is a Quantity in meters.
+        """
+        def wrapper(self, wavelength: Quantity = None, *args, **kwargs):
+            if wavelength is None:
+                if self.wavelength_bound is None:
+                    raise ValueError('Wavelenght must be provided for computation.')
+                wavelength = numpy.linspace(self.wavelength_bound[0].magnitude, self.wavelength_bound[1].magnitude, 100) * self.wavelength_bound.units
+
+            import PyOptik.units
+            if not isinstance(wavelength, PyOptik.units.Quantity):
+                wavelength = wavelength * meter
+            return func(self, wavelength, *args, **kwargs)
+        return wrapper

PyOptik/material/sellmeier_class.py

@@ -0,0 +1,198 @@
+#!/usr/bin/env python
+# -*- coding: utf-8 -*-
+
+import yaml
+import numpy
+import itertools
+from MPSPlots.styles import mps
+import matplotlib.pyplot as plt
+from pydantic.dataclasses import Field, ConfigDict, dataclass
+from typing import Tuple, Optional, Union
+from PyOptik.units import micrometer, Quantity
+from PyOptik.directories import sellmeier_data_path
+from PyOptik.material.base_class import BaseMaterial
+
+
+config_dict = ConfigDict(
+    arbitrary_types_allowed=True
+)
+
+
+@dataclass(config=config_dict, slots=True, eq=False)
+class SellmeierMaterial(BaseMaterial):
+    """
+    Class representing a material with Sellmeier coefficients for refractive index computation.
+
+    Attributes
+    ----------
+    filename : str
+        The name of the YAML file containing material properties.
+    coefficients : numpy.ndarray
+        The Sellmeier coefficients used for calculating the refractive index.
+    wavelength_range : Optional[Tuple[float, float]]
+        The allowable wavelength range for the material in micrometers.
+    reference : Optional[str]
+        Reference information for the material data.
+    formula_type : int
+        The formula type to use for refractive index calculation.
+    """
+    filename: str
+
+    coefficients: Optional[numpy.ndarray] = Field(default=None)
+    wavelength_bound: Optional[Tuple[float, float]] = Field(default=None)
+    reference: Optional[str] = Field(default=None)
+    formula_type: Optional[int] = Field(default=None)
+
+    def __repr__(self) -> str:
+        return self.__str__()
+
+    def __str__(self) -> str:
+        return self.filename
+
+    def __post_init__(self) -> None:
+        """
+        Post-initialization method to load coefficients, wavelength range, formula type, and reference from a YAML file.
+        """
+        self._load_coefficients()
+
+    def _load_coefficients(self) -> None:
+        """
+        Loads the Sellmeier coefficients, wavelength range, formula type, and reference from the specified YAML file.
+        """
+        file_path = sellmeier_data_path / f'{self.filename}'
+
+        if not file_path.with_suffix('.yml').exists():
+            raise FileNotFoundError(f"YAML file {file_path} not found.")
+
+        with file_path.with_suffix('.yml').open('r') as file:
+            parsed_yaml = yaml.safe_load(file)
+
+        # Extract the formula type
+        self.formula_type = int(parsed_yaml['DATA'][0]['type'].split()[-1])
+
+        # Extract coefficients and ensure the list has exactly 7 coefficients by padding with zeros if necessary
+        coefficients_str = parsed_yaml['DATA'][0]['coefficients']
+        coefficients = list(map(float, coefficients_str.split()))
+        if len(coefficients) < 7:
+            coefficients.extend([0.0] * (7 - len(coefficients)))
+        self.coefficients = numpy.array(coefficients)
+
+        # Extract wavelength range
+        if 'wavelength_range' in parsed_yaml['DATA'][0]:
+            data_str = parsed_yaml['DATA'][0]['wavelength_range'].split()
+
+            self.wavelength_bound = numpy.array([float(val) for val in data_str]) * micrometer
+
+        else:
+            self.wavelength_bound = None
+
+        # # Extract reference
+        self.reference = parsed_yaml.get('REFERENCES', None)
+
+    @BaseMaterial.ensure_units
+    def compute_refractive_index(self, wavelength: Union[Quantity]) -> Union[float, numpy.ndarray]:
+        r"""
+        Computes the refractive index n(\u03bb) using the appropriate formula (either Formula 1, 2, 5, or 6).
+
+        Parameters
+        ----------
+        wavelength : Union[Quantity]
+            The wavelength \u03bb in meters, can be a single float or a numpy array.
+
+        Returns
+        -------
+        Union[Quantity]
+            The refractive index n(\u03bb) for the given wavelength or array of wavelengths.
+
+        Raises
+        ------
+        ValueError
+            If the wavelength is outside the specified range or if an unsupported formula type is encountered.
+        """
+        return_as_scalar = numpy.isscalar(wavelength.magnitude)
+
+        wavelength = numpy.atleast_1d(wavelength)
+        self._check_wavelength(wavelength)
+
+        # Compute the refractive index based on the formula type
+        zipped_coefficients = itertools.zip_longest(*[iter(self.coefficients[1:])] * 2)
+
+        match self.formula_type:
+            case 1:  # Formula 1 computation (standard Sellmeier)
+                n_squared = 1.0
+                for (B, C) in zipped_coefficients:
+                    n_squared += (B * wavelength.to(micrometer).magnitude**2) / (wavelength.to(micrometer).magnitude**2 - C**2)
+
+                n = numpy.sqrt(n_squared)
+
+            case 2:  # Formula 2 computation (extended Sellmeier)
+                n_squared = 1 + self.coefficients[0]
+                for (B, C) in zipped_coefficients:
+                    n_squared += (B * wavelength.to(micrometer).magnitude**2) / (wavelength.to(micrometer).magnitude**2 - C)
+                n = numpy.sqrt(n_squared)
+
+            case 5:  # Formula 5 computation (extended Sellmeier)
+                n = 1 + self.coefficients[0]
+                for (B, C) in zipped_coefficients:
+                    n = B * wavelength.to(micrometer).magnitude**C
+
+            case 6:
+                n = 1 + self.coefficients[0]
+                for (B, C) in zipped_coefficients:
+                    n = B / (C - wavelength.to(micrometer).magnitude**-2)
+
+            case _:
+                raise ValueError(f"Unsupported formula type: {self.formula_type}")
+
+        return n[0] if return_as_scalar else n
+
+    @BaseMaterial.ensure_units
+    def plot(self, wavelength: Optional[Quantity] = None) -> None:
+        """
+        Plots the refractive index as a function of wavelength over a specified range.
+
+        Parameters
+        ----------
+        wavelength : Optional[Quantity]
+            The range of wavelengths to plot, in meters. If not provided, the entire allowable wavelength range is used.
+
+        Raises
+        ------
+        ValueError
+            If the wavelength is not a 1D array or list of float values.
+        """
+        if wavelength.ndim != 1:
+            raise ValueError("wavelength must be a 1D array or list of float values.")
+
+        # Calculate the refractive index over the wavelength range
+        refractive_index = self.compute_refractive_index(wavelength)
+
+        with plt.style.context(mps):
+            fig, ax = plt.subplots()
+
+        ax.set(
+            ylabel='Refractive Index',
+            xlabel=r'Wavelength [$\mu$m]',
+            title=f"Refractive Index vs. Wavelength: [{self.filename}]"
+        )
+        ax.plot(wavelength.to(micrometer).magnitude, refractive_index.real, linewidth=2, label='Real Part')
+        ax.legend()
+
+        plt.show()
+
+    def print(self) -> str:
+        """
+        Provides a formal string representation of the Material object, including key attributes.
+
+        Returns
+        -------
+        str
+            Formal representation of the Material object.
+        """
+        return (
+            f"\nMaterial: '{self.filename}',\n"
+            f"coefficients: {self.coefficients},\n"
+            f"wavelength_range: {self.wavelength_bound},\n"
+            f"formula_type: {self.formula_type},\n"
+            f"reference: '{self.reference}')"
+        )