Where You Protonate Matters: Deciphering the Unimolecular Chemistry of Protonated Myrcene and Linalool.

The unimolecular reactions of protonated myrcene and linalool were investigated by collision-induced dissociation and density functional theory calculations. Experiments on a triple quadrupole mass spectrometer showed that protonated myrcene undergoes two major unimolecular reactions losing propene and isobutene, and two minor reactions of ethene and propane loss. In each case, the product ion consists of a substituted five-member ring. Protonation of myrcene was found to form four distinct protomers, three of which can be significantly populated in the ion source. The observed fragmentation reactions were calculated and found to depend on the starting protomer. Each pathway consisted of several hydrogen-migration and ring-forming/opening steps on the way to the observed products. Likewise, protonation of linalool also produces three distinct protomers, with the global minimum being formed by protonation of a central double bond. The major reaction is water loss to form protonated myrcene, but two minor channels were also observed resulting in loss of acetone and isobutene. The calculated minimum energy reaction pathways were found to be consistent with the relative abundances of the ions in the experimental breakdown diagrams.

Application of PLA/GO/ZnO and PLA/GO/Cu2O as sensor.

Polylactic acid modified with graphene oxide (PLA/GO) is proposed to interact with ZnO through 6 different schemes. Density functional theory at B3LYP/LANL2DZ level was utilized to calculate total dipole moment (TDM), HOMO/LUMO energy gap (ΔE) and to map the molecular electrostatic potential (MESP). Results indicated that PLA/GO interacted with ZnO through O-atom forming PLA/GO/OZn composite. This composite interacts with methane, hydrogen sulfide, humidity (H2O), carbon dioxide and ethanol. The same gases were supposed to interact further with PLA/GO/Cu2O. Adsorption energy for the interaction between each composite and the proposed gases were calculated. Both PLA/GO/OZn and PLA/GO/Cu2O composites interacted favorably with H2O. Adsorption energy for interaction of other gases with studied structures are generally low compared to H2O. PLA/GO/OZn have adsorption energy slightly higher than that of PLA/GO/Cu2O. PLA/GO/OZn has higher TDM values than those of PLA/GO/Cu2O, indicating a more polar material. Conversely, PLA/GO/Cu2O exhibited larger ΔE values than those of PLA/GO/OZn. TDM and energy gap results for both studied structures indicated good sensing capabilities. Further insights come from analyzing the calculated density of states (DOS) and partial density of states (PDOS). PLA/GO/Cu2O exhibited high peak for copper in its DOS and PDOS spectra compared to zinc and oxygen in case of PLA/GO/OZn. This means a higher density of available electronic states associated with Cu.

Molecular modeling analysis for functionalized graphene/sodium alginate composite.

This study examined the functionalization of graphene with easily ionizable elements, such as lithium, and subsequently its interaction with the biopolymer sodium alginate (SA), to highlight its potential for biomedical applications. Utilizing Density Functional Theory (DFT), the research comprehensively investigated the structural, electronic, and spectroscopic properties of these graphene-based composites. The electronic properties of functionalized graphene were investigated using DFT at the B3LYP/6-31G(d,p) level. Among the various configurations studied, graphene exhibited weak interaction with two lithium atoms, displaying the highest reactivity in terms of total dipole moment (TDM) at 5.967 Debye and a HOMO/LUMO energy gap (ΔE) of 0.748 eV. Electrostatic potential mapping revealed that graphene when enhanced with lithium and three units of SA, exhibited an augmented potential density on its surface, a finding corroborated by other investigated physical properties. Notably, the configuration of graphene/3SA/Li, with weak interaction occurring at two side carbons, demonstrated the highest reactivity with a TDM of 15.509 Debye and ΔE of 0.280 eV. Additionally, a shift in the spectral characteristics of graphene towards lower wavenumbers was observed as lithium and SA interacted with the graphene substrate. The PDOS plot for Graphene/3SA/Li, showed the highest contribution in the HOMO orbitals was equally from lithium, sodium, hydrogen, and oxygen, while the lowest contribution was from carbon. This computational analysis provides comprehensive insights into the functionalized graphene systems, aiding in their further development and optimization for practical biomedical use.

Mechanistic and kinetic insights of the formation of allene and propyne from the C3H3 reaction with water.

CONTEXT: Allene (H2C = C = CH2) and propyne (CH3-C≡CH) are important compounds in the combustion chemistry. They can be created from the reaction of proparyl radicals with water. In this study, therefore, a computational study into the C3H3 + H2O potential energy landscape has been carefully conducted. The computed results indicate that the reaction paths forming the products (allene: CH2CCH2 + •OH) and (propyne: HCCCH3 + •OH) prevail under the 300-2000 K temperature range, where the latter is much more predominant compared to the former. However, these two products are not easily formed under ambient conditions due to the high energy barriers. In the 300 - 2000 K temperature range, the branching ratio for the propyne + •OH product declines from 100 to 86%, whereas the allene + •OH product shows an increase, reaching 14% at 2000 K. The overall bimolecular rate constant of the title reaction can be presented by the modified Arrhenius expression of ktotal = 1.94 × 10-12 T0.14 exp[(-30.55 kcal.mol-1)/RT] cm3 molecule-1 s-1. The total rate constant at the ambient conditions in this work, 2.37 × 10-34 cm3 molecule-1 s-1, was found to be over five orders of magnitude lower than the total rate constant of the C3H3 + NH3 reaction, 7.98 × 10-29 cm3 molecule-1 s-1, calculated by Hue et al. (Int. J. Chem. Kinet. 2020, 4(2), 84-91). The results in this study contribute to elucidating the mechanism of allene and propylene formation from the C3H3 + H2O reaction, and they can be used for modeling C3H3-related systems under atmospheric and combustion conditions. METHODS: All the geometric structures of the C3H3 + H2O system were optimized by the B3LYP method in conjunction with the 6-311 + + G(3df,2p) basis set. Single-point energies of these species were calculated at the CCSD(T)/6-311 + + G(3df,2p) level of theory. The CCSD(T)/CBS level has also been used to compute single-point energies for the two major reaction channels (C3H3 + H2O → allene + •OH and C3H3 + H2O → propyne + •OH). Rate constants and branching ratios of the key reaction channels were calculated in the 300-2000 K temperature interval using the Chemrate software based on the transition state theory (TST) with Eckart tunneling corrections.

New Insights in the Computational pKb Determination of Primary Amines and Anilines.

Extensive research has already provided reliable methods for the in silico prediction of pKa, while a trustworthy strategy for pKb determination is still being sought. Indeed, the approaches previously exploited for computing pKa have shown their weakness in predicting pKb. In the light of the exceptional reliability demonstrated in the pKa calculation of a wide panel of organic acids, in this work, we exploited our "easy to use methodology", based on the direct approach, to predict the pKb of primary amines. Herein, CAM-B3LYP was compared to WB97XD and B3PW91, exploring the solvation model based on density (SMD) and the polarizable continuum model (PCM), in the presence of two explicit water molecules. Noteworthy, CAM-B3LYP and WB97XD returned completely different solvent accessible surfaces (SAS) and electron potential maps (EPM) for the bases and the conjugated acids, independently from the nature of the substituents. Once again, CAM-B3LYP/SMD/2H2O method confirmed its remarkable reliability, leading to a minimum average error (MAE) lower than 0.3. This outstanding result strengthens the trustworthiness of our method, already successfully applied to predict the pKa of different substituted phenols and carboxylic acids. Thus, our "easy-to-use" process can predict also the pKb of primary ammines and anilines, always ensuring consistent outputs.

Development of cellulose nanowhisker-gallic acid antioxidant bioconjugate via covalent conjugation and supramolecular interactions: A comparative study.

A new supramolecular antioxidant bioconjugate based on cellulose nanowhisker (CNW) and gallic acid (GA) was developed by grafting ß-CD on the surface of CNW and then employing host- guest chemistry to involve GA. Our challenge was to explore the effect of supramolecular conjugation of antioxidant molecules versus their covalent binding on the CNW backbone on the antioxidant activity. The synthesis of these products was confirmed using Fourier transform infrared (FT-IR) and differential scanning calorimetry (DSC) analyses. The antioxidant activity of gallic acid (GA) containing products, both products including its non-covalent interactions with CNW-g-ß-CD and covalent bonding with CNW were experimentally evaluated using DPPH test. Theoretical calculations using Gaussian software and the density functional theory (DFT) method were also performed. The results showed that GA's antioxidant activity increased in non-covalent conjugated form. Hydrogen atom transfer (HAT) was used to predict the antioxidant activity of GA in computational methods. These findings not only expand our understanding of the structure-activity relationships in antioxidant systems but also provide valuable insights that can aid in the design and development of novel biopolymer-based antioxidants with improved properties.

Towards the "Eldorado" of pKa Determination: A Reliable and Rapid DFT Model.

The selection of a "perfect tool" for the theoretical determination of acid-base dissociation constants (Ka) is still puzzling. Recently, we developed a user-friendly model exploiting CAM-B3LYP for determining pKa with impressive reliability. Herein, a new challenge is faced, examining a panel of functionals belonging to different rungs of the "Jacob's ladder" organization, which classifies functionals according to their level of theory. Specifically, meta-generalized gradient approximations (GGAs), hybrid-GGAs, and the more complex range-separated hybrid (RSH)-GGAs were investigated in predicting the pKa of differently substituted carboxylic acids. Therefore, CAM-B3LYP, WB97XD, B3PW91, PBE1PBE, PBEPBE and TPSSTPSS were used, with 6-311G+(d,p) as the basis set and the solvation model based on density (SMD). CAM-B3LYP showed the lowest mean absolute error value (MAE = 0.23) with relatively high processing time. PBE1PBE and B3PW91 provided satisfactory predictions (MAE = 0.34 and 0.38, respectively) with moderate computational time cost, while PBEPBE, TPSSTPSS and WB97XD led to unreliable results (MAE > 1). These findings validate the reliability of our model in predicting carboxylic acids pKa, with MAE well below 0.5 units, using a simplistic theoretical level and a low-cost computational approach.

DFT Studies of Dimethylaminophenyl-Substituted Phthalocyanine and Its Silver Complexes.

The dimethylaminophenyl-substituted silver phthalocyanine [dmaphPcAg] can be used as a UV-vis photoinitiator for in situ preparation of a silver/polymer nanocomposite. To verify early steps of the supposed mechanism of radical polymerization, we performed quantum chemical calculations of m[dmaphPcAg]q complexes with charges q = +1 to -2 in the two lowest spin states m, of a free ligand and its dehydrogenated/deprotonated products m[dmaphPcHn]q, n = 2 to 0, q = 0, -1 or -2, in the lowest spin states m. The calculated electronic structures and electron transitions of all the optimized structures in CHCl3 solutions are compared with experimental EPR and UV-vis spectra, respectively. The unstable 3[dmaphPcAg]+ species deduced only from previous EPR spin trap experiments was identified. In addition to 2[dmaphPcAg]0, our results suggest the coexistence of both reaction intermediates 1[dmaphPcAg]- and 3[dmaphPcAg]- in reaction solutions. Silver nanoparticle formation is a weak point of the supposed reaction mechanism from the energetic, stereochemistry, and electronic structure points of view.

A computational study of H-bonded networks in cyclic water clusters, (H2O)n (n = 3-12).

CONTEXT: We have performed a detailed MM and DFT investigation of neutral water clusters (H2O)n (n = 3-12). Our results show the trend of interaction energies in these clusters as a function of the size of the cluster. They show that the H-bond strength increases with cluster size and that the model of water is better described if two different partial charges are used on the hydrogen, depending on whether hydrogen is H-bonded or not. The average binding enthalpy change due to the formation of H-bonds between water molecules is found to be - 25.9 kJ mol-1 at B3LYP/aug-cc-pVDZ level of theory. We observe the formation of cyclic H-bonded networks through the analysis of frontier orbitals and IR vibrational frequencies spectra. For the water cluster with n = 11, we observe an unusual reduction of the bandgap indicative of a cyclic H-bonded network. METHODS: Calculations were performed with the MMFF94 force field and the B3LYP method using various large basis sets. Molecular orbital diagrams and population analysis were done using standard tools in Gaussian.

Conformational preferences of guanine-containing threose nucleic acid building blocks in B3LYP studies.

In this paper, detailed and systematic gas-phase B3LYP conformational studies of four monomers of threose nucleic acid (TNA) with guanine attached at the C1' atom and bearing different substituents (OH, OP(=O)OH2 and OCH3) in the C2' and C3' positions of the α-l-threofuranose moiety are presented. All exocyclic single-bond (χ, ε and γ) rotations, as well as the ν0-ν4 endocyclic torsion angles, were taken into consideration. Three (threoguanosines TG1-TG3) or two (TG4) energy minima were found for the rotation about the χ torsion angle. The syn orientation (the A rotamer family) is strongly privileged in geometries TG1 and TG2, whereas the anti orientation (the C rotamer family) and the syn orientation are observed to be in equilibrium (with populations of 56% and 44%, respectively) for TG3. In the case of TG4, the high-anti orientation (the B rotamer family) turned out to be by far the most favourable, with the contribution exceeding 90% in equilibrium. Such a preference can be attributed to the inability of H-bonding between sugar and nucleobase and possibly because of the steric strains. The low-energy conformers of TG1-TG4 occupy the northeastern (P âˆ¼ 40°) and/or southern (P âˆ¼ 210°) parts of the pseudorotational wheel, which fits the A- and B-type DNA helices quite well. Additionally, in the case of TG4, some relatively stable geometries have the furanoid ring in conformation lying on the northwestern part of the pseudorotational wheel (P âˆ¼ 288°).

Easy to Use DFT Approach for Computational pKa Determination of Carboxylic Acids.

In pKa computational determination, the challenge in exploring and fostering new methodologies and approaches goes in parallel with the amelioration of computational performances. In this paper a "ready to use methodology" has been compared to other strategies, such as the re-shaping in solvation cavity (Bondi radius re-shaping), wanting to assess its reliability in predicting the pKa of a broad list of carboxylic acids. Thus, the functionals B3LYP and CAM-B3LYP have been selected, using SMD as continuum solvation model. Exploiting our previous results, two water molecules were made explicit on the reaction centre. Data show that our model (CAM-B3LYP/2H2 O) is capable to accurately predict pKa, leading to mean absolute error (MAE) values lower than 0.5. Noteworthy, good results were achieved in computing the pKa of substituents bearing nitro and cyano groups. Focusing on B3LYP, eventually remarkable outputs were obtained only when Bondi correction was applied to the complex with two water molecules. Hence, massive outcomes were obtained in foreseeing the trichloro and trifluoro acetic acid pKa. These findings demonstrated that no complex level of theory nor external factor is required to accurately predict carboxylic acids pKa, with MAE well below 0.5 units.

Band gap, Jahn-Teller deformation, octahedra rotation in transition metal perovskites LaTiO 3 .

The LaTiO 3 perovskite (where Ti is in a d1 state) is investigated by using an all electron Gaussian basis and many functionals, ranging from pure GGA (PBE), to hybrids (full range, B3LYP and PBE0, and range separated, HSE06) to Hartree Fock. Recently, Varignon et al. (Phys. Rev. Res 1, 033131, 2019), showed that, when GGA+U or HSE06 are used, a metallic solution and fractional occupancy of the t 2 g subshell are obtained. Here, it is shown that when a full range hybrid functional is used, an integer occupancy is obtained, as suggested by the Jahn-Teller theorem. When the exact exchange percentage varies from 0 to 100, the system is insulating when it exceeds 20. By reducing progressively the symmetry from cubic down to orthorhombic, the relative importance of the Jahn-Teller deformation and of the rotation of the octahedra is explored.

Review of Systematic Tendencies in (001), (011) and (111) Surfaces Using B3PW as Well as B3LYP Computations of BaTiO3, CaTiO3, PbTiO3, SrTiO3, BaZrO3, CaZrO3, PbZrO3 and SrZrO3 Perovskites.

We performed B3PW and B3LYP computations for BaTiO3 (BTO), CaTiO3 (CTO), PbTiO3 (PTO), SrTiO3 (STO), BaZrO3 (BZO), CaZrO3 (CZO), PbZrO3 (PZO) and SrZrO3 (SZO) perovskite neutral (001) along with polar (011) as well as (111) surfaces. For the neutral AO- as well as BO2-terminated (001) surfaces, in most cases, all upper-layer atoms relax inwards, although the second-layer atoms shift outwards. On the (001) BO2-terminated surface, the second-layer metal atoms, as a rule, exhibit larger atomic relaxations than the second-layer O atoms. For most ABO3 perovskites, the (001) surface rumpling s is bigger for the AO- than BO2-terminated surfaces. In contrast, the surface energies, for both (001) terminations, are practically identical. Conversely, different (011) surface terminations exhibit quite different surface energies for the O-terminated, A-terminated and BO-terminated surfaces. Our computed ABO3 perovskite (111) surface energies are always significantly larger than the neutral (001) as well as polar (011) surface energies. Our computed ABO3 perovskite bulk B-O chemical bond covalency increases near their neutral (001) and especially polar (011) surfaces.

DFT Investigation of the Structural, Electronic, and Optical Properties of AsTi (Bi)-Phase ZnO under Pressure for Optoelectronic Applications.

Pressure-induced phases of ZnO have attracted considerable attention owing to their excellent electronic and optical properties. This study provides a vital insight into the electronic structure, optical characteristics, and structural properties of the AsTi (Bi) phase of ZnO under high pressure via the DFT-based first-principles approach. The phase transformation from BN(Bk) to the Bi phase of ZnO is estimated at 16.1 GPa using local density approximation, whereas the properties are explored precisely by the hybrid functional B3LYP. The electronic structure exploration confirms that the Bi phase is an insulator with a wider direct bandgap, which expands by increasing pressure. The dielectric function evidenced that the Bi phase behaves as a dielectric in the visible region and a metallic material at 18 eV. Optical features such as the refractive index and loss function revealed the transparent nature of the Bi phase in the UV range. Moreover, the considered Bi phase is found to possess a high absorption coefficient in the ultraviolet region. This research provides strong theoretical support for the development of Bi-phase ZnO-based optoelectronic and photovoltaic devices.

Effect of charge transfer on the first hyper-polarizability of N,N-dimethylaniline and julolidine: a DFT based comparative study.

CONTEXT: Quantum mechanical calculations involving electron correlation, frequency dispersion, and solvent effects were carried out to examine the second-order nonlinear optical response of various acceptor, X (-CF3, -CN, -NO2) substituted in N,N-dimethylaniline (DMA) and julolidine(JLD). Here, both DMA and JLD acts as donor and the three substituted groups, X (-CF3, -CN and -NO2) at the para position of both the ring systems as acceptor. The NLO response (ßHRS) of -CF3 and -CN substituted DMA and JLD is relatively lower compared to DMA-NO2 and JLD-NO2. The charge distribution is found higher in case of -NO2 substituted DMA and JLD (±443 and ±449) compared to their -CF3 or -CN substitution. Electronic characteristics such as UV-Vis absorption spectra, crucial excited state parameters and charge transfer contribution to ßHRS have been used to explain the NLO parameter of DMA-X and JLD-X. Variation of the incident optical frequency of light shows fluctuation of ßHRS value and highest values of ßHRS are obtain at the λmax frequency of each compound. Solvent polarity variation study on ßHRS shows that ßHRS varies linearly with the Kirkwood-Onsagar dielectric factor (D). METHODS: All computational studies have been carried out using density functional theory (DFT) based method. Since CAM-B3LYP based hybrid functional improves the asymptotic behavior of the exchange interaction by dividing into short-range and long-range components, first hyperpolarizability values in the present study were computed using DFT/ CAM-B3LYP/ 6-31G+(d,p) level of theory.

Extractive Spectrophotometric Determination and Theoretical Investigations of Two New Vanadium(V) Complexes.

Two new vanadium (V) complexes involving 6-hexyl-4-(2-thiazolylazo)resorcinol (HTAR) and tetrazolium cation were studied. The following commercially available tetrazolium salts were used as the cation source: tetrazolium red (2,3,5-triphenyltetrazol-2-ium;chloride, TTC) and neotetrazolium chloride (2-[4-[4-(3,5-diphenyltetrazol-2-ium-2-yl)phenyl]phenyl]-3,5-diphenyltetrazol-2-ium;dichloride, NTC). The cations (abbreviated as TT+ and NTC+) impart high hydrophobicity to the ternary complexes, allowing vanadium to be easily extracted and preconcentrated in one step. The complexes have different stoichiometry. The V(V)-HTAR-TTC complex dimerizes in the organic phase (chloroform) and can be represented by the formula [(TT+)[VO2(HTAR)]]2. The other complex is monomeric (NTC+)[VO2(HTAR)]. The cation has a +1 charge because one of the two chloride ions remains undissociated: NTC+ = (NT2+Cl-)+. The ground-state equilibrium geometries of the constituent cations and final complexes were optimized at the B3LYP and HF levels of theory. The dimer [(TT+)[VO2(HTAR)]]2 is more suitable for practical applications due to its better extraction characteristics and wider pH interval of formation and extraction. It was used for cheap and reliable extraction-spectrophotometric determination of V(V) traces in real samples. The absorption maximum, molar absorptivity coefficient, limit of detection, and linear working range were 549 nm, 5.2 × 104 L mol-1 cm-1, 4.6 ng mL-1, and 0.015-2.0 µg mL-1, respectively.

Electronic, optical and spectroscopic properties of N-dialkyl-imidazolium hexafluorophosphate (CNMIM.PF6) ionic liquid crystal molecules investigated by computational methods.

In this present work, we calculate the electronic, spectroscopic and nonlinear optical properties (NLO) of N-dialkyl-imidazolium hexafluorophosphate (CNMIM.PF6, where N = 10, 12, 14, 16, 18, 20) ionic liquid crystal molecules under the effect of alkyl chain length variation in cation moiety [CNMIM]+ with fixed anion [PF6]-. CONTEXT: The majority of research on ionic liquid crystal to date has been focused on experiments, while theoretical studies on the optical properties of ionic liquid crystal have been extremely rare. Nonlinear phenomena in optical devices have attracted many researchers. Therefore, results of NLO properties may favor facile synthesis and fabrication of novel-type of materials as well as optoelectronic devices. Spectroscopic studies elucidate further insight into ionic liquid crystal behavior. The results demonstrate that variations in alkyl chain length have an impact on the conformers' electrical, spectroscopic, and NLO properties as well as their stability. The stability of ionic liquid crystal molecules increases with increase in the alkyl chain length and the energy band gap range is 6.64-6.29 eV. Understanding ionic liquid crystal's physical behavior requires an understanding of their dipole moments and NLO features, which are covered in this article. The results of NLO characteristics for all ionic liquid crystal molecules show that their first-order hyperpolarizabilities are higher than the reference molecule (urea). METHODS: The electronic (molecular energy band gap, electrostatic potential map, as well as HOMO-LUMO orbitals) and spectroscopic (IR-RAMAN, UV) properties were evaluated with the help of theoretical model at B3LYP/6-31G(d) while the NLO study has been performed using B3LYP and M06-2X with different basis sets 6-31G(d) and 6-311++G(d,p), as implemented in Gaussian09 software.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory.

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol-1), followed by B3LYP (0.69 kcal mol-1) and the HF theories (0.81-1.0 kcal mol-1). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol-1), ΜΜ3-96 (1.40 kcal mol-1) and the Halgren's MMFF94 force field (1.30 kcal mol-1) and then for the MM2-91 clones MMX (1.77 kcal mol-1) and MM+ (2.01 kcal mol-1) and MM4 (2.05 kcal mol-1). The DREIDING (3.63 kcal mol-1) and UFF (3.77 kcal mol-1) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization.

Tetra-tert-butyl-s-indacene is a Bond-Localized C2h Structure and a Challenge for Computational Chemistry.

Whether tetra-tert-butyl-s-indacene is a symmetric D2h structure or a bond-alternating C2h structure remains a standing puzzle. Close agreement between experimental and computed proton chemical shifts based on minima structures optimized at the M06-2X, ωB97X-D, and M11 levels confirm a bond-localized C2h symmetry, which is consistent with the expected strong antiaromaticity of TtB-s-indacene.

Tailoring the band-gap, optical and fluorescent properties of 3,5-diaminobenzoic acid via functionalization with chemical groups: A DFT study.

3,5-diaminobenzoic acid (3,5-DABA) with chemical formula C7H8N2O2 was functionalized with CH3-, OH-, NH2- and NO2- to obtain: CH3-3,5 DABA, OH-3,5 DABA, NH2-3,5DABA and NO2-3,5DABA. These molecules were built with Gauss view 6.0 and their structural, spectroscopic, optoelectronic and molecular properties were investigated using density functional theory (DFT). B3LYP (Becke's 3-parameter exchange functional with Lee-Yang-Parr correlation energy) functional and 6-311+ G (d, p) basis set were used to understand their reactivity, stability and optical activity. Integral equation formalism polarizable continuum model (IEF - PCM) was used to calculate the absorption wavelength, energy required to excite the molecules and oscillator strength. Our results reveal that the functionalization of 3,5 DABA with the groups caused a decrease of the energy gap from 0.1563 eV, to 0.1461 eV, 0.13818 eV and 0.13811 eV in NO2-3,5DABA, OH-3,5DABA and NH2-3,5DABA respectively. The lowest energy gap of 0.13811 eV for NH2-3,5DABA is in good agreement with its highest reactivity value (global softness of 7.240). The most observed significant donor - acceptor NBO interactions where found to occur between *ΠC16-O17 → *ΠC1-C2, *ΠC3-C4→ *ΠC1-C2, *ΠC1-C2 → *ΠC5-C6, *ΠC3-C4 → *ΠC5-C6, *ΠC2-C3 →*ΠC4-C5 natural bond orbitals having second- order stabilization energies of 101.95 kcal/mol, 368.41 kcal/mol, 174.51 kcal/mol, 255.63 kcal/mol and 235.92 kcal/mol in 3,5-DABA, CH3-3,5-DABA, OH-3,5-DABA, NH2-3,5-DABA and NO2-3,5-DABA respectively. The highest perturbation energy was observed in CH3-3,5DABA while the lowest perturbation energy was observed in 3,5DABA. The absorption band of the compounds were observed in the order: NH2-3,5DABA (404 nm) > N02-3,5DABA (393 nm) > OH-3,5DABA (386 nm) > 3,5DABA (349 nm) > CH3-3,5DABA (347 nm).