A computational study of two promising tweezers

A DFT study was carried out to explore the properties of two nonplanar π-conjugated systems that share dibenzo[a,e]-cyclooctatetraene (DBCOT) as a fundamental element. These systems were presented as molecules with potential use as tweezers and in optoelectronic applications due to their expected nonlinear optical effects. Structure optimizations, TD-DFT calculations, molecular orbitals and topological analysis were performed for C36N2H18O4 and C26H18 structures, formed by one and two 1,3,5,7-cyclooctatetraene rings, respectively, with and without the Na cation. This alkali cation was used as a model to investigate the potential of these molecules as tweezers. The results show a V-shaped structure as a minimum in the ground state for both molecules. This angular structure can trap the Na cation very effectively, especially for one of the cases. In addition, the present analysis opens new studies to explore the use of these systems as tweezers of other cations, anions or aromatic molecules that may exhibit π-stacking. The analysis of the molecular orbitals involved in the main electronic transitions allows us to propose that these systems will have interesting optical properties.


Introduction
In recent years, the study of nonplanar π-conjugated systems has reached a great development.These systems are expected to have important applications in fields, such as molecular electronics, solar cells, optoelectronic applications or sensing materials.
Several different structures have been proposed, ranging from corannulene [1] or cyclooctatetraene [2] to expanded porphyrins [3] or twisted perylene bisimides [4].The size and shape of these structures play a fundamental role in achieving specific optical and electronic properties.Yamaguchi et al. [5] have pointed out the advantages of π-conjugated systems, which combine flexibility and rigidity.Among all these systems, we have had a closer look at two systems that are characterized by a combination of flexibility and rigidity, to be constituted by a flexible central part joining two rigid wings.The flexible unit gives rise to a dynamic conformational change in the excited state, from a nonplanar V-shaped structure to a planar structure.These systems, shown in Fig. 1, can also function as molecular tweezers, picking up other molecules.
The concept of molecular tweezers was first employed by Whitlock in 1978 [6] to name a caffeine-based receptor.A molecular tweezer was defined as a receptor in which two flat arms, usually aromatic, are separated by a rigid spacer.This spacer provides a suitable distance between the flat arms for the inclusion of an aromatic guest molecule through π-π interactions (around 7 Å).
In this study, we worked with two types of tweezers with different nonplanar bending schemes and their planar counterparts.On the one hand, tweezer one (tw1) (C 36 N 2 H 18 O 4 ), synthesized by Yamaguchi et al. [5] is composed of a cyclooctatetraene ring in the center attached by three benzene rings and a pyrrolidine-2,5-dione on each side of the tweezer, as can be observed in Fig. 1A, C. On the other hand, tweezer two (tw2) (C 26 H 18 ), based on proposed in previous works [7,8], is formed by two cyclooctatetraene rings within three benzene rings, whose optimized angular and planar geometries, minimum and transition state in S0 electronic state, respectively, are shown in Fig. 1B, D, respectively.In addition, Fig. 1 includes vertical and front views of the angular and planar structures, respectively.We have conducted a computational study to analyze the suitability of these tweezers for use with cations, taking Na cation as a model.With this objective, not only a structural study of the systems with and without cation was performed, but also a topological analysis and interaction energies, which had not been done previously and which we consider essential to explore the feasibility of these systems as tweezers.

Computational details
All the tweezers' structures were fully optimized at the levels of B3LYP [9], CAM-B3LYP [10], M062X [11], MN12L [12], MN12SX [13], and ωB97X-D [14] density functional theories with two different Pople basis sets, 6-31G(d,p) and 6-31G(2d,2p).Initially for this work we considered the use of fuzzy functions through the 6-31++G(d,p) basis set, but no general improvements in their electronic properties were observed in many of the structures and, in addition, several of them suffered from SCF instabilities that did not allow us to consider them.The optimizations, wave function stability and geometry characterization with frequencies calculations were performed with Gaussian 09 [15].TD-DFT calculations [16][17][18] were performed with Gaussian 16 [19], keeping the same mentioned functionals but using only with the 6-31G(d,p) basis set.
The angle of all nonplanar structures (Ѳ) was measured using a home-made Fortran program designed to calculate the angle between two estimated planes by means of multilinear regression techniques.In this way, we do not consider the dihedral angle of four atoms uniquely but an estimation of the real angle between two planes constituted by a larger set of atoms from two different planes, a more realistic way to study the planarity distortion.
A topological and charge analysis of the electron density, ρ(r), at CAM-B3LYP/6-31G(d,p) was performed using the quantum theory of atoms in molecules (QTAIM) approach [20] for both tweezers with and without Na cation.Moreover, these populations have been compared with other real space partitioning based on the Hirshfeld idea of fuzzy atoms, called fractional occupation Hirshfeld-Iterative (FOHI) [21], where an iterative promolecular density is obtained by atomic SCF calculations optimizing the occupation numbers.The atomic charge of an unspecific atom A was computed as the difference between the nuclear charge Z(A) and the corresponding integrated atomic population N(A): where N(A) corresponds to the integral of a molecular region representing the atom A defined by the following expression: (1)  A, C) and tweezer two (tw2) (B, D), in both angular and planar conformations in the S0 electronic state.All carbon and hydrogen atoms are represented in black and white, respectively, while the oxygen (nitrogen) atoms are colored in red (blue).In the case of the angular structures, the angles were calculated via multilinear regression, as well as the maximum distance in Å between the heavy atoms of the tweezer arms are also included In the case of the QTAIM approach, the weight function, w A , is equal to 1 if w A ∈ Ω A , i.e., if w A belongs to the atomic basin (Ω) of A. In the case of the Hirshfeld-based partitionings, in contrast with the discrete QTAIM definition, w A belongs to the entire molecular space and is expressed as: where NAT is the total number of atoms, and ρ A and ρ n stand for the proatomic densities obtained by the FOHI approach for the atom A and the atom n, respectively.Both intramolecular and intermolecular critical points (CPs), the corresponding bond paths, together with the value of the Laplacian of ρ(r) (∇ 2 ρ b ), the gradient kinetic density (G(r b )), the virial field (V(r b )), and the ratio of the two previous energy densities at the bond critical points (BCPs) were computed using the AIMAll software package [22].The accuracy achieved in the calculation of the QTAIM atomic properties was verified.Thus, the sums over the q(Ω) and the atomic energy values for each molecule reproduce the total charge and molecular electronic energies to within 1.7 × 10 -3 au and 1.70 kJ mol −1 , respectively.All atoms were integrated with an absolute value of the integrated Laplacian function (L(Ω)) less than 4.95 × 10 -3 au.The sum over all absolute values of the atomic L(Ω) was less than 1.28 × 10 -3 au.For the FOHI partitioning, the BRABO package [23] was used in conjunction with the STOCK program [24].
Finally, gas-phase interaction energies were calculated for both tweezers in both angular and planar structures with sodium cation according to Eq. (4).The basis set superposition error (BSSE) was included via the counterpoise (CP) method of Boys and Bernardi [25]: where ΔE int is the interaction energy for the tweezers, E tweezer + Na + is the optimized energy for the tweezer and the cation, E tweezer is the energy from the isolated tweezer, and E Na + is the energy for the sodium cation.
Molecular orbitals and geometries were plotted using the GaussView and Chemcraft visualization software [26,27].

Results and discussion
The two tweezers showed an angular structure at all levels of theory due to the well-known lack of planarity of the cyclooctatetraene.In this sense, we must point out that when a constrained optimization freezing the dihedral angles is performed to obtain a planar structure, a (2) transition state is found for both molecules.The energy differences between the conformations are shown in Tables S1 at different theoretical levels.
V-shaped and planar structures were employed as a starting geometry to perform a full optimization of the complex (with the Na cation), observing that the initial angular structure was in all cases the most stable with all functionals.
Furthermore, sodium cation distance from the tweezers was ~ 2.0 Å in both tw1 and tw2 in the virtual planar geometry, while distance in the angular geometry was ~ 3.3 Å for tw1, and ~ 2.7 Å for tw2.The influence of the DFT functional and basis set used in the calculation on the angle obtained was also analyzed.Angles measurements obtained from multilinear regression in both tweezers were investigated in both basis sets mentioned before in the computational details.Comparing the tweezers with the complexes with sodium, the angles clearly changed in tw1 than tw2 in both basis sets in all the computational levels.Both B3LYP and CAM-B3LYP showed the same change in angle values of 10.4° at the basis set of 6-31G(d,p) and 11.2° at the basis set of 6-31G(2d,2p) for tw1.The most notable change in angles was noticed at tw1 in both basis sets.The most noticeable variation when the cation is included is observed with the functional MN12L for both tweezers in both basis sets.For tw2 the change was lower than tw1 with a value of 1°-5°.The rest of the measurements are summarized in Table 1.On the other hand, Yamaguchi et al. included the dihedral angle of the nonplanar bonds in the cyclooctatetraene [5] from the X-ray structure without the sodium cation (127.2°).According to our results at 6-31G(d,p) basis set, this dihedral angle corresponds to values of 129.4°, 124.7°, 122.7°, 130.9°, 126.9° and 121.7° for B3LYP, CAM-B3LYP, M062X, MN12L, MN12SX and ωB97XD, respectively.Curiously, all these values are systematically lower than those obtained from the multilinear regression analysis, indicating an underestimation of the real molecular angle.Anyway, the closest angles with respect to the experimental one are MN12SX, B3LYP and CAM-B3LYP.
Figure 2 collects the interaction energy data calculated through Eq. 4. It is significant that the interaction is always stronger with the angular structure, being especially very powerful in the second tweezer.This is not surprising if we look at the optimized structures shown in Fig. 3.The tw2 in the angular structure closes around the cation, interacting with the entire structure.However, in tw1 and the planar structures, the interaction is reduced only to the central part.
TDDFT calculations were performed at CAM-B3LYP level with 6-31G(d,p) basis set.The first 20 excited states were considered.For tw2, different functional and basis set were employed as figure S1 shows, concluding the selection 107 Page 4 of 12 of CAM-B3LYP/6-31G(d,p) as theoretical level.This DFT functional is presented as a good selection for TDDFT in different studies [28].Table S2 summarizes the transition energies, main molecular orbitals implied and oscillator strengths (f) for the calculated excitations.Thanks to the computational study, we can assign the two fundamental peaks of the experimental spectrum of tw1 [5] to transits 5 and 11 (we must remember that this is approximate since our spectrum is in gas phase and the experimental in solvent).Figure 4 summarizes the UV-vis spectrum for both tweezers with and without Na cation.For tw1, the introduction of cation does not change the shape of spectrum, it only shifts the band towards a shorter wavelength.In tw2 case, the effect of the cation modifies the shape and height of the spectrum, as well as shifting towards shorter wavelengths.
Also, the greater interaction between cation and tw2 is observed at the orbital level.Thus, as shown in Fig. 5, the HOMO-LUMO gap in tw1 almost does not vary when introducing the cation, while for tw2 there is an increase in the energy difference of almost 1 eV.In Fig. 5 we can also see the shape of the HOMO and LUMO orbitals for the four systems analyzed.
The analysis of the molecular orbitals involved in the main electronic transitions is of great interest when analyzing the system's possibility of presenting important optical properties.Thus, several studies [29] have postulated the relationship between nonlinear optical responses and intramolecular charge transfer, ICT.For this reason, Figures S2 collect the molecular orbitals involved in firsts electronic transitions for both tweezers, showing that some transitions can be involved in process of ICT.As example, we observe that, for the transition of the largest oscillator strength in tw1 (286.6 nm) with 138, 140, 141 and 143 molecular orbitals involved (Table S2), the electron distribution in the occupied molecular orbitals, H-2 (138) and H (140), is mainly localized in the central part of the tweezer and moved to an external part in orbitals L (141) and L + 2 (143).
To understand the interactions between the tweezers and the cation, a QTAIM topological study was performed for both planar and angular tweezers, including sodium as well, to shed light on the nature of this non-covalent interaction.First, we analyzed the electron density critical points (CPs).We show only the results obtained at  CAM-B3LYP/6-31G(d,p), since we did not find significant changes with the other functionals and the larger basis set (QTAIM calculations with all DFT functionals were performed and, as an example, the results with B3LYP, CAM-B3LYP, M062X, MN12L functionals for tw1 are shown in Tables S4).The values of the electron density computed at the ring critical points (RCPs) for tw1 and tw2 are displayed in blue in Figs. 6 and 7, respectively.In all cases, the effect of the cation in modifying the electron density of the RCPs is moderate.In turn, the cation clearly affects the ring(s) closer to the cation for tw1 (cyclooctatetraene), while all the rings for tw2 are slightly modified when the structure is V-shaped, but especially those closer to the Na once again.Therefore, the significant changes are mainly attributed to the effect of the cation-π interaction and an improved confinement and stabilization of the positive charge of the cation on the electronic structure of the tweezers.Regarding the bond critical points (BCPs) and the corresponding associated bond paths (BPs), sodium and tw1 are connected by two BPs, each of which seems to connect the cation to a C=C bond (Fig. 6b, d).However, the direction changes depending on the structure considered.Thereby, in the planar (angular) geometry, the bond path ends at the C1-C1' (C2) double bond and its symmetrically equivalent counterpart of the cyclooctatetraene ring.It is also worth noting the change in direction: while the bond paths in the angular tw1 follow the direction of the molecule, those of the planar counterpart follow a perpendicular direction.
When the sodium cation interacts with the planar tw2 (as shown in Fig. 7b), the corresponding BPs are directly connected to the carbon labeled C1 and its symmetrically equivalent counterpart.In contrast, in the case of tw2 + Na + , the BPs are again attached to double bonds.In addition to the same carbons, as in the angular tw1, the alkali ion is affected by two BPs linked to double bonds from atom C5 and its symmetrically equivalent atoms.This increase in the number of BCPs leads to the appearance of several RCPs and cage critical points (CCPs) of low electron density (between 2 and 4 times less) compared to the typical RCPs observed in the benzene rings of the tweezers, as can be observed in Tables S3 in ESI.The cation tends to approach the carbon atoms with which it can form a stronger interaction, as far as its displacement within the tweezer cavity is allowed.The rest of the BPs and CPs are typical with distances and electron density values very similar to the benzene, cyclooctatetraene and pyrrolidine-2,5-dione molecules.
In Table 2, we summarize the results of some properties of the BCPs involved in the interaction between the cation and the carbons, id est, the atoms of the tweezer implicated in the interaction with Na through a BP.Broadly speaking, the strength of a chemical bond is related to its bond order and is manifested in the electron density at the BCP (ρ b ) [30].Typically, in this molecule, we have observed that C-H and C-C ρ b are between 0.24 and 0.35 au (as seen in the purple values in Figs. 6, 7).This range is far from the values observed in Table 2, where the strength is at least 10 times lower, indicating a weaker interaction than the conventional bonds forming the skeleton of both tweezers.However, it is still possible to highlight some differences between the systems.In the case of both tw1 conformers in Table 2, the two values are similar (~ 0.012 au vs. ~ 0.013 au), in striking contrast to the two tw2 conformations.On the one hand, the planar structure receives the largest value of these, implying a direct strong interaction between C1 and Na atoms.On the other hand, the V-shaped tw2 contains four "anchors" but with weaker BCPs (considering the symmetry, we have only included two in Table 2. Curiously, the sum of both is of the same order as the two BCPs observed in the angular tw1 but slightly larger, suggesting a stronger interaction than those of the two tw1 conformers. The Laplacian of the electron density at the BCP (∇ 2 ρ b ) is a useful tool for classifying bond types.In general, it can indicate whether the interaction is dominated by a covalent (shared) or ionic behavior if the values are negative or positive, respectively [30,31].The negative sign represents an accumulation of electron density on the BCP along the bond path, while values greater than zero denote a depletion of electron density on the BCP toward each nucleus, provoking a closed-shell (ionic) interaction.For example, all C-C and C-H BCPs are less than zero (as seen in Table S7 in the ESI), revealing that these atoms share electrons, while all C-Na + BCPs in Table 2 are clearly positive.This suggests a closed-shell systems interaction, typical of ionic bonds or cation-π interactions.To obtain more information about the nature of the interaction, the value of ∇ 2 ρ b can be decomposed by the local virial theorem [28] into two terms at the BCP, r b , namely the gradient kinetic energy density (G(r b )), which shows the tendency of the system to escape electrons from r b , and the virial field or potential energy density (V(r b )), which represents the ability of the system to accumulate electrons at r b [31]: This expression is defined in au.If we divide it by -G(r b ) and change the second term to the second member of Eq. 5, we obtain the following expression: The value of minus the fraction in the first member of Eq. 6 is positive and depends on the sign of the Laplacian, since the virial field and the gradient kinetic energy density are negative and positive, respectively.If ∇ 2 ρ b is positive (negative), the fraction -V(r b )/G(r b ) is less (greater) than 2 and reflects an ionic or noncovalent (covalent) behavior, since G(r b ) (V(r b )) dominates, and the system prefers to avoid (accumulate) electrons on the BCP.In parallel, the sum of G(r b ) and V(r b ), namely the total energy density, H(r b ), is sometimes combined with this ratio to study the type of interaction formed.However, covalent or closedshell behavior is not absolute, so Espinosa et al. [32] defined a scale consisting of three regions: Region (I) for pure closed-shell interactions whether and region (III) for shared-shell character, where According to Table 2, the values of this ratio are clearly less than one for all the tweezers studied, which means that the type of interaction between the sodium cation and all the tweezers points to a closed-shell behavior.The same observation is certainly true for the total energy density, since all values are greater than zero, even though this long-range interaction is not as strong as that of a typical ionic bond according to the small value obtained.For comparison, we have included the values of H(r b ) (in green) for the rest of the

Tweezers
Bond bonds in Figs. 6 and 7.As can be seen, they are of opposite sign and much stronger than those of C-Na.The previous topological parameters do not provide a readily illustrative tool for analyzing the changes induced by the inclusion of the alkali cation, i.e., they do not give an intuitive idea of which atoms or groups of atoms are mainly involved in the distribution of the positive charge coming from the Na cation.Nevertheless, the QTAIM approach also provides a very useful method for partitioning the molecular space into discrete mononuclear regions identified as atoms.The QTAIM atom definition, with solid and natural foundations, assumes that the boundaries of an atom depend on the surfaces of zero flux in the gradient vector field of the electron density.This partitioning strategy is also compared to another atomic partitioning based on the Hirshfeld idea (Eq.3), where the concept of atom emerges from fuzzy regions defined by a weight function [33].The way the molecular space is separated is completely different for both methods, but together they can provide more insight from different points of view.
We will look at both tweezers without the Na cation.Three types of carbon can be observed, as shown in Figs. 6  and 7. First, those in tw1 (Fig. 6) around the oxygen and nitrogen atoms are strongly charged in both partitionings, even though we will discuss this below.Second, the bridging carbons that are bonded to other C atoms are positively charged (except for C8 in tw1 due to a mesomeric effect).Third, the carbons bonded to the hydrogens (C1, C3 and C5 in tw1, C1, C3, C4, C6 and C7 in tw2, Fig. 7), whose charge is null or slightly negative in the case of QTAIM and, on the contrary, more than -0.122 au when the FOHI scheme is considered.These charges seem to be against the latter partitioning.However, considering the Pauling electronegativity scale, where χ(O) = 3.44 > χ(N) = 3.04 > χ(C) = 2.55 > χ(H) = 2.2, the difference between C and H is large enough to observe this type of electron donation from H atoms.In addition, the carbon atoms are relatively confined in the QTAIM basin compared to the entire molecular space considered in Hirshfeld-based schemes.In any case, the QTAIM basins of these atoms recover the six electrons of the carbon, implying that the FOHI negative charge seems to be caused by the p z delocalization.On the contrary, the bridge carbons (C2, C4 and C6 in tw1, and C2 and C5 in tw2) are positively charged in both partitionings.Again, the FOHI values are larger in absolute value than those from the QTAIM partitioning.Finally, the carbon atoms near the oxygen and nitrogen atoms in tw1 are affected by the strong electronegativity difference.In fact, they are very negatively charged (q(O) = -1.196au and q(N) = -1.259au for the QTAIM partitioning, q(O) ≃ -0.530 au and q(N) = -0.800au for the FOHI partitioning), even though this charge effect is more mitigated in the case of the Hirshfeld-based partitioning scheme.Due to these strong charges coming from O and N atoms, C9, the atom in α, presents a large electron deficit, with 1.510 and 0.805 electrons less for QTAIM and FOHI, respectively.In contrast, C8, the atom in β with respect to the O and N atoms, increases its electron population affected by a mesomeric effect coming from C9, N and O atoms as was observed previously in the indole [31], where the bridging carbon was also negatively charged due to a π-donation in view of the QTAIM analysis at B3LYP/6-311++G(2d,2p) (6.007 au and 1.029 au for the total and π populations, respectively).However, a more detailed study focusing on σ/π populations and 2-delocalization indices to investigate this effect is beyond the scope of this work.
Regarding the hydrogen atoms, as already mentioned, some differences can be observed between the two partitioning schemes.On the one hand, the QTAIM approach keeps the hydrogens slightly charged (H1-H3 and H1-H5 in tw1 in Fig. 6 and in tw2 in Fig. 7, respectively) with values less than 0.015 au in absolute value.However, the hydrogens in tw1 are positively charged (except for H1 in the planar structure) in contrast to tw2 where they are negatively charged.On the other hand, the FOHI hydrogen charges are always positive in any case (between + 0.084 au and + 0.120 au for the hydrogens mentioned above), strictly following the widely accepted electronegativity scale.H4 and H5 in tw1 are a special case, since the former is close to and influenced by the O atom, and the latter is bonded to the N atom, the most negatively charged atom in the molecule, and loses between 0.405 and 0.468 electrons.
When the sodium cation interacts with the tweezers, more than 75% of its initial charge (+ 1 au) is retained in the alkali atom for both tw1 (Fig. 6) and tw2 (Fig. 7).The largest values are observed in the case of the planar structures where the alkali atom uniquely has a ring close to polarize the π electrons, either cyclooctatetraene (tw1) or benzene (tw2), and the charge is recovered at 92.0-92.9% and 91.0-94.4% for tw1 and tw2, respectively.Curiously, planar tw1 provides a very similar QTAIM and FOHI Na charge (+ 0.920 vs + 0.929, respectively), but tw2 presents a relatively important difference, 0.084 au in favor of the QTAIM approach (Na is less charged).Conversely, both angular structures, and especially tw2, counteract the positive charge better because more rings (and therefore more available π electrons) are implied in the interaction with Na.Thus, the alkali atom recovers the 90.9% and 86.8% in tw1, and 89.9% and 77.2% in tw2 for the QTAIM and FOHI splitting, respectively.In this case, the FOHI scheme recovers less than the QTAIM approach, probably due to the difference between the whole molecular space considered with the FOHI partitioning, which gives the opportunity to take more electron density to counteract the positive charge of the sodium than the confined QTAIM basin.The loss of electron density induced by the Na cation reorganizes the charges of both tweezers in the two conformations.By 107 Page 10 of 12 and large, the atoms most affected are the hydrogens, whose population decreases (they are more positively charged), as can be observed in Table 3, where we have included the sum of the population variations of these atoms (∑ΔN(H)) caused by Na + .In this way, H atoms lose between 0.078 au (FOHI value for planar tw1) and 0.195 au (QTAIM value for angular tw2), where the sign indicates the depletion.Furthermore, in both partitioning schemes we observe an important difference between the electron donation from the hydrogens associated with the charge transfer to the sodium cation.The more the hydrogens are depleted, the lower the Na population is in Figs.6B, D and 7B, D. Peculiarly, the same lessening effect happens with the oxygen atom in tw1, i.e., the negative charge decreases by 0.015 au and 0.023 au (0.017 au and 0.028 au) for the planar (angular) structure for QTAIM and FOHI charges, respectively.In the case of the nitrogen, however, the effect is smaller for FOHI charges (Δq = 0.012 au) and virtually zero for QTAIM.
In Table 3, we have collected all the relative populations (ΔN) for all the carbon atoms of tw1 and tw2 interacting with the sodium cation with respect to the isolated molecule.The positive charge coming from the alkali cation attracts the electron density towards the atoms closer to this atom, causing a general depletion of the hydrogen populations, as mentioned before, in view of the electron conservation.This affects the carbons skeleton in various ways.For example, electron shifts provoke several carbon atoms to gain electron density despite the positive Na charge, especially C1 in planar tw1, C2 in angular tw1, and C1 and C2 in tw2.This effect is particularly evident in explaining the bond paths that follow the maximum electron density observed in Figs.6B, D and 7B, D. In the case of planar tw1, C1 increases its population (i.e., it is more negatively charged), and the BP ends in the middle of the C1-C1' bond.On the other hand, in the V-shaped structure, C2 is the atom whose electron density increases more (0.058 au and 0.045 au for QTAIM and FOHI schemes, respectively), and the BP ends in the corresponding bond associated with this atom.However, this is more difficult to distinguish for tw2, where C1 and C2 clearly grow their populations, especially C1 for the planar structure, but the behavior of these atoms in the V-shaped structure is more erratic, with some discrepancies with both partitionings, produced by the four BPs found (Fig. 7D) that share the Na stabilization.In fact, positive variations with the QTAIM approach are observed in two more carbons, C5 and C6, as predicted by the BPs.The FOHI partitioning does not follow the same trend, probably because the grid integration considers the whole molecular space instead of the discrete QTAIM basins.In the planar tw1, the same happens with C3-C8, i.e., while the planar structure presents virtually zero ΔN (< 0.009 au) with the QTAIM approach, the FOHI results are larger in absolute value and C4 and C7 increase their relative population.Strangely, these carbon atoms represent different environments.Thus, C4 is a bridging atom bonded to three carbon atoms, while C7 is one of the atoms bonded to hydrogens and in ɣ position to the O and N atoms (probably affected by resonance effects).A very similar situation is repeated with the angular tw1, but from C5 to C8, because C3 and C4 are closer to the sodium cation than the atoms in the planar conformation and they clearly enlarge their population in both partitionings.In the case of tw2, the effect from C3 to C7 atoms is more atypical.The planar structure loses population at C3 and C4 with the QTAIM approach, while the FOHI scheme shows a different trend, with a significant population increase (decrease) at C3 (C4), since ΔN is greater (less) than zero.Moreover, C5 is unchanged in the QTAIM approach in contrast to the FOHI relative populations, where the value grows.C6 and C7 suffer from electron depletion in both partitionings.On the contrary, the V-shaped tw2 carbons are significantly different.They are closer to the sodium atom and the partial positive charge of the Na causes population accumulation in all atoms except C3 and C4, where the population is slightly modified.

Conclusions
A structural, energetic, orbital and topological analysis of two promising molecules to use as tweezers or in optoelectronic applications were performed.Calculations were carried out with different DFT functional and two basis sets.The formation of complex with Na cation were also considered.
This extensive computational study allows us to verify the potential of the systems analyzed, presenting itself as introductory since it opens new paths of study that we intend to carry out.Thus, on the one hand, we consider necessary an indepth study of the excited states that would allow us to present the potential of these systems in optoelectronic applications.On the other hand, we propose to carry out dynamic studies of the complexation processes with cations or other systems.
Nevertheless, the analysis allows us to conclude important data such as: the preference of the systems for the angular structure, being, in both systems and for all the levels of calculation used, the flat structure a transition state in the ground state.
Quantum Theory of Atoms In Molecules (QTAIM) topology analysis of the electron density clearly indicates a closed shell interaction between the sodium cation and both tweezers.Population analysis, carried out using QTAIM partitioning and a Hirshfeld-based atomic scheme called Fractional Occupation Hirshfeld-I, indicates a small charge transfer from the tweezers to the sodium cation, between 0.034 and 0.250 au, depending on whether these molecules are open or closed, anchored to the cation.Most of these electron distributions come from the hydrogen atoms in both partitioning schemes.The cation causes an increase in the electron populations of the closer atoms whose bonding paths end in the alkali atom.
It is also interesting to note that complexation leads to a major modification on the UV-Vis spectrum for tw2.
In summary, tw1 and tw2 are presented as systems of great interest and new studies are proposed for their use as tweezers, in molecular electronics, sensing or optical applications.

Fig. 1
Fig. 1 Schematic representation of the structures tweezer one (tw1) (A, C) and tweezer two (tw2) (B, D), in both angular and planar conformations in the S0 electronic state.All carbon and hydrogen atoms are represented in black and white, respectively, while the oxygen

Fig. 3
Fig. 3 Schematic representation of optimized angular geometries for tw1 (bottom), tw2 (top) and their interaction with a sodium cation.Three different parameters are presented: angle (Ѳ) in °, the distance

Fig. 6
Fig. 6 Critical points, values and numbering (in black) for tweezer tw1 planar, A without Na cation and B with Na cation, and angular, C without Na cation and D with Na cation, conformations.Values in purple and green at the bond critical points (BCPs), r b , are the Laplacian electron density ∇ 2 ρ(r b ) and the electronic energy density, H(r), respectively.Values in red and orange correspond to the FOHI and QTAIM charges, respectively.Values in blue correspond to the ρ(r b )

Fig. 7 Table 2
Fig. 7 Critical points, values and numbering (in black) for tweezer tw2 planar, A without Na cation and B with Na cation, and angular, C without Na cation and D with Na cation, conformations.Values in purple and green at the bond critical points (BCPs), r b , are the Laplacian electron density ∇ 2 ρ(r b ) and the electronic energy density, H(r), respectively.Values in red and orange correspond to the FOHI and QTAIM charges, respectively.Values in blue correspond to the ρ(r b )

Table 3
QTAIM and FOHI relative populations (ΔN) of both tw1 and tw2 tweezers bonded to the sodium cation with respect to the isolated tweezers P and A denote planar and angular conformations, respectively.∑ΔN(H)representsthesum of all hydrogen atoms, taking symmetry into account.The carbons used in this table follow the numbering from Figs.6 and 7.All values are in au and multiplied by10 3