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.

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

Complex Linear Response Functions for a Multiconfigurational Self-Consistent Field Wave Function in a High Performance Computing Environment

Efficient parallel linear scaling method to get the response density matrix in all-electron real-space density-functional perturbation theory

Self-consistent Quantum Linear Response with a Polarizable Embedding environment

Consistent third-order one-particle transition and excited-state properties within the algebraic-diagrammatic construction scheme for the polarization propagator

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

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

An efficient quantum algorithm for generation of ab initio n-th order susceptibilities for non-linear spectroscopies

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

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

responsefun: Fun with Response Functions in the Algebraic Diagrammatic Construction Framework

Quantum Computation Of Molecular Response Properties

Efficient Spin-Adapted Implementation of Multireference Algebraic Diagrammatic Construction Theory. I. Core-Ionized States and X-Ray Photoelectron Spectra

Third-order algebraic diagrammatic construction theory for electron attachment and ionization energies: Conventional and Green's function implementation

Implementation of extended Lagrangian dynamics in GROMACS for polarizable simulations using the classical Drude oscillator model

Hybrid MPI and OpenMP Parallel Implementation of Large-Scale Linear-Response Time-Dependent Density Functional Theory with Plane-Wave Basis Set

Accurate Spectral Properties within Double-Hybrid Density Functional Theory: A Spin-Scaled Range-Separated Second-Order Algebraic-Diagrammatic Construction-Based Approach

Excitations and spectra from equilibrium real-time Green's functions

Efficient construction of involutory linear combinations of anti-commuting Pauli generators for large-scale iterative qubit coupled cluster calculations

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

Massively Parallel Implementation of Iterative Eigensolvers in Large-Scale Plane-Wave Density Functional Theory