We describe the new version 2.0.0 of the code DIRQFAM that calculates the multipole response of even-even axially symmetric deformed nuclei using the framework of relativistic self-consistent mean-field models. The response is calculated by implementing the finite amplitude method for relativistic quasiparticle random phase approximation. In the new version, we have implemented the following features: (i) meson-exchange interactions, (ii) low-rank method for calculating induced densities and currents, (iii) generalized minimal residual method (GMRES), (iv) evaluation of the nucleon localization function, (v) extraction of the QRPA transition matrix elements and eigenfrequencies, (vi) evaluation of the H11 part of the induced Hamiltonian, (vii) multipole excitation operators up to J≤5 are now included. Program summaryProgram title:DIRQFAM v2.0.0CPC Library link to program files:https://doi.org/10.17632/6gv8nxg4ns.2Licensing provisions: GPLv3Programming language: Fortran 90.External routines/libraries: BLAS/LAPACK, version 3.1. or higher.Journal reference of previous version: A. Bjelčić, T. Nikšić, Comp. Phys. Comm. 253, 107184 (2020).Does the new version supersede the previous version?: YesNature of problem: Multipole response of deformed even-even open-shell nuclei can be calculated using the quasiparticle finite amplitude method (QFAM), based on the relativistic self-consistent mean-field models. The particle-hole channel is described by a zero-range relativistic effective interaction or interaction based on meson exchange, while the particle-particle channel of the effective inter-nucleon interaction is described by a separable finite-range pairing force. The method can be applied to perform systematic studies of collective modes even in heavy deformed nuclei.Solution method: The current implementation computes the multipole response by solving the QFAM equations in a self-consistent iteration scheme. At each iteration the QFAM solutions are updated using the GMRES method. The QFAM amplitudes are expanded in the simplex-y harmonic oscillator basis.Summary of revisions:1.meson-exchange interaction are included in the code;2.calculation of the nucleon localization function;3.extraction of the QRPA transition matrix elements and eigenfrequencies;4.evaluation of the H11 part of the induced Hamiltonian;5.multipoles excitation operators up to J≤5 are included;6.low-rank method for calculating induced densities and currents;7.generalized minimal residual method;8.dynamic memory allocation;9.reduced memory requirements;Additional comments including restrictions and unusual features: Open-shell even-even nuclei with parity conserved axially symmetric ground states are considered. An electric multipole operator with J≤5 is used to calculate the response function.