Brings GPU-accelerated molecular dynamics and N-body simulations to Claude through 15 tools that handle particle system creation, Lennard-Jones potentials, NVT/NVE ensemble runs, and trajectory analysis. You can spin up a 100k particle MD simulation, compute mean squared displacement for diffusion coefficients, or model galaxy collisions with gravitational dynamics. Automatically detects CUDA and falls back to CPU when needed. Part of a larger scicomp suite that includes quantum mechanics and neural network servers, with cross-server data sharing via URIs. Supports async tasks for long-running simulations so you're not blocking on a 10k timestep calculation. Install with uvx or pip, configure in Claude Desktop's JSON, and you've got production molecular dynamics without leaving your chat.
Public tool metadata for what this MCP can expose to an agent.
addAdds two numbers together2 paramsAdds two numbers together
firstNumbernumbersecondNumbernumbersubtractSubtracts the second number from the first number2 paramsSubtracts the second number from the first number
minuendnumbersubtrahendnumbermultiplyMultiplies two numbers together2 paramsMultiplies two numbers together
firstNumbernumberSecondNumbernumberdivisionDivides the first number by the second number2 paramsDivides the first number by the second number
numeratornumberdenominatornumbersumAdds any number of numbers together1 paramsAdds any number of numbers together
numbersarraymoduloDivides two numbers and returns the remainder2 paramsDivides two numbers and returns the remainder
numeratornumberdenominatornumbermeanCalculates the arithmetic mean of a list of numbers1 paramsCalculates the arithmetic mean of a list of numbers
numbersarraymedianCalculates the median of a list of numbers1 paramsCalculates the median of a list of numbers
numbersarraymodeFinds the most common number in a list of numbers1 paramsFinds the most common number in a list of numbers
numbersarrayminFinds the minimum value from a list of numbers1 paramsFinds the minimum value from a list of numbers
numbersarraymaxFinds the maximum value from a list of numbers1 paramsFinds the maximum value from a list of numbers
numbersarrayfloorRounds a number down to the nearest integer1 paramsRounds a number down to the nearest integer
numbernumberceilingRounds a number up to the nearest integer1 paramsRounds a number up to the nearest integer
numbernumberroundRounds a number to the nearest integer1 paramsRounds a number to the nearest integer
numbernumbersinCalculates the sine of a number in radians1 paramsCalculates the sine of a number in radians
numbernumberarcsinCalculates the arcsine of a number in radians1 paramsCalculates the arcsine of a number in radians
numbernumbercosCalculates the cosine of a number in radians1 paramsCalculates the cosine of a number in radians
numbernumberarccosCalculates the arccosine of a number in radians1 paramsCalculates the arccosine of a number in radians
numbernumbertanCalculates the tangent of a number in radians1 paramsCalculates the tangent of a number in radians
numbernumberarctanCalculates the arctangent of a number in radians1 paramsCalculates the arctangent of a number in radians
numbernumberradiansToDegreesConverts a radian value to its equivalent in degrees1 paramsConverts a radian value to its equivalent in degrees
numbernumberdegreesToRadiansConverts a degree value to its equivalent in radians1 paramsConverts a degree value to its equivalent in radians
numbernumberGPU-accelerated Model Context Protocol servers for computational mathematics, physics simulations, and machine learning.
| Guide | Description |
|---|---|
| Installation | Setup instructions for pip, uv, and uvx |
| Configuration | Claude Desktop & Claude Code setup |
| Quick Start | Get running in 5 minutes |
| API Reference | Complete tool documentation |
| Visual Demos | Interactive physics simulations |
This system enables AI assistants to perform real scientific computing — from solving differential equations to running molecular dynamics simulations.
Quantum Wave Mechanics Double-slit interference pattern from solving the time-dependent Schrödinger equation |
N-Body Dynamics Galaxy merger simulation using gravitational N-body calculations |
Crystal Diffraction Bragg scattering from a hexagonal (graphene-like) lattice |
Multi-Slit Interference Complex interference patterns from three coherent sources |
This system provides 4 specialized MCP servers that bring scientific computing capabilities to AI assistants like Claude:
| Server | Description | Tools |
|---|---|---|
| Math MCP | Symbolic algebra (SymPy) + numerical computing | 14 |
| Quantum MCP | Wave mechanics & Schrodinger simulations | 12 |
| Molecular MCP | Classical molecular dynamics | 15 |
| Neural MCP | Neural network training & evaluation | 16 |
Key Features:
Run any MCP server directly without installation:
# Run individual servers
uvx scicomp-math-mcp
uvx scicomp-quantum-mcp
uvx scicomp-molecular-mcp
uvx scicomp-neural-mcp
# Install individual servers
pip install scicomp-math-mcp
pip install scicomp-quantum-mcp
pip install scicomp-molecular-mcp
pip install scicomp-neural-mcp
# Or install all at once
pip install scicomp-math-mcp scicomp-quantum-mcp scicomp-molecular-mcp scicomp-neural-mcp
# With GPU support (requires CUDA)
pip install scicomp-math-mcp[gpu] scicomp-quantum-mcp[gpu] scicomp-molecular-mcp[gpu] scicomp-neural-mcp[gpu]
Add to your Claude Desktop configuration file:
macOS: ~/Library/Application Support/Claude/claude_desktop_config.json
Windows: %APPDATA%\Claude\claude_desktop_config.json
{
"mcpServers": {
"math-mcp": {
"command": "uvx",
"args": ["scicomp-math-mcp"]
},
"quantum-mcp": {
"command": "uvx",
"args": ["scicomp-quantum-mcp"]
},
"molecular-mcp": {
"command": "uvx",
"args": ["scicomp-molecular-mcp"]
},
"neural-mcp": {
"command": "uvx",
"args": ["scicomp-neural-mcp"]
}
}
}
Add to your project's .mcp.json:
{
"mcpServers": {
"math-mcp": {
"command": "uvx",
"args": ["scicomp-math-mcp"]
},
"quantum-mcp": {
"command": "uvx",
"args": ["scicomp-quantum-mcp"]
}
}
}
Or configure globally in ~/.claude/settings.json.
# Solve equations symbolically
symbolic_solve(equations="x**3 - 6*x**2 + 11*x - 6")
# Result: [1, 2, 3]
# Compute derivatives
symbolic_diff(expression="sin(x)*exp(-x**2)", variable="x")
# Result: cos(x)*exp(-x**2) - 2*x*sin(x)*exp(-x**2)
# GPU-accelerated matrix operations
result = matrix_multiply(a=matrix_a, b=matrix_b, use_gpu=True)
# Create a Gaussian wave packet
psi = create_gaussian_wavepacket(
grid_size=[256],
position=[64],
momentum=[2.0],
width=5.0
)
# Solve time-dependent Schrodinger equation
simulation = solve_schrodinger(
potential=barrier_potential,
initial_state=psi,
time_steps=1000,
dt=0.1,
use_gpu=True
)
# Create particle system
system = create_particles(
n_particles=1000,
box_size=[20, 20, 20],
temperature=1.5
)
# Add Lennard-Jones potential
add_potential(system_id=system, potential_type="lennard_jones")
# Run MD simulation
trajectory = run_nvt(system_id=system, n_steps=100000, temperature=1.0)
# Analyze diffusion
msd = compute_msd(trajectory_id=trajectory)
# Define model
model = define_model(architecture="resnet18", num_classes=10, pretrained=True)
# Load dataset
dataset = load_dataset(dataset_name="CIFAR10", split="train")
# Train
experiment = train_model(
model_id=model,
dataset_id=dataset,
epochs=50,
batch_size=128,
use_gpu=True
)
# Export for deployment
export_model(model_id=model, format="onnx", output_path="model.onnx")
# Clone the repository
git clone https://github.com/andylbrummer/math-mcp.git
cd math-mcp
# Install dependencies
uv sync --all-extras
# Install MCP servers in editable mode (required for entry points)
uv pip install --python .venv/bin/python \
-e servers/math-mcp \
-e servers/quantum-mcp \
-e servers/molecular-mcp \
-e servers/neural-mcp
# Run tests
uv run pytest -m "not gpu" # CPU only
uv run pytest # All tests (requires CUDA)
# Run with coverage
uv run pytest --cov=shared --cov=servers
Note: The editable install step is required because
uv syncdoesn't install entry point scripts for workspace packages. After this step, you can run servers directly withuv run scicomp-math-mcp.
See CONTRIBUTING.md for development guidelines.
GPU acceleration provides significant speedups for compute-intensive operations:
| MCP | Operation | CPU | GPU | Speedup |
|---|---|---|---|---|
| Math | Matrix multiply (4096x4096) | 2.1s | 35ms | 60x |
| Quantum | 2D Schrodinger (512x512, 1000 steps) | 2h | 2min | 60x |
| Molecular | MD (100k particles, 10k steps) | 1h | 30s | 120x |
| Neural | ResNet18 training (1 epoch) | 45min | 30s | 90x |
For technical details about the system architecture, see ARCHITECTURE.md.
MIT License - see LICENSE for details.
Contributions are welcome! Please see CONTRIBUTING.md for guidelines.