We report a variational calculation of ground state energies and radii of 4 HeN droplets (3 ≤ N ≤ 40), using the atom-atom interaction HFD-B(HE). The trial wave function has a simple structure, combining two-and three-body correlation functions coming from a translationally invariant configuration-interaction description, and Jastrow-type short-range correlations. The calculated ground state energies differ by around 2% from the diffusion Monte Carlo results. The research on liquid helium clusters has attracted a great interest both experimentally and theoretically [1,2]. This research allows for the analysis of the evolution of various physical properties for increasing size of the system, going from single atoms to the bulk. Helium clusters are expected to remain liquid under all conditions of formation, offering thus the possibility to study finite-size effects in the superfluid state [3]. Moreover, it has been suggested that Bose condensation could be detected by means of helium atom-cluster collisions [4]. The experimental research has faced the difficulties of size detecting the clusters. Recently, molecular beam diffraction from a transmission gratting [5] has proven to be successful to detect even the 4 He dimer [6], giving further impetus to the study of helium clusters.As the atom-atom interaction is well-known and relatively simple, the solution of the Schrödinger equation has been obtained using several microscopic methods, mainly based on Monte Carlo techniques [7,8,9,10,11]. Variational Monte Carlo (VMC) calculations, using a Jastrow-like ansatz for the many-body wave function, are currently used for systems dominated by strong shortrange interactions. The VMC wave functions are the input for diffusion Monte Carlo (DMC), Green function Monte Carlo or path integral techniques which provide essentially exact, within statistical errors, ground state energies of 4 He clusters at zero temperature. Conversely, these calculations constitute a useful benchmark to test other many-body methods.In this work we present a new method of obtaining high-quality variational wave functions to describe small bosonic clusters. The basic idea is to write the trial wave function as the product of three terms, each of them with well defined roles. The first term is the familiar two-body Jastrow correlation factor, which controls the strong atom-atom repulsion at very short distances. The second term is related to a single-particle description of the cluster and is written as the product of N single particle wave functions referred to the center-of-mass coordinate; it provides the required confinement of the constituents and fixes basically the size of the droplet. Finally, the third term corresponds to a special version of the configuration interaction (CI) expansion describing two-and three-particle excitations, its role being to incorporate fine details to the wave function for medium and long ranges, as well as some collective effects. In summary, the trial wave functions we shall consider to describe the ground state...