md-vib-spectra-simulation

Compute vibrational spectra from MD trajectories using MLatom autocorrelation methods.

52|11|Updated Mar 10, 2026
One-click install
npx skills add https://github.com/InternScience/ChemClaw --skill md-vib-spectra-simulation
Or copy as Structured Prompt for Agent
Please help me install this Agent Skill.
Skill: md-vib-spectra-simulation
Source: https://github.com/InternScience/ChemClaw/tree/main/skills/md-vib-spectra-simulation
Command: npx skills add https://github.com/InternScience/ChemClaw --skill md-vib-spectra-simulation

SYSTEM DOCUMENTATION & REQUIREMENTS

💡 This Skill requires mlatom, numpy, matplotlib, h5py, pyh5md, rdkit, xtb, and includes scripts (resource) components.

What problem does it solve?

Compute vibrational spectra from MD trajectories using MLatom to obtain power spectra (velocity autocorrelation) and IR spectra (dipole autocorrelation).

Core Features & Use Cases

  • Load MD data from h5md, or from XYZ/VXYZ, or generate MD trajectories from SMILES using GFN2-xTB.
  • Compute power spectrum, IR spectrum, or both; outputs include PNG plots and numpy arrays.
  • Flexible workflow support for users needing on-demand MD or post-processing of existing trajectories.

Quick Start

Run the script with your trajectory file (e.g., ethanol_traj.h5) to generate the desired spectrum.

Frequently Asked Questions about md-vib-spectra-simulation

High-intent search queries and answers about installing and using this skill.

FAQPage Schema
How do I compute a vibrational power spectrum from an MD trajectory?

You compute a vibrational power spectrum from an MD trajectory by processing velocity data through velocity autocorrelation. This skill uses MLatom to calculate the autocorrelation and outputs the resulting spectrum as PNG plots and numpy arrays.

Can I generate an IR spectrum directly from a SMILES string?

Yes, you can generate an IR spectrum directly from a SMILES string. The skill uses SMILES inputs to run on-demand MD simulations with GFN2-xTB, then calculates the IR spectrum via dipole autocorrelation using MLatom.

What input formats work for calculating vibrational spectra from MD data?

Supported input formats for calculating vibrational spectra include h5md trajectory files and XYZ/VXYZ coordinate files. You can also input SMILES strings to generate MD trajectories on-demand if you do not have existing data.

Do I need rdkit and xtb to calculate IR spectra from existing h5md files?

No, you do not need rdkit and xtb if you are calculating IR spectra from existing h5md files. These dependencies are only required when generating MD trajectories on-demand from SMILES inputs.

What is the difference between a power spectrum and an IR spectrum in molecular dynamics?

A power spectrum derives from velocity autocorrelation, while an IR spectrum derives from dipole autocorrelation. This skill uses MLatom to compute either or both from your MD trajectories, outputting the results as NPY arrays and PNG visualizations.