#!/usr/bin/python3
# Plot the hardnesses of some pure metals on a scatterplot.

# AI slop output from GPT-5 Mini.

# Data purportedly from:
# <https://www.engineeringtoolbox.com/conductors-d_1381.html>
# <http://hyperphysics.phy-astr.gsu.edu/hbase/Tables/rstiv.html>
# <https://en.wikipedia.org/wiki/Hardnesses_of_the_elements_(data_page)>

import matplotlib.pyplot as plt

# Data (approximate values from the table)
elements = ["Ag", "Cu", "Au", "Al", "Ni", "Fe", "Zn", "Pb", "Ti", "Cr", "W", "Mg"]
conductivity = [6.29e7, 5.96e7, 4.52e7, 3.77e7, 1.43e7, 1.03e7, 1.66e7, 4.5e6, 2.2e6, 7.74e6, 1.89e7, 2.26e7]  # S/m
mohs_hardness = [2.8, 2.8, 2.8, 2.75, 4.0, 4.0, 2.8, 1.5, 6.0, 8.5, 7.5, 2.5]

fig, ax = plt.subplots(figsize=(8,6))
sc = ax.scatter(mohs_hardness, conductivity, s=80, color='tab:blue', edgecolor='k')

# Annotate each point
for i, el in enumerate(elements):
    ax.annotate(el, (mohs_hardness[i]+0.08, conductivity[i]*1.02), fontsize=9)

ax.set_yscale('log')  # conductivity spans orders of magnitude
ax.set_xlabel('Mohs Hardness')
ax.set_ylabel('Electrical Conductivity (S/m)')
ax.set_title('Conductivity vs Mohs Hardness for Selected Pure Metals (approx.)')
ax.grid(True, which='both', ls='--', alpha=0.5)

plt.tight_layout()
plt.savefig('conductivity_vs_hardness.png', dpi=300)
plt.show()