planet_generator 0.0.7-pre-alpha

use crate::internal::*;
use crate::prelude::*;

#[derive(
    Copy,
    Clone,
    Eq,
    PartialEq,
    Ord,
    PartialOrd,
    Hash,
    Debug,
    SmartDefault,
    Serialize,
    Deserialize,
    EnumIter,
)]
pub enum ChemicalComponentPresence {
    /// More than 50%
    #[default]
    Dominant,
    /// 50% to 10%
    Significant,
    /// 10% to 1%
    Notable,
    /// 1% to 0.1%
    Minor,
    /// 0.1% or less
    Traces,
}
#[derive(
    Copy,
    Clone,
    Eq,
    PartialEq,
    Ord,
    PartialOrd,
    Hash,
    Debug,
    SmartDefault,
    Serialize,
    Deserialize,
    EnumIter,
)]
pub enum ChemicalComponent {
    #[default]
    Hydrogen,
    Helium,
    Carbon,
    Nitrogen,
    Oxygen,
    Silicon,
    Magnesium,
    Iron,
    Sulfur,
    Sodium,
    Potassium,
    Calcium,
    Aluminum,
    Phosphorus,
    Chlorine,
    Argon,
    Neon,
    Titanium,
    Chromium,
    Manganese,
    Nickel,
    Water,
    CarbonMonoxide,
    CarbonDioxide,
    Methane,
    Ammonia,
    HydrogenSulfide,
    SulfurDioxide,
    SiliconDioxide,
    AluminiumOxide,
    IronOxide,
    IronSulfide,
    CalciumOxide,
    SodiumChloride,
    SodiumOxide,
    PotassiumChloride,
    PotassiumOxide,
    MagnesiumOxide,
    TitaniumDioxide,
    TitaniumTetrachloride,
    PhosphorusPentoxide,
    ChromiumOxide,
    ChromiumChloride,
    ManganeseDioxide,
    ManganeseOxide,
    NickelOxide,
    NickelSulfide,
    NitricOxide,
    NitrogenDioxide,
    NitricAcid,
    SulfuricAcid,
    Methanol,
    Hydroxyl,
    Formaldehyde,
    Ethylene,
    Ethane,
    Acetylene,
    Benzene,
    Acetonitrile,
    Methylamine,
    HydrogenCyanide,
    Glycine,
    Silicates,
    PolycyclicAromaticHydrocarbons,
}

impl ChemicalComponent {
    pub fn molecular_weight_amu(&self) -> f64 {
        match self {
            ChemicalComponent::Hydrogen => 1.008,
            ChemicalComponent::Helium => 4.0026,
            ChemicalComponent::Carbon => 12.01,
            ChemicalComponent::Nitrogen => 14.007,
            ChemicalComponent::Oxygen => 16.00,
            ChemicalComponent::Silicon => 28.0855,
            ChemicalComponent::Magnesium => 24.305,
            ChemicalComponent::Iron => 55.845,
            ChemicalComponent::Sulfur => 32.06,
            ChemicalComponent::Sodium => 22.99,
            ChemicalComponent::Potassium => 39.10,
            ChemicalComponent::Calcium => 40.08,
            ChemicalComponent::Aluminum => 26.98,
            ChemicalComponent::SiliconDioxide => 60.084,
            ChemicalComponent::AluminiumOxide => 101.96,
            ChemicalComponent::IronOxide => 159.69,
            ChemicalComponent::IronSulfide => 87.91,
            ChemicalComponent::CalciumOxide => 56.08,
            ChemicalComponent::SodiumChloride => 58.44,
            ChemicalComponent::SodiumOxide => 61.98,
            ChemicalComponent::PotassiumChloride => 74.55,
            ChemicalComponent::PotassiumOxide => 94.20,
            ChemicalComponent::MagnesiumOxide => 40.30,
            ChemicalComponent::TitaniumDioxide => 79.87,
            ChemicalComponent::TitaniumTetrachloride => 189.71,
            ChemicalComponent::PhosphorusPentoxide => 141.94,
            ChemicalComponent::ChromiumOxide => 151.99,
            ChemicalComponent::ChromiumChloride => 158.36,
            ChemicalComponent::ManganeseDioxide => 86.94,
            ChemicalComponent::ManganeseOxide => 70.94,
            ChemicalComponent::NickelOxide => 74.69,
            ChemicalComponent::NickelSulfide => 90.76,
            ChemicalComponent::NitricOxide => 30.006,
            ChemicalComponent::NitrogenDioxide => 46.0055,
            ChemicalComponent::NitricAcid => 63.01,
            ChemicalComponent::SulfuricAcid => 98.08,
            ChemicalComponent::Phosphorus => 30.97,
            ChemicalComponent::Chlorine => 35.45,
            ChemicalComponent::Argon => 39.95,
            ChemicalComponent::Neon => 20.18,
            ChemicalComponent::Titanium => 47.867,
            ChemicalComponent::Chromium => 51.996,
            ChemicalComponent::Manganese => 54.938,
            ChemicalComponent::Nickel => 58.693,
            ChemicalComponent::Water => 18.015,
            ChemicalComponent::CarbonMonoxide => 28.01,
            ChemicalComponent::CarbonDioxide => 44.01,
            ChemicalComponent::Methane => 16.04,
            ChemicalComponent::Ammonia => 17.031,
            ChemicalComponent::HydrogenSulfide => 34.08,
            ChemicalComponent::SulfurDioxide => 64.066,
            ChemicalComponent::Hydroxyl => 17.007,
            ChemicalComponent::NitricOxide => 30.006,
            ChemicalComponent::NitrogenDioxide => 46.0055,
            ChemicalComponent::Formaldehyde => 30.03,
            ChemicalComponent::Methanol => 32.04,
            ChemicalComponent::Ethylene => 28.05,
            ChemicalComponent::Ethane => 30.07,
            ChemicalComponent::Acetylene => 26.04,
            ChemicalComponent::Benzene => 78.11,
            ChemicalComponent::Acetonitrile => 41.05,
            ChemicalComponent::Methylamine => 31.06,
            ChemicalComponent::HydrogenCyanide => 27.025,
            ChemicalComponent::Glycine => 75.07,
            ChemicalComponent::Silicates => 60.08, // Approximate weight of the silicate group (SiO₄)
            ChemicalComponent::PolycyclicAromaticHydrocarbons => 128.16, // Approximate weight of naphthalene (C₁₀H₈)
        }
    }

    pub fn triple_point(&self) -> Option<(f64, f64)> {
        match self {
            ChemicalComponent::Hydrogen => Some((13.8, 0.070)),
            ChemicalComponent::Helium => Some((2.2, 0.0052)),
            ChemicalComponent::Carbon => None, // Carbon does not have a standard triple point
            ChemicalComponent::Nitrogen => Some((63.15, 0.1235)),
            ChemicalComponent::Oxygen => Some((54.36, 0.0015)),
            ChemicalComponent::Silicon => None, // Silicon does not have a standard triple point
            ChemicalComponent::Magnesium => None, // Magnesium does not have a standard triple point
            ChemicalComponent::Iron => None,    // Iron does not have a standard triple point
            ChemicalComponent::Sulfur => None,  // Sulfur does not have a standard triple point
            ChemicalComponent::Sodium => None,  // Sodium does not have a standard triple point
            ChemicalComponent::Potassium => None, // Potassium does not have a standard triple point
            ChemicalComponent::Calcium => None, // Calcium does not have a standard triple point
            ChemicalComponent::Aluminum => None, // Aluminum does not have a standard triple point
            ChemicalComponent::Phosphorus => None, // Phosphorus does not have a standard triple point
            ChemicalComponent::Chlorine => Some((172.2, 0.4)),
            ChemicalComponent::Argon => Some((83.81, 0.687)),
            ChemicalComponent::Neon => Some((24.56, 0.43)),
            ChemicalComponent::Titanium => None, // Titanium does not have a standard triple point
            ChemicalComponent::Chromium => None, // Chromium does not have a standard triple point
            ChemicalComponent::Manganese => None, // Manganese does not have a standard triple point
            ChemicalComponent::Nickel => None,   // Nickel does not have a standard triple point
            ChemicalComponent::SiliconDioxide => None,
            ChemicalComponent::AluminiumOxide => None,
            ChemicalComponent::IronOxide => None,
            ChemicalComponent::IronSulfide => None,
            ChemicalComponent::CalciumOxide => None,
            ChemicalComponent::SodiumChloride => None,
            ChemicalComponent::SodiumOxide => None,
            ChemicalComponent::PotassiumChloride => None,
            ChemicalComponent::PotassiumOxide => None,
            ChemicalComponent::MagnesiumOxide => None,
            ChemicalComponent::TitaniumDioxide => None,
            ChemicalComponent::TitaniumTetrachloride => None,
            ChemicalComponent::PhosphorusPentoxide => None,
            ChemicalComponent::ChromiumOxide => None,
            ChemicalComponent::ChromiumChloride => None,
            ChemicalComponent::ManganeseDioxide => None,
            ChemicalComponent::ManganeseOxide => None,
            ChemicalComponent::NickelOxide => None,
            ChemicalComponent::NickelSulfide => None,
            ChemicalComponent::NitricOxide => Some((109.5, 0.00015)),
            ChemicalComponent::NitrogenDioxide => Some((261.93, 0.001)),
            ChemicalComponent::NitricAcid => None,
            ChemicalComponent::SulfuricAcid => None,
            ChemicalComponent::Water => Some((273.16, 0.00604)),
            ChemicalComponent::CarbonMonoxide => Some((68.15, 0.00015)),
            ChemicalComponent::CarbonDioxide => Some((216.55, 5.11)),
            ChemicalComponent::Methane => Some((90.67, 0.117)),
            ChemicalComponent::Ammonia => Some((195.4, 0.060)),
            ChemicalComponent::HydrogenSulfide => Some((187.61, 0.0276)),
            ChemicalComponent::SulfurDioxide => Some((197.67, 0.0169)),
            ChemicalComponent::Hydroxyl => None, // Hydroxyl radical does not have a standard triple point
            ChemicalComponent::NitricOxide => Some((109.5, 0.00015)),
            ChemicalComponent::NitrogenDioxide => Some((261.93, 0.001)),
            ChemicalComponent::Formaldehyde => Some((155.2, 0.016)),
            ChemicalComponent::Methanol => Some((175.47, 0.08)),
            ChemicalComponent::Ethylene => Some((104.0, 0.00033)),
            ChemicalComponent::Ethane => Some((89.89, 0.00014)),
            ChemicalComponent::Acetylene => Some((192.34, 0.0127)),
            ChemicalComponent::Benzene => Some((278.68, 0.048)),
            ChemicalComponent::Acetonitrile => Some((229.4, 0.042)),
            ChemicalComponent::Methylamine => Some((175.8, 0.0014)),
            ChemicalComponent::HydrogenCyanide => Some((260.8, 0.02)),
            ChemicalComponent::Glycine => None, // Glycine does not have a standard triple point
            ChemicalComponent::Silicates => None, // Silicates as a group do not have a standard triple point
            ChemicalComponent::PolycyclicAromaticHydrocarbons => None, // PAHs do not have a standard triple point
        }
    }

    pub fn boiling_point(&self) -> Option<f64> {
        match self {
            ChemicalComponent::Hydrogen => Some(20.28),
            ChemicalComponent::Helium => Some(4.22),
            ChemicalComponent::Carbon => None, // Carbon does not have a standard boiling point
            ChemicalComponent::Nitrogen => Some(77.36),
            ChemicalComponent::Oxygen => Some(90.20),
            ChemicalComponent::Silicon => None, // Silicon does not have a standard boiling point
            ChemicalComponent::Magnesium => None, // Magnesium does not have a standard boiling point
            ChemicalComponent::Iron => None,      // Iron does not have a standard boiling point
            ChemicalComponent::Sulfur => None,    // Sulfur does not have a standard boiling point
            ChemicalComponent::Sodium => None,    // Sodium does not have a standard boiling point
            ChemicalComponent::Potassium => None, // Potassium does not have a standard boiling point
            ChemicalComponent::Calcium => None,   // Calcium does not have a standard boiling point
            ChemicalComponent::Aluminum => None,  // Aluminum does not have a standard boiling point
            ChemicalComponent::Phosphorus => None, // Phosphorus does not have a standard boiling point
            ChemicalComponent::Chlorine => Some(239.11),
            ChemicalComponent::Argon => Some(87.3),
            ChemicalComponent::Titanium => None, // Titanium does not have a standard boiling point
            ChemicalComponent::Chromium => None, // Chromium does not have a standard boiling point
            ChemicalComponent::Manganese => None, // Manganese does not have a standard boiling point
            ChemicalComponent::Nickel => None,    // Nickel does not have a standard boiling point
            ChemicalComponent::Water => Some(373.16),
            ChemicalComponent::CarbonMonoxide => Some(82.9),
            ChemicalComponent::CarbonDioxide => Some(304.25),
            ChemicalComponent::SiliconDioxide => Some(2950.0),
            ChemicalComponent::AluminiumOxide => Some(3250.0),
            ChemicalComponent::IronOxide => Some(3200.0),
            ChemicalComponent::IronSulfide => None, // Decomposes before boiling
            ChemicalComponent::CalciumOxide => Some(3500.0),
            ChemicalComponent::SodiumChloride => Some(1413.0),
            ChemicalComponent::SodiumOxide => Some(1950.0),
            ChemicalComponent::PotassiumChloride => Some(1420.0),
            ChemicalComponent::PotassiumOxide => None, // Decomposes before boiling
            ChemicalComponent::MagnesiumOxide => Some(3600.0),
            ChemicalComponent::TitaniumDioxide => Some(3200.0),
            ChemicalComponent::TitaniumTetrachloride => Some(409.0),
            ChemicalComponent::PhosphorusPentoxide => Some(613.0),
            ChemicalComponent::ChromiumOxide => Some(2973.0),
            ChemicalComponent::ChromiumChloride => Some(1307.0),
            ChemicalComponent::ManganeseDioxide => None, // Decomposes before boiling
            ChemicalComponent::ManganeseOxide => Some(2700.0),
            ChemicalComponent::NickelOxide => Some(2973.0),
            ChemicalComponent::NickelSulfide => None, // Decomposes before boiling
            ChemicalComponent::NitricOxide => Some(121.36),
            ChemicalComponent::NitrogenDioxide => None, // Decomposes before boiling
            ChemicalComponent::NitricAcid => Some(359.0),
            ChemicalComponent::SulfuricAcid => Some(611.0),
            ChemicalComponent::Methane => Some(111.66),
            ChemicalComponent::Ammonia => Some(239.81),
            ChemicalComponent::Neon => Some(27.07),
            ChemicalComponent::HydrogenSulfide => Some(212.9),
            ChemicalComponent::SulfurDioxide => Some(263.05),
            ChemicalComponent::Hydroxyl => None, // Hydroxyl radical does not have a standard boiling point
            ChemicalComponent::NitricOxide => Some(121.36),
            ChemicalComponent::NitrogenDioxide => None, // NO₂ decomposes before boiling
            ChemicalComponent::Formaldehyde => Some(252.2),
            ChemicalComponent::Methanol => Some(337.85),
            ChemicalComponent::Ethylene => Some(169.42),
            ChemicalComponent::Ethane => Some(184.55),
            ChemicalComponent::Acetylene => Some(189.34),
            ChemicalComponent::Benzene => Some(353.25),
            ChemicalComponent::Acetonitrile => Some(354.6),
            ChemicalComponent::Methylamine => Some(266.1),
            ChemicalComponent::HydrogenCyanide => Some(299.2),
            ChemicalComponent::Glycine => None, // Glycine does not have a standard boiling point
            ChemicalComponent::Silicates => None, // Silicates as a group do not have a standard boiling point
            ChemicalComponent::PolycyclicAromaticHydrocarbons => None, // PAHs do not have a standard boiling point
        }
    }

    pub fn is_chemically_stable(&self, temperature: u32, pressure: f32) -> bool {
        match self {
            ChemicalComponent::Hydrogen => temperature < 100 && pressure < 10.0,
            ChemicalComponent::Helium => true, // Helium is generally stable under all conditions
            ChemicalComponent::Carbon => true, // Carbon is stable in most conditions
            ChemicalComponent::Nitrogen => temperature < 500 && pressure < 50.0,
            ChemicalComponent::Oxygen => temperature < 500 && pressure < 50.0,
            ChemicalComponent::Silicon => true, // Silicon is stable in most conditions
            ChemicalComponent::Magnesium => temperature < 650 && pressure < 50.0,
            ChemicalComponent::Iron => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::Sulfur => temperature < 200 && pressure < 20.0,
            ChemicalComponent::Sodium => temperature < 1000 && pressure < 50.0,
            ChemicalComponent::Potassium => temperature < 800 && pressure < 50.0,
            ChemicalComponent::Calcium => temperature < 1000 && pressure < 50.0,
            ChemicalComponent::Aluminum => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::Phosphorus => temperature < 400 && pressure < 20.0,
            ChemicalComponent::Chlorine => temperature < 400 && pressure < 50.0,
            ChemicalComponent::Argon => true, // Argon is stable under most conditions
            ChemicalComponent::Neon => true,  // Neon is stable under most conditions
            ChemicalComponent::Titanium => temperature < 2000 && pressure < 100.0,
            ChemicalComponent::Chromium => temperature < 2000 && pressure < 100.0,
            ChemicalComponent::Manganese => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::Nickel => temperature < 2000 && pressure < 100.0,
            ChemicalComponent::Water => temperature < 100 && pressure < 10.0,
            ChemicalComponent::CarbonMonoxide => temperature < 500 && pressure < 50.0,
            ChemicalComponent::CarbonDioxide => temperature < 500 && pressure < 50.0,
            ChemicalComponent::Methane => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Ammonia => temperature < 400 && pressure < 50.0,
            ChemicalComponent::SiliconDioxide => true,
            ChemicalComponent::AluminiumOxide => true,
            ChemicalComponent::IronOxide => true,
            ChemicalComponent::IronSulfide => temperature < 1200 && pressure < 100.0,
            ChemicalComponent::CalciumOxide => true,
            ChemicalComponent::SodiumChloride => true,
            ChemicalComponent::SodiumOxide => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::PotassiumChloride => true,
            ChemicalComponent::PotassiumOxide => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::MagnesiumOxide => true,
            ChemicalComponent::TitaniumDioxide => true,
            ChemicalComponent::TitaniumTetrachloride => temperature < 500 && pressure < 100.0,
            ChemicalComponent::PhosphorusPentoxide => true,
            ChemicalComponent::ChromiumOxide => true,
            ChemicalComponent::ChromiumChloride => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::ManganeseDioxide => temperature < 1200 && pressure < 100.0,
            ChemicalComponent::ManganeseOxide => temperature < 1500 && pressure < 100.0,
            ChemicalComponent::NickelOxide => true,
            ChemicalComponent::NickelSulfide => temperature < 1200 && pressure < 100.0,
            ChemicalComponent::NitricOxide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::NitrogenDioxide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::NitricAcid => temperature < 400 && pressure < 50.0,
            ChemicalComponent::SulfuricAcid => temperature < 500 && pressure < 50.0,
            ChemicalComponent::HydrogenSulfide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::SulfurDioxide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Hydroxyl => temperature < 100 && pressure < 10.0,
            ChemicalComponent::NitricOxide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::NitrogenDioxide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Formaldehyde => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Methanol => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Ethylene => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Ethane => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Acetylene => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Benzene => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Acetonitrile => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Methylamine => temperature < 300 && pressure < 20.0,
            ChemicalComponent::HydrogenCyanide => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Glycine => temperature < 300 && pressure < 20.0,
            ChemicalComponent::Silicates => true, // Silicates are stable under most conditions
            ChemicalComponent::PolycyclicAromaticHydrocarbons => {
                temperature < 300 && pressure < 20.0
            }
        }
    }

    pub fn can_exist_in_atmosphere(&self, magnetic_field: MagneticFieldStrength) -> bool {
        match (self, magnetic_field) {
            // None
            (_, MagneticFieldStrength::None) => match self {
                ChemicalComponent::Water
                | ChemicalComponent::Methane
                | ChemicalComponent::Ammonia
                | ChemicalComponent::HydrogenSulfide
                | ChemicalComponent::CarbonMonoxide
                | ChemicalComponent::Formaldehyde
                | ChemicalComponent::Methanol
                | ChemicalComponent::Ethylene
                | ChemicalComponent::Ethane
                | ChemicalComponent::Acetylene
                | ChemicalComponent::Methylamine
                | ChemicalComponent::HydrogenCyanide
                | ChemicalComponent::NitricOxide
                | ChemicalComponent::NitrogenDioxide
                | ChemicalComponent::NitricAcid
                | ChemicalComponent::SulfurDioxide
                | ChemicalComponent::SulfuricAcid => false,
                _ => true,
            },
            // Weak
            (_, MagneticFieldStrength::Weak) => match self {
                ChemicalComponent::Methane
                | ChemicalComponent::Ammonia
                | ChemicalComponent::HydrogenSulfide
                | ChemicalComponent::Formaldehyde
                | ChemicalComponent::Ethylene
                | ChemicalComponent::Ethane
                | ChemicalComponent::Acetylene
                | ChemicalComponent::Methylamine
                | ChemicalComponent::HydrogenCyanide
                | ChemicalComponent::NitricOxide
                | ChemicalComponent::NitrogenDioxide
                | ChemicalComponent::NitricAcid
                | ChemicalComponent::SulfurDioxide
                | ChemicalComponent::SulfuricAcid => false,
                _ => true,
            },
            // Moderate
            (_, MagneticFieldStrength::Moderate) => match self {
                ChemicalComponent::Ammonia
                | ChemicalComponent::HydrogenSulfide
                | ChemicalComponent::Formaldehyde
                | ChemicalComponent::Ethylene
                | ChemicalComponent::Acetylene
                | ChemicalComponent::HydrogenCyanide
                | ChemicalComponent::NitricOxide
                | ChemicalComponent::NitrogenDioxide
                | ChemicalComponent::NitricAcid
                | ChemicalComponent::SulfuricAcid => false,
                _ => true,
            },
            // Strong
            (_, MagneticFieldStrength::Strong) => match self {
                ChemicalComponent::HydrogenSulfide
                | ChemicalComponent::Formaldehyde
                | ChemicalComponent::Acetylene
                | ChemicalComponent::NitricOxide
                | ChemicalComponent::NitrogenDioxide
                | ChemicalComponent::SulfuricAcid => false,
                _ => true,
            },
            // Very Strong
            (_, MagneticFieldStrength::VeryStrong) => match self {
                ChemicalComponent::Formaldehyde
                | ChemicalComponent::NitricOxide
                | ChemicalComponent::NitrogenDioxide => false,
                _ => true,
            },
            // Extreme
            (_, MagneticFieldStrength::Extreme) => match self {
                ChemicalComponent::Formaldehyde => false,
                _ => true,
            },
        }
    }

    pub fn is_related_element(&self, other: &ChemicalComponent) -> bool {
        let related_pairs = vec![
            // Hydrogen compounds
            (ChemicalComponent::Water, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Water, ChemicalComponent::Oxygen),
            (ChemicalComponent::Methane, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Methane, ChemicalComponent::Carbon),
            (ChemicalComponent::Ammonia, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Ammonia, ChemicalComponent::Nitrogen),
            (
                ChemicalComponent::HydrogenSulfide,
                ChemicalComponent::Hydrogen,
            ),
            (
                ChemicalComponent::HydrogenSulfide,
                ChemicalComponent::Sulfur,
            ),
            (ChemicalComponent::Methanol, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Methanol, ChemicalComponent::Carbon),
            (ChemicalComponent::Methanol, ChemicalComponent::Oxygen),
            (ChemicalComponent::Formaldehyde, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Formaldehyde, ChemicalComponent::Carbon),
            (ChemicalComponent::Formaldehyde, ChemicalComponent::Oxygen),
            (ChemicalComponent::Ethylene, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Ethylene, ChemicalComponent::Carbon),
            (ChemicalComponent::Ethane, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Ethane, ChemicalComponent::Carbon),
            (ChemicalComponent::Acetylene, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Acetylene, ChemicalComponent::Carbon),
            (ChemicalComponent::Benzene, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Benzene, ChemicalComponent::Carbon),
            (ChemicalComponent::Acetonitrile, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Acetonitrile, ChemicalComponent::Carbon),
            (ChemicalComponent::Acetonitrile, ChemicalComponent::Nitrogen),
            (ChemicalComponent::Methylamine, ChemicalComponent::Hydrogen),
            (ChemicalComponent::Methylamine, ChemicalComponent::Carbon),
            (ChemicalComponent::Methylamine, ChemicalComponent::Nitrogen),
            (
                ChemicalComponent::HydrogenCyanide,
                ChemicalComponent::Hydrogen,
            ),
            (
                ChemicalComponent::HydrogenCyanide,
                ChemicalComponent::Carbon,
            ),
            (
                ChemicalComponent::HydrogenCyanide,
                ChemicalComponent::Nitrogen,
            ),
            // Carbon compounds
            (ChemicalComponent::CarbonDioxide, ChemicalComponent::Carbon),
            (ChemicalComponent::CarbonDioxide, ChemicalComponent::Oxygen),
            (ChemicalComponent::CarbonMonoxide, ChemicalComponent::Carbon),
            (ChemicalComponent::CarbonMonoxide, ChemicalComponent::Oxygen),
            (
                ChemicalComponent::PolycyclicAromaticHydrocarbons,
                ChemicalComponent::Carbon,
            ),
            (
                ChemicalComponent::PolycyclicAromaticHydrocarbons,
                ChemicalComponent::Hydrogen,
            ),
            // Sulfur compounds
            (ChemicalComponent::SulfurDioxide, ChemicalComponent::Sulfur),
            (ChemicalComponent::SulfurDioxide, ChemicalComponent::Oxygen),
            (ChemicalComponent::SulfuricAcid, ChemicalComponent::Sulfur),
            (ChemicalComponent::SulfuricAcid, ChemicalComponent::Oxygen),
            (ChemicalComponent::SulfuricAcid, ChemicalComponent::Water),
            (ChemicalComponent::SulfuricAcid, ChemicalComponent::Hydrogen),
            (
                ChemicalComponent::HydrogenSulfide,
                ChemicalComponent::Hydrogen,
            ),
            (
                ChemicalComponent::HydrogenSulfide,
                ChemicalComponent::Sulfur,
            ),
            // Nitrogen compounds
            (ChemicalComponent::NitricOxide, ChemicalComponent::Nitrogen),
            (ChemicalComponent::NitricOxide, ChemicalComponent::Oxygen),
            (
                ChemicalComponent::NitrogenDioxide,
                ChemicalComponent::Nitrogen,
            ),
            (
                ChemicalComponent::NitrogenDioxide,
                ChemicalComponent::Oxygen,
            ),
            (ChemicalComponent::NitricAcid, ChemicalComponent::Nitrogen),
            (ChemicalComponent::NitricAcid, ChemicalComponent::Oxygen),
            (ChemicalComponent::NitricAcid, ChemicalComponent::Hydrogen),
            // Other compound relationships
            (ChemicalComponent::CalciumOxide, ChemicalComponent::Calcium),
            (ChemicalComponent::CalciumOxide, ChemicalComponent::Oxygen),
            (ChemicalComponent::SodiumChloride, ChemicalComponent::Sodium),
            (
                ChemicalComponent::SodiumChloride,
                ChemicalComponent::Chlorine,
            ),
            (ChemicalComponent::SodiumOxide, ChemicalComponent::Sodium),
            (ChemicalComponent::SodiumOxide, ChemicalComponent::Oxygen),
            (
                ChemicalComponent::PotassiumChloride,
                ChemicalComponent::Potassium,
            ),
            (
                ChemicalComponent::PotassiumChloride,
                ChemicalComponent::Chlorine,
            ),
            (
                ChemicalComponent::PotassiumOxide,
                ChemicalComponent::Potassium,
            ),
            (ChemicalComponent::PotassiumOxide, ChemicalComponent::Oxygen),
            (
                ChemicalComponent::MagnesiumOxide,
                ChemicalComponent::Magnesium,
            ),
            (ChemicalComponent::MagnesiumOxide, ChemicalComponent::Oxygen),
            (
                ChemicalComponent::TitaniumDioxide,
                ChemicalComponent::Titanium,
            ),
            (
                ChemicalComponent::TitaniumDioxide,
                ChemicalComponent::Oxygen,
            ),
            (
                ChemicalComponent::TitaniumTetrachloride,
                ChemicalComponent::Titanium,
            ),
            (
                ChemicalComponent::TitaniumTetrachloride,
                ChemicalComponent::Chlorine,
            ),
            (
                ChemicalComponent::PhosphorusPentoxide,
                ChemicalComponent::Phosphorus,
            ),
            (
                ChemicalComponent::PhosphorusPentoxide,
                ChemicalComponent::Oxygen,
            ),
            (
                ChemicalComponent::ChromiumOxide,
                ChemicalComponent::Chromium,
            ),
            (ChemicalComponent::ChromiumOxide, ChemicalComponent::Oxygen),
            (
                ChemicalComponent::ChromiumChloride,
                ChemicalComponent::Chromium,
            ),
            (
                ChemicalComponent::ChromiumChloride,
                ChemicalComponent::Chlorine,
            ),
            (
                ChemicalComponent::ManganeseDioxide,
                ChemicalComponent::Manganese,
            ),
            (
                ChemicalComponent::ManganeseDioxide,
                ChemicalComponent::Oxygen,
            ),
            (
                ChemicalComponent::ManganeseOxide,
                ChemicalComponent::Manganese,
            ),
            (ChemicalComponent::ManganeseOxide, ChemicalComponent::Oxygen),
            (ChemicalComponent::NickelOxide, ChemicalComponent::Nickel),
            (ChemicalComponent::NickelOxide, ChemicalComponent::Oxygen),
            (ChemicalComponent::NickelSulfide, ChemicalComponent::Nickel),
            (ChemicalComponent::NickelSulfide, ChemicalComponent::Sulfur),
        ];

        related_pairs.contains(&(*self, *other)) || related_pairs.contains(&(*other, *self))
    }

    pub fn thermal_decomposition_temperature(&self) -> Option<u32> {
        match self {
            ChemicalComponent::Water => Some(2000),
            ChemicalComponent::CarbonMonoxide => Some(1200),
            ChemicalComponent::CarbonDioxide => Some(2000),
            ChemicalComponent::SiliconDioxide => Some(2500),
            ChemicalComponent::AluminiumOxide => Some(2980),
            ChemicalComponent::IronOxide => Some(3000),
            ChemicalComponent::IronSulfide => Some(1100),
            ChemicalComponent::CalciumOxide => Some(2850),
            ChemicalComponent::SodiumChloride => Some(1465),
            ChemicalComponent::SodiumOxide => Some(1250),
            ChemicalComponent::PotassiumChloride => Some(1500),
            ChemicalComponent::PotassiumOxide => Some(1100),
            ChemicalComponent::MagnesiumOxide => Some(2800),
            ChemicalComponent::TitaniumDioxide => Some(2900),
            ChemicalComponent::TitaniumTetrachloride => Some(200),
            ChemicalComponent::PhosphorusPentoxide => Some(550),
            ChemicalComponent::ChromiumOxide => Some(3000),
            ChemicalComponent::ChromiumChloride => Some(1800),
            ChemicalComponent::ManganeseDioxide => Some(535),
            ChemicalComponent::ManganeseOxide => Some(1780),
            ChemicalComponent::NickelOxide => Some(2100),
            ChemicalComponent::NickelSulfide => Some(1000),
            ChemicalComponent::NitricOxide => Some(800),
            ChemicalComponent::NitrogenDioxide => Some(700),
            ChemicalComponent::NitricAcid => Some(500),
            ChemicalComponent::SulfuricAcid => Some(500),
            ChemicalComponent::Methane => Some(1000),
            ChemicalComponent::Ammonia => Some(700),
            ChemicalComponent::HydrogenSulfide => Some(600),
            ChemicalComponent::SulfurDioxide => Some(600),
            ChemicalComponent::Hydroxyl => Some(500),
            ChemicalComponent::Formaldehyde => Some(600),
            ChemicalComponent::Methanol => Some(700),
            ChemicalComponent::Ethylene => Some(800),
            ChemicalComponent::Ethane => Some(600),
            ChemicalComponent::Acetylene => Some(700),
            ChemicalComponent::Benzene => Some(1100),
            ChemicalComponent::Acetonitrile => Some(850),
            ChemicalComponent::Methylamine => Some(500),
            ChemicalComponent::HydrogenCyanide => Some(800),
            ChemicalComponent::Glycine => Some(900),
            ChemicalComponent::Silicates => Some(1700),
            ChemicalComponent::PolycyclicAromaticHydrocarbons => Some(1200),
            _ => None,
        }
    }

    pub fn can_stay_stable_at_temperature(&self, temperature: u32) -> bool {
        if let Some(decomposition_temp) = self.thermal_decomposition_temperature() {
            temperature < decomposition_temp
        } else {
            true // Assume stable if no decomposition temperature is listed
        }
    }

    pub fn dissociation_products(&self) -> Vec<ChemicalComponent> {
        match self {
            ChemicalComponent::Water => {
                vec![ChemicalComponent::Hydrogen, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::Methane => {
                vec![ChemicalComponent::Carbon, ChemicalComponent::Hydrogen]
            }
            ChemicalComponent::Ammonia => {
                vec![ChemicalComponent::Nitrogen, ChemicalComponent::Hydrogen]
            }
            ChemicalComponent::HydrogenSulfide => {
                vec![ChemicalComponent::Sulfur, ChemicalComponent::Hydrogen]
            }
            ChemicalComponent::CarbonMonoxide => {
                vec![ChemicalComponent::Carbon, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::Methanol => {
                vec![ChemicalComponent::Methane, ChemicalComponent::Water]
            }
            ChemicalComponent::Formaldehyde => vec![
                ChemicalComponent::CarbonMonoxide,
                ChemicalComponent::Hydrogen,
            ],
            ChemicalComponent::Ethylene => {
                vec![ChemicalComponent::Carbon, ChemicalComponent::Hydrogen]
            }
            ChemicalComponent::Ethane => {
                vec![ChemicalComponent::Ethylene, ChemicalComponent::Hydrogen]
            }
            ChemicalComponent::Acetylene => {
                vec![ChemicalComponent::Carbon, ChemicalComponent::Hydrogen]
            }
            ChemicalComponent::Methylamine => {
                vec![ChemicalComponent::Methane, ChemicalComponent::Ammonia]
            }
            ChemicalComponent::HydrogenCyanide => vec![
                ChemicalComponent::Carbon,
                ChemicalComponent::Nitrogen,
                ChemicalComponent::Hydrogen,
            ],
            ChemicalComponent::NitricOxide => {
                vec![ChemicalComponent::Nitrogen, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::NitrogenDioxide => {
                vec![ChemicalComponent::NitricOxide, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::NitricAcid => {
                vec![ChemicalComponent::NitrogenDioxide, ChemicalComponent::Water]
            }
            ChemicalComponent::SulfurDioxide => {
                vec![ChemicalComponent::Sulfur, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::SulfuricAcid => {
                vec![ChemicalComponent::SulfurDioxide, ChemicalComponent::Water]
            }
            ChemicalComponent::Hydroxyl => {
                vec![ChemicalComponent::Hydrogen, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::CarbonDioxide => {
                vec![ChemicalComponent::Carbon, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::IronOxide => {
                vec![ChemicalComponent::Iron, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::IronSulfide => {
                vec![ChemicalComponent::Iron, ChemicalComponent::Sulfur]
            }
            ChemicalComponent::CalciumOxide => {
                vec![ChemicalComponent::Calcium, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::SodiumChloride => {
                vec![ChemicalComponent::Sodium, ChemicalComponent::Chlorine]
            }
            ChemicalComponent::SodiumOxide => {
                vec![ChemicalComponent::Sodium, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::PotassiumChloride => {
                vec![ChemicalComponent::Potassium, ChemicalComponent::Chlorine]
            }
            ChemicalComponent::PotassiumOxide => {
                vec![ChemicalComponent::Potassium, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::MagnesiumOxide => {
                vec![ChemicalComponent::Magnesium, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::TitaniumDioxide => {
                vec![ChemicalComponent::Titanium, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::TitaniumTetrachloride => {
                vec![ChemicalComponent::Titanium, ChemicalComponent::Chlorine]
            }
            ChemicalComponent::PhosphorusPentoxide => {
                vec![ChemicalComponent::Phosphorus, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::ChromiumOxide => {
                vec![ChemicalComponent::Chromium, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::ChromiumChloride => {
                vec![ChemicalComponent::Chromium, ChemicalComponent::Chlorine]
            }
            ChemicalComponent::ManganeseDioxide => {
                vec![ChemicalComponent::Manganese, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::ManganeseOxide => {
                vec![ChemicalComponent::Manganese, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::NickelOxide => {
                vec![ChemicalComponent::Nickel, ChemicalComponent::Oxygen]
            }
            ChemicalComponent::NickelSulfide => {
                vec![ChemicalComponent::Nickel, ChemicalComponent::Sulfur]
            }
            _ => vec![],
        }
    }

    /// Determines if a substance can exist in liquid state at given temperature and pressure.
    ///
    /// # Parameters
    /// - `substance`: The substance to check.
    /// - `temperature`: The temperature in Kelvin.
    /// - `pressure`: The pressure in atm.
    ///
    /// # Returns
    /// `true` if the substance can exist in liquid state, `false` otherwise.
    ///
    /// # Example
    /// ```
    /// use crate::planet_generator::prelude::*;
    ///
    /// let water = ChemicalComponent::Water;
    /// let can_exist = water.can_exist_as_liquid(280, 1.0);
    /// println!("Can water exist as liquid at 280 K and 1 atm? {}", can_exist);
    /// ```
    pub fn can_exist_as_liquid(&self, temperature: u32, pressure: f32) -> bool {
        let temperature = temperature as f64;
        let pressure = pressure as f64;
        if let Some((triple_point_temp, triple_point_pressure)) = self.triple_point() {
            if let Some(boiling_point_temp) = self.boiling_point() {
                return temperature > triple_point_temp
                    && temperature < boiling_point_temp
                    && pressure > triple_point_pressure;
            }
        }
        false
    }

    /// Determines if a substance can exist in gas state at given temperature and pressure.
    ///
    /// # Parameters
    /// - `temperature`: The temperature in Kelvin.
    /// - `pressure`: The pressure in atm.
    ///
    /// # Returns
    /// `true` if the substance can exist in gas state, `false` otherwise.
    ///
    /// # Example
    /// ```
    /// use crate::planet_generator::prelude::*;
    ///
    /// let water = ChemicalComponent::Water;
    /// let can_exist = water.can_exist_as_gas(400, 1.0);
    /// println!("Can water exist as gas at 400 K and 1 atm? {}", can_exist);
    /// ```
    pub fn can_exist_as_gas(&self, temperature: u32, pressure: f32) -> bool {
        let temperature = temperature as f64;
        let pressure = pressure as f64;
        if let Some((triple_point_temp, triple_point_pressure)) = self.triple_point() {
            if let Some(boiling_point_temp) = self.boiling_point() {
                // Check if the temperature is above the boiling point or above the triple point temperature
                if temperature > boiling_point_temp as f64 {
                    return true;
                } else if temperature > triple_point_temp && pressure < triple_point_pressure {
                    return true;
                }
            } else if temperature > triple_point_temp {
                // For substances that sublimate, check if the temperature is above the triple point temperature
                return true;
            }
        }
        false
    }

    pub fn can_be_retained_as_atmospheric_gas(
        &self,
        temperature: u32,
        pressure: f32,
        magnetic_field: MagneticFieldStrength,
    ) -> bool {
        self.can_exist_as_gas(temperature, pressure)
            && self.can_exist_in_atmosphere(magnetic_field)
            && self.can_stay_stable_at_temperature(temperature)
    }

    pub fn components_liquid_at(temperature: u32, pressure: f32) -> Option<Vec<ChemicalComponent>> {
        let mut possible_liquids = Vec::new();
        let mut possible_lowest_pressure = None;

        if (ChemicalComponent::Water.can_exist_as_liquid(temperature, pressure)) {
            possible_liquids.push(ChemicalComponent::Water);
        } else {
            for component in ChemicalComponent::iter() {
                if component.can_exist_as_liquid(temperature, pressure) {
                    possible_liquids.push(component);
                } else if let Some((triple_point_temp, triple_point_pressure)) =
                    component.triple_point()
                {
                    if component.can_exist_as_liquid(temperature, triple_point_pressure as f32) {
                        if let Some((_, lowest_pressure)) = possible_lowest_pressure {
                            if triple_point_pressure < lowest_pressure {
                                possible_lowest_pressure = Some((component, triple_point_pressure));
                            }
                        } else {
                            possible_lowest_pressure = Some((component, triple_point_pressure));
                        }
                    }
                }
            }
        }

        if !possible_liquids.is_empty() {
            return Some(possible_liquids);
        } else if let Some((component, _)) = possible_lowest_pressure {
            return Some(vec![component]);
        }

        None
    }

    pub fn molecular_weight_kg(&self) -> f64 {
        const AMU_TO_KG: f64 = 1.66053906660e-27;
        self.molecular_weight_amu() * AMU_TO_KG
    }
}

impl Display for ChemicalComponent {
    fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
        match self {
            ChemicalComponent::Hydrogen => write!(f, "Hydrogen"),
            ChemicalComponent::Helium => write!(f, "Helium"),
            ChemicalComponent::Carbon => write!(f, "Carbon"),
            ChemicalComponent::Nitrogen => write!(f, "Nitrogen"),
            ChemicalComponent::Oxygen => write!(f, "Oxygen"),
            ChemicalComponent::Silicon => write!(f, "Silicon"),
            ChemicalComponent::Magnesium => write!(f, "Magnesium"),
            ChemicalComponent::Iron => write!(f, "Iron"),
            ChemicalComponent::Sulfur => write!(f, "Sulfur"),
            ChemicalComponent::Sodium => write!(f, "Sodium"),
            ChemicalComponent::Potassium => write!(f, "Potassium"),
            ChemicalComponent::Calcium => write!(f, "Calcium"),
            ChemicalComponent::Aluminum => write!(f, "Aluminum"),
            ChemicalComponent::Phosphorus => write!(f, "Phosphorus"),
            ChemicalComponent::Chlorine => write!(f, "Chlorine"),
            ChemicalComponent::Argon => write!(f, "Argon"),
            ChemicalComponent::Titanium => write!(f, "Titanium"),
            ChemicalComponent::Neon => write!(f, "Neon"),
            ChemicalComponent::Chromium => write!(f, "Chromium"),
            ChemicalComponent::Manganese => write!(f, "Manganese"),
            ChemicalComponent::Nickel => write!(f, "Nickel"),
            ChemicalComponent::Water => write!(f, "Water"),
            ChemicalComponent::CarbonMonoxide => write!(f, "Carbon Monoxide"),
            ChemicalComponent::CarbonDioxide => write!(f, "Carbon Dioxide"),
            ChemicalComponent::Methane => write!(f, "Methane"),
            ChemicalComponent::Ammonia => write!(f, "Ammonia"),
            ChemicalComponent::HydrogenSulfide => write!(f, "Hydrogen Sulfide"),
            ChemicalComponent::SulfurDioxide => write!(f, "Sulfur Dioxide"),
            ChemicalComponent::Hydroxyl => write!(f, "Hydroxyl"),
            ChemicalComponent::NitricOxide => write!(f, "Nitric Oxide"),
            ChemicalComponent::NitrogenDioxide => write!(f, "Nitrogen Dioxide"),
            ChemicalComponent::Formaldehyde => write!(f, "Formaldehyde"),
            ChemicalComponent::Methanol => write!(f, "Methanol"),
            ChemicalComponent::Ethylene => write!(f, "Ethylene"),
            ChemicalComponent::Ethane => write!(f, "Ethane"),
            ChemicalComponent::Acetylene => write!(f, "Acetylene"),
            ChemicalComponent::Benzene => write!(f, "Benzene"),
            ChemicalComponent::Acetonitrile => write!(f, "Acetonitrile"),
            ChemicalComponent::Methylamine => write!(f, "Methylamine"),
            ChemicalComponent::HydrogenCyanide => write!(f, "Hydrogen Cyanide"),
            ChemicalComponent::Glycine => write!(f, "Glycine"),
            ChemicalComponent::Silicates => write!(f, "Silicates"),
            ChemicalComponent::PolycyclicAromaticHydrocarbons => {
                write!(f, "Polycyclic Aromatic Hydrocarbons")
            }
            ChemicalComponent::SiliconDioxide => write!(f, "Silicon Dioxide"),
            ChemicalComponent::AluminiumOxide => write!(f, "Aluminium Oxide"),
            ChemicalComponent::IronOxide => write!(f, "Iron Oxide"),
            ChemicalComponent::IronSulfide => write!(f, "Iron Sulfide"),
            ChemicalComponent::CalciumOxide => write!(f, "Calcium Oxide"),
            ChemicalComponent::SodiumChloride => write!(f, "Sodium Chloride"),
            ChemicalComponent::SodiumOxide => write!(f, "Sodium Oxide"),
            ChemicalComponent::PotassiumChloride => write!(f, "Potassium Chloride"),
            ChemicalComponent::PotassiumOxide => write!(f, "Potassium Oxide"),
            ChemicalComponent::MagnesiumOxide => write!(f, "Magnesium Oxide"),
            ChemicalComponent::TitaniumDioxide => write!(f, "Titanium Dioxide"),
            ChemicalComponent::TitaniumTetrachloride => write!(f, "Titanium Tetrachloride"),
            ChemicalComponent::PhosphorusPentoxide => write!(f, "Phosphorus Pentoxide"),
            ChemicalComponent::ChromiumOxide => write!(f, "Chromium Oxide"),
            ChemicalComponent::ChromiumChloride => write!(f, "Chromium Chloride"),
            ChemicalComponent::ManganeseDioxide => write!(f, "Manganese Dioxide"),
            ChemicalComponent::ManganeseOxide => write!(f, "Manganese Oxide"),
            ChemicalComponent::NickelOxide => write!(f, "Nickel Oxide"),
            ChemicalComponent::NickelSulfide => write!(f, "Nickel Sulfide"),
            ChemicalComponent::NitricAcid => write!(f, "Nitric Acid"),
            ChemicalComponent::SulfuricAcid => write!(f, "Sulfuric Acid"),
        }
    }
}

pub const NON_METALS_ELEMENTS: [ChemicalComponent; 2] =
    [ChemicalComponent::Hydrogen, ChemicalComponent::Helium];

pub const MOST_COMMON_ELEMENTS: [ChemicalComponent; 24] = [
    ChemicalComponent::Hydrogen,
    ChemicalComponent::Helium,
    ChemicalComponent::Carbon,
    ChemicalComponent::Nitrogen,
    ChemicalComponent::Oxygen,
    ChemicalComponent::Silicon,
    ChemicalComponent::Magnesium,
    ChemicalComponent::Iron,
    ChemicalComponent::Sulfur,
    ChemicalComponent::Sodium,
    ChemicalComponent::Potassium,
    ChemicalComponent::Calcium,
    ChemicalComponent::Aluminum,
    ChemicalComponent::Neon,
    ChemicalComponent::Phosphorus,
    ChemicalComponent::Chlorine,
    ChemicalComponent::Argon,
    ChemicalComponent::Water,
    ChemicalComponent::CarbonMonoxide,
    ChemicalComponent::CarbonDioxide,
    ChemicalComponent::Methane,
    ChemicalComponent::Ammonia,
    ChemicalComponent::HydrogenSulfide,
    ChemicalComponent::SulfurDioxide,
];

pub const ALL_ELEMENTS: [ChemicalComponent; 64] = [
    ChemicalComponent::Hydrogen,
    ChemicalComponent::Helium,
    ChemicalComponent::Carbon,
    ChemicalComponent::Nitrogen,
    ChemicalComponent::Oxygen,
    ChemicalComponent::Silicon,
    ChemicalComponent::Magnesium,
    ChemicalComponent::Iron,
    ChemicalComponent::Sulfur,
    ChemicalComponent::Sodium,
    ChemicalComponent::Potassium,
    ChemicalComponent::Calcium,
    ChemicalComponent::Aluminum,
    ChemicalComponent::Phosphorus,
    ChemicalComponent::Chlorine,
    ChemicalComponent::Argon,
    ChemicalComponent::Titanium,
    ChemicalComponent::Chromium,
    ChemicalComponent::Manganese,
    ChemicalComponent::Nickel,
    ChemicalComponent::Water,
    ChemicalComponent::CarbonMonoxide,
    ChemicalComponent::CarbonDioxide,
    ChemicalComponent::Methane,
    ChemicalComponent::Ammonia,
    ChemicalComponent::Neon,
    ChemicalComponent::HydrogenSulfide,
    ChemicalComponent::SulfurDioxide,
    ChemicalComponent::Hydroxyl,
    ChemicalComponent::NitricOxide,
    ChemicalComponent::NitrogenDioxide,
    ChemicalComponent::Formaldehyde,
    ChemicalComponent::Methanol,
    ChemicalComponent::Ethylene,
    ChemicalComponent::Ethane,
    ChemicalComponent::Acetylene,
    ChemicalComponent::Benzene,
    ChemicalComponent::Acetonitrile,
    ChemicalComponent::Methylamine,
    ChemicalComponent::HydrogenCyanide,
    ChemicalComponent::Glycine,
    ChemicalComponent::Silicates,
    ChemicalComponent::PolycyclicAromaticHydrocarbons,
    ChemicalComponent::SiliconDioxide,
    ChemicalComponent::AluminiumOxide,
    ChemicalComponent::IronOxide,
    ChemicalComponent::IronSulfide,
    ChemicalComponent::CalciumOxide,
    ChemicalComponent::SodiumChloride,
    ChemicalComponent::SodiumOxide,
    ChemicalComponent::PotassiumChloride,
    ChemicalComponent::PotassiumOxide,
    ChemicalComponent::MagnesiumOxide,
    ChemicalComponent::TitaniumDioxide,
    ChemicalComponent::TitaniumTetrachloride,
    ChemicalComponent::PhosphorusPentoxide,
    ChemicalComponent::ChromiumOxide,
    ChemicalComponent::ChromiumChloride,
    ChemicalComponent::ManganeseDioxide,
    ChemicalComponent::ManganeseOxide,
    ChemicalComponent::NickelOxide,
    ChemicalComponent::NickelSulfide,
    ChemicalComponent::NitricAcid,
    ChemicalComponent::SulfuricAcid,
];

#[derive(
    Copy, Clone, Eq, PartialEq, Ord, PartialOrd, Debug, Hash, SmartDefault, Serialize, Deserialize,
)]
pub enum ElementPresenceOccurrence {
    Absence,
    VeryLow,
    Low,
    #[default]
    Normal,
    High,
    VeryHigh,
    Omnipresence,
}

impl Display for ElementPresenceOccurrence {
    fn fmt(&self, f: &mut std::fmt::Formatter<'_>) -> std::fmt::Result {
        match self {
            ElementPresenceOccurrence::Absence => write!(f, "Absence"),
            ElementPresenceOccurrence::VeryLow => write!(f, "Very Low Occurrence"),
            ElementPresenceOccurrence::Low => write!(f, "Low Occurrence"),
            ElementPresenceOccurrence::Normal => write!(f, "Normal Occurrence"),
            ElementPresenceOccurrence::High => write!(f, "High Occurrence"),
            ElementPresenceOccurrence::VeryHigh => write!(f, "Very High Occurrence"),
            ElementPresenceOccurrence::Omnipresence => write!(f, "Omnipresence"),
        }
    }
}

pub(crate) fn generate_random_non_metal_element(rng: &mut SeededDiceRoller) -> ChemicalComponent {
    NON_METALS_ELEMENTS[rng.gen_range(0..NON_METALS_ELEMENTS.len())]
}

pub(crate) fn generate_random_common_element(rng: &mut SeededDiceRoller) -> ChemicalComponent {
    MOST_COMMON_ELEMENTS[rng.gen_range(0..MOST_COMMON_ELEMENTS.len())]
}

pub(crate) fn generate_random_element(rng: &mut SeededDiceRoller) -> ChemicalComponent {
    ALL_ELEMENTS[rng.gen_range(0..ALL_ELEMENTS.len())]
}

pub(crate) fn liquid_majority_composition_likelihood(
    component: ChemicalComponent,
    peculiarities: &[StarPeculiarity],
) -> f64 {
    let base_likelihood = match component {
        ChemicalComponent::Water => 1.0,
        ChemicalComponent::Methane => 0.8,
        ChemicalComponent::Ammonia => 0.6,
        ChemicalComponent::CarbonDioxide => 0.5,
        ChemicalComponent::SulfurDioxide => 0.4,
        ChemicalComponent::Formaldehyde => 0.2,
        ChemicalComponent::Hydrogen => 0.3,
        ChemicalComponent::Helium => 0.1,
        ChemicalComponent::Carbon => 0.2,
        ChemicalComponent::Nitrogen => 0.4,
        ChemicalComponent::Oxygen => 0.3,
        ChemicalComponent::Silicon => 0.1,
        ChemicalComponent::Magnesium => 0.1,
        ChemicalComponent::Iron => 0.1,
        ChemicalComponent::Sulfur => 0.3,
        ChemicalComponent::Sodium => 0.1,
        ChemicalComponent::Potassium => 0.1,
        ChemicalComponent::Calcium => 0.1,
        ChemicalComponent::Aluminum => 0.1,
        ChemicalComponent::Phosphorus => 0.2,
        ChemicalComponent::Chlorine => 0.3,
        ChemicalComponent::Argon => 0.1,
        ChemicalComponent::Titanium => 0.1,
        ChemicalComponent::Chromium => 0.1,
        ChemicalComponent::Manganese => 0.1,
        ChemicalComponent::Nickel => 0.1,
        ChemicalComponent::HydrogenSulfide => 0.3,
        ChemicalComponent::Hydroxyl => 0.1,
        ChemicalComponent::NitricOxide => 0.2,
        ChemicalComponent::NitrogenDioxide => 0.2,
        ChemicalComponent::Methanol => 0.4,
        ChemicalComponent::Ethylene => 0.3,
        ChemicalComponent::Ethane => 0.4,
        ChemicalComponent::Neon => 0.2,
        ChemicalComponent::Acetylene => 0.2,
        ChemicalComponent::Benzene => 0.3,
        ChemicalComponent::Acetonitrile => 0.2,
        ChemicalComponent::Methylamine => 0.2,
        ChemicalComponent::HydrogenCyanide => 0.2,
        ChemicalComponent::Glycine => 0.1,
        ChemicalComponent::Silicates => 0.1,
        ChemicalComponent::PolycyclicAromaticHydrocarbons => 0.2,
        ChemicalComponent::CarbonMonoxide => 0.4,
        ChemicalComponent::SiliconDioxide => 0.1,
        ChemicalComponent::AluminiumOxide => 0.1,
        ChemicalComponent::IronOxide => 0.1,
        ChemicalComponent::IronSulfide => 0.1,
        ChemicalComponent::CalciumOxide => 0.1,
        ChemicalComponent::SodiumChloride => 0.1,
        ChemicalComponent::SodiumOxide => 0.1,
        ChemicalComponent::PotassiumChloride => 0.1,
        ChemicalComponent::PotassiumOxide => 0.1,
        ChemicalComponent::MagnesiumOxide => 0.1,
        ChemicalComponent::TitaniumDioxide => 0.1,
        ChemicalComponent::TitaniumTetrachloride => 0.2,
        ChemicalComponent::PhosphorusPentoxide => 0.1,
        ChemicalComponent::ChromiumOxide => 0.1,
        ChemicalComponent::ChromiumChloride => 0.1,
        ChemicalComponent::ManganeseDioxide => 0.1,
        ChemicalComponent::ManganeseOxide => 0.1,
        ChemicalComponent::NickelOxide => 0.1,
        ChemicalComponent::NickelSulfide => 0.1,
        ChemicalComponent::NitricAcid => 0.2,
        ChemicalComponent::SulfuricAcid => 0.2,
    };

    let stability = match component {
        ChemicalComponent::Hydrogen => 0.7,
        ChemicalComponent::Helium => 1.0,
        ChemicalComponent::Carbon => 0.9,
        ChemicalComponent::Nitrogen => 0.9,
        ChemicalComponent::Oxygen => 0.8,
        ChemicalComponent::Silicon => 0.8,
        ChemicalComponent::Magnesium => 0.7,
        ChemicalComponent::Iron => 0.7,
        ChemicalComponent::Sulfur => 0.6,
        ChemicalComponent::Sodium => 0.5,
        ChemicalComponent::Potassium => 0.5,
        ChemicalComponent::Calcium => 0.6,
        ChemicalComponent::Aluminum => 0.7,
        ChemicalComponent::Phosphorus => 0.6,
        ChemicalComponent::Neon => 1.0,
        ChemicalComponent::Chlorine => 0.6,
        ChemicalComponent::Argon => 1.0,
        ChemicalComponent::Titanium => 0.8,
        ChemicalComponent::Chromium => 0.7,
        ChemicalComponent::Manganese => 0.7,
        ChemicalComponent::Nickel => 0.7,
        ChemicalComponent::Water => 1.0,
        ChemicalComponent::CarbonMonoxide => 0.5,
        ChemicalComponent::CarbonDioxide => 0.9,
        ChemicalComponent::Methane => 0.8,
        ChemicalComponent::Ammonia => 0.6,
        ChemicalComponent::HydrogenSulfide => 0.4,
        ChemicalComponent::SulfurDioxide => 0.5,
        ChemicalComponent::Hydroxyl => 0.3,
        ChemicalComponent::NitricOxide => 0.4,
        ChemicalComponent::NitrogenDioxide => 0.4,
        ChemicalComponent::Formaldehyde => 0.2,
        ChemicalComponent::Methanol => 0.6,
        ChemicalComponent::Ethylene => 0.6,
        ChemicalComponent::Ethane => 0.7,
        ChemicalComponent::Acetylene => 0.5,
        ChemicalComponent::Benzene => 0.7,
        ChemicalComponent::Acetonitrile => 0.5,
        ChemicalComponent::Methylamine => 0.4,
        ChemicalComponent::HydrogenCyanide => 0.3,
        ChemicalComponent::Glycine => 0.8,
        ChemicalComponent::Silicates => 0.9,
        ChemicalComponent::PolycyclicAromaticHydrocarbons => 0.6,
        ChemicalComponent::SiliconDioxide => 0.1,
        ChemicalComponent::AluminiumOxide => 0.1,
        ChemicalComponent::IronOxide => 0.1,
        ChemicalComponent::IronSulfide => 0.1,
        ChemicalComponent::CalciumOxide => 0.1,
        ChemicalComponent::SodiumChloride => 0.1,
        ChemicalComponent::SodiumOxide => 0.1,
        ChemicalComponent::PotassiumChloride => 0.1,
        ChemicalComponent::PotassiumOxide => 0.1,
        ChemicalComponent::MagnesiumOxide => 0.1,
        ChemicalComponent::TitaniumDioxide => 0.1,
        ChemicalComponent::TitaniumTetrachloride => 0.2,
        ChemicalComponent::PhosphorusPentoxide => 0.1,
        ChemicalComponent::ChromiumOxide => 0.1,
        ChemicalComponent::ChromiumChloride => 0.1,
        ChemicalComponent::ManganeseDioxide => 0.1,
        ChemicalComponent::ManganeseOxide => 0.1,
        ChemicalComponent::NickelOxide => 0.1,
        ChemicalComponent::NickelSulfide => 0.1,
        ChemicalComponent::NitricAcid => 0.2,
        ChemicalComponent::SulfuricAcid => 0.2,
    };

    let mut adjusted_likelihood = base_likelihood * stability;

    for peculiarity in peculiarities {
        if let StarPeculiarity::UnusualElementPresence((peculiar_component, occurrence)) =
            peculiarity
        {
            if *peculiar_component == component {
                adjusted_likelihood *= get_occurence_adjustment_factor(*occurrence);
            }
        }
    }

    adjusted_likelihood
}

fn get_occurence_adjustment_factor(occurrence: ElementPresenceOccurrence) -> f64 {
    match occurrence {
        ElementPresenceOccurrence::Absence => 0.0,
        ElementPresenceOccurrence::VeryLow => 0.2,
        ElementPresenceOccurrence::Low => 0.5,
        ElementPresenceOccurrence::Normal => 1.0,
        ElementPresenceOccurrence::High => 1.5,
        ElementPresenceOccurrence::VeryHigh => 2.0,
        ElementPresenceOccurrence::Omnipresence => 3.0,
    }
}

#[cfg(test)]
mod tests {
    use super::*;

    #[test]
    fn test_molecular_weight_amu() {
        assert_eq!(ChemicalComponent::Hydrogen.molecular_weight_amu(), 1.008);
        assert_eq!(ChemicalComponent::Helium.molecular_weight_amu(), 4.0026);
        assert_eq!(ChemicalComponent::Oxygen.molecular_weight_amu(), 16.00);
        assert_eq!(ChemicalComponent::Water.molecular_weight_amu(), 18.015);
        assert_eq!(
            ChemicalComponent::CarbonDioxide.molecular_weight_amu(),
            44.01
        );
    }

    #[test]
    fn test_triple_point() {
        assert_eq!(
            ChemicalComponent::Hydrogen.triple_point(),
            Some((13.8, 0.070))
        );
        assert_eq!(
            ChemicalComponent::Helium.triple_point(),
            Some((2.2, 0.0052))
        );
        assert_eq!(
            ChemicalComponent::Oxygen.triple_point(),
            Some((54.36, 0.0015))
        );
        assert_eq!(
            ChemicalComponent::Water.triple_point(),
            Some((273.16, 0.00604))
        );
        assert_eq!(
            ChemicalComponent::CarbonDioxide.triple_point(),
            Some((216.55, 5.11))
        );
        assert_eq!(ChemicalComponent::Silicon.triple_point(), None);
    }

    #[test]
    fn test_boiling_point() {
        assert_eq!(ChemicalComponent::Hydrogen.boiling_point(), Some(20.28));
        assert_eq!(ChemicalComponent::Helium.boiling_point(), Some(4.22));
        assert_eq!(ChemicalComponent::Oxygen.boiling_point(), Some(90.20));
        assert_eq!(ChemicalComponent::Water.boiling_point(), Some(373.16));
        assert_eq!(ChemicalComponent::Silicon.boiling_point(), None);
    }

    #[test]
    fn test_can_exist_as_liquid() {
        assert!(ChemicalComponent::Water.can_exist_as_liquid(280, 1.0));
        assert!(!ChemicalComponent::Water.can_exist_as_liquid(280, 0.001));
        assert!(ChemicalComponent::Oxygen.can_exist_as_liquid(60, 0.002));
        assert!(!ChemicalComponent::Oxygen.can_exist_as_liquid(60, 0.0005));
        assert!(ChemicalComponent::CarbonDioxide.can_exist_as_liquid(220, 6.0));
        assert!(!ChemicalComponent::CarbonDioxide.can_exist_as_liquid(220, 4.0));
    }

    #[test]
    fn test_can_exist_as_gas() {
        assert!(ChemicalComponent::Water.can_exist_as_gas(400, 1.0));
        assert!(ChemicalComponent::Water.can_exist_as_gas(300, 0.005));
        assert!(ChemicalComponent::Oxygen.can_exist_as_gas(100, 0.001));
        assert!(!ChemicalComponent::Oxygen.can_exist_as_gas(50, 1.0));
        assert!(ChemicalComponent::CarbonDioxide.can_exist_as_gas(220, 5.0));
        assert!(!ChemicalComponent::CarbonDioxide.can_exist_as_gas(220, 6.0));
    }

    #[test]
    fn test_components_liquid_at() {
        let components = ChemicalComponent::components_liquid_at(280, 1.0);
        assert!(components.is_some());
        let components = components.unwrap();
        assert!(components.contains(&ChemicalComponent::Water));
        assert!(!components.contains(&ChemicalComponent::Oxygen));

        let components = ChemicalComponent::components_liquid_at(60, 0.002);
        assert!(components.is_some());
        let components = components.unwrap();
        assert!(components.contains(&ChemicalComponent::Oxygen));
        assert!(!components.contains(&ChemicalComponent::Water));

        let components = ChemicalComponent::components_liquid_at(220, 6.0);
        assert!(components.is_some());
        let components = components.unwrap();
        assert!(components.contains(&ChemicalComponent::CarbonDioxide));
    }

    #[test]
    fn test_molecular_weight_kg() {
        let amu_to_kg = |amu: f64| amu * 1.66053906660e-27;
        assert_eq!(
            ChemicalComponent::Hydrogen.molecular_weight_kg(),
            amu_to_kg(1.008)
        );
        assert_eq!(
            ChemicalComponent::Helium.molecular_weight_kg(),
            amu_to_kg(4.0026)
        );
        assert_eq!(
            ChemicalComponent::Oxygen.molecular_weight_kg(),
            amu_to_kg(16.00)
        );
        assert_eq!(
            ChemicalComponent::Water.molecular_weight_kg(),
            amu_to_kg(18.015)
        );
        assert_eq!(
            ChemicalComponent::CarbonDioxide.molecular_weight_kg(),
            amu_to_kg(44.01)
        );
    }
}