In this paper, we propose a replica-based method called Vendi Sampling where multiple copies of a molecular system are simulated in parallel. To enhance the conformational space exploration, the replicas are coupled via a global statistical measure, the Vendi score . The Vendi score reflects the instantaneous and internal diversity of the replicas and is a function of the eigenvalues of a Gram matrix computed via a pre-specified kernel. The Gram matrix is constructed by evaluating the kernel between all pairs of simulation replicas. As such, the proposed Vendi sampling method does not rely on modulating thermodynamic variables, Hamiltonians, or defining collective variables to enhance sampling. Instead, an extended ensemble is defined, which can drive sampling in infinite dimensional feature spaces, without explicit calculation of collective variables.