Structural characterization of liquid pyrrole has been obtained by performing classical molecular dynamics simulations with a new parameterization of electrostatic interactions. Despite the relatively simple molecular structure of pyrrole, a correct and accurate representation of its intermolecular interactions in bulk phase is a challenging task, since these are affected at short range by the quadrupole-quadrupole term. This new parameterization permits not only to correctly describe the liquid structure but also to obtain macroscopic properties in excellent agreement with experiments.bulk properties, force field, molecular dynamics, quadrupole moment 1 | I N TR ODU C TI ON Pyrrole (C 4 H 5 N) is one of the simplest and most representative heteroaromatic molecules. Its five-membered ring is the building block of several molecules of biological and pharmacological relevance, such as amino acids, indole alkaloids, [1] and porphyrins. [2] Some molecules containing the pyrrole moiety find application also as therapeutic compounds like fungicides, antibiotics, anti-inflammatory drugs, [3] cholesterol reducing drugs, [4] and antitumor agents. [5] Pyrrole has been used in technological applications like, for example, polymers with remarkable electric conductivity [6] and in dye-sensitized solar cells. [7,8] The structure of pyrrole (see Figure 1) features an NAH group and a p system, which permit effective self-aggregation into clusters through the formation of a weak NAHÁ Á Á p hydrogen bonds. This behavior has been confirmed by experimental and theoretical studies which have shown the tendency to form T-shaped aggregates, [9][10][11][12][13] which are typical structures found in the so called "quadrupolar crystals". In particular, it has been observed in microwave spectroscopy experiments [14] that the angle formed by two monomers approaches to $55 with a distance between the centers of mass of $4:1 Å.The nature of the NAHÁ Á Á p hydrogen bond and the relevance of high multipole contributions to the liquid structure are not yet completely understood, especially in the liquid state. In fact, although several studies about small clusters are available, only few experimental and theoretical studies have been performed on liquid pyrrole. [13,[15][16][17] Hence, in this work, we present a computational study on liquid pyrrole with the purpose of accurately determining its structural features and bulk properties.The structure of liquid pyrrole and its evolution over time has been described by classical molecular dynamics (MD), which allows to estimate macroscopic properties from microscopic models. To obtain results comparable with experiments, it is mandatory to employ an accurate force field (FF), which describes the interactions at atomic and molecular level using a set of functional forms based on a limited number of parameters. Some properties are very sensitive to the electrostatic part of the force field, which in the most widely used FFs is described using partial atomic charges.Int J Quantum Chem. 2018;118:e25554.