In this work, we examine nonadiabatic population dynamics for 1 1 B 1 and 1 1 A 2 states of ozone molecule (O 3 ). In O 3 , two lowest singlet excited states, 1 A 2 and 1 B 1 , can be coupled. Thus, population transfer between them occurs through the seam involving these two states. At any point of the seam (conical intersection), the Born-Oppenheimer approximation breaks down, and it is necessary to investigate nonadiabatic dynamics. We consider a linear vibronic coupling Hamiltonian model and evaluate vibronic coupling constant, diabatic frequencies for three modes of O 3 , bilinear and quadratic coupling constants for diabatic potentials, displacements, and Huang-Rhys coupling constants using ab initio calculations. The electronic structure calculations have been performed at the multireference configuration interaction and complete active space with second-order perturbation theory with a full-valence complete active space self-consistent field methods and augmented Dunning's standard correlation-consistent-polarized quadruple zeta basis set to determine ab initio potential energy surfaces for the ground state and first two excited states of O 3 , respectively. We have chosen active space comprising 18 electrons distributed over 12 active orbitals. Our calculations predict the linear vibronic coupling constant 0.123 eV. We have obtained the population on the 1 1 B 1 and 1 1 A 2 excited electronic states for the first 500 fs after photoexcitation.