In this thesis, we investigate the entanglement and magnetic properties of frustrated spin systems and correlated electronic properties of conjugated carbon systems. In chapter 1, we present different approaches to solve the time-independent, nonrelativistic Schr¨odinger equation for a many-body system. We start with the full non-relativistic Hamiltonian of a multi nuclear system to describe the Born - Oppenheimer approximation which allows the study of electronic Hamiltonian which treats nuclear positions parametrically. We then also describe ab initio techniques such as the Hartree-Fock Method and density functional theories. We then introduce model Hamiltonians for strongly correlated systems such as the Hubbard, Pariser-Parr-Pople and Heisenberg models, and show how they result from the noninteracting one-band tight-binding model. In chapter 2, we discuss various numerical techniques like the exact diagonalization methods and density matrix renormalization group (DMRG) method. We also discuss quantum entanglement and the success of DMRG which can be attributed to the area law of entanglement entropy. In chapter 3, we study here different regions in phase diagrams of the spin-1/2, spin-1 and spin-3/2 one-dimensional antiferromagnetic Heisenberg systems with nearest-neighbor (J1) and next-nearest-neighbor (J2) interactions and dimerization (d ). Frustration arises for specific relative signs of the interactions J1 and J2. In particular, we analyze the behavior of the bipartite entanglement entropy and fidelity at the gapless to gapped phase transitions and across the lines separating different phases in the J2−d plane. All the calculations in this work are based on exact diagonalizations of finite systems. In chapter 4, we study Heisenberg spin-1/2 and spin-1 chains with alternating ferromagnetic (JF 1 ) and antiferromagnetic (JA 1 ) nearest-neighbor interactions and a ferromagnetic next-nearest-neighbor interaction (JF 2 ). In this model frustration is present due to non-zero JF 2 . The model with site spin s behaves like a Haldane spin chain with site spin 2s in the limit of vanishing JF 2 and large JF 1 /JA 1 . We show that the exact ground state of the model can be found along a line in the parameter space. For fixed JF 1 , the phase diagram in the space of JA 1 −JF 2 is determined using numerical techniques complemented by analytical calculations. A number of quantities, including the structure factor, energy gap, entanglement entropy and zero temperature magnetization, are studied to understand the complete phase diagram. An interesting and potentially important feature of this model is that it can exhibit a macroscopic magnetization jump in the presence of a magnetic field; we study this using an effective Hamiltonian. In chapter 5, we study correlated electronic properties of zigzag and armchair fused naphthalenes and polyperylene systems in the presence of long-range electronelectron interactions. We find that the ground state of zigzag fused naphthalene system is a higher spin state, while the ground state of armchair fused naphthalene is a singlet. The spin gap of polyperylene is unusually small and the ground state is a singlet. Our calculations of optical gap and two-photon gap suggest that polyperylene should exhibit fluorescence. From the charge gap calculation, we predict that in zigzag fused naphthalene and polyperylene systems, excitons are weakly binding. Peierls type of distortion is negligible in zigzag fused naphthalene and polyperylene systems, however, in armchair fused naphthalene system, interior bonds have tendency to distort in low-lying excited states. In chapter 6, we study the ground state spin of the Heisenberg spin-1/2 nearestneighboring antiferromagnetic exchange models of systems with fused odd member rings. In particular, we compute the ground state spin of fused three and five membered rings as well as fused five membered rings. In the thermodynamic limit, the ground state of the fused three and five membered system is a higher spin state, while fused five membered system shows a singlet ground state, for all system sizes.

The present thesis deals with a theoretical study of electronic structures in -conjugated molecular materials with focus on their application in organic elec-tronics. We also discuss a modified and efficient symmetrized DMRG algorithm for studying excited states in these systems. In recent times, organic conjugated systems have emerged as potential candidates in a wide range of fascinating fields by virtue of their tunable electronic properties, easy processability and low cost. Tunability in the electronic and optical properties primarily are centered on the or-dering and nature of the low-lying excited states. Probing these important excited states also demands development of efficient and adaptable techniques. Chapter 1 provides a basic overview of conjugated organic polymers which have been utilized over decades in diverse fields as in organic light emitting diodes (OLED), organic solar cells (OSC) and non-linear optical (NLO) devices. These systems also contribute significantly to theoretical understanding as they pro vide important insights of one and quasi-one dimensional systems. In this chapter, we have given basic description of the electronic processes in OLED and OSC along with a brief theoretical description of -conjugated organic systems. Chapter 2 gives an account of the numerical techniques which are necessary for the study of low-dimensional strongly correlated systems like -conjugated sys-tems. For this purpose, effective low-energy model Hamiltonians viz. Huckel,¨ Hubbard and Pariser-Parr-Pople Hamiltonians are discussed. Exact diagonalization technique within the diagrammatic valence bond (DVB) basis and density matrix renormalization group (DMRG) technique are discussed in details. We have also given brief accounts of the methods employed to study real-time dynamics. A short description of different computational techniques for the study of NLO properties in -conjugated systems is also provided. Engineering the position of the lowest triplet state (T1) relative to the first excited singlet state (S1) is of great importance in improving the efficiencies of organic light emitting diodes and organic photovoltaic cells. In chapter 3, we have carried out model exact calculations of substituted polyene chains to understand the fac-tors that affect the energy gap between S1 and T1. The factors studied are backbone dimerization, different donor-acceptor substitutions and twisted backbone geome-try. The largest system studied is an eighteen carbon polyene which spans a Hilbert space of about 991 million in the triplet subspace. We show that for reverse inter-system crossing (RISC) process, the best choice involves substituting all carbon sites on one half of the polyene with donors and the other half with acceptors. Singlet fission (SF) is a potential pathway for significant enhancement of efficiency in OSC. In chapter 4, we study singlet fission in a pair of polyene molecules in two different stacking arrangements employing exact many-body wave packet dy-namics. In the non-interacting model, SF is absent. The individual molecules are treated within Hubbard and Pariser-Parr-Pople (PPP) models and the interac-tion between them involves transfer terms, intersite electron repulsions and site-charge—bond-charge repulsion terms. Initial wave packet is construc ted from ex-cited singlet state of one molecule and ground state of the other. Time develop-ment of this wave packet under the influence of intermolecular interactions is fol-lowed within the Schrodinger¨ picture by an efficient predictor-corrector scheme. In unsubstituted Hubbard and PPP chains, 21A state leads to significant SF yield while the 11B state gives negligible fission yield. On substitution by donor-acceptor groups of moderate strength, the lowest excited state will have sufficient 2 1A char-acter and hence gives significant SF yield. Because of rapid internal c onversion, the nature of the lowest excited singlet will determine the SF contribution to OSC effi - ciency. Furthermore, we find the fission yield depends considerably on th e stacking arrangement of the polyene molecules. In chapter 5, we have given an account of a new modified algorithm for symmetry adaptation within symmetrized density matrix renormalization group (SDMRG) technique. SDMRG technique has been an efficient method for studying low-lying eigenstates in one and quasi-one dimensional electronic systems. However, SDMRG method until now, had bottlenecks involving construction of linearly in-dependent symmetry adapted basis states as the symmetry matrices in the DMRG basis were not sparse. Our modified algorithm overcomes this bottleneck. T he new method incorporates end-to-end interchange symmetry (C2), electron-hole symmetry (J) and parity or spin-flip symmetry (P) in these calculations. The one-to-one correspondence between direct-product basis states in the DMRG Hilbert space for these symmetry operations renders the symmetry matrices in the new ba-sis with maximum sparseness, just one non-zero matrix element per row. Using methods similar to those employed in exact diagonalization technique for Pariser-Parr-Pople (PPP) models, developed in the eighties, it is possible to construct or-thogonal SDMRG basis states while bypassing the slow step of Gram-Schmidt orthonormalization procedure. The method together with the PPP model which incorporates long-range electronic correlations is employed to study the correlated excited states of 1,12-benzoperylene. In chapter 6, we have studied the correlated excited states of coronene and ova-lene within Pariser-Parr-Pople (PPP) model employing symmetry adapted density matrix renormalization group technique. These polynuclear aromatic hydrocar-bons can be considered as graphene nanoflakes and study of their ele ctronic struc-tures will shed light on the electron correlation effects in these finite-size gr aphene analogues. The electron correlation effect usually diminishes on going from one-dimensional to higher-dimensional systems, yet, it is significant within these fin ite-size graphene derivatives where it depends on the molecular topology. We have characterized these low-lying energy states by calculating bond orders, spin den-sities in the lowest triplet state and two-photon absorption cross-sections for low-lying two-photon states. vi