Synthesis and Molecular Modeling Study of Two Bromo-Dimethoxybenzaldehydes

This work describes the synthesis, molecular structure, and packaging of the compounds 6-bromo-2,3-dimethoxybenzaldehyde (6-BRB) and 5-bromo-2,3-dimethoxybenzaldehyde (5-BRB). Characterization in the solid-state was carried out by X-ray diffraction, and the analysis of the interactions was described by the Hirshfeld surface, which helps in understanding the effect of replacing the bromine position in the aromatic ring. Both compounds (6-BRB and 5-BRB) crystallized in the space group P21/c and are stabilized with C–HO interactions, but only the 6-BRB has a halogen-type interaction. Theoretical calculations, carried out by the density functional theory at the CAM-B3LYP/6-311++G(d,p) level of theory in gas phase, provided information on their electronic properties. The geometric properties were compared to experimental data, and the analyses of the frontier molecular orbitals and the molecular electrostatic potential (MEP) maps were obtained to predict the physical-chemical properties. The supramolecular arrangements were analyzed throughout complexation interactions and corrected by counterpoise method. Hyperconjugation energies were employed to examine the stability of the intermolecular interactions of the compounds. The interaction intensities were determined through the topological parameters obtained through the quantum theory of atoms in molecules model.


Introduction
The brominated methoxyphenyl moiety is widely encountered in natural alkaloids, such as convolutamines, amathamides, and volutamides. 1 This important class of compounds exhibits an array of biological properties, including antibacterial, antifungal, and antitumor activity. 2,3 In this context, it has been found that the biological activity depends on the bromination of the methoxyphenyl unit. 4,5 In particular, we are interested in the use of methoxyaldehydes brominated in different positions as building blocks for the synthesis of these natural alkaloids and derivatives. Out of a series of compounds, two brominated 2,3-dimethoxybenzaldehydes were suitable for study in the solid state.
With the increased understanding of their structure and structure-activity relationships, it has become possible to create new designs and fine-tune molecular functional structures by choosing the nature of the molecules involved, increasingly motivating the synthetic production of compounds. 6,7 The understanding of noncovalent interactions between molecules is advantageous for the design of new crystalline materials, in which hydrogen bonding provides an additional stabilizing effect for the crystal structure. [8][9][10] Among these compounds, those which present a halogen bond have their potential for useful applications described in the literature in the areas of crystal engineering, material science, nonlinear optics, biochemical chemistry, and medicine. 6,[11][12][13][14] This class of interactions plays a fundamental role, for example, by establishing a halogen bond, acting as protein-ligand complexes. 13,15 Halogen bonds are directional interactions between a halogen atom (chlorine, bromine, and iodine) and a free electron pair from a Lewis base, where the distance between the two atoms is less than the sum of their van der Waals radii. They can be classified as a noncovalent interaction, which, due to its specificity, can influence the crystal packing. 6,13,16 Although these are specific intermolecular noncovalent interactions, they considerably influence the stabilization of crystalline states. Density functional theory (DFT) is an approach to quantum mechanics that provides important descriptors for understanding molecular structure based on the electron density of molecular systems. 17,18 Through quantum mechanics calculations, it is possible to predict important information about the molecular and electronic structure of chemical compounds. 19 Furthermore, the molecular structure, interactions, and supramolecular arrangement were characterized by geometric and topological parameters through Hirshfeld surfaces (HS), to understand crystal packaging and help future correlations and structure-activity applications. This theoretical study contains considerable information regarding the replacement of bromine, which can be very important for the design and synthesis of new materials containing halogenated compounds.

Crystallographic analysis
The single crystal X-ray diffraction (XRD) data of 6-BRB and 5-BRB were obtained using an AtlasS2 diffractometer (Agilent SuperNova, Tijuana, B.C. Mexico). The data collection, cell refinements and reduction were carried out via CrysAlis Pro software. 21 The molecular structure was solved on SHELXT and refined on SHELXL through the OLEX2 platform. [22][23][24] The hydrogen atoms were refined by the riding model with the individual isotropic displacement parameters or U iso (H) = 1.5 U eq , according to the C-H model (0.93 Å for aromatic groups and 0.96 Å for methyl groups). The crystal-structure ellipsoid diagrams were generated from an Oak Ridge thermal ellipsoid plot (ORTEP) with thermal-motion probability ellipsoids on the average atomic sites, derived from the atomic displacement matrix, and all supramolecular representations were generated using Mercury. [25][26][27] Furthermore, the hydrogen bonds and molecular interactions were checked by PLATON. 28 Crystallographic information files were deposited in the Cambridge structural database. 29 CrystalExplorer 17 30 was used to plot the HS and their fingerprints. [31][32][33] These surfaces are located in a region of the crystal space around the molecule, defined by a continuous weight function for a molecule in a crystal, with a specified weight function W(r) = 0.5, with the space in the crystal being divided into specific regions, and where the distribution of electrons summed over all the atoms in the molecule (ρ imolecule ) dominates the corresponding sum over the crystal (ρ icrystal ), which is located in a region of the crystalline space around the molecule. 31,34 The HS describes intermolecular interactions through color-scaled mapping that permits the identification and description of surface properties. 31 On a HS, the d norm (normalized distance) is seen in red, white, and blue, where red indicates regions with short intermolecular contacts, white indicates contacts around the separation of the van der Waals radii interface, and blue indicates long contact distances. 35 The normalized distance is constructed by the surface distance to the closest exterior atom (d e ), where the molecules act as the strongest intermolecular contact receptors, by the surface distance to the closest interior atom (d i ), where the molecule acts as the strongest intermolecular donor contacts, and by the van der Waals radii of the internal and external atoms (r i vdW and r e vdW ). 35 The shape index maps on HS can be used to identify the interactions involving the π system (π⋯π interactions) and the complementary hollows (red), where two molecular surfaces touch one another, and bumps (blue), indicating possible interactions, can be described using the shape index. This allows us to identify characteristic packaging modes and flat arrangements, defined as those in which the interval of S is [-1 +1], K 1 and K 2 are eigenvalues of the diagonalization of a hessian matrix and are related to the main curvatures of the surface. In addition, an analysis of intermolecular contacts and their contributions to the packaging of crystals, based on the combination of d e and d i , can generate unique fingerprints for each molecule in the form of 2D plots, which summarize the percentage contribution to the nature and type of intermolecular interaction present in the molecule. 33

Molecular modeling calculations
Theoretical calculations at CAM-B3LYP 36 / 6-311G++(d,p) level of theory, as implemented in the Gaussian09, 36 were carried out for compounds 6-BRB and 5-BRB in order to assess the electronic properties due to the change of the bromine atom position in the compound 2,3-dimethoxybenzaldehyde. The input of gas phase geometry optimization was taken from the X-ray diffraction data, and all the calculations were carried out without any constraint on atomic coordinates. The wave functions were optimized with the help of exchange and correlation functionals CAM-B3LYP and population balance equation (PBE), respectively. The geometric properties were compared with the crystallographic data.
Frontier molecular orbital (FMO) analysis was used to describe the reactivity and kinetics of the compounds. 37,38 Molecular electrostatic potential (MEP) maps were obtained to verify the chemical properties throughout the electronic isodensities surfaces. 36,39 Quantum theory of atoms in molecules (QTAIM), [40][41][42] carried out by Multiwfn software, 43 was used to analyze the types of interactions in the supramolecular arrangement of 6-BRB and 5-BRB. Natural bond orbital (NBO) analyses were carried out to estimate the intermolecular interactions stabilities in the crystalline solid-state. 44,45 The hyperconjugation energies from an electron delocalization between donor (Lewis type) and acceptor orbitals (non-Lewis type) were estimated using the second-order perturbation.

Solid state description
The ORTEP diagram of 5-BRB and 6-BRB compounds ( Figure 1) indicates that both molecular conformations are almost planar. Their crystallographic parameters are presented in Table 1, in which both compounds crystallized on the centrosymmetric monoclinic space group P2 1 /n, with Z = 4 (one independent molecule in the asymmetric unit). In addition, other details such as bond distance, angles and dihedral angles are shown in Table 2, where the molecular coplanarity is indicated by C 6 -C 1 -C 2 -C 3 and C 7 -C 1 -C 2 -C 3 , with the respective dihedral angles of -0.5º and 178.2º for 6-BRB, and of -0.1º and 176.9º for 5-BRB. Furthermore, we have a noticeable change in the dihedral angles of C 3 -C 4 -C 5 -Br, with values of 179º for 6-BRB and -180º for 5-BRB, due to the intermolecular interaction of bromine, which shows the experimental and theoretical results.
The molecular structure skeletons of 6-BRB and 5-BRB with their asymmetric units consist of one aromatic ring with a bromine radical, two methyl groups and one carbonyl. The crystal structure 6-BRB is stabilized by the intra and intermolecular interaction C-H⋯O and the halogen bond Br⋯Br, while in structure 5-BRB there are the intermolecular interactions of type C-H⋯O. All interaction parameters are shown in Tables 3 and 4.
115.578 119.629  (8), with a long carbonyl group along the a and b axes, respectively. These intermolecular interactions contribute to supramolecular arrangements, with a bifurcated interaction involving O1 in intermolecular interactions, and C 9 -H 9C ⋯O 1 and C 4 -H 4 ⋯O 1 forming a tetramer. They also increase in the a and c axis directions, which can be described by the ring R 4 4 (24). The HS contributes to understanding the crystalline packing, analyzing the intermolecular interactions, and recognizing predominant non-covalent interactions, identified by the red spots displayed by the d norm , thus identifying the regions. Mapped over d norm (ranging from -0.0862 to 0.9988 Å for 6-BRB and from -0.1709 to 1.3021 Å for 5-BRB) shows in Figure 3, that the structure of 6-BRB C 8 -H 8B …O 2 (1a), and the halogen bond Br⋯Br(2a) contribute to supramolecular packaging, while the compound 5-BRB is stabilized by hydrogen bonds, C 8 -H 8B ⋯O 2 (1b), C 9 -H 9C ⋯O 1 (2b), C 4 -H 4 ⋯O 1 (3b).
The shape index assists in the interpretation of the interactions that occur in the π system. The π⋯π interaction, represented by the red and blue triangles (bow tie), is where two aromatic rings overlap. For this interaction, the centroid of the aromatic rings for the compounds was calculated, with the distance between the two rings. For the 6-BRB structure (Figure 4a), the overlap appears at a distance of 4.125 Å, showing the contribution of these interactions in packaging and molecular stability, while in the 5-BRB (Figure 4b) we have a distance of 3.961 Å.
The relative surface area resulting from the contributions of each interaction that is exposed in the supramolecular arrangement is quantified by a fingerprint analysis, generated from the graph of (d e vs. d i ). The interactions are highlighted in the fingerprint graphs for structure 6-BRB ( Figure 5a) and for structure 5-BRB (Figure 5b), in which most contributions are related to the H⋯H interaction -35.3 and 36.7%, respectively, and the H⋯O interactions, 26.4 and 21.5%, respectively. They are characterized by peaks at the bottom of the fingerprint graphs and indicate the C-H⋯O interactions of the total HS. Another characteristic is the presence of C⋯C interactions: 7.3 and 6.2% present in the center, in addition to a noticeable difference in Br-Br, of 2.5 and 0.9%, due to the interaction between  them in structure 6-BRB. Furthermore, the surfaces provide information about each donor-acceptor pair and allow us to measure how much the shape divides the surfaces into their scale of colors and stains.
Molecular modeling analysis Table 2 shows the optimized geometric parameters in the gas phase for 6-BRB and 5-BRB at CAM-B3LYP/6-311++G(d,p) level of theory. The mean absolute percentage deviations (MAPD) for both compounds were calculated using the equation: where χ XRD stands for the geometric parameters taken from the XRD data and χ CAM-B3LYP represents the geometric parameters obtained from the theoretical calculations.
The MAPD values obtained for the bond length in 6-BRB were 0.76 and 1.08% for the functionals CAM-B3LYP and PBE (a ratio of 1:1.4), respectively; for 5-BRB, the MAPD values were 0.40 and 0.97% (a ratio of 1:2.4), respectively. These results show that both functionals describe molecular systems similarly; however, the functional CAM-B3LYP resulted in more accurate values for the geometry of both compounds, when compared to the experimental values. The coefficients of determination show that the bond lengths in 6-BRB are described with 99.29% accuracy against 99.24%, respectively for the functionals CAM-B3LYP and PBE; for 5-BRB, the bond lengths are described as 99.83 and  The bromine atom bonded to the C 6 carbon in the 6-BRB molecule exerts a repulsive force on the O 1 atom from the carbonyl group, increasing the C 1 -C 7 -O 1 bond angle slightly when compared to the 5-BRB. The dihedral angles show that the substituent groups on the aromatic ring are in the same plane, except for the methyl groups attached to the O 2 atom, which are outside that plane. The electron density ρ(r) is the fundamental variable in DFT, 17,18 and it allows the definition of some important chemical concepts, such as chemical reactivity of molecules. 37 The highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) enable us to obtain valuable chemical descriptors, such as chemical potential,  hardness, and the global electrophilicity index. These descriptor indices allow us to understand the chemical reactivity and kinetic stability of molecules. 19 Figures 8 and 9 show the HOMO and LUMO energies obtained at the CAM-B3LYP/6-311++G(d,p) level of theory for the compounds 6-BRB and 5-BRB, respectively. HOMO energy is associated with a molecule's ability to donate electrons, while LUMO energy is related to the molecule's ability to receive electrons. The larger LUMO-HOMO gap energy (∆E HOMO-LUMO ) always refers to higher kinetic stability and lower chemical reactivity. The calculated gap values for both compounds are -161. 40 and -162.90 kcal mol -1 , respectively. Therefore, we can conclude that both molecules are chemically stable. The quantum descriptors electronegativity (χ), chemical potential (µ = -χ), hardness (η), softness (s = 1/η), and the global electrophilicity index 46 (ω) can be obtained using the expressions where E is the energy of the system, N is the number of electrons, υ(r) is the external potential, I ≅ -E HOMO is the ionization potential and A ≅ -E LUMO is the electron affinity. The transfer of electrons between two molecules in a chemical process is facilitated if their chemical potentials are different. An electron migrates from a molecule with a larger μ to another with a smaller μ. In chemical processes, a molecule with µ < 0 is stable and does not decompose spontaneously. The calculation results show that compounds 6-BRB and 5-BRB are electronically stable (see Table 5). Furthermore, we can state that both molecules are hard and not polarizable, i.e., the molecules are resistant to small deformations of their electronic clouds. Figure 10 shows the MEP map for 6-BRB. The molecular electrostatic potential map is often used to predict molecular reactivity, and it provides information on a wide range of interactions. 47 The red regions represent regions with higher electronic charge densities. In these   regions, the MEP is negative. The blue areas represent depleted electronic charge density and, therefore, these regions show positive electrostatic potential. Finally, the green spots are associated with null electrostatic potential. Table S2 (SI section) shows the calculated molecular electrostatic potential values for the oxygen and bromine atoms. The bromine atom bonded to the C 6 from 6-BRB has much higher electrostatic potential when compared to the compound 5-BRB (see Figures 10c and 10d). The repulsion in the region between Br and O 1 explains the slight increase observed in the C 1 -C 7 -O 1 bond angle.  Figure 11 shows the HOMO and LUMO plots for two dimers taken directly from the crystal structure. These dimers were chosen because they represent the most important interactions observed for the compound 6-BRB in the solid state.
The C 6 -Br⋯Br-C 6 interaction is the main interaction observed for the first dimer (Figure 11a), and the dominant interaction for the second dimer is C 8 -H⋯O 2 -C 8 (Figure 11b). The first interaction is interesting since it occurs frontally between the bromine atoms of the molecules. The MEP map also shows that the Br⋯Br interaction induces the polarization of the electronic clouds of the molecules. The MEP map also shows that the O 2 and O 3 atoms of the aromatic ring methoxy groups form an electronically charged surface on one side of the molecule. In contrast, on the other side, the electronic population is depleted by the presence of the off-axis methyl group of the molecule, justifying the interaction.
The calculated C 6 -Br⋯Br-C 6 and C 8 -H⋯O 2 -C 8 interactions are respectively -0.47 and -1.51 kcal mol -1 ,  corrected by the counterpoise method. 51 Therefore, the later interaction is the main driving force for the molecular assembling in the crystal. Table 6 shows the natural bond orbital 44,45 analysis for the two 6-BRB selected dimers. The hyperconjugation energy due to electron delocalization between occupied bonding orbitals (donor) and unoccupied antibonding orbitals (acceptor) helps to stabilize the dimer interaction. The hyperconjugation energies were obtained by the second order perturbation energy equation. 45 (5) where 〈s|F|s〉 2 or F ij 2 is the Fock matrix element between the i and j natural bond orbitals, ε s* is the energy of the s * antibonding orbital and ε s is the energy of the bonding orbital s, n s stands for the population occupation of the donor orbital. For the C 6 -Br⋯Br-C 6 interaction, we observed that the s bonding orbital of C 6 -Br from one molecular unit hyperconjugates with the s * antibonding orbital of C 6 -Br of the other molecular unit with a second-order perturbative energy of 0.06 kcal mol -1 . The η 1 (Br), η 2 (Br), and η 3 (Br) lone pair orbitals of the Br from unit 1 hyperconjugates with the s * antibonding orbital of C 6 -Br of unit 2, with hyperconjugation energies 0.44, 0.26, and 0.17 kcal mol -1 , respectively. The same hyperconjugation behavior was observed from unit 2 to unit 1 for C 6 -Br⋯Br-C 6 interaction. For the C 8 -H⋯O 2 -C 8 interaction, the most important contribution to the hyperconjugation from unit 1 to unit 2 is π(C 1 -C 2 ) → π * (C 3 -C 4 ) and π(C 5 -C 6 ) → π * (C 3 -C 4 ) with hyperconjugation energies of 0.16 and 0.13 kcal mol -1 , respectively. From unit 2 to unit 1, the most significant contribution is the hyperconjugation η 2 (O 3 ) → s * (C 8 -H) with second-order perturbation energy of 0.1 kcal mol -1 . In general, interactions between NBO in intramolecular interactions present high E 2 i→j* values, while in intermolecular interactions, E 2 i→j* values are very small. This means that the latter are less stable than to the former. [52][53][54] The Laplacian sign of the electronic density ∇ 2 ρ(r) at bond critical point (BCP) is directly correlated with the electronic density between two nuclear attractors in a bond. If ∇ 2 ρ(r) > 0 and ρ(r) < 0.1 at the BCP, where nuclear attractors support the entire charge concentration; then the interaction is classified as a closed-shell interaction (ionic, van der Waals interactions, or hydrogen bonds). 55 If ∇ 2 ρ(r) < 0 and ρ(r) > 0.2, the electronic charge is accumulated at the BCP; then the interaction is classified as covalent. Finally, if 0.1 < ρ(r) < 0.2 and ∇ 2 ρ(r) < 0, the interaction is classified as partially covalent; but if ∇ 2 ρ(r) > 0 at the BCP, then we have a hydrogen bond and di-hydrogen bond interactions.
The QTAIM analysis shows that the intermolecular interactions for the compound 6-BRB in the crystals have low intensity. Table 7 shows the topological parameters ρ BCP and ∇ 2 ρ BCP calculated at the bond critical points. Figure 12 shows the molecular graph for 6-BRB with the bond paths and the bond critical points (BCP). The electronic density ρ at the bond critical point for Br⋯Br interaction is 0.006 a.u., with a positive Laplacian, and the distance between the two halide attractors is 3.54 Å. The C 8 -H⋯O 2 -C 8 axial interaction electronic densities are 0.011 a.u. and the distance between the two attractors is 2.43 Å. In addition, two more BCPs are observed: C 6 ⋯C 3 and C 4 ⋯C 6 . Both contacts have low ρ values in the intranuclear region. Furthermore, the low total energy of QTAIM value H(r) and the value of the G/|H| > 1.0 indicate that the electronic flux in the internuclear region is very low, resulting in a weak interaction, closed-shell type. 56 Thus, 6-BRB crystals are formed in such a way that the molecules interact through van der Waals forces.
The C 4 -H⋯O 1 -C 7 , C 7 -O 1 ⋯H-C 9, and C 8 -H⋯O 2 -C 8 interactions stabilize the supramolecular arrangement of 5-BRB. In the first interaction (Figure 13a), the high electron density on the O 1 atom of the carbonyl group of one molecule (negative electrostatic potential region) attracts another molecule that has depleted electronic density (positive electrostatic potential region). The energy observed for this interaction is -2.91 kcal mol -1 . For the C 7 -O 1 ⋯H-C 9 (Figure 13b) interaction, the negative electrostatic potential on the O 1 atom of the carbonyl group Table 6. Second-order perturbation theory analysis in the NBO basis obtained at the CAM-B3LYP/6-311G(d,p) level of theory for 6-BRB Donor (i) Occupancy (i) Acceptor (j) Occupancy (j) E 2 / (kcal mol -1 ) E j -E i / a.u. F(i,j) / a.u.  attracts the positive electrostatic potential region on the methyl group (C 9 ) of the second molecule, as shown by the MEP map ( Figure 10). The energy for this interaction corrected by the counterpoise method is -1.95 kcal mol -1 . The C 8 -H⋯O 2 -C 8 (Figure 13c) axial interaction occurs in the region with negative electrostatic potential on the O 2 atom, with a positive electrostatic potential region on the methyl group from another molecular. The calculated complexation energy is -1.30 kcal mol -1 .
The NBOs analysis show that the hyperconjugation energies E 2 i→j are very small for these interactions (see Table 8); that is, the donor ligand orbitals (Lewis type) hyperconjugate very weakly with the acceptor antibonding orbitals (non-Lewis type). Therefore, the interactions that occur in the crystals of the 5-BRB are weak. For the dimer shown in Figure 13a, two hyperconjugations were observed between the lone pair orbitals from the O 1 atom with the s * antibonding orbital of the C 4 -H bond; that is, η 1 (O 1 ) → s * (C 4 -H) and η 2 (O 1 ) → s * (C 4 -H) with hyperconjugation energies of 0.09 and 0.77 kcal mol -1 , respectively. Lone pairs from the O 1 atom also hyperconjugate with the C 9 -H s * antibonding orbitals: η 1 (O 1 ) → s * (C 9 -H) and η 2 (O 1 ) → s * (C 9 -H) with hyperconjugation energies of 0.89 and 0.10 kcal mol -1 . For the dimer depicted in Figure 13b, the most important hyperconjugation occurs between the C 7 -O 1 π bonding orbital with the s * orbital antibonding of C 9 -H with 1.01 kcal mol -1 . For the dimer shown in Figure 13c, we   found two significant hyperconjugations between the lone pair orbitals from O 2 atom and the s * antibonding orbital of C 12 -H, η 2 (O 2 ) → s * (C 12 H) with hyperconjugation energy of 0.99 kcal mol -1 , and with the π* antibonding orbital of C 8 -O 2 , η 1 (O 2 ) → π * (C 8 O 2 ) with hyperconjugation energy of 0.06 kcal mol -1 .
The topological analysis of the molecular system of the compound 5-BRB (Table 9) shows that the molecular electronic densities calculated at the bond critical points (BCP) have small values with ∇ 2 ρ BCP > 0 for all critical points analyzed. The values obtained for H(r) ca. 0, indicating low intensity interactions between nuclear attractors. These results indicate that the interactions between the molecules in the crystal are all of the closedshell type, so the molecules interact by van der Waals forces to form the tetramer. Figure 14a shows the bond paths and the bond critical points for a tetramer taken from the crystal, and Figure 14b shows a dimer from the second packaging layer. The electronic density ρ at the BCP between H and O atoms for the C 4 -H⋯O 1 -C 7 interaction is 0.029 a.u., and the distance between the attractors is 2.52 Å. The electronic density at the BCP for the C 7 -O 1 ⋯H-C 9 interaction is 0.027 a.u., and the distance between the attractors is 2.38 Å. These interactions are bifurcated, with the O 1 atom from the carbonyl group interacting with two hydrogens from other molecular units. Figure 14b shows the structure of the second layer of the packaging of the tetramer structure of 5-BRB. The electronic density at BCP for the C 8 -H⋯O 2 -C 8 interaction shown in Figure 14b is 0.010 a.u., and the distance between the attractors is 2.47 Å. Interestingly, the Br atom from one molecular unit interacts in a bifurcated way with a second molecular unit. The first BCP appears in the Br⋯H-C 7 intranuclear region with an electronic density of 0.032 a.u. and with internuclear attractors' distance of 3.22 Å. The second BCP appears in the intranuclear region of Br⋯H-C 8 interaction with an electronic density of 0.031 a.u. and the attractors' distance of 3.41 Å. In the layered structure, two BCPs for the Br atom appear: the first one appears in the intranuclear Br⋯Br region with an electronic density of 0.005 a.u. The second one is located on the Br⋯C 4 path with electronic density of 0.004 a.u.
The molecular topology analysis shows that the aromatic rings of the 5-BRB molecules in the crystals interact weakly. This interaction is characterized by a BCP on the C 4 ⋯C 4 path with an electron density of 0.005 a.u. A critical cage point (CCP) is also observed in tetramer layers with an electronic density of 0.002 a.u.

Conclusions
The compound 6-BRB is crystallized in the P2 1 /n special group, with four molecules in the unit cell. Its crystalline state is stabilized by C-H⋯O and Br⋯Br interactions and π⋯π stacking. The compound 5-BRB also crystallized in the P2 1 /n special group. With four molecules in the unit cell, however, its crystalline state is stabilized by C-H⋯O interactions and π⋯π stacking. Both interactions were confirmed by their HS and shape index. Theoretical calculations showed that the molecules of the compounds are electronically stable, and that the carbonyl groups have electrophilic activity in both compounds. In their respective crystals, the molecules interact through low-energy interactions, thus being closed-shell, interactions and it is noteworthy that the values obtained for these energies are only an estimate for the crystal system. NBO analyses showed that the donor orbitals interact weakly with the acceptor orbitals in these interactions. Finally, molecular topology showed that electronic charge densities are low in the internuclear regions of the interaction sites, in order to present a van der Waals character.