Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2311.11314 · Journal of the Optical Society of America B

Simulation of matrix product states to unveil the initial state dependency of non-Gaussian dynamics of Kerr nonlinearity

Souvik Agasti

Published Apr 4, 2024• 6 citations
arXivHTMLPDF

Abstract

We simulate a free dissipative and coherent-driven Kerr nonlinear system using a time-evolving block decimation (TEBD) algorithm to study the impact of the initial state on the exact quantum dynamics of the system. The superposition of two coherent branches results in non-classical time dynamics. The Wigner state representation confirms that the system ends up saturating to two different branches, through evolving different trajectories, resulting in de-Gaussification throughout evolution. Furthermore, we also see that the time evolution suffers the residual effect of the initial state.

DOI
10.1364/josab.519373

Related papers

Simulation of kerr nonlinearity: revealing initial state dependency

Souvik Agasti

Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using time-evolving block decimation (TEBD) algorithm and time propagation on the Heisenberg equation of motion using Euler’s method to study how the numerical results are analogous to classical bistability . The system evolves through different trajectories to stabilize different branches for different external drives and initial conditions. The Wigner state reprentation confirms the system to suffer a residual effect of initial state throughout the non-classical dynamical evolution and the metastable states of the system . Furthermore, we also see the numerically simulated spectral density remains significantly different from analytical counterparts when initial states do not lie to the same branch of the final state.

arXiv:2308.10667•Physica Scripta•Sep 25, 2023•8 citations

Numerical simulation of Kerr nonlinear systems; analyzing non-classical dynamics

Souvik Agasti

Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using Euler’s method by solving Heisenberg equation of motion and time evolving block decimation (TEBD) algorithm, and demonstrate how the numerical results are analogous to classical bistability. The comparison with analytics show that the TEBD numerics follow the quantum mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation . Comparing between two different numerical techniques, we see that the semi-classical Euler’s method gives the dynamics of the system field of one among two coherent branches, whereas TEBD numerics generate the superposition of both of them. Therefore, the time dynamics determined by TEBD numerical method undergoes through a non-classical state which is also shown by determining second order correlation function.

Journal of Physics Communications•Sep 20, 2019•10 citations

Dynamical simulations of classical stochastic systems using matrix product states

T. H. Johnson, Stephen R. L. Clark, Dieter Jaksch

We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of one-dimensional quantum systems, to simulate the time evolution of nonequilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time evolution is efficiently simulated by repeatedly updating and approximately refactorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations, we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size.

arXiv:1006.2639•Physical Review E•Sep 16, 2010•44 citations

Coherent-State Propagation: A Computational Framework for Simulating Bosonic Quantum Systems

Nikita Guseynov, Zoë Holmes, Armando Angrisani

We introduce coherent-state propagation, a computational framework for simulating bosonic systems. We focus on bosonic circuits composed of displaced linear optics augmented by Kerr nonlinearities, a universal model of bosonic quantum computation that is also physically motivated by driven Bose-Hubbard dynamics. The method works in the Schrödinger picture representing the evolving state as a sparse superposition of coherent states. We develop approximation strategies that keep the simulation cost tractable in physically relevant regimes, notably when the number of Kerr gates is small or the Kerr nonlinearities are weak, and prove rigorous guarantees for both observable estimation and sampling. In particular, bosonic circuits with logarithmically many Kerr gates admit quasi-polynomial-time classical simulation at exponentially small error in trace distance. We further identify a weak-nonlinearity regime in which the runtime is polynomial for arbitrarily small constant precision. We complement these results with numerical benchmarks on the Bose-Hubbard model with all-to-all connectivity. The method reproduces Fock-basis and matrix-product-state reference data, suggesting that it offers a useful route to the classical simulation of bosonic systems.

quant-pharXiv:2604.19625•Apr 22, 2026•0 citations

Simulation of Matrix Product States For Dissipation and Thermalization Dynamics of Open Quantum Systems

Souvik Agasti

Abstract We transform the system/reservoir coupling model into a one-dimensional semi-infinite discrete chain through unitary transformation to simulate the open quantum system numerically with the help of time evolving block decimation (TEBD) algorithm. We apply the method to study the dynamics of dissipative systems. We also generate the thermal state of a multimode bath using minimally entangled typical thermal state (METTS) algorithm, and investigate the impact of the thermal bath on an empty system. For both cases, we give an extensive analysis of the impact of the modeling and simulation parameters, and compare the numerics with the analytics.

Journal of Physics Communications•Dec 12, 2019•8 citations

Nonlinear effects on the dynamics of quantum harmonic modes coupled through angular momentum

N. Canosa, R. Rossignoli, Javier Gonzalo Garcia +2 more

Abstract We investigate nonlinear effects on the dynamics of entanglement and other quantum observables in a system of two harmonic modes coupled through angular momentum. The nonlinearity arises from a Kerr-type anharmonic term in each mode. The emergence and evolution of entanglement, non-Gaussianity, photon number, photon antibunching and squeezing are examined for different initial coherent product states and couplings, through exact diagonalization in a truncated basis. It is shown that the anharmonic terms, even if weak, can lead to very significant effects for such initial states, considerably enhancing and stabilizing entanglement and leading to a non negligible non-Gaussianity of the evolved states. They also affect other observables, stabilizing the dynamics after an initial transient regime, for not too small initial average populations of each mode. Analytic short-time approximate expressions are also provided.

arXiv:2004.11945•Journal of Physics B Atomic Molecular and Optical Physics•Aug 28, 2020•6 citations
We simulate a free dissipative and coherent-driven Kerr nonlinear system using a time-evolving block decimation (TEBD) algorithm to study the impact of the initial state on the exact quantum dynamics of the system. The superposition of two coherent branches results in non-classical time dynamics. The Wigner state representation confirms that the system ends up saturating to two different branches, through evolving different trajectories, resulting in de-Gaussification throughout evolution. Furthermore, we also see that the time evolution suffers the residual effect of the initial state.
Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using time-evolving block decimation (TEBD) algorithm and time propagation on the Heisenberg equation of motion using Euler’s method to study how the numerical results are analogous to classical bistability . The system evolves through different trajectories to stabilize different branches for different external drives and initial conditions. The Wigner state reprentation confirms the system to suffer a residual effect of initial state throughout the non-classical dynamical evolution and the metastable states of the system . Furthermore, we also see the numerically simulated spectral density remains significantly different from analytical counterparts when initial states do not lie to the same branch of the final state.
Abstract We simulate coherent driven free dissipative Kerr nonlinear system numerically using Euler’s method by solving Heisenberg equation of motion and time evolving block decimation (TEBD) algorithm, and demonstrate how the numerical results are analogous to classical bistability. The comparison with analytics show that the TEBD numerics follow the quantum mechanical exact solution obtained by mapping the equation of motion of the density matrix of the system to a Fokker-Plank equation . Comparing between two different numerical techniques, we see that the semi-classical Euler’s method gives the dynamics of the system field of one among two coherent branches, whereas TEBD numerics generate the superposition of both of them. Therefore, the time dynamics determined by TEBD numerical method undergoes through a non-classical state which is also shown by determining second order correlation function.
We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of one-dimensional quantum systems, to simulate the time evolution of nonequilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time evolution is efficiently simulated by repeatedly updating and approximately refactorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations, we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size.
We introduce coherent-state propagation, a computational framework for simulating bosonic systems. We focus on bosonic circuits composed of displaced linear optics augmented by Kerr nonlinearities, a universal model of bosonic quantum computation that is also physically motivated by driven Bose-Hubbard dynamics. The method works in the Schrödinger picture representing the evolving state as a sparse superposition of coherent states. We develop approximation strategies that keep the simulation cost tractable in physically relevant regimes, notably when the number of Kerr gates is small or the Kerr nonlinearities are weak, and prove rigorous guarantees for both observable estimation and sampling. In particular, bosonic circuits with logarithmically many Kerr gates admit quasi-polynomial-time classical simulation at exponentially small error in trace distance. We further identify a weak-nonlinearity regime in which the runtime is polynomial for arbitrarily small constant precision. We complement these results with numerical benchmarks on the Bose-Hubbard model with all-to-all connectivity. The method reproduces Fock-basis and matrix-product-state reference data, suggesting that it offers a useful route to the classical simulation of bosonic systems.
Abstract We transform the system/reservoir coupling model into a one-dimensional semi-infinite discrete chain through unitary transformation to simulate the open quantum system numerically with the help of time evolving block decimation (TEBD) algorithm. We apply the method to study the dynamics of dissipative systems. We also generate the thermal state of a multimode bath using minimally entangled typical thermal state (METTS) algorithm, and investigate the impact of the thermal bath on an empty system. For both cases, we give an extensive analysis of the impact of the modeling and simulation parameters, and compare the numerics with the analytics.
Abstract We investigate nonlinear effects on the dynamics of entanglement and other quantum observables in a system of two harmonic modes coupled through angular momentum. The nonlinearity arises from a Kerr-type anharmonic term in each mode. The emergence and evolution of entanglement, non-Gaussianity, photon number, photon antibunching and squeezing are examined for different initial coherent product states and couplings, through exact diagonalization in a truncated basis. It is shown that the anharmonic terms, even if weak, can lead to very significant effects for such initial states, considerably enhancing and stabilizing entanglement and leading to a non negligible non-Gaussianity of the evolved states. They also affect other observables, stabilizing the dynamics after an initial transient regime, for not too small initial average populations of each mode. Analytic short-time approximate expressions are also provided.