  We examine the non-Markovian dynamics in a multipartite system of two
initially correlated atomic qubits, each located in a single-mode leaky cavity
and interacting with its own bosonic reservoir. We show the dominance of
non-Markovian features, as quantified by the difference in fidelity of the
evolved system with its density matrix at an earlier time, in three specific
two-qubit partitions associated with the cavity-cavity and atom-reservoir
density matrices within the same subsystem, and the cavity-reservoir reduced
matrix across the two subsystems. The non-Markovianity in the cavity-cavity
subsystem is seen to be optimized in the vicinity of the exceptional point. The
CHSH-Bell inequality computed for various two-qubit partitions show that high
non-locality present in a specific subsystem appears in conjunction with
enhanced non-Markovian dynamics in adjacent subsystems. This is in contrast to
the matching existence of non-locality and quantum correlations in regions
spanned by time t and the cavity decay rate for select partitions. We discuss
the applicability of these results to photosynthetic systems.
