Development

Convert molecular file formats using Open Babel. Handles SMILES, mol2, sdf, pdb, xyz, cif, and 100+ other formats. Also performs 3D coordinate generation and hydrogen addition.

Use when the user asks about OER (oxygen evolution reaction) overpotential, water oxidation catalysis, or the 4-electron water splitting pathway on a surface catalyst.

Use when the user asks about NRR (nitrogen reduction reaction), ammonia synthesis, N2 fixation, or the electrochemical reduction of N2 to NH3 on a catalyst surface.

Route LAMMPS molecular dynamics requests to force-field-specific sub-skills. Use when the user requests LAMMPS with a specific potential type (DeePMD, ReaxFF).

Run LAMMPS molecular dynamics with DeePMD-kit machine learning potentials. Use when the user wants MD simulations driven by a trained DP model.

Automatically logs new discoveries to CLAUDE.md files. Triggers when fixing bugs, discovering pitfalls, finding performance issues, or learning new patterns in the codebase.

Use when the user asks about HER (hydrogen evolution reaction), hydrogen adsorption free energy, or volcano plot descriptor for HER catalysts.

Generate and manage GPAW Python-based DFT calculations. Use when the user requests GPAW, Python DFT, real-space grid DFT, or LCAO-DFT with ASE integration.

Use when the user asks for Gibbs free energy, zero-point energy (ZPE), thermal corrections, or thermodynamic properties from DFT + frequency data.

Generate and manage Gaussian calculations. Use when the user requests Gaussian, G16, GJF files, or needs hybrid functionals (B3LYP), MP2, CCSD(T), or molecular quantum chemistry with Gaussian basis sets.

Use when the user asks to generate a reaction energy diagram, free energy profile, potential energy surface plot, or pathway comparison diagram for catalysis or reaction mechanism studies.

Convert between computational chemistry data formats using dpdata. Handles VASP, QE, CP2K, Gaussian, LAMMPS, and DeePMD formats. Essential for preparing ML potential training data.

Use when the user asks about density of states (DOS), projected DOS (PDOS), d-band center, spin-resolved DOS, or electronic structure analysis from completed DFT calculations.

Generate and manage DFTB+ calculations. Use when the user requests DFTB+, tight-binding DFT, SCC-DFTB, or needs fast approximate DFT for large systems or MD.

Run DeePMD-kit inference to predict energies, forces, and stresses using a trained DP model. Also covers model evaluation and testing.

CP2K single point energy calculation. Uses GTH pseudopotentials, Gaussian-plane-wave (GPW) method, and DZVP basis sets for periodic DFT.

CP2K geometry optimization. Handles bulk, slab, and molecular systems with GPW method. Efficient for large systems (200+ atoms).

Use when the user asks to test ENCUT convergence, k-point convergence, or any parameter sweep to determine converged computational settings.

Use when the user asks about CO2 reduction reaction (CO2RR), CO2 electroreduction intermediates, Faradaic efficiency, or selectivity toward CO, methanol, methane, formic acid, etc.

Use when the user asks about Bader charge analysis, charge transfer, oxidation states from DFT, or electron density partitioning.

Use when the user asks for adsorption energy, binding energy, or wants to compare how strongly a molecule binds to a surface.

Generate and manage ABINIT DFT calculations. Use when the user requests ABINIT, or needs DFPT phonons, GW calculations, or BSE optical spectra.

Recover from stale RFDB nodes that survive `commit_batch` cleanup despite having their `file` field listed in `changed_files`. Use when: (1) re-running an enricher/orchestrator step but old data with the same file field stays in the graph, (2) `commit_batch(changed_files=[X], ...)` reports success but a query later still returns nodes with `file = X`, (3) you re-pointed structural nodes (e.g., DIRECTORY/FILE) from a synthetic file path to real file paths and the old synthetic-path nodes won't...

Fix silent failures when parsing Grafema semantic IDs that are in URI format instead of legacy arrow format. Use when: (1) code splits semantic IDs by "->" but gets the whole string back because IDs are grafema:// URIs, (2) file path extraction from semantic IDs returns empty or wrong values, (3) derived edges (DEPENDS_ON, etc.) produce 0 results despite source edges existing, (4) any code that processes semantic IDs after the analysis pipeline's to_uri_format() has run. The grafema:// URI fo...