Cloud-based quantum chemistry platform with Python API. Preferred for computational chemistry workflows including pKa prediction, geometry optimization, conformer searching, molecular property calculations, protein-ligand docking (AutoDock Vina), and AI protein cofolding (Chai-1, Boltz-1/2). Use when tasks involve quantum chemistry calculations, molecular property prediction, DFT or semiempirical methods, neural network potentials (AIMNet2), protein-ligand binding predictions, or automated computational chemistry pipelines. Provides cloud compute resources with no local setup required.
npx skills add https://github.com/foryourhealth111-pixel/Vibe-Skills --skill rowan
Use this skill only for explicit Rowan context, including Rowan, rowan-python, labs.rowansci.com, the Rowan API, or a Rowan-specific chemistry workflow. Do not use it for generic chemistry, RDKit, PubChem, ChEMBL, docking, pKa, conformer search, quantum chemistry, molecular machine learning, Boltz, or Chai-1 unless explicit Rowan context is present.
Rowan is a cloud-based computational chemistry platform that provides programmatic access to quantum chemistry workflows through a Python API. It enables automation of complex molecular simulations without requiring local computational resources or expertise in multiple quantum chemistry packages.
Key Capabilities:
Why Rowan:
uv pip install rowan-python
Generate an API key at labs.rowansci.com/account/api-keys.
Option 1: Direct assignment
import rowan
rowan.api_key = "your_api_key_here"
Option 2: Environment variable (recommended)
export ROWAN_API_KEY="your_api_key_here"
The API key is automatically read from ROWAN_API_KEY on module import.
import rowan
# Check authentication
user = rowan.whoami()
print(f"Logged in as: {user.username}")
print(f"Credits available: {user.credits}")
Calculate the acid dissociation constant for molecules:
import rowan
import stjames
# Create molecule from SMILES
mol = stjames.Molecule.from_smiles("c1ccccc1O") # Phenol
# Submit pKa workflow
workflow = rowan.submit_pka_workflow(
initial_molecule=mol,
name="phenol pKa calculation"
)
# Wait for completion
workflow.wait_for_result()
workflow.fetch_latest(in_place=True)
# Access results
print(f"Strongest acid pKa: {workflow.data['strongest_acid']}") # ~10.17
Generate and optimize molecular conformers:
import rowan
import stjames
mol = stjames.Molecule.from_smiles("CCCC") # Butane
workflow = rowan.submit_conformer_search_workflow(
initial_molecule=mol,
name="butane conformer search"
)
workflow.wait_for_result()
workflow.fetch_latest(in_place=True)
# Access conformer ensemble
conformers = workflow.data['conformers']
for i, conf in enumerate(conformers):
print(f"Conformer {i}: Energy = {conf['energy']:.4f} Hartree")
Optimize molecular geometry to minimum energy structure:
import rowan
import stjames
mol = stjames.Molecule.from_smiles("CC(=O)O") # Acetic acid
workflow = rowan.submit_basic_calculation_workflow(
initial_molecule=mol,
name="acetic acid optimization",
workflow_type="optimization"
)
workflow.wait_for_result()
workflow.fetch_latest(in_place=True)
# Get optimized structure
optimized_mol = workflow.data['final_molecule']
print(f"Final energy: {optimized_mol.energy} Hartree")
Dock small molecules to protein targets:
import rowan
# First, upload or create protein
protein = rowan.create_protein_from_pdb_id(
name="EGFR kinase",
code="1M17"
)
# Define binding pocket (from crystal structure or manual)
pocket = {
"center": [10.0, 20.0, 30.0],
"size": [20.0, 20.0, 20.0]
}
# Submit docking
workflow = rowan.submit_docking_workflow(
protein=protein.uuid,
pocket=pocket,
initial_molecule=stjames.Molecule.from_smiles("Cc1ccc(NC(=O)c2ccc(CN3CCN(C)CC3)cc2)cc1"),
name="EGFR docking"
)
workflow.wait_for_result()
workflow.fetch_latest(in_place=True)
# Access docking results
docking_score = workflow.data['docking_score']
print(f"Docking score: {docking_score}")
Predict protein-ligand complex structures using AI models:
import rowan
# Protein sequence
protein_seq = "MENFQKVEKIGEGTYGVVYKARNKLTGEVVALKKIRLDTETEGVPSTAIREISLLKELNHPNIVKLLDVIHTENKLYLVFEFLHQDLKKFMDASALTGIPLPLIKSYLFQLLQGLAFCHSHRVLHRDLKPQNLLINTEGAIKLADFGLARAFGVPVRTYTHEVVTLWYRAPEILLGCKYYSTAVDIWSLGCIFAEMVTRRALFPGDSEIDQLFRIFRTLGTPDEVVWPGVTSMPDYKPSFPKWARQDFSKVVPPLDEDGRSLLSQMLHYDPNKRISAKAALAHPFFQDVTKPVPHLRL"
# Ligand SMILES
ligand = "CCC(C)CN=C1NCC2(CCCOC2)CN1"
# Submit cofolding with Chai-1
workflow = rowan.submit_protein_cofolding_workflow(
initial_protein_sequences=[protein_seq],
initial_smiles_list=[ligand],
name="kinase-ligand cofolding",
model="chai_1r" # or "boltz_1x", "boltz_2"
)
workflow.wait_for_result()
workflow.fetch_latest(in_place=True)
# Access structure predictions
print(f"Predicted TM Score: {workflow.data['ptm_score']}")
print(f"Interface pTM: {workflow.data['interface_ptm']}")
For users working with RDKit molecules, Rowan provides a simplified interface:
import rowan
from rdkit import Chem
# Create RDKit molecule
mol = Chem.MolFromSmiles("c1ccccc1O")
# Compute pKa directly
pka_result = rowan.run_pka(mol)
print(f"pKa: {pka_result.strongest_acid}")
# Batch processing
mols = [Chem.MolFromSmiles(smi) for smi in ["CCO", "CC(=O)O", "c1ccccc1O"]]
results = rowan.batch_pka(mols)
for mol, result in zip(mols, results):
print(f"{Chem.MolToSmiles(mol)}: pKa = {result.strongest_acid}")
Available RDKit-native functions:
run_pka, batch_pka - pKa calculationsrun_tautomers, batch_tautomers - Tautomer enumerationrun_conformers, batch_conformers - Conformer generationrun_energy, batch_energy - Single-point energiesrun_optimization, batch_optimization - Geometry optimizationSee references/rdkit_native.md for complete documentation.
# List recent workflows
workflows = rowan.list_workflows(size=10)
for wf in workflows:
print(f"{wf.name}: {wf.status}")
# Filter by status
pending = rowan.list_workflows(status="running")
# Retrieve specific workflow
workflow = rowan.retrieve_workflow("workflow-uuid")
# Submit multiple workflows
workflows = rowan.batch_submit_workflow(
molecules=[mol1, mol2, mol3],
workflow_type="pka",
workflow_data={}
)
# Poll status of multiple workflows
statuses = rowan.batch_poll_status([wf.uuid for wf in workflows])
# Create folder for project
folder = rowan.create_folder(name="Drug Discovery Project")
# Submit workflow to folder
workflow = rowan.submit_pka_workflow(
initial_molecule=mol,
name="compound pKa",
folder_uuid=folder.uuid
)
# List workflows in folder
folder_workflows = rowan.list_workflows(folder_uuid=folder.uuid)
Rowan supports multiple levels of theory:
Neural Network Potentials:
Semiempirical:
DFT:
Methods are automatically selected based on workflow type, or can be specified explicitly in workflow parameters.
For detailed API documentation, consult these reference files:
references/api_reference.md: Complete API documentation - Workflow class, submission functions, retrieval methodsreferences/workflow_types.md: All 30+ workflow types with parameters - pKa, docking, cofolding, etc.references/rdkit_native.md: RDKit-native API functions for seamless cheminformatics integrationreferences/molecule_handling.md: stjames.Molecule class - creating molecules from SMILES, XYZ, RDKitreferences/proteins_and_organization.md: Protein upload, folder management, project organizationreferences/results_interpretation.md: Understanding workflow outputs, confidence scores, validationimport rowan
import stjames
smiles_list = ["CCO", "c1ccccc1O", "CC(=O)O"]
# Submit all pKa calculations
workflows = []
for smi in smiles_list:
mol = stjames.Molecule.from_smiles(smi)
wf = rowan.submit_pka_workflow(
initial_molecule=mol,
name=f"pKa: {smi}"
)
workflows.append(wf)
# Wait for all to complete
for wf in workflows:
wf.wait_for_result()
wf.fetch_latest(in_place=True)
print(f"{wf.name}: pKa = {wf.data['strongest_acid']}")
import rowan
# Upload protein once
protein = rowan.upload_protein("target.pdb", name="Drug Target")
protein.sanitize() # Clean structure
# Define pocket
pocket = {"center": [x, y, z], "size": [20, 20, 20]}
# Screen compound library
for smiles in compound_library:
mol = stjames.Molecule.from_smiles(smiles)
workflow = rowan.submit_docking_workflow(
protein=protein.uuid,
pocket=pocket,
initial_molecule=mol,
name=f"Dock: {smiles[:20]}"
)
import rowan
import stjames
mol = stjames.Molecule.from_smiles("complex_molecule_smiles")
# Generate conformers
conf_wf = rowan.submit_conformer_search_workflow(
initial_molecule=mol,
name="conformer search"
)
conf_wf.wait_for_result()
conf_wf.fetch_latest(in_place=True)
# Analyze lowest energy conformers
conformers = sorted(conf_wf.data['conformers'], key=lambda x: x['energy'])
print(f"Found {len(conformers)} unique conformers")
print(f"Energy range: {conformers[0]['energy']:.4f} to {conformers[-1]['energy']:.4f} Hartree")
rowan.whoami().credits to check balanceimport rowan
try:
workflow = rowan.submit_pka_workflow(
initial_molecule=mol,
name="calculation"
)
workflow.wait_for_result(timeout=3600) # 1 hour timeout
if workflow.status == "completed":
workflow.fetch_latest(in_place=True)
print(workflow.data)
elif workflow.status == "failed":
print(f"Workflow failed: {workflow.error_message}")
except rowan.RowanAPIError as e:
print(f"API error: {e}")
except TimeoutError:
print("Workflow timed out")
Convert PyTorch AT_DISPATCH macros to AT_DISPATCH_V2 format in ATen C++ code. Use when porting AT_DISPATCH_ALL_TYPES_AND*, AT_DISPATCH_FLOATING_TYPES*, or other dispatch macros to the new v2 API. For ATen kernel files, CUDA kernels, and native operator implementations.
Write docstrings for PyTorch functions and methods following PyTorch conventions. Use when writing or updating docstrings in PyTorch code.
Statistical models library for Python. Use when you need specific model classes (OLS, GLM, mixed models, ARIMA) with detailed diagnostics, residuals, and inference. Best for econometrics, time series, rigorous inference with coefficient tables. For guided statistical test selection with APA reporting use statistical-analysis.
Answer questions about the AI SDK and help build AI-powered features. Use when developers: (1) Ask about AI SDK functions like generateText, streamText, ToolLoopAgent, embed, or tools, (2) Want to build AI agents, chatbots, RAG systems, or text generation features, (3) Have questions about AI providers (OpenAI, Anthropic, Google, etc.), streaming, tool calling, structured output, or embeddings, (4) Use React hooks like useChat or useCompletion. Triggers on: "AI SDK", "Vercel AI SDK", "generateText", "streamText", "add AI to my app", "build an agent", "tool calling", "structured output", "useChat".
Create an llms.txt file from scratch based on repository structure following the llms.txt specification at https://llmstxt.org/
Use when working directly with the `esm` Python SDK, ESM3 or ESMC model IDs, Forge/Biohub inference clients, or ESMFold2 folding workflows.
Modal is a serverless cloud platform for running Python on demand, including on-demand GPUs. Use when deploying or serving AI/ML models, running GPU-accelerated workloads (training, fine-tuning, inference), serving web endpoints, scheduling batch jobs, or scaling Python code to cloud containers with the Modal SDK.
Use Therapeutics Data Commons through the PyTDC Python package for registry discovery, approved dataset access, task-aware splits, evaluator metrics, benchmark groups, and bounded molecular-oracle workflows.
Take foryourhealth111-pixel/rowan from the repository into ~/.claude/skills for personal
use, or into .claude/skills inside a project.
The agent identifies a skill by the name field in its header. Two skills with the
same name cannot sit side by side — one of them will be ignored.
The instructions reference pip, uv.
Without those the skill loads but fails at the first command.