Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

# 2-Nitro- and 4-fluorocinnamaldehyde based receptors as naked-eye chemosensors to potential molecular keypad lock

## Abstract

New-generation chemosensors desire small organic molecules that are easy to synthesise and cost-effective. As a new interdisciplinary area of research, the integration of these chemosensors into keypad locks or other advanced communication protocols is becoming increasingly popular. Our lab has developed new chemosensor probes that contain 2-nitro- (1–3) and 4-fluoro-cinnamaldehyde (4–6) and applied them to the anion recognition and sensing process. Probes 1–6 are colorimetric sensors for naked-eye detection of AcO/CN/F, while probes 4–6 could differentiate between F and AcO/CN anions in acetonitrile. Using the density functional theory (DFT), it was found that probes 1–6 acted as effective chemosensors. By using Probe 5 as a chemosensor, we explored colorimetric recognition of multiple anions in more detail. Probe 5 was tested in combination with a combinatorial approach to demonstrate pattern-generation capability and its ability to distinguish among chemical inputs based on concentration. After pattern discrimination using principal component analysis (PCA), we examined anion selectivity using DFT computation. In our study, probe 5 demonstrates excellent performance as a chemosensor and shows promise as a future molecular-level keypad lock system.

## Introduction

Rationally designed molecular systems that create new functional materials are of fundamental interest for application in molecular electronics1,2,3,4,5,6,7, opto-electronics7,8,9,10, chemo-sensing11,12, and information security systems13,14. Many organic molecules with light-emitting and light-absorbing properties have been used to build chemical logic devices15,16. The advantage of chemical logic devices over conventional devices are their sensitivity and specificity to chemical structure and their ability to respond to unconventional input signals17. A chemical logic system could be activated only by a specific chemical input. The encryption protects them from detection, tampering, and cracking. Chemical logic gate research has made strides ever since de Silva and coworkers made their major breakthrough18,19,20,21,22. In the field of chemosensors, there has been a rapid evolution from very simple chemical applications to molecular logic computation systems. Chemistry-based logic sensors are now capable of processing chemical information in the same way as electronic logic devices19. Based on chemical systems, a number of molecular switches, logic gates, and logic circuits have been developed in the past decade14,15.

The idea of selective detection, low production costs, and ease of processing has always been central to the success of the research. Organic molecules provide these qualities; however, they display some limitations such as poor solubility in aqueous media, non-specific fluorescence quenching caused by electron transfer mechanisms, heavy transition metals, and spin–orbit coupling, which need to be eliminated23. Small organic molecules have been successfully employed in colorimetric sensing of medically and environmentally important ions and neutral species. It offers onsite and short-time detection, simplicity, selectivity, sensitivity, and reversibility along with minimal or no sample preparation and manual expertise. These chemosensors are increasingly getting popular in chemical logic devices17.

Recent reports describe the application of thiosemicarbazone-based sensors in host–guest interactions and Molecular Logic Gates24. Hydrazone-based sensor molecules have been used in optical molecular switches and sensors for biologically important ions and molecules25,26. There are also pyrrole/indole and phenyl/hydroxyl based sensors which have been proposed for molecular keypad locks and as a molecular switch27. Moreover, calyx-4-pyrrole-based sensors have been reported for pattern generation and recognition sensing for application as potential molecular keypad locks15. These system has added advantage of miniaturization and potential for functional group diversification. There is scope for exploring the pattern-generating molecules to personalized medicine for specific targeted diagnosis and therapy as in the logic-based therapy and they are not confined to sensing and molecular computing28. Moreover, the integration of these anion sensor molecules into molecular level keypad lock devices has become a dynamic area of research29,30,31.

Our interest in substituted thiosemicarbazones is because chemosensors designed with C=S and N–H functional groups have been shown to be efficient30. As well as their low-cost, unique sensing potential for fluoride and cyanide32, and the ability to tune their conjugation via functional group diversification. Recently, we reported improvements in detection limits and binding constant values by optimizing the substituents on the thiosemicarbazides33. Two new chemosensor systems were developed keeping in view the unique properties of nitro- and fluoro- substituted molecules: 2-nitro-thiosemicarbazides (1–3) and 4-fluoro-based-thiosemicarbazides (4–6). These molecular sensors confer an advantage in terms of synthesis ease, modularity, and unique structural properties. Probes 1 through 3 and probes 4 through 6 showed an intense change from colorless to dark red and from colorless to yellow, respectively, in response to fluoride ions, which could be seen by naked eye.

The present work describes the synthesis, characterization, anion sensing properties and binding energies of probes 1–6. In our investigation, we found that 5 was able to distinguish between AcO- and F- in a particular sequential addition, as previously reported for calyx-4-pyrrole based sensors15. These observations led us to investigate a potential molecular keypad lock. In order to simulate molecular level keypad lock, we expect that probe 5 would be able to recognize chemical input patterns and distinguish among the concentrations of the inputs15,34. We tested several combinations of inputs and created a library of data. We performed principal component analysis (PCA) to distinguish the patterns generated in the experiment, and we applied DFT to investigate the chemistry. Therefore, this way of using a molecular system provides a new insight into the detection of anions, as well as novel possibilities for their integration into multifunctional logic devices, which would be of future interest for molecular computing.

## Results and discussion

The chemosensors designed with the presence of –NH and C=S functional groups on probes to increase the reactivity of the probes. The two step synthesis yielded compound 1–6 as probe (receptors) in good yield (76–88%). The structures of all the probes were characterized by with different spectroscopic techniques including FTIR, 1HNMR and 13C NMR, Mass spectrometry. The IR spectrum of probes 1–6 showed characteristics C=S stretching frequencies at 1186–1198 cm−1, for C=N at 1580–1590 cm−1 and for N–H at 3119–3219 and 3318–3334 cm−1. The 1H NMR showed two characteristics singlet peaks at 11.67–10.03 and 9.13 ppm for N–H protons.

### Colorimetric analysis and UV–visible spectral studies

The UV–visible spectra of probes 1–6 were recorded without the anion. The chromophore effect could be seen among the two groups, i.e. nitro-derivatives receptors 1–3 and fluoro-derivatives receptor 4–6. Specifically; upon comparing the absorption maxima of 1 (318 nm) and 4 (342 nm), where the receptor unit were the same however the chromosphere unit were different. The difference calculated due to the chromophore effect was about ~ 24 nm which was attributed to π–π* and/or n–π* transitions in the two separate range groups as expected depending on the electron-withdrawing effect of substituents (Scheme 1). Initially, probes 1–6 were evaluated for anion sensing ability by the naked eye. We found that the reaction of the probe 1–6 with the anions solutions (F, Cl, Br, I, ClO4, AcO, HSO4, CN, and SCN) caused a change in the color which could be seen by naked eyes. The intense color change was observed with fluoride, while less prominent with cyanide and acetate ion solutions. However, other anions did not show a visible change in color (Fig. 1). Probe 1–3 which were nitro-derivatives showed a change from colorless to dark red while the fluoro-derivatives 4–6 changed from colorless to yellow.

The absorption spectrum of probes 1–6 with the anions (F, Cl, Br, I, ClO4, AcO, HSO4, CN, and SCN) reveal that only fluoride and cyanide could cause obvious changes in the spectrum while some change by acetate anion. As shown in Fig. 2 new broad absorption bands between 350 and 500 nm appeared upon addition of fluoride and cyanide anions for probes 1–3 while new absorption band between 350 and 450 nm appeared after the addition of the same anions in probes 4–6. The addition of acetate anion to probes 1–6 showed a decrease in the intensity of the ICT band to a broad shoulder band and concomitant formation of new bands similar to fluoride and cyanide anion addition, however, subtle differences were noted. These overall changes indicate anion binding, through hydrogen bonding at the electron-rich thiourea, attached to electron-withdrawing functional groups which create a push–pull effect and consequently shift the absorption wavelength towards red.

The ratiometric response of the probes 1–6 was studied by incremental addition of fluoride, acetate, and cyanide ions (1 × 10−2 M) to the fixed volume of receptor probes (3 ml, 3 × 10−5 M). During titration of receptor 1–3 with fluoride and cyanide anions the absorption maxima at ~ 320 nm progressively collapsed and the new band at ~ 245 nm and red-shift broad absorption bands around 350 and ~ 450 nm developed with isosbestic points observed at 275, 367 nm (1), 274, 370 nm (2) and 272, 380 nm (3) respectively, indicating a single component was produced. A large ∆λmax of 122, 123, and 139 nm were observed with the change in color of the solutions from colorless to red. These observations were attributed to experimental more favorable towards ICT process (π–π* and/or n–π*) in 1–3 due to charge propagation to the electron-deficient (–NO2) moiety, and a significant intense color change due to extended conjugation35. During titration of receptors with acetate anion, we observed specific and unique changes. We noted that absorption maxima for receptor 1–3 at around 320 nm decreased slightly and new broadband appeared around 350 nm and 460 nm with significant ‘tailing off’. Thus we notice a subtle difference in the spectra of titration with acetate anion.

Titration of receptor 4–6 with fluoride and cyanide revealed that absorption maxima for 4, 5 and 6, at 342, 343, and 339 nm collapsed and a new band at around 246 nm and a new red-shift sharp band at 400, 401 and 402 nm appeared and developed with isosbestic points at 287, 363 nm (4) 296, 362 nm (5), and 292, 359 nm (6) respectively with a change in color of the solutions from colorless to yellow. The ∆λmax of 63, 64, and 67 nm were observed. The new red-shift band induced by fluoride anions match well with recent reports of anion sensors36, suggesting that it could be the response towards deprotonating of receptors 1–6 as also predicted by our theoretical calculation. For receptor 4–6 upon titration with acetate anion as noted above the absorption maxima reduced only slightly and new broad bands appear around 248 nm and at around 406 nm giving a slight yellow color solution. The spectra for 1–6 upon titration with acetate revealed only 10–20% reduction in the absorption maxima. Among all the receptors the least changes in the absorption maxima observed for receptors 3 and 5 (Fig. 2c–e).

A comparison of the fluoride and acetate anions titration plot with model receptor 5 are shown in Fig. 3. It could be seen that the absorption maxima of 5 collapsed upon titration with fluoride while the absorption maxima reduced only 20% upon titration with acetate. Hence we propose that these changes are due to the direct interaction of the acetate anion with the receptor probes in a specific and unique manner and not by a common mechanism such as deprotonation of the thiourea protons, as such a mechanism would have led to the formation of the same product for all the anions and hence, identical absorption spectrum.

The detection limit is an essential function of the anion sensor. All the receptors 1–6 were efficient sensors for fluoride anion with the detection limits in the range of 4.49 × 10–7 M to 2.57 × 10–6 M. The detection limit of 4.49 × 10–7 M for fluoride was much lower compared to some recent reports on chemosensors37. The cyanide anion detection limits was in the range of 4.78 × 10–7 M to 1.98 × 10–6 M. The detection limit of 1.98 × 10–6 M was much less that phenylthiourea derivatives38, and azo-benzyldenethiourea reported previously According to the World Health Organization, concentrations lower than 5.3 mM for fluoride are acceptable in drinking water39, which means that probe 1–6 are valuable tools for the detection of fluoride in safe drinking. The stoichiometry of the receptor to anion ratio was determined by taking equimolar combinations of receptors and tetra-butyl-ammonium salts of anions in CH3CN. The calculations were done using Job’s Plot. All the compounds show receptor: anion ratio as 1:2 for 1–6 (Fig. S2a–f). From the above titration results, the anions indeed interacted with the (N–H) protons of the semicarbazone40. The Benesi-Hildebrand equation was applied to calculate the association constant41. It was found that nitro-based receptors show a lower association constant value compared to fluoro-based receptors. Specifically, comparing 1 and 4, their associations constant were 1.64 × 104 and 1.63 × 105 respectively where the receptor unit are the same only the chromophoric units are different (Table 1). Similar substituent effect have been reported recently42. Among the nitro-based receptors, compound 3 has the highest association constant which could be attributed to the push effect of the tertiary amine group on the receptor. Whereas, among the fluoro-based receptors, probes and 6 had high association constants which are attributed to the para- and ortho-para position of the fluorine substituent on the receptors respectively. It was very predictable that the highest association constant would be for 6 with ortho-para difluoride substituted receptor unit which provides a greater electron-withdrawing effect. Our result are inline with our previous report of Coumarin-thiosemicarbazones and Acridine based thiosemicarbazone receptor compounds as shown in the Table 1, however, the association constant values for compound 6 are higher in the current study compared to the earlier reports. The nature of interaction between thiourea protons and the anion plays a critical role. Most likely the sensing for acetate was due to hydrogen bonding interaction with thiourea protons, however for fluoride initial hydrogen bonding interaction followed by ion-pair complex formation (bifluoride ion) HF2, due to its deprotonating ability43.

Benesi–Hildebrand equation was used to calculate binding constants;

$$\frac{{\varvec{b}}}{\Delta {\varvec{A}}}=\frac{1}{{{\varvec{S}}}_{{\varvec{t}}}{{\varvec{K}}}_{{\varvec{a}}}\Delta{\varvec{\varepsilon}}}\times \frac{1}{\left[{\varvec{L}}\right]}+\frac{1}{{\varvec{S}}{\varvec{t}}\Delta{\varvec{\varepsilon}}}.$$
(1)

ΔA = Asubstrate + anionAsubstrate, St = Total concentration of substrate, Ka = Association constant or Binding Constant, Δε = εsubstrate+anion − εsubstrate − εanion.

Dividing the intercept with slope from BH Plot, Ka was calculated. Compound 1–3 did not show fluorescence. This could be due release of absorbed energy in a non-radiative process such as transferring excess energy to surrounding molecules, decay to lower but close energy states or reacting with surrounding molecule as has been observed for non-fluorescent compounds19. Compound 4, 5, and 6 give fluorescence as expected reason being imine derivatives of 4-fluorocinnamaldehyde having fluoro-substituent attached to phenyl rings on both sides of thiourea moiety. They were further investigated for fluorescence studies. It was observed that a new band has appeared at 525 nm when excited at 405 nm. Incremental addition of TBAF (10 equiv) in a fixed volume of compounds 4–6 (1 × 10–4 mol L−1) was observed to increase the intensity at 525 nm as shown in Fig. 4 below.

### Spectroscopic pattern analysis

All compound (1–6) exhibit capability for sensing for F, AcO, and CN but derivatives 3 and 5 showed sensing responses to AcO anions with distinct spectral patterns. Some recent studies have reported fluoro-substituted phenyl giving a curved conformation to the compounds42. A curved conformations where a cavity is formed limits the receptor/anions interaction by allowing the anions to approach for the –N–H proton based on the anion size44. Among probes 4–6 we selected probe 5 with a fluoro-substituent at meta-position (Fig. 5a) as it has intermediate selectivity for anions but distinct spectral response (Fig. 5b). Hence we anticipated that the fluoro-derivative 5 due to the difference in the interaction of the receptor and anions (AcO/F/CN) could lead us to optimize sequential chemical input and generating distinct spectral responses.

We first analyzed the ability of 5 to distinguish between ratios of the same input anions. We noticed that 5 was able to distinguish between the addition of anions (F or C), followed by the second addition of the same anion (FF and CC), indicating that 5 could lead to the generation of a distinguishable absorption pattern (Fig. 6a). According to this, changing the concentration affects the ratio between them. Secondly, we analyzed the ability of 5 to distinguish between different input anions. We used changing the anion i.e. input to CF, which revealed altered interaction of the probe 5 to give changed spectra (Fig. 6b). Upon reversing the sequence the spectra showed minor differences (Fig. 6c). However, probe 5 could recognize AcO and F at different spectral patterns. Probe 5 on the addition of AcO followed by F could induce 348 nm and 401 nm absorption bands; however, only 401 nm absorption band was observed when reversed the addition sequence. The absorption response at four different wavelengths 254, 317, 340, and 402 nm were recorded to see the change in the spectra between different sequential input systems. The plot in Fig. 6d shows a significant change between the different sequences analyzed. This indicated that the sequential addition of anions favorably interacts with the same sites. This tendency are higher in multivalent receptors which exhibit binding and conformation dynamics which allows the compounds to discriminate among different sequential inputs46. Thus we found that our system could produce spectral patterns of distinctive nature and could also distinguish among “passwords” containing different inputs and different ratios of inputs.

The principal component analysis (PCA, Fig. 7) of the complete spectral data shows the probe 5 could discriminate between all possible permutations of anions (F, C, A, FF, CC, FC, CF, etc.) akin to an equivalent 2-digit electronic keypad device. The PCA score plot in Fig. 7 shows the complete segregation of keypad passwords of different combinations. However, the complete separation among the different combinations of passwords are clearly seen along PC-1. The score plot of the PCA model also shows that PC-1 contributes more to separate among sequential combinations as compare to factor-2. PC-1 carries 57% of the spectral variation and PC2, 32% respectively.

Since probe 5 shows absorption maxima at 348 nm and 401 nm for AcO and F respectively hence sequential chemical input AcO and Fcould be utilized as password of the molecular keypad lock for molecular traffic signal with 2-chemical inputs. These observations inspire us to explore probe 5 to develop molecular keypad locks or chemical passwords. To develop a simple molecular keypad lock, chemical input-1 and input-2 were assigned at “A” and “B” while, the outputs signals at 348 nm and 401 nm were assigned at “M” and “E”. Key “R” was assigned for the “OPEN” state of the lock. Now when the input “A” was added first giving 348 nm followed by “B” giving 401 nm absorption band, when detected by the system, this would result in “OPEN” state of the lock. This sequence of chemical inputs would open the lock and the password would be “AMBER” (Fig. 8). All the other sequences would fail to open the lock.

The changes in the anion (A vs. B) and /or its concentration (i.e. A vs. AA or B vs. BB) could lead to distinct patterns for passwords A, B, AA, and BB. Thus this system could generate unique optical “fingerprints”. Upon addition of OAc, F, and CN anion, as distinct code entries denoted by A, B, or C a more advanced molecular security system that has the capability of the 3-digit passwords consisting of different permutations of 3-input keys could be developed. Although this possess the challenge of differentiating between the 27 combination as shown in Table 2 below, however, many could be readily be differentiated by the molecular device.

We conducted the sequential additions and the final 3-digit molecular password system was created after testing a wide range of permutation and eliminating any group of sequences that generated overlapping spectra. The non-overlapping patterns generated by different passwords consisting of acetate (A), fluoride (B), and cyanide (C) as input signals are presented in Fig. 9a. The resulting PCA plot was applied to differentiate between the unique chemical passwords (i.e., BBB, CCC, CCB, BCA, BCC, CBB, and AAA) and the remaining combinations that resulted in overlapping patterns i.e. Groups A–C. The pattern analysis of the combinations reveals that seven 3-digit passwords could be authorized by the unimolecular security systems (Fig. 9b). The overlapping spectra formed in three groups. Group A formed from overlapping BBA, CCA, ACC, BAA, AAC, CAA, CAC, ACA while Group B formed from overlapping ABA, BAB, CAB, ABC, ABB, CBA, and group C had BBC, ACB, AAB, BCB, CBC, BAC.

Thus our attempted to explore the ability of this system to process a diverse set of chemical inputs shows that probe possesses the ability to generate patterns based on the chemical inputs, secondly the pattern generated by the UV–Vis system was efficient in distinguishing among different concentrations of the chemical inputs. This system could authorize multiple users without changing the “lock” or “the key”. Each user could change his password by altering the chemical input. The chemical anions are all colorless and transparent in the visible region and hence their structural and concentration level cannot be determined visually which gives an added security to this system. This work presents an insight to design anion sensors for future information security systems and multifunctional logic devices.

### Theoretical analysis

The scheme of optimized geometries with B3LYP-D347,48, /6–311 + G(d,p)49,50,51, with water as solvent by applying the implicit solvent model SMD52, are shown in the Fig. 10 and Figs. S9 and S10. These calculations were used to obtain the thermodynamic data to calculate the equilibrium constants, pKa and Hammet constants of sensors. These geometries were after utilized to perform the TD-DFT calculations with ωB97X-D353, functional and basis set ma-def2-TZVP44,54, with water as solvent by applying the SMD model.

Initially it was investigated the reactivity of the chemosensors 16 through frontier molecular orbital (FMO) analysis (Table 3). Though the FMO analysis was verified that overall, the chemosensors showed FMO energies remarkably close with an even larger similarity between the chemosensors derived from nitrocinnamaldehyde, particularly 1 and 2 and between derived from fluorocinnamaldehyde, 46. The compounds 1 and 2 showed lower HOMO and LUMO values and with lower ΔEHOMO-LUMO energy than 4–6 (fluorocinnamaldehyde derivatives). Compound 3 showed an equivalent LUMO energy to 2, however with higher HOMO energy than 1 and 2 thus presenting the lowest ΔEHOMO-LUMO between them. This behavior could be explained by the donor nature of the amine group. The compounds derived of fluorocinnamaldehyde (46) showed a behavior remarkably analogous to compound 0 (non-substituted) demonstrating a low influence of –F groups on the FMO of chemosensors in relation to the –NO2 group. Although the compound 6 (with two –F substituent groups) demonstrated lowest HOMO and LUMO energies than 4 and 5 and higher ΔEHOMO-LUMO.

The analysis of the surfaces of HOMO and LUMO orbitals of compound 0 and 16 showed that the HOMO and LUMO orbitals are remarkably similar for all compounds (Table S1), however, the chemosensors derived from nitrocinnamaldehyde (13) have their LUMO surface with a large concentration in the –NO2 group while the derived of fluorocinnamaldehyde have low participation of –F in the LUMO orbital. Demonstrating the higher influence of –NO2 group on the FMO of chemosensors. Figure 11 demonstrates this behavior. Besides that, it was verified that the HOMO of compound 3 (Fig. 12) was concentrated near of amine group, which explains its slightly higher energy than other chemosensors as discussed before.

The GRD indices of chemosensors are shown in the Table 4 and demonstrated that the compounds 16 have a significant analogous reactivity which was expected from the FMO analysis. The high ΔEHOMO-LUMO showed by the compounds 16 demonstrate high stability in which the 46 indicates slightly higher chemical stability than 13.

The compound 2 present the highest ionization potential (IP) values which indicate higher chemical stability and inertness. Overall, the compounds 13 showed higher values of electronic affinity (EA) than 46 which demonstrate that compounds 13 can easily take electron than 46; this behavior should be related to the –NO2 group. The sensors 13 also showed slightly higher values of electronegativity (χ) and higher negative values of chemical potential than 46, these two indices associated indicate that compounds 13 have a higher tendency to attract electrons than 46. However, the compounds 46 showed higher hardness (η) values than 13; this behavior should be related to the presence of –F groups. The compounds 13 also showed higher values of electrophilicity (ω) than 46. Thus, the chemosensors 13 should attract the anions more than 46. The sensors 13 also presented higher values of softness (σ) than 45.

To further, investigate the excitations of chemosensors were performed the hole-electron analysis of chemosensors. The first ten excitations (states 1–10) of chemosensors 1–6 were calculated (Tables S2–S8). However, to facilitate our analysis just the two main excitations of chemosensors (16) were further evaluated. Thus, the two excitations with the highest fosc were chosen. These excitations are arranged in Table 5 (largest fosc) and Table 6 (second largest fosc). Were observed that for chemosensors 13 (nitrocinnamaldehyde derivatives) the excitations with larger fosc are the excitations S0 → S4 (state 4) and for chemosensors 46 (fluorocinnamaldehyde derivatives), are the excitation S0 → S2 (state 2). Overall, the nitrocinnamaldehyde showed lower wavelengths and higher excitation energies (for its equivalent excitations) than fluorocinnamaldehyde derivatives.

As for FMO and GRD analysis, the hole-electron analysis showed indices values remarkably similar for compounds 16. However, the nitrocinnamaldehyde derivatives (13) presented some slight differences in relation to fluorocinnamaldehyde derivatives (46). Though the analysis of hole and electron distribution were verified that the D indices for 13 present medium values, however, the compounds 46 presented low D indices values which indicate that the S0 → S2 excitations of 46 have a lower distance between hole and electron distribution than S0 → S4 excitations of 13 compounds, and consequently higher ECoul than 46. The Sr index of compounds 46 was also higher than 13. However, the H index of compounds 13 was higher than 46, which indicates that a higher extension of the hole and electron distribution of compounds 13. The t index was substantially < 0 indicating that the main excitation of compounds 16 are a local excitation (LE). Finally, the HDI and EDI indices of sensors 13 showed slightly lower than 46 indicating a higher spatial extension of hole and electron of 13. The hole–electron analysis indices demonstrated that the main excitation of compounds 1-6 is a local excitation (LE)of π–π* type, as verified by experimental results.

The analysis of hole + electron surfaces (Tables S9–S15) demonstrates that the excitation with highest fosc along 16 are again similar. However, it could be observed in Fig. 13 that 1 respective excitation has a higher hole extension and electron spatial distribution than 4 besides a higher separation between hole and electron. This behavior should be related to the presence of –NO2 group in the sensors 13 and demonstrates a higher CT character in 13 than 46.

For the excitations of sensors 16 with the second highest fosc, although the values of hole-electron analysis indices are remarkably similar, the trend between 13 and 46 were exchanged. The sensors 13 showed higher wavelength and lower excitation energies than 46. For this excitation, the sensors 46 showed medium values of D indices as for 13 sensors, although also lower than 13. The Sr index values of this excitation are lower than the excitation with the highest fosc which indicates a possibility of n–π* excitation type. The t index remained with a value < 0, but with lower values, which indicates an increase of CT in the excitation. The HDI and EDI also are higher than the previous excitation and indicate an increase of CT in the excitation. However, this excitation were still a local excitation.

It was verified, through the analysis of hole + electron surface (Fig. 14) and the analysis of atomic contributions to the hole, that the excitation S0 → S1 of compound 1 have a contribution of 11.93% of S atom to the hole and the excitation S0 → S3 o sensor 4 have a contribution of 50.45% of S atom to the hole (Tables S17, S20, respectively). Therefore, along with the hole–electron indices analysis, the hole + electron surfaces indicate that the excitation S0 → S1 of compound 1 was an excitation of π–π* type while the S0 → S3 of compound 4 was an excitation of n–π* type which demonstrate the importance of the S atom of thiosemicarbazides to the excitations of the molecule which should be reflected in the UV–vis spectrum of sensors.

Besides that, other important excitation with lower fosc beyond low excitation energies are excitation of n-π* type with the n orbitals of S atom. This excitation was present in all sensors and has in general the lowest excitation energy with an exception for compounds 1 and 3.

The plot of the hole and electron distribution of this excitation were shown in the Fig. 15, for compound 1 this excitation was the excitation S0 → S2, but for compounds 0 (reference, non-substituted), 2, 4, 5 and 6 was the excitation S0 → S1, and finally for sensor 3 was the excitation S0 → S3 as could be verified in Tables S9S15. This excitation demonstrates the importance of the S atom as a donor in the sensors. In the compounds 46 the influence of the S atom were even more highlighted representing on average 80% of contribution to the hole of the first excitation (Tables S16S22).

Is also, important highlight, that the hole and electron of compounds 16 are composed of more than one molecular orbital (MO) with different percentages of contributions to the formation of holes and electrons. In the supplementary material are showed the contributions of each MO with the percentage of contribution above 5% to the hole or electron (Tables S16S22). The main MO that contributes to the composition of hole and electron of main excitations of compounds 1–6 are shown in Tables S23S29.

In order to further understand the nucleophilic regions of the sensors were performed a charge analysis of compounds 16. Through the analysis was verified that the bond N–H of thiosemicarbazides were the more nucleophilic bonds and the charges of H atoms and atoms that are neighbored to it were shown in Table 7. The GRD analysis indicates that sensors 13 showed slightly high nucleophilicity and the CHelpG charge distribution of H atoms bonded to N atom from thiosemicarbazides agree with this behavior. The charge analysis demonstrates that H11 and H14 atoms of compounds 13 were more positively charged than 46. Besides that, the H14 atom was more positively charged than H11 for all compounds along with the N14 atom (bonded to H14) was negatively charged against the N11 atom (bonded to H11) was positively charged, which indicates that H14 should be the H atom first dissociated in an acid–base equilibrium. Thus, the H14 atom was considered more acid for the dissociation of the sensor in the acid–base equilibrium to the calculation of acid–base equilibrium and pKa of 16. The optimized structures of sensor anions are shown in the Fig. S1 and the geometries for compound 0 and its anion are shown in the Fig. S2.

To study the acidity of chemosensors and the influence of substituents its $$p{K}_{a}$$ and theoretical Hammet constants were calculated. To validate our results were performed a brief test of methodology by calculating the $$p{K}_{a}$$ and $${\sigma }_{X}$$ of the benzoic acid and some and substituted benzoic acids (R-bza, R = m-NO2–, p-NO2–, m-F–, or p-F) and compare with experimental results55. The results are shown in Tables S30 and S31. It was observed a large error of the calculated $${p{K}_{a}}_{calc}$$ and $${{\sigma }_{X}}_{calc}$$ with implicit solvent SMD in relation to experimental values, as verified by Thapa and Schlegel56. Although, the average error for $$p{K}_{a}$$ was about 34%, the acidity trends and influence of substituents in the reactivity of the compounds was maintained. Thus we considered this method suitable to the study of chemosensors.

The calculated results for receptors 16 are shown in the Table 8 and the thermochemical data used to calculate it are shown in Table S32. The analysis demonstrated that in general compounds 13 was slightly more acids than 46, with an exception for compound 6 that showed a second higher acidity followed by 3. This behavior demonstrates that the push–pull effect in 3 (due to the presence of the donor amine group and the electron-withdrawing substituent –NO2) besides the inductive electron-withdrawing effect of two –F substituents in 6make the receptors more capable to interact with anions as verified by experimental results of association constants. However, the compound 4 showed the lowest acidity along with all receptors, including the reference 0 without substituents, which do not agree with experimental association constants. This behavior demonstrates that other effects should be associated with receptors reactivity beyond the widespread substituent effects.

The calculated Hammet constants $${\sigma }_{X}$$ reveal that all sensors showed positive value of $${\sigma }_{X}$$, which demonstrate that the substituents favored the proton dissociation, except for the compound 4 that showed a $${\sigma }_{X}$$ =  − 0.30. It was observed by analyzing the Hammet constants, that the sensors which most favor the proton dissociation are the compounds 3 and 6, with the highest $${\sigma }_{X}$$ values (3.74 and 2.79, respectively), in agreement with experimental constant dissociations.

The equilibrium constants of proton dissociation equilibrium between sensors and anions (F, AcO, Br, Cl, HSO4, ClO4, CN and SCN) were also calculated and are shown in the Table S34. It was observed through our analysis that just the anions F, AcO and CN presented equilibrium constants above the order of 10–10, demonstrating equilibrium with these three anions are favored, in agreement with experimental results. It was observed through our analysis that just the anions F, AcO and CN presented equilibrium constants above the order of 10–10, demonstrating a equilibrium with this three anions are favored, in agreement with experimental results. Furthermore, the $${K}_{eq}$$ for F, AcO and CN anions showed the follow trend: CN > F > AcO (Table 9), in agreement with experimental trend for binding constants, except for compounds 3 and 6. Is important to highlight that, although the experimental trend was retained, the absolute $${K}_{eq}$$ values are far from binding constants experimental values because of implicit solvent model limitations for the calculation of solvation energy (ΔGsol) for such chemical systems. However, as discussed before, this model was considered suitable for a qualitative evaluation of dissociation equilibrium between sensors and anions.

Is expected that anions interact with sensors though H-bonds with thiourea protons before its deprotonation though a excited state proton transfer (ESPT)57,58. In order to evaluate this anion-sensor interactions were performed geometry optimization calculations of receptors interacting with anions through H-bonds with thiourea protons. The relative Gibbs free energies and enthalpy energies for the receptor-anion interaction were calculated and their values are shown in the Tables S35S36. It was verified that the Gibbs free energies (Table S35) presented negative values for receptor-anion interactions for the complexes with F anion. Meanwile the receptor complexes with CN and AcO anions showed endergonic, except for compound 2 that showed slightly negative values of interactions with CN and AcO (− 0.03 and − 0.04 kcal/mol, respectively). However, the enthalpy energies (Table S36) showed negative values and demonstrated that the strength of interaction between receptors and anions follows the trend F > AcO > CN. This trend is different of the binding energies trend that showed higher binding constants for CN than AcO, but this behavior is associated to the H-bonds with thiourea protons and neglect the equilibrium and interactions between anions and solvent as discussed before.

To further investigate the H-bonds between receptors and anions F, CN and AcO were performed AIM analysis of these complexes. The AIM molecular graphs for the receptor-anion complexes are shown in Fig. 16 (Compound 1 interacting with F, CN and AcO anions) and Fig. S11 (compounds 16 interacting with F anion), Fig. S12 (compounds 16 interacting with CN anion), and Fig. S13 (compounds 16 interacting with AcO anion). The AIM properties for BCPs a and b with are relate to the H-bonds with thiourea hydrogen are shown in Table S37. It was observed, through the analysis of molecular graphs, the presence of bond paths, BPs (orange lines), between H atoms (of thiosemicarbazide group) and anions which proved the existence of the H-bond between sensor and anions. The bond critical points, BCPs (tiny orange spheres), of the interaction were labeled as a and b.

The analysis of AIM properties (Table S37) demonstrate that H-bonds (BCP a and b) are non-covalent interactions due the positive values of Laplacian of the electronic density, 2ρ(r), a long with electronic density values, ρ(r) < 0.1 a.u. In order to simplify the analysis of H-bonds strength through the AIM properties were performed calculations of binding energy of the H-bonds through the equation presented by Espinosa59, that allows to calculate the BE by applying the density of potential energy, V(r), at the BCP. The BE values are shown in the Table 10 and demonstrate a similar behavior to the enthalpy values for receptor-anions interaction. It was observed the same trend of H-bond strength (H-bond stronger between sensor and F followed by AcO and CN, respectively) than verified trough enthalpy energy analysis.

Moreover, other important observation is that compound 4 has higher experimental association constant values than 1. The theoretical calculations (through GRD, hole-electron and Hammet constant analysis) demonstrate that 1 have more electrophilic character than 4, due to the electron-withdrawing effect of –NO2 group that exerts a withdrawal effect by resonance, which was verified through the hole–electron analyses of excitations with great participation of MO at –NO2 group. Although, this electrophilic behaviour for 1, the AIM analysis demonstrate that the BCP a, related to the H-bond between anion and the hydrogen that should be deprotonated, has higher BE for 4 (BE =  − 13.68 kcal mol−1) than 1 (BE =  − 12.86 kcal mol−1), that should related with higher experimental association constant values for 4 than 1.

This trend observed for compounds 1 and 4, through the theoretical calculations agrees with the experimental55, trends presented by Hammett constants($${{\sigma }_{m-N{O}_{2}}}_{exp}$$= 0.71 and $${{\sigma }_{p-F}}_{exp}$$= 0.06) which demonstrate a higher electron withdrawing of m-NO2 group (substituent in 1) than p-F (substituent in 4). Thus, the 4 experimental association constants higher than 1 should be related with other effects instead of widespread p-F electron-withdrawing. Overall, the theoretical results indicate that compounds 3 and 6 showed the higher electrophilic character and capacity to interact in the anions which agree with experimental association constants, although all the sensors presented in this work showed a remarkably similar reactivity behavior.

## Conclusion

In Summary, we have developed new types of cinnamaldehyde (2-nitro & 4-fluoro) sensors 1–3 and 4–6, respectively, for anion sensing based on substituted thiosemicarbazides. These have been characterized by UV–Vis, IR, 1H NMR, 13C NMR, and mass and theoretical calculations. The theoretical analysis demonstrates that all receptors 16 showed a similar reactivity, although the GRD beyond the analysis of calculated pKa and Hammet constants showed that compounds 3 and 6 showed higher electrophilicity in agreement with experimental association constants. The analysis of proton dissociation equilibrium between 16 and anions showed that higher equilibrium constants for equilibrium between 16 and anions F, CN and AcO. The electronic structures analysis of substituent factors demonstrates that –NO2 substituent showed a high withdrawing effect than –F substituent which reflected in the excitations of compounds and their nucleophilicity, pKa, and Hammet constants. The presence of –NO2 demonstrates highlighted influence when the amino group were also a substituent in receptor due to the push–pull effect which increases the electrophilicity of 3. It was observed that the main excitations in 16 was π–π* type, however, it was verified that overall n–π* excitations showed lower excitation energies with a large participation of S atom to the excitation. Besides that, compounds 46 showed higher participation of the S atom in the excitations due to the influence of –F substituent. It was observed that all receptors 1–3 could recognize fluoride anion in ACN solvent and receptor 4–6 could recognize AcO- and F- in ACN medium. Further, receptor 5 with AcO and F has been used to construct a password for molecular keypad lock with AcO and F as inputs. Combinatorial sensing was demonstrated as a powerful technique for password protection at the molecular scale. Eleven 3-digit passwords could be authorized by the uni-molecular security system. This work provides new insight into the design of anions sensors for future multifunctional logic devices.

## Materials and methods

All the reactions were carried out at room temperature. All chemicals and solvents (HPLC grade) were procured from Merck and were used without further purification. Tetra n-butyl ammonium salts for anions employed for colorimetric analysis were obtained from Sigma Aldrich and stored in a fridge or desiccators until use. Progress of the reaction was followed by thin layer Chromatography. Shimadzu 840/Shimadzu prestige-21 and Bruker alpha FT-IR were used for the measurement of FT-IR spectra of the samples. 1HNMR and 13C NMR spectra were recorded using Bruker (Rhenistetten-Forchheim, Germany) AM 300 MHz and 75 MHz spectrometers. MestNova software was used for integration and calculating coupling constants. Deuterated Dimethyl Sulfoxide (DMSO) and Chloroform (Chloroform-d) were used as standards with specific peaks. UV–Visible and fluorescence spectra were recorded using Shimadzu UV-1800 and spectrofluorophotometer Shimadzu RF-6000 respectively. UV absorbance data were recorded in the range from 200 to 800 nm. Absorbance was recorded in quartz cell (1-cm Width).

### Synthesis of chemosensor compounds (1–6)

#### General synthesis

The functionalization of receptors compounds 1–6 were carried out in two steps following our previously reported protocol. Briefly, the key thiosemicarbazide based receptor-ligands were prepared by reflux reaction in ethanol and addition of hydrazine hydrate (18 mmol) drop wise in reaction vial containing isothiocyanate (15 mmol) at 80 °C for 1–2 h. The reaction was followed by TLC until completion. Precipitates were collected by placing the reaction mixture in an ice bath and then filtered and washed repeatedly using cold ethanol. In the next step reflux condensation reaction of 2-nitrocinnamaldehyde (1–3) and 4-fluorocinnamaldehyde (4–6) with corresponding thiosemicarbazides in a 1:1 stoichiometry, at 80 °C in ethanol in presence of Hydrochloric acid (0.1 mmol) as catalyst afforded the pure products (Scheme 1).

#### General method for the synthesis of receptor compounds (1–6)

The receptor compounds (1–6) were prepared by reacting 2-Nitro (3.28 mmol, 0.581 g) or 4-Fluorocinnamaldehyde (3.28 mmol, 0.492 g) with corresponding thiosemicarbazide (3.31 mmol) at 80 °C in ethanol. Hydrochloric acid (0.1 mmol) was used as catalyst. Solid obtained was separated by the process of filtration and then wash with cold ethanol. Product was dried in desiccators and recrystallized in Chloroform.

#### 4-(4-Fluorophenyl)-1-((E)-3-(2-nitrophenyl) allylidene) thiosemicarbazide (1)

Yield (0.959 g, 85%). M.P 146–148. ESI–MS: m/z 345.077 (M+); Calculated: 344.363, IR (CH2Cl2): 1186 cm−1 ν(C=S), 1585 cm−1 ν(C=N), 3119 & 3318 cm−1 ν(N–H). 1HNMR (600 MHz, Chloroform-d) 11.54 (s, 1H, NH, H-13), 9.13 (s, 1H, NH, H-11), 8.24 (bs, 1H, H3), 7.81 (d, J = 7.6 Hz, 2H, H16,18), 7.59 (d, J = 7.3 Hz, 1H, H6), 7.48 (bs, 1H, H7), 7.26–7.32 (m, 2H, H4,5), 6.97–7.0 (m, 3H, H15,19, H9), 6.74 (dd, J = 15.3 Hz, 9.5 Hz, 1H, H8), 13CNMR (151 MHz, Chloroform-d) 175.3, 155.4, 153.7, 147.3, 143.8, 133.4, 132.8, 130.8, 129.5, 128.7, 127.7, 126.1, 125.9, 125.7, 125.7, 125.1, 124.3, 123.3, 123.3, 114.8, 114.7.

#### 4-(4-Cyanophenyl)-1-((E)-3-(2-nitrophenyl) allylidene) thiosemicarbazide (2)

Yield (1.014 g, 88%). M.P 115–153. ESI–MS: m/z 352.0484 (M+), Calculated: 351.382, IR (CH2Cl2): 1189 cm−1 ν(C=S), 1584 cm−1 ν(C = N), 3138 & 3337 cm−1 ν(N–H). 1HNMR (600 MHz, Chloroform-d) 11.67 (s, 1H, NH, H-13), 9.24 (s, 1H, NH, H-11), 7.94 (bs, 1H, H3), 7.72 (dd, 2H, J = 7.7 Hz, H16,18), 7.67 (bs, 1H, H6),7.46 (bs, 1H, H7), 7.39 (bs, 1H, H4), 7.23–7.15 (m, 4H, H5, 9, 15, 19), 6.51 (dd, J = 13.26, 10.74 Hz, 1H, H8). 13CNMR (151 MHz, Chloroform-d) 187.3, 175.1, 147.2, 144.2, 138.8, 133.6, 132.9, 130.7, 129.3, 128.9, 128.8, 128.1, 127.7, 127.6, 126.4, 124.4, 117.9, 111.6, 77.3, 77.1, 76.8, 39.9, 39.8, 39.6, 39.5, 39.4, 39.2, 39.1.

#### 4-(4-(Dimethylamino)phenyl)-1-((E)-3-(2-nitrophenyl) allylidene) thiosemicarbazide (3)

Yield (0. 958 g, 79%). M.P 165–167. ESI–MS: m/z 370.1462 (M+), Calculated: 369.440. IR (CH2Cl2): 1198 cm−1 ν(C=S), 1590 cm−1 ν(C=N), 3132 & 3334 cm−1 ν(N–H). 1HNMR (600 MHz, Chloroform-d) 10.03 (s, 1H, NH, H-13), 8.96 (s, 1H, NH, H-11), 8.02 (d, J = 8.0 Hz, 1H, H3), 7.77 (d, J = 9.1 Hz, 1H, H6), 7.72 (d, J = 7.7 Hz, 1H, H7), 7.64 (t, J = 7.5 Hz, 1H, H5), 7.50–7.46 (m, 2H, H4, 9), 7.44–7.41 (m, 2H, H15, 19), 6.88 (dd, J = 15.7 Hz, 9.2, 1H, H8), 6.77 (d, J = 8.3 Hz, 2H, H16, 18), 2.99 (s, 6H, H-20, 21). 13CNMR (151 MHz, Chloroform-d) 176.4, 149.3, 147.9, 143.0, 134.3, 133.4, 131.5, 129.4, 128.3, 126.4, 125.2, 112.5, 77.37, 77.2, 76.9, 40.8.

#### (E)-4-(4-fluorophenyl)-1-((E)-3-(4-fluorophenyl) allylidene) thiosemicarbazide (4)

Yield (0. 794 g, 76%). M.P 148–149. ESI–MS: m/z 318.0918 (M+), Calculated: 317.356. IR (CH2Cl2): 1188 cm−1 ν(C=S), 1586 cm−1 ν(C=N), 3219 & 3328 cm−1 ν(N–H). 1HNMR (600 MHz, Chloroform-d) 10.88 (s, 1H, NH, H-13), 9.19 (s, 1H, NH, H-11), 8.37 (bs, 1H, H19), 7.78 (d, J = 9.2 Hz, 1H, H9), 7.40–7.38 (m, 2H, H16, 18), 7.09–7.06 (m, 3H, H3, 5, 15), 6.99 (t, 2H, J = 8.3 Hz, H2, 6), 6.86 (d, 1H, J = 16.0 Hz, H7), 6.74 (dd, 1H, J = 15.9 Hz, 9.3 Hz, 8H). 13CNMR (151 MHz, Chloroform-d) 175.5, 164.0, 162.3, 155.8, 154.2, 145.2, 139.2, 132.2, 128.8, 128.8, 126.5, 126.2, 126.2, 125.5, 124.5, 123.8, 116.0, 115.9, 115.3, 115.2.

#### (E)-4-(3-fluorophenyl)-1-((E)-3-(4-fluorophenyl)allylidene) thiosemicarbazide (5)

Yield (0.835 g, 80%). M.P 156–157. ESI–MS: m/z 318.0900 (M+), Calculated: 317.356. IR (CH2Cl2): 1189 cm−1 ν(C=S), 1580 cm−1 ν(C=N), 3212 & 3337 cm−1 ν(N–H). 1HNMR (600 MHz, Chloroform-d) 9.89 (s, 1H, NH, H-13), 9.16 (s, 1H, NH, H-11), 7.67 (t, 2H, J = 9.24 Hz, H3,5), 7.42–7.44 (m, 2H, H15,17), 7.35 (d, 1H, J = 8.0 Hz, H7), 7.30 (dd, J = 15.0 Hz, 1H), 7.05 (t, 2H, J = 8.3 Hz, H2,6), 6.87–6.92 (m, 2H, H15, 19), 6.77 (dd, J = 16.0 Hz, 9.2 Hz, H8). 13CNMR (151 MHz, Chloroform-d) 175.4, 164.3, 163.6, 162.6, 161.9, 144.7, 139.8, 139.6, 132.1, 129.9, 129.8, 129.0, 128.9, 123.9, 119.1, 119.1, 116.3, 116.1, 112.7, 112.6, 111.1, 110.9, 77.4, 77.2, 76.9.

#### (E)-4-(2,4-difluorophenyl)-1-((E)-3-(4-fluorophenyl) allylidene) thiosemicarbazide (6)

Yield (0. 970 g, 88%). M.P 138–140. ESI–MS: m/z 336.0817 (M+), Calculated: 335.346. IR (CH2Cl2): 1188 cm−1 ν(C=S), 1585 cm−1 ν(C=N), 3125 & 3334 cm−1 ν(N–H). 1HNMR (600 MHz, Chloroform-d) 10.33 (s, 1H, NH, H-13), 9.01 (s, 1H, NH, H-11), 8.17 (m, 1H, H7), 7.74 (d, 1H, J = 9.2 Hz, H9), 7.42–7.44 (m, 2H, H3,5), 7.04 (t, 2H, J = 8.4 Hz, H, 18,19), 6.92–6.86 (m, 3H, H2,6,16), 6.76 (dd, J = 16.02 Hz, 9.21 Hz, H8). 13CNMR (151 MHz, Chloroform-d) 176.1, 164.3, 162.6, 161.3, 161.2, 159.6, 159.5, 156.5, 156.5, 154.9, 154.8, 145.2, 139.9, 132.1, 132.0, 129.0, 128.9, 127.6, 127.6, 124.0, 122.6, 116.2, 116.1, 111.1, 110.9, 104.3, 104.1, 103.9, 77.4, 77.2, 76.9.

### Spectrophotometric experiments (UV/visible and florescence studies)

UV–visible absorption titrations were performed in water at room temperature. Freshly prepared solutions of sensor 1–6 (5 × 10−5 mol L−1) were used. Aliquots of aqueous stock solutions of anions were added fluoride (F), chloride (Cl), bromide (Br), iodide (I), perchlorate (ClO4), acetate (CH3COO), bisulphate (HSO4), cyanide (CN)-, and thiocyanate (SCN) diluted with distilled water. The solutions were prepared with 5 × 10−5 mol L−1 of sensor in distilled water with varying volumes of anions. For Job’s plot equimolar (1 × 10−4 M) solutions of Tetra-n-butylammonium fluoride (TBAF) and receptors were prepared.

The absorption data were fitted to the equation

$$-\mathrm{log}\left[{A}^{-}\right]=\mathrm{ log}{\beta }_{\mathrm{A}-} +\mathrm{ log}\left[\frac{{A}_{\mathrm{max}}-A}{A-{A}_{\mathrm{min}}}\right],$$

where log [A] is the logarithm molar anion concentration at that point, log βA− is the dissociation constant, Amax is the maximum absorbance at the selected wavelength, Amin is the minimum absorbance at the given wavelength and A is the observed absorbance at that specific wavelength.

Plotting log [(Amax − A)/(A − Amin)] versus the − log [A], the log βA− was derived from the slope of the resulting plot. Titrations were repeated at least twice until the results were reproducible.

Fresh stock solutions were prepared for ligands of strength (5 × 10–5 mol L−1) and anions of strength (5 × 10–4 mol L−1). All compound 1–6 were screened for chemo sensing properties. UV scanning was done between 200 and 800 nm upon addition of 30 equiv. of anion solution in ligands solution. Variation in absorption in a fixed volume of ligands 3 mL (3 × 10–5 mol L−1) was plotted against incremental addition of fluoride (1 × 10–2 mol L−1) Association constants was calculated using Benesi–Hildebrand Eq. (1);

Florescence studies were carried out using emission spectrum of (4–6) molecules having fluorophore nuclei on both sides of the thiourea moiety.

### UV–visible and fluorescence titration

Freshly prepared solutions of probes 1–6 (5 × 10−5 mol L−1) in acetonitrile were used. Similarly the stock solution of (1 × 10−4 mol L−1) of different anions like fluoride (F), chloride (Cl), bromide (Br), iodide (I), perchlorate (ClO4 ), acetate (AcO), bisulfate (HSO4), cyanide (CN), and thiocyanate (SCN)from tetra-n-butylammonium salts were also prepared in acetonitrile. The sensing and recognition abilities were evaluated by measuring absorption and emission during the titration with anions.

### Measurements and pattern analysis

The dependence of the absorption pattern on the sequence of addition was demonstrated by adding 10 equivalents of two different anion to acetonitrile solution of 5 (30 µM, 3 mL) in different orders. The mixture was allowed to equilibrate for about 1 min after each anion addition. Solutions of acetate anion (10 mM) and fluoride (10 mM); Fluoride (10 mM) and Cyanide (10 mM), were used for obtaining the various 2-input chemical passwords. Absorption experiments were performed in three replicates and principal component analysis (PCA) was applied to distinguish between the absorption patterns generated.

As a first step toward developing 3-input chemical passwords, solutions of various Anion were prepared (75 µL, 10 mM) and each solution was added separately to a solution of 5 (30 µM, 3 mL) in Acetonitrile. Following a 1 min equilibration absorbance spectra were recorded in three replicates and PCA was applied to identify chemical inputs that could potentially be used for obtaining 3-digit chemical passwords, namely, the anions that induced the most distinguishable changes to the absorption patterns.

### Principle component analysis

PCA is a non-supervised multivariate exploratory data analysis tool used to reveal the hidden structure within large data sets, to reduce the dimensionality of a complex data set without much loss of information, to extract the most important information from the data, to identify noise and outlier in the data set and to visualize the pattern of grouping based on similarities and dis-similarities in the data. It is also known as a projection method, because it takes information carried by the original variables and projects them onto a smaller number of latent variables called Principal Components (PC). Each PC explains a certain amount of the total information contained in the original data and the first PC contains the greatest source of variation in the data set.

The procedure of PCA is based on to convert a set of correlated variables into a new set of uncorrelated variables called principal components. PCA redistributes the total variance of the data set in such a way that the first principal component has maximum variance, followed by second component and so on. The covariance of any of the principal component with any other principal component is zero (uncorrelated) and they are orthogonal to each other. By plotting PCs important sample and variable interrelationships could be revealed, leading to the interpretation of certain sample groupings, similarities or differences.

In this study PCA was applied to classify sensing behavior of a molecular level keypad lock stimulated by two sequential chemical inputs (AcO- and F-), which has potential for application in security devices. In order to build the PCA models UV–Vis absorption spectral data were used data to build the PCA model and model was internally validated, using a leave one out full cross validation procedure using Singular Value Decomposition (SVD) algorithm with 67 segments and total of 7 components.

### Computational details

All the compounds 16 were optimized with Gaussian 0960, package by applying the DFT method, with B3LYP-D347,48, functional and 6–311 + G(d,p)49,50,51, basis set with water as solvent by applying the implicit continuum solvent model SMD52, B3LYP-D3/6–311 + G(d,p)/SMD. Were also performed the geometrical optimization of a non-substituted thiosemicarbazide derivatives (compound 0) to be used as reference.

In order to explore the non-covalent interactions between receptors and anions it was performed AIM (Atom in Molecules) analysis61,62,63, by applying the Multiwfn package64. The AIM data was calculated with B3LYP-D3/6–311 + G(d,p) method with solvent model SMD (water). In order to obtain a more accurate analysis of the excitations and frontier molecular orbitals (FMO) were performed TD-DFT calculations from the optimized structures of 0 and 16 with ωB97X-D353, functional and ma-def2-TZVP44,54, basis set with continuum solvent SMD (water) in Orca package (version 4.2.1)65,66, were performed and analysis of the global reactivity descriptors (GRD) which were obtained by applying the equations shown below. The GDR analysis is based in FMO that is in the highest occupied molecular orbital (HOMO) and lowest unoccupied molecular orbital (LUMO) energies (in eV) of compounds 0 and 16.

The ionization potential (IP) and electron affinity (EA) were given by the Eqs. (2) and (3), respectively67.

$$IP=-{\mathrm{E}}_{\mathrm{HOMO}},$$
(2)
$$EA=-{\mathrm{E}}_{\mathrm{LUMO}}.$$
(3)

The Koopman’s theorem was used to calculate the global hardness (η)68, the electronegativity (χ)69, and the chemical potential (μ)70, by applying the equations below71:

$$\chi =\frac{\left[IP+EA\right]}{2},$$
(4)
$$\eta =\frac{\left[IP-EA\right]}{2},$$
(5)
$$\mu =\frac{{\mathrm{E}}_{\mathrm{HOMO}}{+\mathrm{E}}_{\mathrm{LUMO}}}{2}.$$
(6)

The electrophilicity index (ω)72, was calculated through the Eq. (7) reported by Parr et al.:

$$\omega =\frac{{\mu }^{2}}{2\eta }.$$
(7)

The global softness (σ)73, was calculate by the Eq. (8).

$$\sigma =\frac{1}{2\eta }.$$
(8)

To study the excitations in the chemosensors was applied the hole–electron theory74,75,76. In this theory the electron leaves the hole and goes to the electron in an excitation. For an excitation that is perfect described as a HOMO → LUMO transition, the HOMO could be described as the hole and LUMO as the electron. However, as for Natural Transition Orbitals (NTOs), commonly the hole and electron are composed by multiple molecular orbitals (MO) contributions with different participation percentages to its representations.

The hole and electron distribution could be quantitative characterized through indexes that measure the extension of hole and electron distributions and the separation between them which could be used to evaluate the type of excitation (charge transference, CT, or local excitation, LE). The hole delocalization index (HDI) and electron delocalization index (EDI) measures the breadth of spatial distribution of hole and electron. The HDI and EDI indices showed larger values for concentrate spatial distributions of hole and electron and lesser HDI and EDI indices for delocalized spatial distributions. In a local excitation (LE) the spatial distribution of hole and electron are close. The charge-transfer excitation (CT) has a separation between the hole and electron spatial distribution. The Sr index is a measure of the extension of the overlap between hole and electron distribution. A large Sr index indicates an excitation of type π–π* while lower Sr index indicates a n–π* type of excitation. The D index is the magnitude of charge transference (CT) length that is a measure of the distance between centroid of hole and electron, thus high values of D index reveals a CT excitation. The H index is a measure of spatial extension of hole and electro distribution. The t index measures the separation degree between hole and electron in CT direction. A t index < 0 indicates a small separation between hole and electron due to CT while t index > 0 reveals a substantially separation between hole and electron distribution that indicate a CT excitation. Finally, the Coulomb attraction energy (Ecoul, also called excitation binding energy between hole and electron) is also related with the excitation type along with D index in which higher D index leads to lower Ecoul values.

To evaluate the de acidity through $$p{K}_{a}$$ values and the influence of –R substituent groups in thiosemicarbazide derivatives were calculate the Hammet constants77, in terms of calculated $$p{K}_{a}$$ values, beyond that these results were compared with calculated aromatic indices. All this indices were calculated from thermochemical data obtained through the geometrical optimizations with B3LYP/6–311 + G(d,p)/SMD method.

The $$p{K}_{a}$$ values were calculated using as reference the work ofThapa and Schlegel56, by applying the Eq. (9) below:

$$p{K}_{a}=\frac{\Delta {G}_{aq}^{*}}{2.303RT}.$$
(9)

The compounds 0 and 16 was considered as a protonated sensor ($$SH$$) in a proton dissociation equilibrium in water (simulated through continuum solvent SMD), thus the equilibrium can be written as $$SH \leftrightharpoons {S}^{-}+{H}^{+}$$. Therefore, the by applying the Eq. (10):

$$\Delta {G}_{aq}^{*} = {G}_{aq}^{*}{(S}^{-})+{G}_{aq}^{*}\left({H}^{+}\right)-{G}_{aq}^{*}\left(SH\right).$$
(10)

In which, $${G}_{aq}^{*}{(S}^{-})$$ and $${G}_{aq}^{*}\left(SH\right)$$ are the Gibbs free energy of deprotonated and protonated, respectively, chemosensor in water (SMD) and $${G}_{aq}^{*}\left({H}^{+}\right)$$ were calculated as demonstrated in the Thapa and Schlegel56, work and presented a value of − 270.296808 kcal/mol. The theoretical Hammet constants77, were calculated from an adaptation of the Eq. (11) presented in the work of Hansch et al.78.

$${\sigma }_{X}=\mathrm{log}{K}_{X}-\mathrm{log}{K}_{H}.$$
(11)

Thus, adapting the Eq. (4) we can write the Hammet constant as the Eq. (12):

$${\sigma }_{X}={p{K}_{a}}_{H}-{p{K}_{a}}_{X}.$$
(12)

In our approach in the adapted Hammet Eq. (5), the $${p{K}_{a}}_{H}$$ is the $$p{K}_{a}$$ of non-substituted thiosemicarbazide derivative (compound 0) and the $${p{K}_{a}}_{X}$$ is the $$p{K}_{a}$$ of chemosensors 16.

To evaluate the chemosensor acid–base equilibrium between sensors ($$SH$$) and anions ($${A}^{-}$$), given by $$SH+{A}^{-}\leftrightharpoons {S}^{-}+AH$$, were calculated its equilibrium constants through the Eq. (13):

$${\Delta G}_{eq} = {[G}_{aq}^{*}{(S}^{-})+{G}_{aq}^{*}\left(AH\right)]-\left[{G}_{aq}^{*}\left(SH\right)+{G}_{aq}^{*}\left({A}^{-}\right)\right] .$$
(13)

The equilibrium constant was calculated through the Eq. (14):

$${K}_{eq}={e}^{\left(\frac{-{\Delta G}_{eq}}{RT}\right).}$$
(14)

The electron excitation analysis, FMO, hole-electron analysis and the AIM analysis were performed by applying the Multiwfn program64, (version 3.8) and the visualization of analysis was performed with VMD software79, and the visualization of optimized structures were performed with CYLview software80.

## References

1. 1.

Chen, H. & Fraser Stoddart, J. From molecular to supramolecular electronics. Nat. Rev. Mater. https://doi.org/10.1038/s41578-021-00302-2 (2021).

2. 2.

Schenning, A. P. & Meijer, E. Supramolecular electronics; nanowires from self-assembled π-conjugated systems. Chem. Commun. 26, 3245–3258 (2005).

3. 3.

Munakata, M., Wu, L. P. & Kuroda-Sowa, T. Toward the construction of functional solid-state supramolecular metal complexes containing copper(I) and silver(I). In Advances in Inorganic Chemistry Vol. 46 (ed. Sykes, A. G.) 173–303 (Academic Press, 1998).

4. 4.

Shah, H. H. et al. New multi-ferrocenyl- and multi-ferricenyl- materials via coordination-driven self-assembly and via charge-driven electro-crystallization. Inorg. Chem. 52, 12012–12022. https://doi.org/10.1021/ic401803p (2013).

5. 5.

Shah, H. H. et al. Long-range intramolecular electronic communication in bis(ferrocenylethynyl) complexes incorporating conjugated heterocyclic spacers: Synthesis, crystallography, and electrochemistry. Inorg. Chem. 52, 4898–4908. https://doi.org/10.1021/ic3024887 (2013).

6. 6.

Shah, H. H. et al. Estimation of energy levels of self-assembled ferrocenyls and investigation of charge-driven electro-crystallization of ferricenyl materials. Energy Procedia 100, 149–154. https://doi.org/10.1016/j.egypro.2016.10.157 (2016).

7. 7.

Shah, H. H., Al-Belushi, R. A. & Khan, M. S. Ferrocenylethynyl—A versatile platform for new molecules to novel materials. SQU J. Sci. 19, 15 (2014).

8. 8.

Shah, H. H. et al. New di-ferrocenyl-ethynylpyridinyl triphenylphosphine copper halide complexes and related di-ferricenyl electro-crystallized materials. Dalton Trans. 43, 9497–9507. https://doi.org/10.1039/C3DT52914B (2014).

9. 9.

Khan, M. S. et al. Synthesis and characterization of platinum(ii) di-ynes and poly-ynes incorporating ethylenedioxythiophene (EDOT) spacers in the backbone. Dalton Trans. 40, 10174–10183. https://doi.org/10.1039/C1DT11010A (2011).

10. 10.

Haque, A., Al-Balushi, R. A., Al-Busaidi, I. J., Khan, M. S. & Raithby, P. R. Rise of conjugated poly-ynes and poly (metalla-ynes): From design through synthesis to structure-property relationships and applications. Chem. Rev. 118, 8474–8597 (2018).

11. 11.

Resendez, A. & Malhotra, S. V. Boronic acid appended naphthyl-pyridinium receptors as chemosensors for sugars. Sci. Rep. 9, 6651. https://doi.org/10.1038/s41598-019-42812-8 (2019).

12. 12.

Babazadeh, S., Moghaddam, P. A., Keshipour, S. & Mollazade, K. Colorimetric sensing of imidacloprid in cucumber fruits using a graphene quantum dot/Au (III) chemosensor. Sci. Rep. 10, 14327. https://doi.org/10.1038/s41598-020-71349-4 (2020).

13. 13.

Jiang, S. et al. Information security scheme based on computational temporal ghost imaging. Sci. Rep. 7, 1–8 (2017).

14. 14.

Andréasson, J. & Pischel, U. Molecules for security measures: From keypad locks to advanced communication protocols. Chem. Soc. Rev. 47, 2266–2279 (2018).

15. 15.

Li, Q., Wang, Z., Xu, J., Yue, Y. & Shao, S. Recognition and sensing of AcO− and F− using a calix[4]pyrrole-derived hydrazone: A potential molecular keypad lock. RSC Adv. 4, 34470–34473. https://doi.org/10.1039/C4RA05499G (2014).

16. 16.

Qu, Y., Hua, J. & Tian, H. Colorimetric and ratiometric red fluorescent chemosensor for fluoride ion based on diketopyrrolopyrrole. Org. Lett. 12, 3320–3323 (2010).

17. 17.

De Silva, A. P. Molecular Logic-Based Computation (Royal Society of Chemistry, 2016).

18. 18.

de Silva, P. A., Gunaratne, N. H. & McCoy, C. P. A molecular photoionic AND gate based on fluorescent signalling. Nature 364, 42–44 (1993).

19. 19.

Wu, L. et al. Förster resonance energy transfer (FRET)-based small-molecule sensors and imaging agents. Chem. Soc. Rev. 49, 5110–5139 (2020).

20. 20.

Daly, B. et al. Molecular memory with downstream logic processing exemplified by switchable and self-indicating guest capture and release. Nat. Commun. 10, 49. https://doi.org/10.1038/s41467-018-07902-7 (2019).

21. 21.

de Silva, A. P. Sense and versatility. Nature 445, 718–719 (2007).

22. 22.

Ling, J., Daly, B., Silverson, V. A. & De Silva, A. P. Taking baby steps in molecular logic-based computation. Chem. Commun. 51, 8403–8409 (2015).

23. 23.

Kaur, B., Kaur, N. & Kumar, S. Colorimetric metal ion sensors—A comprehensive review of the years 2011–2016. Coord. Chem. Rev. 358, 13–69 (2018).

24. 24.

Basheer, S. M., Bhuvanesh, N. S. P. & Sreekanth, A. Analytical and computational investigation on host-guest interaction of cyclohexyl based thiosemicarbazones: Construction of molecular logic gates using multi-ion detection. Mater. Sci. Eng. C 105, 110127. https://doi.org/10.1016/j.msec.2019.110127 (2019).

25. 25.

Hossain, S. M., Prakash, V., Mamidi, P., Chattopadhyay, S. & Singh, A. K. Pyrene-appended bipyridine hydrazone ligand as a turn-on sensor for Cu 2+ and its bioimaging application. RSC Adv. 10, 3646–3658 (2020).

26. 26.

Ray, D., Foy, J. T., Hughes, R. P. & Aprahamian, I. A switching cascade of hydrazone-based rotary switches through coordination-coupled proton relays. Nat. Chem. 4, 757–762 (2012).

27. 27.

Li, Q., Yue, Y., Guo, Y. & Shao, S. Fluoride anions triggered “OFF–ON” fluorescent sensor for hydrogen sulfate anions based on a BODIPY scaffold that works as a molecular keypad lock. Sens. Actuators B Chem. 173, 797–801 (2012).

28. 28.

Peri-Naor, R., Ilani, T., Motiei, L. & Margulies, D. Protein–protein communication and enzyme activation mediated by a synthetic chemical transducer. J. Am. Chem. Soc. 137, 9507–9510 (2015).

29. 29.

Xu, Y. et al. A near infrared fluorescent dye for trivalent ions sensing and working as a molecular keypad lock. RSC Adv. 4, 11634–11639 (2014).

30. 30.

Sheykhi-Estalkhjani, A., Mahmoodi, N. O., Yahyazadeh, A. & Nadamani, M. P. Synthesis of new bis-benzylidene-hydrazides as a sensitive chromogenic sensor for naked-eye detection of CN and AcO ions. Tetrahedron 74, 4868–4874. https://doi.org/10.1016/j.tet.2018.07.050 (2018).

31. 31.

Srivastava, P., Verma, M., Sivakumar, S. & Patra, A. K. A smart FRET probe exhibiting a molecular keypad lock device based on rapid detection of nitric oxide mediated by Cu2+ ion. Sens. Actuators B Chem. 291, 478–484 (2019).

32. 32.

Santos-Figueroa, L. E. et al. Synthesis and evaluation of thiosemicarbazones functionalized with furyl moieties as new chemosensors for anion recognition. Org. Biomol. Chem. 10, 7418–7428. https://doi.org/10.1039/C2OB26200B (2012).

33. 33.

Kaur, N. & Jindal, G. “Switch on” fluorescent sensor for the detection of fluoride ions in solution and commercial tooth paste. Spectrochim. Acta A Mol. Biomol. Spectrosc. 223, 117361 (2019).

34. 34.

Margulies, D., Felder, C. E., Melman, G. & Shanzer, A. A molecular keypad lock: A photochemical device capable of authorizing password entries. J. Am. Chem. Soc. 129, 347–354 (2007).

35. 35.

Singh, A., Mohan, M. & Trivedi, D. R. Design and synthesis new colorimetric receptors for naked-eye detection of biologically important fluoride and acetate anions in organic and arsenite in aqueous medium based on ICT mechanism: DFT study and test strip application. Spectrochim. Acta A Mol. Biomol. Spectrosc. 225, 117522. https://doi.org/10.1016/j.saa.2019.117522 (2020).

36. 36.

Ashraf, A. et al. Synthesis of sensitive novel dual signaling pyridopyrimidine-based fluorescent “turn off” chemosensors for anions determination. Measurement 151, 107267. https://doi.org/10.1016/j.measurement.2019.107267 (2020).

37. 37.

Arabahmadi, R. A reversible fluorescence “ON–OFF–ON” sensor for sequential detection of F− and Cu2+ ions and its application as a molecular-scale logic device and security keypad lock. J. Coord. Chem. https://doi.org/10.1080/00958972.2019.1595604 (2019).

38. 38.

Lin, Y.-S., Zheng, J.-X., Tsui, Y.-K. & Yen, Y.-P. Colorimetric detection of cyanide with phenyl thiourea derivatives. Spectrochim. Acta A Mol. Biomol. Spectrosc. 79, 1552–1558 (2011).

39. 39.

WHO. Guidelines for Drinking-Water Quality Vol. 216, 303–304 (World Health Organization, 2011).

40. 40.

Pandurangan, K., Kitchen, J. A. & Gunnlaugsson, T. Colorimetric ‘naked eye’ sensing of anions using a thiosemicarbazide receptor: A case study of recognition through hydrogen bonding versus deprotonation. Tetrahedron Lett. 54, 2770–2775. https://doi.org/10.1016/j.tetlet.2013.02.107 (2013).

41. 41.

Bayindir, S. & Yararli, K. The easy synthesis of new N-substituted 5-oxindoline-rhodanines and their sensing ability: The recognition of acetate ions in aqueous solution. New J. Chem. 43, 8168–8178. https://doi.org/10.1039/C9NJ01732A (2019).

42. 42.

Islam, M. et al. Receptor-spacer-fluorophore based coumarin-thiosemicarbazones as anion chemosensors with “turn on” response: Spectroscopic and computational (DFT) studies. ChemistrySelect 3, 7633–7642. https://doi.org/10.1002/slct.201801035 (2018).

43. 43.

Ali, H. D. P., Kruger, P. E. & Gunnlaugsson, T. Colorimetric ‘naked-eye’ and fluorescent sensors for anions based on amidourea functionalised 1,8-naphthalimide structures: Anion recognition via either deprotonation or hydrogen bonding in DMSO. New J. Chem. 32, 1153–1161. https://doi.org/10.1039/B715533F (2008).

44. 44.

Weigend, F. & Ahlrichs, R. Balanced basis sets of split valence, triple zeta valence and quadruple zeta valence quality for H to Rn: Design and assessment of accuracy. PCCP 7, 3297–3305 (2005).

45. 45.

Isaac, I. O. et al. Novel acridine-based thiosemicarbazones as ‘turn-on’ chemosensors for selective recognition of fluoride anion: A spectroscopic and theoretical study. R. Soc. Open Sci. 5, 180646. https://doi.org/10.1098/rsos.180646 (2018).

46. 46.

Fasting, C. et al. Multivalency as a chemical organization and action principle. Angew. Chem. Int. Ed. 51, 10472–10498 (2012).

47. 47.

Becke, A. D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 38, 3098 (1988).

48. 48.

Becke, A. D. Density-functional thermochemistry. III. The role of exact exchange. J. Chem. Phys. 98, 5648–5652. https://doi.org/10.1063/1.464913 (1993).

49. 49.

Binning, R. Jr. & Curtiss, L. Compact contracted basis sets for third-row atoms: Ga–Kr. J. Comput. Chem. 11, 1206–1216 (1990).

50. 50.

McGrath, M. P. & Radom, L. Extension of Gaussian-1 (G1) theory to bromine-containing molecules. J. Chem. Phys. 94, 511–516 (1991).

51. 51.

Curtiss, L. A. et al. Extension of Gaussian-2 theory to molecules containing third-row atoms Ga–Kr. J. Chem. Phys. 103, 6104–6113 (1995).

52. 52.

Marenich, A. V., Cramer, C. J. & Truhlar, D. G. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J. Phys. Chem. B 113, 6378–6396. https://doi.org/10.1021/jp810292n (2009).

53. 53.

Lin, Y.-S., Li, G.-D., Mao, S.-P. & Chai, J.-D. Long-range corrected hybrid density functionals with improved dispersion corrections. J. Chem. Theory Comput. 9, 263–272 (2013).

54. 54.

Zheng, J., Xu, X. & Truhlar, D. G. Minimally augmented Karlsruhe basis sets. Theor. Chem. Acc. 128, 295–305 (2011).

55. 55.

Hollingsworth, C. A., Seybold, P. G. & Hadad, C. M. Substituent effects on the electronic structure and pKa of benzoic acid. Int. J. Quantum Chem. 90, 1396–1403 (2002).

56. 56.

Thapa, B. & Schlegel, H. B. Density functional theory calculation of p K a’s of thiols in aqueous solution using explicit water molecules and the polarizable continuum model. J. Phys. Chem. A 120, 5726–5735 (2016).

57. 57.

Basheer, S. M., Haribabu, J., Bhuvanesh, N. S., Karvembu, R. & Sreekanth, A. Naphthalenyl appended semicarbazone as “turn on” fluorescent chemosensor for selective recognition of fluoride ion. J. Mol. Struct. 1145, 347–355 (2017).

58. 58.

Chen, J.-S., Zhou, P.-W., Zhao, L. & Chu, T.-S. A DFT/TDDFT study of the excited state intramolecular proton transfer based sensing mechanism for the aqueous fluoride chemosensor BTTPB. RSC Adv. 4, 254–259 (2014).

59. 59.

Espinosa, E., Molins, E. & Lecomte, C. Hydrogen bond strengths revealed by topological analyses of experimentally observed electron densities. Chem. Phys. Lett. 285, 170–173 (1998).

60. 60.

Frisch, M. et al. Gaussian 09, Revision D 01, 2009 (Gaussian Inc., 2009).

61. 61.

Bader, R. F. W., Atoms in Molecules: A Quantum Theory. International Series of Monographs on Chemistry. (Clarendon Press, 1994).

62. 62.

Bader, R. F. A quantum theory of molecular structure and its applications. Chem. Rev. 91, 893–928 (1991).

63. 63.

Kumar, P. S. V., Raghavendra, V. & Subramanian, V. Bader’s theory of atoms in molecules (AIM) and its applications to chemical bonding. J. Chem. Sci. 128, 1527–1536 (2016).

64. 64.

Lu, T. & Chen, F. Multiwfn: A multifunctional wavefunction analyzer. J. Comput. Chem. 33, 580–592. https://doi.org/10.1002/jcc.22885 (2012).

65. 65.

Neese, F. Software Update: the ORCA program system, version 4.0. Wiley Interdiscipl. Rev. Comput. Mol. Sci. 8, e1327 (2018).

66. 66.

Neese, F. The ORCA program system. WIREs Comput. Mol. Sci. 2, 73–78. https://doi.org/10.1002/wcms.81 (2012).

67. 67.

Parr, R. G., Szentpály, L. V. & Liu, S. Electrophilicity index. J. Am. Chem. Soc. 121, 1922–1924 (1999).

68. 68.

Chattaraj, P. & Roy, D. Perennial review: Update. Chem. Rev. 107, 2065 (2007).

69. 69.

Parthasarathi, R., Padmanabhan, J., Elango, M., Subramanian, V. & Chattaraj, P. Intermolecular reactivity through the generalized philicity concept. Chem. Phys. Lett. 394, 225–230 (2004).

70. 70.

Saravanan, S. & Balachandran, V. Quantum mechanical study and spectroscopic (FT-IR, FT-Raman, UV-Visible) study, potential energy surface scan, Fukui function analysis and HOMO-LUMO analysis of 3-tert-butyl-4-methoxyphenol by DFT methods. Spectrochim. Acta A Mol. Biomol. Spectrosc. 130, 604–620. https://doi.org/10.1016/j.saa.2014.04.058 (2014).

71. 71.

Parr, R. G. & Pearson, R. G. Absolute hardness: Companion parameter to absolute electronegativity. J. Am. Chem. Soc. 105, 7512–7516. https://doi.org/10.1021/ja00364a005 (1983).

72. 72.

Parr, R. G., Donnelly, R. A., Levy, M. & Palke, W. E. Electronegativity: The density functional viewpoint. J. Chem. Phys. 68, 3801–3807. https://doi.org/10.1063/1.436185 (1978).

73. 73.

Koopmans, T. Über die Zuordnung von Wellenfunktionen und Eigenwerten zu den Einzelnen Elektronen Eines atoms. Physica 1, 104–113. https://doi.org/10.1016/S0031-8914(34)90011-2 (1934).

74. 74.

Liu, Z., Lu, T. & Chen, Q. An sp-hybridized all-carboatomic ring, cyclo [18] carbon: Electronic structure, electronic spectrum, and optical nonlinearity. Carbon 165, 461–467 (2020).

75. 75.

Yuan, J. et al. Photoswitchable boronic acid derived salicylidenehydrazone enabled by photochromic spirooxazine and fulgide moieties: Multiple responses of optical absorption, fluorescence emission, and quadratic nonlinear optics. J. Phys. Chem. C 123, 29838–29855 (2019).

76. 76.

Jahantigh, F., Ghorashi, S. B. & Belverdi, A. R. A first principle study of benzimidazobenzophenanthrolin and tetraphenyldibenzoperiflanthene to design and construct novel organic solar cells. Physica B 542, 32–36 (2018).

77. 77.

Hammett, L. P. The effect of structure upon the reactions of organic compounds. Benzene derivatives. J. Am. Chem. Soc. 59, 96–103. https://doi.org/10.1021/ja01280a022 (1937).

78. 78.

Hansch, C., Leo, A. & Taft, R. W. A survey of Hammett substituent constants and resonance and field parameters. Chem. Rev. 91, 165–195. https://doi.org/10.1021/cr00002a004 (1991).

79. 79.

Humphrey, W., Dalke, A. & Schulten, K. VMD: Visual molecular dynamics. J. Mol. Graph. 14, 33–38. https://doi.org/10.1016/0263-7855(96)00018-5 (1996).

80. 80.

Legault, C. Y. CYLview, 1.0b. Université de Sherbrooke, 2009 (http://www.cylview.org)

## Acknowledgements

The authors extend their appreciation to the Researchers Supporting Project number (RSP-2021/201), King Saud University, Riyadh, Saudi Arabia. A.A.C.B. (Grants 2011/07895-8, 2015/01491-3, and 2014/25770-6) is thankful to Fundação de Amparo à Pesquisa do Estado de São Paulo for financial support. A.A.C.B. (Grant 309715/2017-2) also thanks the Brazilian National Research Council (CNPq) for financial support and fellowships. This study was also financed in part by the CAPES—Finance Code 001.

## Author information

Authors

### Contributions

M.I. and J.H. carried out synthesis of probes and their characterization. Z.S., A.A.-H. and H.H.S. designed the idea and wrote manuscript. F.M. and V.S. designed and studied the concept of molecular level keypad lock system. M.K., S.F.A.M., A.A.C.B. and M.U.K. did computational studies by DFT. N.M. and N.S.S.Z. carried out analytical applications and prepared/revised figures and graphs.

### Corresponding authors

Correspondence to Zahid Shafiq or Muhammad Khalid.

## Ethics declarations

### Competing interests

The authors declare no competing interests.

### Publisher's note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Rights and permissions

Reprints and Permissions

Islam, M., Shafiq, Z., Mabood, F. et al. 2-Nitro- and 4-fluorocinnamaldehyde based receptors as naked-eye chemosensors to potential molecular keypad lock. Sci Rep 11, 20847 (2021). https://doi.org/10.1038/s41598-021-99599-w