Schiff bases have many chemical and biological applications in medicine and pharmaceuticals due to the presence of an imine group (?C=N?). These bases are used in many different fields of technology, and in photochemistry because of their photochromic properties. Here, the structural and electronic properties of the Schiff base formed by tacrine and saccharin (TacSac) were explored using density functional theory with the B3LYP, M06-2X, M06L, and ωB97XD functionals in combination with the 6-311++G(d,p) basis set. The time-dependent formalism was used at the B3LYP/6-311++G(d,p) level to obtain electronic transitions. The calculations were repeated in an implicit solvent model mimicking water, using the polarizable continuum model in conjunction with a solvation model based on a density approach. The results indicate that TacSac cannot form spontaneously, but can be obtained in mild reactions. However, the resulting Schiff base displays different characteristics to its monomers. It also has the potential for use in photochemical intramolecular charge-transfer systems.
Quantum chemical computations were used for prediction of the structure and color of alizarin complex with alkali metal hydroxides in methanolic solutions. The color prediction relying on the single Gaussian-like band once again proved the usefulness of the PBE0 density functional due to the observed smallest color difference between computed and experimentally derived values. It was found that the alkali metal hydroxide molecules can bind to the two oxygen atoms of both hydroxyl groups of alizarin or to one of these atoms and the oxygen atom from the keto group in a complex with three methanol molecules. This means that two electronic transitions need to be taken into account when considering the spectra of the studied complexes. The resulting bond lengths and angles are correlated with the properties of the alkali metal atoms. The molar mass, the atomic radius, and the Pauling electronegativity of studied metals are quite accurate predictors of the geometric properties of hydroxide complexes with alizarin in methanol solution.
Graphical abstract The spectra of the neutral and monoanionic form of alizarin together with color changes resulting from addition of different metal hydroxides and represented in CIE color space
Herein we report a study of the switchable [3]rotaxane reported by Huang et al. (Appl Phys Lett 85(22):5391–5393, 1) that can be mounted to a surface to form a nanomechanical, linear, molecular motor. We demonstrate the application of semiempirical electronic structure theory to predict the average and instantaneous force generated by redox-induced ring shuttling. Detailed analysis of the geometric and electronic structure of the system reveals technical considerations essential to success of the approach. The force is found to be in the 100–200 pN range, consistent with published experimental estimates.
The structures and energies for the Huisgen 1,3-dipolar cycloaddition reactions of methyl and ethyl azides with some cyclooctynes and dibenzocyclooctynes were computed at the B3LYP/6-311++G(d,p) level. The activation strain model (ASM) and quantitative molecular orbital (MO) theory were used to investigate the reactivity and regiochemistry in these reactions. The energy decomposition analysis (EDA) was used to identify the intrinsic electronic factor that lead to the preferential formation of 1,7-regiochemistry products. The reactivity order agrees with formation of more synchronous transition states and lower distortion energies. For the reaction of N3Met with azacyclooctyne, the 1,7-regiochemistry preference is attributed to a lower FMO gap and a higher contribution of the polarization term of the interaction energy than for the 1,8-transition state. For the reaction with aza-dibenzocyclooctyne, the 1,7-preference is due to a lower strain energy and a more pronounced contribution of the exchange term of the interaction energy.
Graphical Abstract In the reactions between methyl and ethyl azides with azacyclooctynes the regiochemistry is governed by the intrinsic electronic factors
Unknown force-field parameters for metal organic beryllium complexes used in emitting and electron transporting layers of OLED structures are determined. These parameters can be used for the predictive atomistic simulations of the structure and properties of amorphous organic layers containing beryllium complexes. The parameters are found for the AMBER force field using a relaxed scan procedure and quantum-mechanical DFT calculations of potential energy curves for specific internal (angular) coordinates in a series of three Be complexes (Bebq2; Be(4-mpp)2; Bepp2). The obtained parameters are verified in calculations of some molecular and crystal structures available from either quantum-mechanical DFT calculations or experimental data.
The electronic structures and spectroscopic properties of two complexes [M(pic)3] (M = Ir, Rh) containing picolinate as bidentate ligands have been calculated by means density functional theory (DFT) and time-dependent DFT/TD-DFT using three hybrid functionals B3LYP, PBE0 and mPW1PW91. The PBE0 and mPW1PW91 functionals, which have the same HF exchange fraction (25%), give similar results and do not differ drastically from B3LYP results. Calculated geometric parameters of the complexes are in good agreement with the available experimental data. The UV absorptions observed in acetonitrile were assigned on the basis of singlet state transitions. The most intense band observed in the UV-C region corresponds to ligand-to-ligand charge transfer states (LLCT) in both complexes. The theoretical spectrum of the rhodium complex is characterized by a large degree of mixing between metal-to-ligand-charge-transfer (MLCT), LLCT and metal centered (MC) states in the UV-A region. The presence of low-lying excited states with MC character affects the absorption spectrum under spin-orbit coupling (SOC) effects and play important roles in the photochemical properties.
The factors that explain the competition between intramolecular NO linkage photoisomerization and NO photorelease in five ruthenium nitrosyl complexes were investigated. By applying DFT-based methods, it was possible to characterize the ground states and lowest triplet potential energy surfaces of these species, and to establish that both photoisomerization and photorelease processes can occur in the lowest triplet state of each species. This work highlights the crucial role of the sideways-bonded isomer, a metastable state also known as the MS2 isomer, in the photochemical loss of NO, while the results obtained also indicate that the population of the triplet state of this isomer is compulsory for both processes and show how photoisomerization and photorelease interfere.
Theoretical calculations for the first tri-iron-based extended metal atom chain (EMAC) molecule are reported. The studied triple-high-spin (S?=?6) complex exhibits ferromagnetic ordering (according to Ising and spin-projection approximations), which renders it unique among all previously prepared and theoretically calculated EMAC compounds. This ordering originates from the prevailing ferromagnetic nearest-neighbor interactions, while the magnetic superexchange between terminal Fe2+ sites is weaker and antiferromagnetic. Calculations indicate that this linear chain system based on a tri-iron core shows potential for the development of spin-frustrated behavior, which could be achieved through rational modification of the equatorial and axial ligands.
The conversion of 2-phenylbenzimidazole using o-phenylenediamine and benzaldehyde can be improved significantly under β-cyclodextrin (β-CD). The density functional theory (DFT) method was applied to study the whole process. According to energy parameters (binding energy, deformation energy) and structural deformation, entry models and the reaction process can be pinpointed, with o-phenylenediamine embedding β-CD from a wide rim, and then benzaldehyde passing into the inclusion from the narrow rim. Subsequently, natural bonding orbital (NBO), Mulliken charge, frontier orbital, FuKui function and nuclear magnetic resonance (NMR) methods were employed to reveal the mechanism of electron transfer. The results illustrate that β-CD plays a catalytic role in synthesis reaction mechanism on the secondary side, improving the reactivity and selectivity of the process.
Graphical Abstract Density functional theory study of the effects of β-cyclodextrin in synthesis of 2-phenylbenzimidazole via benzaldehyde and o-phenylenediamine
For the first time, the structures, stabilities and electronic properties of alkaline-earth metal doped B44 fullerenes were investigated by means of density functional theory calculations. Our results reveal that M@B44 (M = Ca, Sr, Ba) possess endohedral configurations as their lowest energy structures, whereas the exohedral form is favored when metal is Be or Mg. The large binding energies and sizable HOMO–LUMO gap energies of Ca@B44, Sr@B44 and Ba@B44 suggest the considerable possibility to achieve these novel endohedral borofullerenes experimentally. Born-Oppenheimer molecular dynamics (BO-MD) simulations at various temperatures further confirmed the extreme dynamic stabilities of these endohedral complexes. Their bonding patterns were also analyzed in detail. Finally, we simulated their infrared absorption spectra and 11B nuclear magnetic resonance spectra to help future structural characterization.
The main absorption peaks were obtained for 1-(1-hydroxy-4-methyl-2-phenylazo)-2-naphthol-4-sulfonic acid. Generalized gradient approximation, hybrid, semi-empirical, and Coulomb attenuating methods were utilized to compare theoretical electronic transitions and experimental absorption spectra at different pH. The main peaks and shoulders observed in experimental spectra were assigned to its correct conformer. In order to find the most populated conformer, thermal effects on stability calculations were investigated to obtain molar fractions of possible isomers present at room and higher temperature. Theoretical electronic transitions at distinct pH could be obtained varying the protonation a deprotonation degree. It was found that generalized gradient approximation performs very well the first transition peak at neutral pH. For higher pH, all methodologies got a bathochromic shift in agreement with experiment and finally, from these theoretical results, it was obtained that this azo dye is hardly protonated in experiments since results presented here, predict a variation of absorption spectra for all proposed methodologies when the molecule is protonated, which is different to experimental results.
Seven models that related the features of molecular surface electrostatic potentials (ESPs) above the bond midpoints and rings, statistical parameters of ESPs to the experimental impact sensitivities h50 of eight strained cyclic explosives with the C–NO2 bonds were theoretically predicted at the DFT-B3LYP/6-311++G** level. One of the models was used to investigate the changes of h50 for the nitrocyclohydrocarbon frameworks in the H-bonded complexes of HF with nitrocyclopropane, nitrocyclobutane, nitrocyclopentane, and nitrocyclohexane. The results show that the correlation coefficients of the obtained models are small. When adding the effect of ring strain, the value of correlation coefficient is increased. According to the calculated h50, the sensitivities in the frameworks are increased after hydrogen bonding. As a global feature of molecules, surface electrostatic potential is more available to judge the sensitivity change than the trigger bond dissociation energy or ring strain energy in H-bonded complex.
Graphical Abstract A theoretical prediction of the relationships between the impact sensitivity and electrostatic potential in strained cyclic explosive and application to H-bonded complex of nitrocyclohydrocarbon?
Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C?=?C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae. BLYP/6-311++G** calculated ν(C?=?C) and ν(C–C) frequencies for all-trans and all-cis polyenes containing up to 14 carbon–carbon double bonds were used to estimate these values for very long chains. Correction of raw, unscaled vibrational data was performed by comparing theoretical and experimental wavenumbers for polyenes chains containing 3 to 12 conjugated C?=?C units with both ends substituted by tert-butyl groups. The corrected ν(C?=?C) and ν(C–C) wavenumbers for all-trans molecules were used to estimate the presence of 9 – 12 C?=?C units in all-trans polyene pigment in red coral.
Graphical abstract Detailed density functional theory (DFT) calculations on the structure and harmonic frequencies of model all-trans and all-cis polyenes were undertaken. For the first time, we report on the convergence of selected B3LYP/6-311++G** and BLYP/6-311++G** calculated structural parameters resulting from a systematic increase in polyene size (chains containing 2 to 14 C=C units). The limiting values of the structural parameters for very long chains were estimated using simple three-parameter empirical formulae.
A theoretical 1H NMR spectroscopy and thermodynamic analysis of the host–guest inclusion process involving the norfloxacin (NFX) into β-cyclodextrin (β-CD) was carried out. DFT structure and stabilization energies were obtained in both gas and aqueous phases. We could establish that the complex formation is enthalpy driven, and the hydrogen bonds established between NFX and β-CD play a major role in the complex stabilization. Besides, a theoretical 1H NMR analysis has shown to be a supplementary proceeding to predict appropriately the inclusion mode of norfloxacin molecule into the β-CD. In this work, a theoretical study of the NFX@β-CD complex is reported for the first time, seeking a deep understanding of topology and thermodynamics of the inclusion complex formation.
Quercetin is a prototypical antioxidant and prominent member of flavonoids, a large group of natural polyphenols. The oxidation of quercetin may lead to its dimerization, which is a paradigm of the more general polyphenol oligomerization. There exist two opposing mechanisms to describe the dimerization process, namely radical-coupling or Diels-Alder reactions. This work presents a comprehensive rationalization of this dimerization process, acquired from density functional theory (DFT) calculations. It is found that the two-step radical-coupling pathway is thermodynamically and kinetically preferred over the Diels-Alder reaction. This is in agreement with the experimental results showing the formation of only one isomer, whereas the Diels-Alder mechanism would yield two isomers. The evolution in bonding, occurring during these two processes, is investigated using the atoms in molecules (AIM) and electron localization function (ELF) topological approaches. It is shown that some electron density is accumulated between the fragments in the transition state of the radical-coupling reaction, but not in the transition state of the Diels-Alder process.
Graphical Abstract Quantum chemistry calculations of the dimerization process of quercetin show that a radical coupling approach is preferred to a Diels-Alder type reaction, in agreement with experimental results. Analysis of the bonding evolution highlights the reaction mechanism.
The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in methanol, water, and chloroform in the molecularly imprinted recognition were investigated systematically. The M06-2X results revealed that: 1) the hydroxyl groups in polar solvents such as methanol and water may markedly influence the weak interactions, and then alter the adsorption and emission spectra; 2) the electronic excitation in absorption spectra of dichlorvos is dominated by the configuration HOMO?→?LUMO, but in the most stable dichlorvos-MAA it becomes the ππ* excitation of HOMO?→?LUMO?+?1; 3) Mulliken charges reveal that dichlorvos almost dissociates to Cl- and a cation in its S1 excitation state; 4) the phosphorescence spectra of dichlorvos-MAA are relatively weak.
Graphical Abstract The absorption and emission spectra of dichlorvos and the dichlorvos-MAA complex in the molecularly imprinted recognition of dichlorvos were investigated systematically in methanol, water, and chloroform as solvents.
Determination of electrophilic and nucleophilic sites of a molecule is the primary task to find the active sites of the lead molecule. In the present study, the active sites of busulfan have been predicted by molecular electrostatic potential surface and Fukui function analysis with the help of dispersion corrected density functional theory. Similarly, the identification of active binding sites of the proteins against lead compound plays a vital role in the field of drug discovery. Rigid and flexible molecular docking approaches are used for this purpose. For rigid docking, Hex 8.0.0 software employing fast Fourier transform (FFT) algorithm has been used. The partial flexible blind docking simulations have been performed with AutoDock 4.2 software; where a Lamarckian genetic algorithm is employed. The results showed that the most electrophilic atoms of busulfan bind with the targets. It is clear from the docking studies that busulfan has inhibition capability toward the targets 12CA and 1BZM.
Multilayer-shaped compression and slide models were employed to investigate the complex sensitive mechanisms of cocrystal explosives in response to external mechanical stimuli. Here, density functional theory (DFT) calculations implementing the generalized gradient approximation (GGA) of Perdew-Burke-Ernzerhof (PBE) with the Tkatchenko-Scheffler (TS) dispersion correction were applied to a series of cocrystal explosives: diacetone diperoxide (DADP)/1,3,5-trichloro-2,4,6-trinitrobenzene (TCTNB), DADP/1,3,5-tribromo-2,4,6-trinitrobenzene (TBTNB) and DADP/1,3,5-triiodo-2,4,6-trinitrobenzene (TITNB). The results show that the GGA-PBE-TS method is suitable for calculating these cocrystal systems. Compression and slide models illustrate well the sensitive mechanism of layer-shaped cocrystals of DADP/TCTNB and DADP/TITNB, in accordance with the results from electrostatic potentials and free space per molecule in cocrystal lattice analyses. DADP/TCTNB and DADP/TBTNB prefer sliding along a diagonal direction on the a?c face and generating strong intermolecular repulsions, compared to DADP/TITNB, which slides parallel to the b?c face. The impact sensitivity of DADP/TBTNB is predicted to be the same as that of DADP/TCTNB, and the impact sensitivity of DADP/TBTNB may be slightly more insensitive than that of DADP and much more sensitive than that of TBTNB.
