chem-ts-optimization

A chemistry workflow for improving a guessed transition-state structure—the arrangement of atoms at a reaction's energy barrier—and checking that it is a first-order saddle point.

In plain words
What is it for?
Optimise non-periodic molecular transition-state guesses with Sella, then use vibrational calculations to check for the required imaginary mode.
Why use it?
It helps distinguish a valid reaction transition state from a structure that only looks plausible.

Skill for Claude CodeCodex

Install

Getting it into your agent

One page per mod, every tool's command on it. A separate URL per tool would split the same page into five that compete with each other.

agentmods
npx agentmods add skills/learningmatter-mit/atomisticskills/chem-ts-optimization
Any agent
npx skills add learningmatter-mit/AtomisticSkills --skill chem-ts-optimization
Clone the repo
git clone --depth 1 https://github.com/learningmatter-mit/AtomisticSkills

Made for: Claude Code, Codex.

Per session 23 Skills are progressive disclosure: only the name and description are preloaded; the body loads when the skill is used.
When invoked 917 The whole file, excluding the scripts and references it only reads on demand.
Security scan A 0 findings. Scan, not verified.
Origin original No closer match found in the catalogue.
Token cost

What it costs to keep this loaded

Counted locally with the o200k_base tokenizer, which is exact for GPT models; Claude uses its own tokenizer and its counts differ. Treat this as one consistent yardstick across the catalogue rather than a bill. Prices are per million input tokens.

ModelPer sessionOnce invoked
Fable 5 $0.00023 $0.00917
Opus 5 $0.00012 $0.00458
Sonnet 5 $0.00005 $0.00183
Haiku 4.5 $0.00002 $0.00092

Measured 2d ago against content hash 88904f151386, method: parsed. Prices are Anthropic first-party input rates as of 2026-08-30, from the pricing page.

Security

Grade A, and why

chem-ts-optimization scanned grade A with 0 findings against 26 rules in 11 categories — prompt injection, anti-refusal, data exfiltration, privilege escalation, supply chain, agent snooping, system-prompt leakage, SSRF and excessive agency — measured 2d ago.

The scan reads SKILL.md. This mod also ships 2 executable files (examples/acetonitrile/run_example.sh, scripts/optimize_ts_sella.py), listed below but not scanned — reading those needs a real analyzer, not pattern matching.

A static scan of the body, not an audit. Every finding is printed with the line that produced it so you can judge whether it matters here. A mod is markdown that instructs an agent; that is exactly why what it instructs is worth reading.

Nothing flagged

None of the 26 patterns this scan looks for appear in this file: no shell pipes, no recursive deletes, no credential paths, no hidden text, no instruction-override or anti-refusal phrasing, no agent-config snooping. That is not a guarantee, it is the absence of the things that are checkable.

.agents/skills/chem-ts-optimization/SKILL.md · 108 lines

How it starts

The opening of the file, as written. The whole thing — 108 lines — stays where its author put it; the contents beside it link to each section on GitHub.

TS Optimization with Sella

Optimize a transition-state guess and check whether it is a first-order saddle point.

Scope

  • Domain: molecular chemistry only (non-periodic systems).
  • Trigger: user has a TS guess and needs TS optimization plus frequency validation.
  • Exclusions: periodic diffusion/path workflows (use chem-neb-barrier instead).

Tool

optimize_ts_sella.py

Runs Sella TS optimization followed by finite-difference vibrations.

Use with MACE

# Env: mace-agent
python .agents/skills/chem-ts-optimization/scripts/optimize_ts_sella.py \
  --ts_guess ts_guess.xyz \
  --model_type mace \
  --model_name MACE-OFF23-small \
  --fmax 0.02 \
  --steps 500 \
  --imag_cutoff_cm1 -50.0 \
  --output_dir results/ts_opt

Use with FAIRChem (UMA)

# Env: fairchem-agent
python .agents/skills/chem-ts-optimization/scripts/optimize_ts_sella.py \
  --ts_guess ts_guess.xyz \
  --model_type fairchem \
  --model_name uma-s-1p1 \
  --task_name omol \
  --fmax 0.02 \
  --steps 500 \
  --imag_cutoff_cm1 -50.0 \
  --output_dir results/ts_opt

Arguments

  • --ts_guess: required TS guess geometry (XYZ supported by ASE I/O).
  • --model_type: required backend (mace or fairchem).
  • --model_name: optional model identifier/checkpoint.
  • --task_name: optional model head/task (for UMA molecular runs use omol).
  • --device: auto|cpu|cuda (default auto).
  • --fmax: Sella convergence threshold in eV/A (default 0.02).
  • --steps: maximum TS optimization steps (default 500).
  • --vib_delta: finite-difference displacement in A (default 0.01).
  • --vib_nfree: finite-difference stencil size (2 or 4, default 2).
  • --imag_cutoff_cm1: imaginary mode cutoff in cm^-1 (default -50.0).
  • --keep_vib_cache: optional flag to keep vibration cache files in output_dir/vib.
  • --output_dir: required output directory.

Outputs

  • ts_optimized.xyz: optimized TS geometry.
  • ts_opt.traj: TS optimization trajectory.
  • ts_opt.log: optimizer log.
  • ts_optimization_results.json: run summary and pass/fail decision.
  • vib/ cache files only when --keep_vib_cache is set.

Read the full file on GitHub · 108 lines

Changes

What this file has done since we first saw it

Hashed on every crawl. A supply-chain change to an agent config is a question of when, not whether, so the history is kept rather than the latest state alone.

  1. 2d ago First seen · 108 lines · 23 tokens per session scan A 88904f151386

Subscribe to this mod's changes

chem-ts-optimization is a skill published in the GitHub repository learningmatter-mit/AtomisticSkills (154 stars, last pushed 7d ago), licensed MIT. It adds 23 tokens to every session and 917 once invoked, about $0.0001 per session on Opus 5. A static security scan graded it A with 0 findings. No closer match exists in the catalogue, so it is treated as the original; first seen 2026-08-30.

Related

Other skills, from other repositories

datamol

Pythonic wrapper around RDKit with simplified interface and sensible defaults. Preferred for standard drug discovery including SMILES parsing, standardization, descriptors, fingerprints, clustering, 3D conformers, parallel processing. Returns native rdkit.Chem.Mol objects. For advanced control or custom parameters…

synthetic-sciences/openscience · 67 tokens

deepchem

Molecular ML with diverse featurizers and pre-built datasets. Use for property prediction (ADMET, toxicity) with traditional ML or GNNs when you want extensive featurization options and MoleculeNet benchmarks. Best for quick experiments with pre-trained models, diverse molecular representations. For graph-first…

synthetic-sciences/openscience · 78 tokens

hypogenic

Automated LLM-driven hypothesis generation and testing on tabular datasets. Use when you want to systematically explore hypotheses about patterns in empirical data (e.g., deception detection, content analysis). Combines literature insights with data-driven hypothesis testing. For manual hypothesis formulation use…

synthetic-sciences/openscience · 69 tokens

binding-affinity

Hybrid ML + physics binding affinity prediction. Empirical scoring, MM/GBSA rescoring, multi-method consensus, and batch virtual screening for protein-ligand complexes.

synthetic-sciences/openscience · 38 tokens

drug-design

End-to-end drug discovery pipeline orchestration. Deterministic Python script that auto-chains structure prediction, pocket detection, de novo design, docking, scoring, and ADMET filtering into reproducible workflows.

synthetic-sciences/openscience · 44 tokens

medchem

Medicinal chemistry filters. Apply drug-likeness rules (Lipinski, Veber), PAINS filters, structural alerts, complexity metrics, for compound prioritization and library filtering.

synthetic-sciences/openscience · 39 tokens