Hamiltonian non-Hermicity: accurate dynamics with the multiple Davydov D$_2$ Ans\"atze

We examine the applicability of the numerically accurate method of time dependent variation with multiple Davydov Ansatze (mDA) to non-Hermitian systems. Three systems of interest includes: a non-Hermitian system of dissipative Landau-Zener transitions, a non-Hermitian, multimode Jaynes-Cummings model, and a dissipative Holstein-Tavis-Cummings model, where complex many-body dynamics are accurately captured by the mDA method. Our findings highlight the versatility of the mDA as a powerful numerical tool for investigating complex many-body non-Hermitian systems, which can be extended to explore diverse phenomena such as skin effects, excited-state dynamics, and spectral topology in the non-Hermitian field.