“…As in ref , we decouple the valence and core-excited states, , yielding the following matrix representation of the molecular Hamiltonian H normalm normalo normall = ( lefttrue boldH normalv 0 0 boldH normalc ) where H v and H c are sub-Hamiltonians acting on the manifolds of the valence and core-excited electronic states, respectively. In both subspaces, a vibronic coupling model up to second order ,, is employed, leading to coupled potential energy surfaces (PESs) in a diabatic representation. , In this framework, each submatrix in eq is expressed as H x = H false( 0 false) + W x , boldx ∈ { v , c } where the zeroth-order Hamiltonian is constituted by the ground-state Hamiltonian in the harmonic approximation H false( 0 false) = prefix∑ i ω normali 2 ( − ∂ 2 ∂ Q i 2 + Q i 2 ) bold1 with ω i representing the frequency of mode Q i .…”