“…However there exist other VQE approaches to directly target these correlation functions in either the (imaginary) time or frequency domains [580,601,605]. These include VQE-based variational approaches to directly solve the linear equations for the response of a system at a given frequency [606], which can be cast as a modified cost function, with similarly parameterized VQE ansatz [397,607,608]. While these approaches can describe the correlation functions to high accuracy over the whole energy range without restricting to a low-energy subspace, their challenge arises chiefly from the substantially more difficult optimization problem for the ansatz, originating from the larger condition number of the cost function, as well as the necessity for Hadamard tests to compute the transition amplitudes between the excited (or response) and ground states (for more details, see Appendix.…”