materials.h

Go to the documentation of this file.

 
#ifndef QPMS_MATERIALS_H
#define QPMS_MATERIALS_H
#include "qpms_types.h"
#include <gsl/gsl_spline.h>
 
#ifndef SPEED_OF_LIGHT
#define SPEED_OF_LIGHT (2.99792458e8)
#endif
 
 
typedef struct qpms_epsmu_generator_t {
 
    qpms_epsmu_t (*function) (complex double omega, const void *params);
    const void *params;
} qpms_epsmu_generator_t;
 
static inline qpms_epsmu_t qpms_epsmu_generator_eval(
        qpms_epsmu_generator_t gen, complex double omega) {
    return gen.function(omega, gen.params);
}
 
qpms_epsmu_t qpms_epsmu_const_g(complex double omega, 
    const void *epsmu 
    );
 
 
static inline complex double qpms_refindex(qpms_epsmu_t em) {
    return csqrt(em.eps * em.mu);
}
 
 
static inline complex double qpms_wavenumber(complex double omega, qpms_epsmu_t em) {
    return qpms_refindex(em)*omega/SPEED_OF_LIGHT;
}
 
 
static inline complex double qpms_waveimpedance(qpms_epsmu_t em) {
    return csqrt(em.mu / em.eps);
}
 
typedef struct qpms_ldparams_triple_t {
    double f;
    double omega;
    double gamma;
} qpms_ldparams_triple_t;
 
 
typedef struct qpms_ldparams_t {
    complex double eps_inf; 
    double omega_p; 
    size_t n; 
    qpms_ldparams_triple_t data[]; 
} qpms_ldparams_t;
 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_AG; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_AU; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_CU; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_AL; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_CR; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_TI; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_BE; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_NI; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_PD; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_PT; 
extern const qpms_ldparams_t *const QPMS_LDPARAMS_W ; 
 
complex double qpms_lorentzdrude_eps(complex double omega, const qpms_ldparams_t *);
 
qpms_epsmu_t qpms_lorentzdrude_epsmu(complex double omega, const qpms_ldparams_t *);
 
qpms_epsmu_t qpms_lorentzdrude_epsmu_g(
        complex double omega, 
        const void *ldparams 
        );
 
 
// TODO use gsl_interp instead of gsl_spline.
typedef struct qpms_permittivity_interpolator_t {
        double *wavelength_m; 
        double *n; 
        double *k; 
        gsl_interp *interp_n; 
        gsl_interp *interp_k; 
        size_t size; 
        // I could add gsl_interp_accel, but that is not necessary.
} qpms_permittivity_interpolator_t;
 
qpms_permittivity_interpolator_t *qpms_permittivity_interpolator_create(const size_t incount,
                const double *wavelength_m, 
                const double *n, 
                const double *k, 
                const gsl_interp_type *iptype 
                );
 
qpms_permittivity_interpolator_t *qpms_permittivity_interpolator_from_yml(
                const char *path, 
                const gsl_interp_type *iptype 
                );
 
complex double qpms_permittivity_interpolator_eps_at_omega(
                const qpms_permittivity_interpolator_t *interp, double omega_SI);
 
 
qpms_epsmu_t qpms_permittivity_interpolator_epsmu_g(
        complex double omega, 
                const void * interpolator 
        );
 
double qpms_permittivity_interpolator_omega_min(
                const qpms_permittivity_interpolator_t *ip);
 
double qpms_permittivity_interpolator_omega_max(
                const qpms_permittivity_interpolator_t *ip);
 
void qpms_permittivity_interpolator_free(qpms_permittivity_interpolator_t *interp);
 
static inline complex double qpms_drude_epsilon(
                complex double eps_inf, 
                complex double omega_p, 
                complex double gamma_p, 
                complex double omega 
                ) {
        return eps_inf - omega_p*omega_p/(omega*(omega+I*gamma_p));
}
 
 
#endif //QPMS_MATERIALS_H