In our present study, the observed linear polarization data of comet Hyakutake are studied at wavelengths λ = 0.365µm, λ = 0.485µm and 0.684µ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 ∼ 1.56 -1.70, turned out to be most suitable which provides the best fits to the observed dust scattering properties at three wavelengths λ = 0.365µm, 0.485µm and 0.684µm. Thus the effective radius of the aggregate (r) lies in the range 0.45µm ≤ r ≤ 0.49µm at λ = 0.365µm; 0.60µm ≤ r ≤ 0.66µm at λ = 0.485µm and 0.88µm ≤ r ≤ 0.94µm at λ = 0.684µm. Now using superposition T-MATRIX code and the power-law size distribution, n(r) ∼ 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 λ = 0.365 µm, (1.743, 0.100) at λ = 0.485 µm and (1.695, 0.100) at λ = 0.684 µ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 in situ measurement of comets by different spacecraft.