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.

VeloxChem: GPU-accelerated Fock matrix construction enabling complex polarization propagator simulations of circular dichroism spectra of G-quadruplexes

Efficient and Parallel Implementation of Real and Complex Response Functions Employing the Second-Order Algebraic-Diagrammatic Construction Scheme for the Polarization Propagator

Accelerating Coupled-Cluster Calculations with GPUs: An Implementation of the Density-Fitted CCSD(T) Approach for Heterogeneous Computing Architectures Using OpenMP Directives

GPUQT: an Efficient Linear-Scaling Quantum Transport Code Fully Implemented on Graphics Processing Units.

Quantum Chemistry for Solvated Molecules on Graphical Processing Units (GPUs)using Polarizable Continuum Models

Efficient Linear-Scaling Quantum Transport Calculations on Graphics Processing Units and Applications on Electron Transport in Graphene

Computing linear optical spectra in the presence of nonadiabatic effects on Graphics Processing Units using molecular dynamics and tensor-network approaches

Overcoming the bottleneck for quantum computations of complex nanophotonic structures: Purcell and FRET calculations using a rigorous mode hybridization method

Distributed Memory, GPU Accelerated Fock Construction for Hybrid, Gaussian Basis Density Functional Theory

Polariton Spectra under the Collective Coupling Regime. I. Efficient Simulation of Linear Spectra and Quantum Dynamics

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

Advanced Techniques for High-Performance Fock Matrix Construction on GPU Clusters

GFCCLib: Scalable and Efficient Coupled-Cluster Green's Function Library for Accurately Tackling Many Body Electronic Structure Problems

KeldyshQFT: A C++ codebase for real-frequency multiloop functional renormalization group and parquet computations of the single-impurity Anderson model

A black-box, general purpose quadratic self-consistent field code with and without Cholesky Decomposition of the two-electron integrals

Formulation and Implementation of Frequency-Dependent Linear Response Properties with Relativistic Coupled Cluster Theory for GPU-accelerated Computer Architectures

A Framework for Lattice QCD Calculations on GPUs

ELEQTRONeX: A GPU-Accelerated Exascale Framework for Non-Equilibrium Quantum Transport in Nanomaterials

Benchmark Force Fields for the Molecular Dynamic Simulation of G-Quadruplexes

An alternative GPU acceleration for a pseudopotential plane-waves density functional theory code with applications to metallic systems