In this thesis, I have investigated gas phase infrared spectroscopy of environmentally as well as astrophysical important large organic molecules such as naphthalene, methy-lated naphthalene, fluorine, methyalted fluorine etc. which are commonly known as polycyclic aromatic hydrocarbons (PAHs). Depending upon the molecular weight these organic molecules can exist both in gaseous as well as in the particulate state at room temperature hence they are the major environmental pollutants. They are also responsible for the unidentified infrared emission bands in the interstellar medium. Chapter 1 provides a brief introduction to my thesis work. A detailed literature survey on the importance, abundance of the PAHs in the environment as well as various spectroscopic techniques useful for identifying the PAHs has been done. Since the objective of my thesis work is to assign the observed fundamental infrared bands of large organic molecules with the help of high level quantum mechanical calculations, a brief introduction to the various high level quantum mechanical techniques that I have used in assigning the bands have been described in this chapter. In Chapter 2 I have presented the experimental and the theoretical methodologies in details. The chapter begins with a detailed description of the experimental procedure used for recording the infrared spectrum of these molecules followed by the theoretical methodologies used for the assignment of the observed infrared bands as well as for identifying the Fermi resonances. In Chapters 3 and 4, of this thesis I have recorded infrared spectrum of 1-and 2-methylnaphthalene (1-and 2-MN), fluorine (FL), 1-methylfluorene (1-MFL) and 1,8-dimethylfluorene (1,8-DMF) in the gas phase. The observed bands were assigned with the help of scaled harmonic frequency, scaled quantum mechanical harmonic force field (SQMFF) and enharmonic frequency calculations. The first two methods are based on the harmonic approximation, whereas the enharmonic frequency calculation is based on the standard second order perturbation theory. All these calculations gave me a partial fit to the fundamental bands in both aromatic and aliphatic C-H stretching as well as in the non C-H stretching region. At the end of both the chapters an error analysis in fitting the spectrum from all the three different calculations have been presented. Evidently the non linear least square fitting method employed in SQMFF calculation gives much better agreement between the experiment and theory than the other two methods. It has been observed in the experimental spectrum of methylated naphthalene that the band structure near the C-H stretch around 3000 cm−1 is very complicated and many bands and shoulders remain unassigned by the methods described in Chapters 3 and 4. Fermi resonance is one of the potential reason for the complicated band structure in this region. In Chapter 5, I have taken naphthalene and have investigated the Fermi resonance around the C-H stretching region using an effective vibrational hamiltonian (EVH) approach. In this method I have constructed an EVH consisting of 8 C-H stretches and 8 H-C-C in-plane bend overtones and 28 H-C-C in-plane bend combination modes as the basis. Both type 1 (stretch overtone) and type 2 (stretch combination) Fermi resonances were investigated. Calculated frequencies belonging to B1u and B2u irreducible representation were compared with the observed bands. Many bands and shoulders have been assigned as the overtone and combination modes of low frequency H-C-C bend motion obtained from the EVH approach. How-ever some bands remain unassigned in this method. This is perhaps due to the neglect of the carbon framework motion in the construction of the EVH. To improve upon the results obtained from the EVH formalism I included the carbon frame degrees of freedom and have carried out a full variation treatment in curvilinear coordinates. I have considered the 8 C-H stretches and 8 H-C-C in-plane bends of naphthalene as local mode oscillators and 17 coordinates belonging to the carbon framework motion as curvilinear normal mode oscillators. A quartic hamiltonian in a mixed local mode -normal mode basis was constructed including up to three body terms in both kinetic and potential energy part. The hamiltonian was subsequently recast into the ladder operator form and diagonal zed in a symmetry adapted basis with polyad constraints. Frequencies so obtained were compared to the experiment All these findings have been presented in Chapter 6 of this thesis. The concluding remark of the thesis and the future direction is presented in Chapter 7