scambio.constant

 1r"""Physical constants."""
 2
 3from scipy.constants import au, c, e, h, k, sigma
 4
 5
 6class ConstantMeta(type):
 7    r"""Metaclass of Constant."""
 8
 9    def __repr__(cls) -> str:
10        r"""Represent Constant class."""
11        return "\n".join(
12            [
13                f"Physical constants:",
14                f"-    au: {cls.au:.2e} cm",
15                f"-     c: {cls.c:.2e} cm·s¯¹",
16                f"-     h: {cls.h:.2e} eV·s",
17                f"-     k: {cls.k:.2e} eV·K¯¹",
18                f"-     q: {cls.q:.2e} mC",
19                f"-     σ: {cls.sigma:.2e} mW·cm¯²·K¯⁴",
20                f"- f_Sun: {cls.fsun:.2e}",
21                f"- T_Sun: {cls.tsun:.2e} K",
22                f"- µ_Sun: {cls.usun:.2e} eV",
23                f"- P_Sun: {cls.psun:.2e} mW·cm¯²",
24            ]
25        )
26
27    @property
28    def au(cls) -> float:
29        r"""Astronomical unit [$\mathrm{cm}$]."""
30        return 1e2 * au
31
32    @property
33    def c(cls) -> float:
34        r"""Speed of light in vacuum [$\mathrm{cm \, s^{-1}}$]."""
35        return 1e2 * c
36
37    @property
38    def h(cls) -> float:
39        r"""Plack's constant [$\mathrm{eV \, s}$]."""
40        return h / e
41
42    @property
43    def k(cls) -> float:
44        r"""Boltzmann's constant [$\mathrm{eV \, K^{-1}}$]."""
45        return k / e
46
47    @property
48    def q(cls) -> float:
49        r"""Elementary charge [$\mathrm{mC}$]."""
50        return 1e3 * e
51
52    @property
53    def sigma(cls) -> float:
54        r"""Stefan-Boltzmann's constant [$\mathrm{mW \, cm^{-2} \, K^{-4}}$]."""
55        return 1e-1 * sigma
56
57    @property
58    def fsun(cls) -> float:
59        r"""Sun-Earth sunlight attenuantion factor."""
60        return (cls.rsun / cls.au) ** 2
61
62    @property
63    def rsun(cls) -> float:
64        r"""Sun radius [$\mathrm{cm}$]."""
65        return 6.957e10
66
67    @property
68    def tsun(cls) -> float:
69        r"""Sun temperature [$\mathrm{K}$]."""
70        return 6000.0
71
72    @property
73    def usun(cls) -> float:
74        r"""Sun chemical potential [$\mathrm{eV}$]."""
75        return 0.0
76
77    @property
78    def psun(cls) -> float:
79        r"""Sunlight power density [$\mathrm{mW \, cm^{-2}}$]."""
80        return cls.fsun * cls.sigma * cls.tsun**4
81
82
83class Constant(metaclass=ConstantMeta):
84    r"""Physical constants."""