Abstract:Nonlinear Optical Spectroscopy is a well-developed field with theoretical and experimental advances that have benefited multiple disciplines, including chemistry, biology, and physics. However, for the accurate interpretation of the corresponding multi-dimensional spectra, there is a need for precise quantum dynamical simulations based on model Hamiltonians. In this article, we present the initial release of our code, QuDPy (Quantum Dynamics in Python), which provides a robust numerical platform for performing quantum dynamics simulations based on model systems, including open quantum systems. A distinguishing feature of our approach is the ability to specify various high-order optical response pathways in the form of double-sided Feynman diagrams through a straightforward input syntax. This syntax outlines the time-ordering of ket-sided or bra-sided optical interactions acting on the time-evolving density matrix of the system. We utilize the quantum dynamics capabilities of QuTip to simulate the spectral response of complex systems, allowing us to compute virtually any n -th order optical response of the model system. To illustrate the utility of our approach, we provide a series of example calculations. Program summary Program Title: QuDPy CPC Library link to program files: https://doi.org/10.17632/5xm9pm24cz.1 Developer's repository link: https://github.com/sa-shah/QuDPy Licensing provisions: MIT License Programming language: Python (v3.7) Supplementary material: Available as Google Colab Files. Example 1: https://tinyurl.com/y3j5jmmr Example 2: https://tinyurl.com/37vwntn5 External packages: • QuTip (v.4.7) and dependencies i.e. Numpy, Matplotlib ( https://qutip.org/ ) • UFSS Automatic Diagram Generator ( https://github.com/peterarose/ufss ) Nature of problem: Accurate quantum simulations of complex systems are required in order to understand and interpret multi-dimensional ultrafast spectroscopic signals. This code provides an open-source/multi-platform method that facilitates the generation of higher-order non-linear optical responses for an arbitrary molecular or material system given a model input Hamiltonian and bath model. Solution method: We use the double-sided Feynman diagram method [1, 2] to generate (symbolically) a set of response functions corresponding to the nth order non-linear response of the system to a series of laser pulses using the UFSS package [3] We then perform a series of accurate quantum dynamics calculations using the QuTip package [4] to generate the numerical response and spectra which correspond to specific experimental conditions. References [1] S. Mukamel, Principles of Nonlinear Optics and Spectroscopy, Oxford University Press, 1995. [2] P. Hamm, M. Zanni, Concepts and Methods of 2D Infrared Spectroscopy, Cambridge University Press, 2011. [3] Peter A. Rose, Jacob J. Krich, J. Chem. Phys. 154 (2021) 034109. [4] J. Johansson, P. Nation, F. Nori, Comput. Phys. Commun. 184 (4) (2013) 1234–1240.

FeynmanPAQS: a Graphical Interface Program for Photonic Analog Quantum Computing

Piquasso: A Photonic Quantum Computer Simulation Software Platform

A manufacturable platform for photonic quantum computing

GraphiQ: Quantum circuit design for photonic graph states

QuDPy: A Python-based tool for computing ultrafast non-linear optical responses

QuFlow: Visualizing Parameter Flow in Quantum Circuits for Understanding Quantum Computation

A Dynamically Programmable Quantum Photonic Microprocessor for Graph Computation

Strawberry Fields: A Software Platform for Photonic Quantum Computing

Quantum circuits with many photons on a programmable nanophotonic chip

Applications of Near-Term Photonic Quantum Computers: Software and Algorithms

SimuQ: A Framework for Programming Quantum Hamiltonian Simulation with Analog Compilation

A versatile single-photon-based quantum computing platform

Very-large-scale integrated quantum graph photonics

Realization of fractional quantum Hall state with interacting photons

Digital Discovery of 100 diverse Quantum Experiments with PyTheus

PauliForest: Connectivity-Aware Synthesis and Pauli-Oriented Qubit Mapping for Near Term Quantum Simulation

Reconfiguration of Quantum Photonic Integrated Circuits Using Auxiliary Fields

Q2Chemistry: A Quantum Computation Platform for Quantum Chemistry

Photonic quantum simulation of ground state configurations of Heisenberg square and checkerboard lattice spin systems