SKILL.md:
name: quantum-sim description: Simulate quantum circuits using a statevector engine. Supports up to 20 qubits with gates H, X, Y, Z, S, T, Rx, Ry, Rz, P, CX, CZ, SWAP. Built-in presets for Bell state, GHZ state, QFT, Grover search, and quantum teleportation. Accepts inline QASM-like syntax or a circuit file. Use when the user asks to simulate a quantum circuit, demonstrate quantum entanglement, run Grover search, apply quantum gates, show measurement probabilities, or explore quantum algorithms. Trigger phrases include "simulate quantum", "quantum circuit", "Bell state", "GHZ", "Grover", "QFT", "quantum gates", "qubit", "statevector", "superposition". metadata: {"openclaw":{"emoji":":atom_symbol:","requires":{"bins":["python3"]}}}
Quantum
Circuit Simulator
Statevector-based quantum circuit simulator in pure Python 3 stdlib. No Qiskit, no pip install. Optional: if numpy is installed, matrix ops run faster for larger qubit counts.
Simulates up to 20 qubits. Provides statevector output, measurement probability histograms, and shot-based measurement sampling.
Quick Start
python3 {baseDir}/scripts/quantum_sim.py --list-presets python3 {baseDir}/scripts/quantum_sim.py --preset bell python3 {baseDir}/scripts/quantum_sim.py --preset ghz --qubits 5 python3 {baseDir}/scripts/quantum_sim.py --preset grover python3 {baseDir}/scripts/quantum_sim.py --qasm "qubits 3; h 0; cx 0 1; cx 0 2" python3 {baseDir}/scripts/quantum_sim.py --preset bell --json
Built-in Presets
Preset | Qubits
| Description | ||||
|---|---|---|---|---|
| bell | 2 | Bell state ( | 00>+ | 11>)/sqrt(2) - max entanglement |
| ghz | N | GHZ state - N-qubit entanglement (use --qubits) | ||
| qft | N | Quantum Fourier Transform (use --qubits) | ||
| grover | 2 | Grover search, marks | 11>, 100% success rate | |
| teleportation | 3 | Full quantum teleportation protocol |
Supported Gates
Single-qubit: H, X, Y, Z, S, T, Sdg, Tdg, Rx(theta), Ry(theta), Rz(theta), P(lambda) Two-qubit: CX (CNOT), CZ, SWAP
QASM-like Syntax
Write one instruction per line (or semicolon-separated inline):
qubits 3 h 0 cx 0 1 cx 0 2 rx 1.5708 0 rz 3.1416 1 measure 2048
Rules:
- First line must be: qu
bits N
- Gate names are case-insensitive
- Single-qubit gates: gate_name qubit_index
- Rotation gates: rx/ry/rz/p theta qubit_index (theta in radians)
- Two-qubit gates: cx/cz/swap qubit1 qubit2
- Comments start with #
Save to a file and run: python3 {baseDir}/scripts/quantum_sim.py --qasm-file my_circuit.qasm
All Flags
| Flag | Effect |
|---|---|
| --preset NAME | Run a built-in preset circuit |
| --qubits N | Override qubit count for ghz/qft presets |
| --qasm "..." | Inline QASM (semicolons separate instructions) |
| --qasm-file PATH | Load circuit from a .qasm text file |
| --shots N | Measurement shots (default 1024) |
| --json | Output statevector + counts as JSON |
--list-presets | Show all presets with descriptions --statevector-only| Skip measurement simulation
Output Format
Each run prints:
- Circuit summary (qubit count, gate sequence)
- Statevector - complex amplitudes for all basis states with prob > 0.001
- Measurement histogram - shot counts with ASCII bar chart
Example output for Bell state: === Bell state |Phi+>: maximally entangled 2-qubit state === Qubits: 2 Dim: 4 Gates applied: 2 Circuit: H(0) -> CX(0,1)
Statevector: |00> +0.7071+0.0000j p=0.5000 [######### ] |11> +0.7071+0.0000j p=0.5000 [######### ]
Measurement (1024 shots): |11> [############ ] 532 ( 52.0%) |00> [########### ] 492 ( 48.0%)
JSON Output
Use --json for machine-readable output (safe to pipe, no ANSI): python3 {baseDir}/scripts/quantum_sim.py --preset bell --json | python3 -c \ "import json,sys; d=json.load(sys.stdin); print(d['probabilities'])"
JSON keys: label, n_qubits, gates, statevector (re/im per state), probabilities, counts, shots
Physics Notes
- Qubit ordering: qubit 0 is least significant bit. |01> means q1=0, q0=1.
- Statevector size: 2^n complex numbers. Memory: 16 * 2^n bytes. 20 qubits = 16MB.
- Grover preset: 1 iteration on 2 qubits achieves 100% success for marked state |11>.
- QFT: includes bit-reversal swaps. Output is frequency-domain representation.
- Teleportation: simulates all 3 qubits unitarily (no mid-circuit measurement collapse).
- All rotation angles are in radians. pi/2 = 1.5708, pi
= 3.1416.