Magnetic circular dichroism (MCD) spectroscopy is a technique based on the measurement of the differential absorption of left and right circular polarized light passing through a sample immersed in a (relatively strong) magnetic field parallel to the direction of propagation of the light. The frequency of the light is typically in a region of electronic absorption of the sample, that is UV/vis, near IR or X-ray, but magnetic circular dichroism of vibrational transitions can also be recorded. A recently proposed experimental protocol utilizes pulsed vector beams, where the magnetic field generated by the vector beams replaces the external static field applied in conventional MCD experiments. MCD has been applied with success to molecular systems to gain insight into structural, electronic, and magnetic properties, and experiments are moving into the investigation of nanosystems as well. Traditionally, the MCD signal is analyzed in terms of three characteristic band shapes, whose underlying electronic-structure intensity descriptors are known as the Faraday A, B, and C terms. The C term relates to the Zeeman splitting of degenerate ground states due to the external field, and its absorption-shaped band is temperature dependent. Only paramagnetic compounds exhibit C term signals. A and B terms have different physical origins and are both independent of temperature. While the A term arises from the magnetic field-induced Zeeman splitting of degenerate electronic states, the B term is the result of mixing of non-degenerate states under the influence of the magnetic field. Thus, all molecules exhibit B-type signals while only molecules with a three-fold or higher symmetry axis will have nonzero A terms. A derivative line shape indicates the presence of an A term and, thus, excited-state degeneracy, while a transition with only a B term contribution gives rise to a Gaussian-shaped absorption feature. When present, the A term dominates the spectrum. MCD experiments can help to deconvolute bands observed in a regular absorption spectrum. Since the MCD signals taking both positive and negative values, transitions can be better distinguished and, as mentioned, (near-)degenerate states can be identified from derivative lineshapes. In the interpretation of experimental MCD spectra, computational spectroscopy can provide valuable support. From quantum chemical calculation, it is possible to gain insight into the underlying electronic transitions, which in a spectrum are convoluted in a band of varying resolution. Among the most recent developments of theoretical approaches targeting MCD are implementations based on response theory within the algebraic diagrammatic construction (ADC) scheme, coupled cluster singles and doubles (CCSD), and resolution-of-identity coupled cluster singles and approximate doubles (RI-CC2) , as well as the linear-response (LR-) time-dependent (TD-) density functional theory (DFT) approach, in which excited states are calculated in a finite magnetic field. The calculation of the individual Faraday terms requires an explicit solution of the final excited states. This approach is limited by the convergence of the excited state eigenvectors and their derivatives which becomes increasingly difficult in regions with a high density of states. Damped response theory and the complex polarization propagator framework offer an alternative approach where the spectrum is computed directly as an interpolation of crosssection grid points in any frequency region. Implementations of this type of approach have been presented for various electronic-structure methods for the computation of linear spectroscopies such as one-photon absorption and electronic circular dichroism (ECD), of C 6 dispersion coefficients, and of the nonlinear spectroscopies such as MCD, 9,24 nuclear spin-induced circular dichroism (NSCD), magneto-chiral dichroism (MChD) and magnetochiral birefringence dispersion, as well as resonant inelastic X-ray scattering (RIXS) and two-photon absorption. In this work we extend the RI-CC2 framework for MCD spectroscopy 10 in TURBO-MOLE to include the calculation of the A term. Moreover, we address the shortcomings of the explicit calculation of the Faraday terms in regions where many close-lying electronic states can be expected by presenting an RI-CC2 implementation of damped quadratic response. The RI-CC2 method allows for applications to molecular systems of sizes that are not amenable for CCSD and higher-order methods, while still producing reasonably accurate excitation energies and, thus, spectra.