Rational design of disordered molecular aggregates and solids for optoelectronic applications relies on our ability to predict the properties of such materials using theoretical and computational methods. However, large molecular systems where disorder is too significant to be considered in the perturbative limit cannot be described using either first principles quantum chemistry or band theory. Multiscale modeling is a promising approach to understanding and optimizing the optoelectronic properties of such systems. It uses first-principles quantum chemical methods to calculate the properties of individual molecules, then constructs model Hamiltonians of molecular aggregates or bulk materials based on these calculations. In this paper, we present a protocol for constructing a tight-binding Hamiltonian that represents the excited states of a molecular material in the basis of Frenckel excitons: electron-hole pairs that are localized on individual molecules that make up the material. The Hamiltonian parametrization proposed here accounts for excitonic couplings between molecules, as well as for electrostatic polarization of the electron density on a molecule by the charge distribution on surrounding molecules. Such model Hamiltonians can be used to calculate optical absorption spectra and other optoelectronic properties of molecular aggregates and solids. Video Link The video component of this article can be found at https://www.jove.com/video/60598/ 13,14 and time-dependent density functional theory (TD-DFT) 15,16. Organic molecules with applications in optoelectronics typically have relatively large π-conjugated systems; many also have donor and acceptor groups. Capturing the correct charge-transfer behavior in such molecules is essential to calculating their optoelectronic properties, but it can only be accomplished using long-range corrected hybrid functionals in TD-DFT 17,18,19,20. Calculations that use such functionals scale super linearly with the size of the system and, at present, they are only practical for modeling the optoelectronic properties of individual organic molecules or small molecular aggregates that can be described using no more thañ 10 4 atomic basis functions. A simulation method that could describe disordered materials that consist of large numbers of chromophores would be very useful for modeling these systems.