Atomic and molecular clusters have been a topic of great interest for last few decades, mainly because of their unusual properties, tunability and vast technological applications. In this thesis, we have explored the optical properties of few clusters using first principles calculations.We have performed systematic large-scale all-electron correlated calculations on boron B n , aluminum Al n and magnesium Mg n clusters (n=2-5), to study their linear optical absorption spectra. Several possible isomers of each cluster were considered, and their geometries were optimized at the coupled-cluster singles doubles (CCSD) level of theory. Using the optimized ground-state geometries, excited states of different clusters were computed using the multi-reference singles-doubles configurationinteraction (MRSDCI) approach, which includes electron correlation effects at a sophisticated level. These CI wavefunctions were used to compute the transition dipole matrix elements connecting the ground and various excited states of different clusters, eventually leading to their linear absorption spectra. The convergence of our results with respect to the basis sets, and the size of the CI expansion was carefully examined.Isomers of a given cluster show a distinct signature spectrum, indicating a strongstructure property relationship. This fact can be used in experiments to distinguish between different isomers of a cluster. Owing to the sophistication of our calculations, our results can be used for benchmarking of the absorption spectra and be used to design superior time-dependent density functional theoretical (TDDFT) approaches.The contribution of configurations to many-body wavefunction of various excited states suggests that in most cases optical excitations involved are collective, and plasmonic in nature.In addition, we have calculated the optical absorption in various isomers of neutral boron B 6 and cationic boron B + 6 clusters using computationally less expensive configuration interaction singles (CIS) approach, and benchmarked these results against more sophisticated equation-of-motion (EOM) CCSD based approach. In all closed shell systems, a complete agreement on the nature of configurations involved is observed in both methods. On the other hand, for open-shell systems, minor contribution from double excitations are observed, which are not captured in the CIS method.Optical absorption in planar boron clusters in wheel shape, B 7 , B 8 and B 9 computed using EOM-CCSD approach, have been compared to the results obtained from TDDFT approach with a number of functionals. This benchmarking reveals that range-separated functionals such as ωB97xD and CAM-B3LYP give qualitatively as well as quantitatively the same results as that of EOM-CCSD.