Recently, the result obtained from `Stardust' mission suggests that the overall ratio of compact to aggregate particles is 65:35 (or 13:7) for Comet 81P/Wild 2 (Burchell et al. 2008). In the present work, we propose a model which considers cometary dust as a mixture of compact and aggregate particles, with composition of silicate and organic. We consider compact particles as spheroidal particles and aggregates as BCCA and BAM2 aggregate with some size distribution. For modeling Comet 1P/ Halley, the power-law size distribution n(a)= a^{-2.6}, for both compact and aggregate particles is taken. We take a mixture of BAM2 and BCCA aggregates with a lower and upper cutoff size around 0.20$\mu m$ and 1$\mu m$. We also take a mixture of prolate, spherical and oblate compact particles with axial ratio (E) from 0.8 to 1.2 where a lower and upper cutoff size around 0.1$\mu m$ and 10$\mu m$ are taken. Using T-matrix code, the average simulated polarization curves are generated which can best fit the observed polarization data at the four wavelengths $\lambda$ = 0.365$\mu m$, 0.485$\mu m$, 0.670$\mu m$ and 0.684$\mu m$. The suitable mixing percentage of aggregates emerging out from the present modeling corresponds to 50% BAM2 and 50% BCCA particles and silicate to organic mixing percentage corresponds to 78% silicate and 22% organic in terms of volume. The present model successfully reproduces the observed polarization data, especially the negative branch, more effectively as compared to other work done in the past. It is found that among the aggregates, the BAM2 aggregate plays a major role, in deciding the cross-over angle and depth of negative polarization branch.

In our present study, the observed linear polarization data of comet Hyakutake are studied at wavelengths $\lambda=0.365\mu m$, $\lambda=0.485\mu m$ and 0.684$\mu m$ through simulations using Ballistic Particle-Cluster Aggregate and Ballistic Cluster-Cluster Aggregate aggregates of 128 spherical monomers. We first investigated that the size parameter of the monomer, $x$ $\sim$ 1.56 -- 1.70, turned out to be most suitable which provides the best fits to the observed dust scattering properties at three wavelengths $\lambda = 0.365$$\mu m$, 0.485$\mu m$ and 0.684$\mu m$. Thus the effective radius of the aggregate (r) lies in the range $0.45 \mu m \le r \le 0.49 \mu m$ at $\lambda = 0.365$$\mu m$; $ 0.60 \mu m \le r \le 0.66 \mu m$ at $\lambda = 0.485$$\mu m$ and $0.88 \mu m \le r \le 0.94 \mu m$ at $\lambda = 0.684$$\mu m$. Now using superposition \textsc{t-matrix} code and the power-law size distribution, $n(r) \sim r^{-3}$, the best-fitting values of complex refractive indices are calculated which can best fit the observed polarization data at the above three wavelengths. The best-fitting complex refractive indices $(n,k)$ are found to be (1.745, 0.095) at $\lambda = 0.365 $ $\mu m$, (1.743, 0.100) at $\lambda = 0.485 $ $\mu m$ and (1.695, 0.100) at $\lambda = 0.684$ $\mu m$. The refractive indices coming out from the present analysis correspond to mixture of both silicates and organics, which are in good agreement with the \textit{in situ} measurement of comets by different spacecraft.