feat(pl): phonon, pourbaix, neb, rdf_msd, wulff and chempot renderers

[Starlitnightly]

Six plots the analysis namespaces produced data for and mv.pl could not draw

Agents fell back to raw matplotlib for these. Each renderer is @registered with a requires contract the probe battery checks by deletion, takes md / level / ax, returns the Axes, puts units on the axes and the level of theory in the title, uses Okabe–Ito for series (viridis for the one surface), and raises a ValueError naming the producer when its slot is absent.

function	reads
mv.pl.phonon	obsm['phonon_dos_{level}'], uns['grids']['phonon_dos'] (method, supercell), obs['n_imaginary_modes_{level}']; dispersion= takes the mv.prop.dispersion bands object on a shared frequency axis
mv.pl.pourbaix	uns['pourbaix'] → maps[name], ph_grid, potential_grid, ph, potential; obs['pourbaix_decomposition']
mv.pl.neb	obsm['neb_profile_{level}'] + grid neb_profile, obs['neb_converged_{level}']
mv.pl.rdf_msd	obsm['rdf_md_{rdf_level}'], obsm['coordination_md_{rdf_level}'], obs['first_shell_{rdf_level}'], obsm['md_msd_trace_{level}'] + grid, obs['diffusivity_{level}']
mv.pl.wulff	uns['wulff'][level][name] → vertices, face_index, face_miller, area_fractions, anisotropy
mv.pl.chempot	uns['chempot_diagram'] (elements, domains[*].vertices); kind='window' reads uns['chempot_limits']

Three producers deposited nothing drawable, extended minimally

• mv.md.run keeps the MSD trace it already fitted D to, as obsm['md_msd_trace_<level>'] on the same picosecond grid as the temperature trace (replaced together on a rerun of a different length). Both traces are now claimed in produces.
• mv.surf.wulff stores the polyhedron as three h5ad-writable arrays (vertices (n,3), face_index (n,), face_miller) plus per-family area_fractions. The bulk.uns['wulff'] claim the README recorded as "never written" is restored — it now is, and the plot depends on it. README row updated.
• mv.thermo.pourbaix evaluates the decomposition energy over a pH/E mesh (ph_range=, potential_range=, n_grid=) on the diagram it already fetched per material and stores it in uns['pourbaix']['maps']. Still network-gated, still in UNPROBEABLE.

Tests

• tests/test_pl_{phonon,pourbaix,neb,rdf_msd,wulff,chempot}.py — real EMT producers where they run offline; the Pourbaix map is synthetic in the producer's own shape.
• Probe cases added to tests/_contract_cases.py (aqueous, barred, diffusing, chempot_mapped, wulffed), so test_no_claim_goes_unprobed stays green.
• One tutorial code cell per function (nb_getting_started, nb_screening, nb_defects_and_diffusion — replaces the mv.pl.spectra NEB workaround —, nb_dynamics, nb_surfaces_and_adsorption), so test_notebook_coverage stays green.
• Full suite: 946 passed, 144 skipped, 0 failed (baseline on main: 914 passed, 143 skipped; +32 new tests, +1 mp-api-gated skip).

3D collections bypass matplotlib ≥ 3.10's autolim, which produced infinite limits on the NaN-padded ragged face array under pytest; both callers set limits themselves.

🤖 Generated with Claude Code