Material Properties of Glass

Scott Prahl

Sept 2023

This replicates many plots found in Ghatak in chapter 6 to illustrate the changes in material properties with wavelength.

# Jupyterlite support for ofiber
try:
    import micropip
    await micropip.install("ofiber")
except ModuleNotFoundError:
    pass

import numpy as np
import matplotlib.pyplot as plt
import ofiber

# to make graphs a bit better
%config InlineBackend.figure_format='retina'

Ghatak figure 6.1

λ = np.linspace(600,1600,50)*1e-9

glass = ofiber.glass(0)
name = ofiber.glass_name(0)
n=ofiber.n(glass,λ)
plt.plot(λ*1e9,n)
plt.title(name+' (compare with Ghatak fig 6.1)')
plt.xlabel('Wavelength [nm]')
plt.ylabel('Index of Refraction')
plt.show()

Ghatak figure 6.2

λ = np.linspace(600,1600,50)*1e-9

glass = ofiber.glass(0)
name = ofiber.glass_name(0)

dn=ofiber.dn(glass,λ)*1e-6
plt.plot(λ*1e9,dn)
plt.title(name+' (compare with Ghatak fig 6.2)')
plt.xlabel('Wavelength [nm]')
plt.ylabel(r'$dn/d\lambda$ [1/$\mu$m]')
plt.show()

Ghatak figure 6.3

λ = np.linspace(600,1600,50)*1e-9

glass = ofiber.glass(0)
name = ofiber.glass_name(0)

d2n=ofiber.d2n(glass,λ)*1e-12
plt.plot(λ*1e9,d2n)
plt.plot([600,1600],[0,0],':k')
plt.title(name+' (compare with Ghatak fig 6.3)')
plt.xlabel('Wavelength [nm]')
plt.ylabel(r'$d^2n/d\lambda^2$ [1/$\mu$m$^2$]')
plt.show()

Ghatak figure 6.4

λ = np.linspace(1000,1600,50)*1e-9

for x in [0,0.063,0.193]:
    glass = ofiber.doped_glass(x)
    name = ofiber.doped_glass_name(x)
    d2n=ofiber.d2n(glass,λ)*1e-12
    plt.plot(λ*1e9,d2n,label=name)

plt.plot([1000,1600],[0,0],':k')
plt.title('Ghatak fig 6.4')
plt.xlabel('Wavelength [nm]')
plt.ylabel(r'$d^2n/d\lambda^2$ [1/$\mu$m$^2$]')
plt.legend()
plt.show()

Ghatak figure 6.5

λ = np.linspace(1500,3500,50)*1e-9
for i in range(6,11):
    glass = ofiber.glass(i)
    name = ofiber.glass_name(i)
    n = ofiber.n(glass,λ)
    plt.plot(λ*1e9,n,label=name)

plt.xlim(1500,3800)
plt.xlabel("Wavelength [nm]")
plt.ylabel('Index of Refraction')
plt.title('figure 6.5 in Ghatak')
plt.legend()
plt.show()

Ghatak figure 6.6

λ = np.linspace(1500,3500,50)*1e-9
for i in range(6,11):
    glass = ofiber.glass(i)
    name = ofiber.glass_name(i)
    d2n = ofiber.d2n(glass,λ)*1e-12
    plt.plot(λ*1e9,d2n,label=name)


plt.axhline(0,color='black', ls=':')
plt.xlabel("Wavelength [nm]")
plt.ylabel(r'$d^2n/d\lambda^2$ [$\mu m^2$]')
plt.legend(loc='upper right')
plt.title('Figure 6.6 in Ghatak')
plt.show()