spdcalc 2.0.1

SPDCalc, the Spontaneous Parametric Downconversion Calculator
Documentation 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53

//! # BBO1 CrystalType
//!
//! [More Information](http://www.newlightphotonics.com/v1/bbo-properties.html)
//!
//! ## Example
//! ```
//! use spdcalc::{crystal::*, dim::ucum, na::Vector3, utils::*};
//! let nm = spdcalc::dim::f64prefixes::NANO * ucum::M;
//! let indices = CrystalType::BBO_1.get_indices(1550.0 * nm, from_celsius_to_kelvin(20.));
//! let expected = ucum::Unitless::new(Vector3::new(
//!   1.6465863905897784,
//!   1.6465863905897784,
//!   1.531015476762964
//! ));
//! assert_eq!(indices, expected)
//! ```

use super::*;
use crate::utils::*;
use dim::{
  f64prefixes::MICRO,
  ucum::{Kelvin, K, M},
};

pub const META: CrystalMeta = CrystalMeta {
  id: "BBO_1",
  name: "BBO ref 1",
  reference_url: "http://www.newlightphotonics.com/v1/bbo-properties.html",
  axis_type: OpticAxisType::NegativeUniaxial,
  point_group: PointGroup::HM_3m,
  transmission_range: Some(ValidWavelengthRange(189e-9, 3_500e-9)),
  temperature_dependence_known: true,
};

// from Newlight Photonics
const DNO: f64 = -9.3e-6;
const DNE: f64 = -16.6e-6;

/// Get refractive Indices
#[allow(clippy::unreadable_literal)]
pub fn get_indices(wavelength: Wavelength, temperature: Kelvin<f64>) -> Indices {
  let l_sq = (wavelength / (MICRO * M)).powi(2); // Convert for Sellmeier Coefficients

  let mut no = (2.7359 + 0.01878 / (l_sq - 0.01822) - 0.01354 * l_sq).sqrt();
  let mut ne = (2.3753 + 0.01224 / (l_sq - 0.01667) - 0.01516 * l_sq).sqrt();

  let f = *((temperature - from_celsius_to_kelvin(20.0)) / K);

  no += f * DNO;
  ne += f * DNE;

  Indices::new(na::Vector3::new(no, no, ne))
}