Code for simulating nonlinear optical spectroscopies of closed and open systems
Updated on August 6, 2024 to version 0.2.0. Significant structural chagnes to the code base. Added composite diagrams. Unified UF2 and RKE code where possible, unified UF2 and RKE user interface where possible. Added helper classes for simulating common spectra. Higher-order signals are now easier for the user to calculate. Old code structure available in version 0.1.3.
pip install ufss
(Note: ufss is only written for python 3, so you may need to run pip3 install ufss if pip points to python 2 on your machine)
Visualizing Feynman diagrams requires a TeX distribution
To try ufss without installing or downloading the repository,
follow this link to see examples using the Diagram Generator:
https://colab.research.google.com/github/peterarose/ufss/blob/master/examples/DiagramGeneratorExample_Colab.ipynb
(Note: Google's Coloaboratory gives a warning message about running Jupyter notebooks not authored by Google. When prompted by the warning, select "RUN ANYWAY", and then click "YES" when it asks you if you would like to reset all runtimes)
Contains 4 separate modules:
- Diagram Generator (DG) - tool for automatically generating all Feynman diagrams for arbitrary order nonlinear optical spectroscopies
- class DiagramGenerator in ufss/diagram_automation.py
- Ultrafast Ultrafast (UF2) - fast algorithm for calculating individual Feynman diagrams including arbitrary pulse shapes
- for closed systems: class Wavepackets in ufss/UF2_core.py
- for open systems: class DensityMatrices in ufss/UF2_open_core.py
- Runga-Kutta-Euler (RKE) - alternative algorithm for calculating individual Feynman diagrams including arbitrary pulse shapes
- for closed systems: to be included in this repository eventually. For a working version see https://github.com/peterarose/ultrafastultrafast
- for open systems: class RKE_DensityMatrices in ufss/RKE_open_core.py
- Hamiltonian/Liouvillian Generator (HLG) - tool for generating vibronic Hamiltonians and (optionally) Liouvillians from simple parameter inputs
- contained in vibronic_eigenstates sub-directory
See README file in the examples subdirectory for more information. Most of the examples are jupyter notebooks.
This code depends upon the following packages:
numpy, matplotlib, scipy>=1, pyyaml, pyx
Note: as of version 0.2.2, pyfftw is an optional package. If you would like to use it, you must install fftw first, and then use pip to install pyfftw.
If you use UFSS, we request that you cite the following articles, which describe how its components work.
[1] Peter A. Rose and Jacob J. Krich, "Automatic Feynman diagram generation for nonlinear optical spectroscopies and application to fifth-order spectroscopy with pulse overlaps", J. Chem. Phys. 154, 034109 (2021)
[2] Peter A. Rose and Jacob J. Krich, "Efficient numerical method for predicting nonlinear optical spectroscopies of open systems", J. Chem. Phys. 154, 034108 (2021)