Tutorials
Below (and in the examples
folder on gitlab) you can find serveral tutorials and examples that demonstate how to use dynasor.
In setting up these examples, suitable parameters were chosen for the MD simulations that provide physically sound (albeit not necessarily highly converged) results. In practice, it is strongly recommended to carefully test the effect of e.g., system size, sampling time, spacing of snapshots, possibly interatomic potential/force field, and dynasor averaging parameters.
Getting started
To begin with, there are two tutorials about dynamics in aluminum. Here you will learn the basics, like setting up a trajectory, setting up the q-points to sample (along a path in the solid case and spherically averaged in the liquid case), running a dynasor calculation, understanding the contents of a dynasor-sample, and some basic visualization of the results.
Advanced applications
These tutorials cover some more advanced topics, such as splitting structure factors by species, some convergence considerations, and spectral energy density calculations.
In this tutorial you will learn how to handle multi-component systems by computing full and partial dynamic structure factors and current correlation for ordered and disordered NiAl.
In this tutorial you will learn how to compute the partial static structure factors for the different components of a halide perovskite and study the convergence of the results with system size and trajectory length.
Here, you will learn how to compute the spectral energy density from a molecular dynamics simulation of a crystal with dynasor, which captures the full anharmonicity, and compare this to a phonopy dispersion, which is based on 0 K lattice dynamics.
Post-processing
In these tutorials you will learn about ways to post-process the dynasor results, for example weighting of structure factors to compare with different scattering experiments, or extracting phonon frequencies and lifetimes using the damped harmonic oscillator model.
Mode projection
It can be useful to analyze the dynamics of the atoms in terms of normal modes.
The normal modes behaves roughly as independent harmonic oscillators for most materials.
It is thus possible to reason about these degrees of freedom as each having a single frequency and lifetime.
Even for very anharmonic materials normal modes can serve as a intuitive basis for further analysis.
In the mode projection tutorial the basic functionality of dynasor.ModeProjector
is demonstrated.
Other topics
This tutorial provides some general insights into auto-correlation functions, Fourier transforms, and numerical tools for decreasing noise.
This example serves to demonstrate the typical convergence of a correlation function with total simulation time.