Theoretical investigation on switchable second-order nonlinear optical (NLO) properties of novel cyclopentadienylcobalt linear [4]phenylene complexes

As a kind of novel organometallic complexes, the cyclopentadienylcobalt (CpCo) linear [4]phenylene complexes (4 = number of benzene rings) display efficient switchable nonlinear optical (NLO) response when CpCo reversibly migrates along the linear [4]phenylene triggered by heating or lighting. In this paper, the second-order NLO properties for CpCo linear [4]phenylene complexes were calculated by using the density functional theory (DFT) methods with four functionals. All of the functionals yield the same order of β _tot values: 1<2<4<3. The effect of solvent on second-order NLO properties has been studied using polarized continuum model (PCM) in the tetrahydrofuran (THF) solution. The solvent leads to a slight enhancement of the NLO responses for the studied complexes relevant to their NLO responses in vacuo. The electronic absorption spectra were investigated by the TDDFT methods. The TDDFT calculations indicate that the maximum absorption peaks of complexes 2–4 in the near-infrared spectrum area show the bathochromic shift together with a decreasing intensity compared to complex 1. We have also found that the cobalt (Co) atom acts as a donor in all the organometallic complexes and the d → π* and π → π* charge transfer (CT) transitions contribute to the enhancement of second-order NLO response. Furthermore, two experimentally existing complexes 1 and 3 are found to have a large difference in β _tot values. It is our expectation that this difference may stimulate the search for a new type of switchable NLO material based on CpCo linear [4]phenylene complexes.

Assessing the quantum mechanical level of theory for prediction of linear and nonlinear optical properties of push-pull organic molecules

In this paper, we assessed the quantum mechanical level of theory for prediction of linear and nonlinear optical (NLO) properties of push-pull organic molecules. The electric dipole moment (μ), mean polarizability (〈α〉) and total static first hyperpolarizability (β_t) were calculated for a set of benzene, styrene, biphenyl and stilbene derivatives using HF, MP2 and DFT (31 different functionals) levels and over 71 distinct basis sets. In addition, we propose two new basis sets, NLO-V and aNLO-V, for NLO properties calculations. As the main outcomes it is shown that long-range corrected DFT functionals such as M062X, ωB97, cam-B3LYP, LC-BLYP and LC-ωPBE work satisfactorily for NLO properties when appropriate basis sets such as those proposed here (NLO-V or aNLO-V) are used. For most molecules with β ranging from 0 to 190 esu, the average absolute deviation was 13.2 esu for NLO-V basis sets, compared to 27.2 esu for the standard 6-31 G(2d) basis set. Therefore, we conclude that the new basis sets proposed here (NLO-V and aNLO-V), together with the cam-B3LYP functional, make an affordable calculation scheme to predict NLO properties of large organic molecules.

CO2 adsorption on polar surfaces of ZnO

Physical and chemical adsorption of CO₂ on ZnO surfaces were studied by means of two different implementations of periodic density functional theory. Adsorption energies were computed and compared to values in the literature. In particular, it was found that the calculated equilibrium structure and internuclear distances are in agreement with previous work. CO₂ adsorption was analyzed by inspection of the density of states and electron localization function. Valence bands, band gap and final states of adsorbed CO₂ were investigated and the effect of atomic displacements analyzed. The partial density of states (PDOS) of chemical adsorption of CO₂ on the ZnO(0001) surface show that the p orbitals of CO₂ were mixed with the ZnO valence band state appearing at the top of the valence band and in regions of low-energy conduction band.

Prodrugs of fumarate esters for the treatment of psoriasis and multiple sclerosis—a computational approach

Density functional theory (DFT) calculations at B3LYP/6-31 G (d,p) and B3LYP/6-311?+?G(d,p) levels for the substituted pyridine-catalyzed isomerization of monomethyl maleate revealed that isomerization proceeds via four steps, with the rate-limiting step being proton transfer from the substituted pyridinium ion to the C=C double bond in INT1. In addition, it was found that the isomerization rate (maleate to fumarate) is solvent dependent. Polar solvents, such as water, tend to accelerate the isomerization rate, whereas apolar solvents, such as chloroform, act to slow down the reaction. A linear correlation was obtained between the isomerization activation energy and the dielectric constant of the solvent. Furthermore, linearity was achieved when the activation energy was plotted against the pK _a value of the catalyst. Substituted-pyridine derivatives with high pK _a values were able to catalyze isomerization more efficiently than those with low pK _a values. The calculated relative rates for prodrugs 1–6 were: 1 (406.7), 2 (7.6?×?10⁶), 3 (1.0), 4 (20.7), 5 (13.5) and 6 (2.2?×?10³). This result indicates that isomerizations of prodrugs 1 and 3–5 are expected to be slow and that of prodrugs 2 and 6 are expected to be relatively fast. Hence, prodrugs 2 and 3–5 have the potential to be utilized as prodrugs for the slow release of monomethylfumarate in the treatment of psoriasis and multiple sclerosis.

Theoretical study of geometrical and nonlinear optical properties of pyridinum N-phenolate betaine dyes

This paper presents an ab initio quantum chemical investigation of the geometrical structures and the non-linear optical properties (NLO) of three structural isomers of pyridinium N-phenolate betaine dye. The ground state geometrical parameters and the first-order hyperpolarizabilities were calculated using the Hartree-Fock (HF) as well as the second-order perturbation Møller-Pleset (MP2) method with the 6–31G, 6–31G(d), 6–31G(d,p), 6–31+G(d), 6–31++G(d,p), 6–311+G(d), aug-cc-PVDZ and the recently developed Z3PolX basis sets. Moreover, the first-order hyperpolarizability was calculated at the coupled cluster singles and doubles (CCSD/6–31+G(d)) level of theory. The analysis of the results of calculations for the investigated isomers indicates that there are important differences in their NLO activities. Additionally, it was shown that Z3PolX basis set works reasonable well for betaine dyes.

B30H8, B39H9 2−, B42H10, B48H10, and B72H12: polycyclic aromatic snub hydroboron clusters analogous to polycyclic aromatic hydrocarbons

Calculations performed at the ab initio level using the recently reported planar concentric π-aromatic B₁₈H₆ ²⁺(1) [Chen Q et al. (2011) Phys Chem Chem Phys 13:20620] as a building block suggest the possible existence of a new class of B_3n H _m polycyclic aromatic hydroboron (PAHB) clusters—B₃₀H₈(2), B₃₉H₉ ^2?(3), B₄₂H₁₀(4/5), B₄₈H₁₀(6), and B₇₂H₁₂(7)—which appear to be the inorganic analogs of the corresponding C _n H _m polycyclic aromatic hydrocarbon (PAHC) molecules naphthalene C₁₀H₈, phenalenyl anion C₁₃H₉ ^?, phenanthrene/anthracene C₁₄H₁₀, pyrene C₁₆H₁₀, and coronene C₂₄H₁₂, respectively, in a universal atomic ratio of B:C?=?3:1. Detailed canonical molecular orbital (CMO), adaptive natural density partitioning (AdNDP), and electron localization function (ELF) analyses indicate that, as they are hydrogenated fragments of a boron snub sheet [Zope RR, Baruah T (2010) Chem Phys Lett 501:193], these PAHB clusters are aromatic in nature, and exhibit the formation of islands of both σ- and π-aromaticity. The predicted ionization potentials of PAHB neutrals and electron detachment energies of small PAHB monoanions should permit them to be characterized experimentally in the future. The results obtained in this work expand the domain of planar boron-based clusters to a region well beyond B₂₀, and experimental syntheses of these snub B_3n H _m clusters through partial hydrogenation of the corresponding bare B_3n may open up a new area of boron chemistry parallel to that of PAHCs in carbon chemistry.

Computational and experimental studies of the electronic excitation spectra of EDTA and DTPA substituted tetraphenylporphyrins and their Lu complexes

Ethylendiaminetetraacetic acid (EDTA) substituted and diethylenetriaminopentaacetic acid (DTPA) substituted aminated free-base tetraphenylporphyrins (H₂ATPP) and the corresponding lutetium(III) complexes have been studied computationally at the density functional theory (DFT) and second-order algebraic diagrammatic construction (ADC(2)) levels using triple-ξ basis sets augmented with polarization functions. The molecular structures were optimized using Becke's three-parameter hybrid functional (B3LYP). The electronic excitation spectra in the range of 400–700 nm were calculated using the ADC(2) and the linear-response time-dependent DFT methods. The calculated spectra are compared to those measured in ethanol solution. The calculated excitation energies agree well with those deduced from the experimental spectra. The excitation energies for the Q_x band calculated at the B3LYP and ADC(2) level are 0.20-0.25 eV larger than the experimental values. The excitation energies for the Q_y band calculated at the B3LYP level are 0.10-0.20 eV smaller than the ADC(2) ones and are thus in good agreement with experiment. The calculated excitation energies corresponding to the B_x and B_y bands are 0.10-0.30 eV larger than the experimental values. The excitation energies of the B_x and B_y bands calculated at the B3LYP level are in somewhat better agreement with experiment than the ADC(2) ones. The calculated and measured band strengths largely agree.

Theoretical study on cooperative effects between X⋯N and X⋯Carbene halogen bonds (X = F,Cl,Br and I)

Quantum chemical calculations are performed to study the interplay between halogen?nitrogen and halogen?carbene interactions in NCX?NCX?CH₂ complexes, where X?=?F, Cl, Br and I. Molecular geometries and interaction energies of dyads and triads are investigated at the MP2/aug-cc-pVTZ level of theory. It is found that the X?N and X?C_carbene interaction energies in the triads are larger than those in the dyads, indicating that both the halogen bonding interactions are enhanced. The estimated values of cooperative energy E _coop are all negative with much larger E _coop in absolute value for the systems including iodine. The nature of halogen bond interactions of the complexes is analyzed using parameters derived from the quantum theory atoms in molecules methodology and energy decomposition analysis.

Design of Lewis acid–base complex: enhancing the stability and first hyperpolarizability of large excess electron compound

In the present paper, a new type of Lewis acid–base complex BX₃???Li@Calix[4]pyrrole (X = H and F) was designed and assembled based on electride molecule Li@calix[4]pyrrole (as a Lewis base) and the electron deficient molecule BX₃ (as a Lewis acid) by employing quantum mechanical calculation. The new Lewis acid–base complex offers an interesting push-excess electron-pull (P-e-P) framework to enhance the stability and nonlinear optical (NLO) response. To measure the nonlinear optical response, static first hyperpolarizabilities (β ₀) are exhibited. Significantly, point-face assembled Lewis acid–base complex BF₃???Li@Calix[4]pyrrole (II) has considerable first hyperpolarizabilities (β ₀) value (1.4?×?10⁶ a.u.), which is about 117 times larger than reported 11,721 a.u. of electride Li@Calix[4]pyrrole. Further investigations show that, in BX₃???Li@Calix[4]pyrrole with P-e-P framework, a strong charge-transfer transition from the ground state to the excited state contributes to the enhancement of first hyperpolarizability. Theory calculation of enthalpies of reaction (Δ_rH⁰) at 298 K demonstrates that it is feasible to synthetize the complexes BX₃???Li@Calix[4]pyrrole. In addition, compared with Li@Calix[4]pyrrole, the vertical ionization potential (VIP) and HOMO–LUMO gap of BX₃???Li@Calix[4]pyrrole have obviously increased, due to the introduction of the Lewis acid molecule BX₃. The novel Lewis acid–base NLO complex possesses not only a large nonlinear optical response but also higher stability.

Density functional studies on photophysical properties and chemical reactivities of the triarylboranes: effect of the constraint of planarity

The geometric and electronic structures, absorption spectra, transporting properties, chemical reactivity indices and electrostatic potentials of the planar three-coordinate organoboron compounds 1-2 and twisted reference compound Mes ₃ B, have been investigated by employing density functional theory (DFT) and conceptual DFT methods to shed light on the planarity effects on the photophysical properties and the chemical reactivity. The results show that the planar compounds 1-2 exhibit significantly lower HOMO level than Mes ₃ B, owing to the stronger electronic induction effect of boron centers. This feature conspicuously induces a blue shifted absorption for 1, although 1 seemingly possesses more extended conjugation framework than Mes ₃ B. Importantly, the reactivity strength of the boron atoms in 1-2 is much lower than that in Mes ₃ B, despite the fact that the tri-coordinate boron centers of 1-2 are completely naked. The interesting and abnormal phenomenon is caused by the strong p-π electronic interactions, that is, the empty p-orbital of boron center is partly filled by π-electron of the neighbor carbon atoms in 1-2, which are confirmed by the analysis of Laplacian of the electron density and natural bond orbitals. Furthermore, the negative electrostatic potentials of the boron centers in 1-2 also interpret that they are not the most preferred sites for incoming nucleophiles. Moreover, it is also found that the planar compounds 1-2 can act as promising electron transporting materials since the internal reorganization energies for electron are really small.

Electron correlation effects and density analysis of the first-order hyperpolarizability of neutral guanine tautomers

Dipole moments (μ), charge distributions, and static electronic first-order hyperpolarizabilities (β _μ ) of the two lowest-energy keto tautomers of guanine (7H and 9H) were determined in the gas phase using Hartree–Fock, Møller–Plesset perturbation theory (MP2 and MP4), and DFT (PBE1PBE, B97-1, B3LYP, CAM-B3LYP) methods with Dunning’s correlation-consistent aug-cc-pVDZ and d-aug-cc-pVDZ basis sets. The most stable isomer 7H exhibits a μ value smaller than that of the 9H form by a factor of ca. 3.5. The β _μ value of the 9H tautomer is strongly dependent on the computational method employed, as it dramatically influences the β _μ (9H)/β _μ (7H) ratio, which at the highest correlated MP4/aug-cc-pVDZ level is predicted to be ca. 5. The Coulomb-attenuating hybrid exchange-correlation CAM-B3LYP method is superior to the conventional PBE1PBE, B3LYP, and B97-1 functionals in predicting the β _μ values. Differences between the largest diagonal hyperpolarizability components were clarified through hyperpolarizability density analyses. Dipole moment and first-order hyperpolarizability are molecular properties that are potentially useful for distinguishing the 7H from the 9H tautomer.

Theoretical study on aluminum carbide endohedral fullerene-Al4C@C80

The possibility of a new endohedral fullerene with a trapped aluminum carbide cluster, Al₄C @C₈₀-I _h , was theoretical investigated. The geometries and electronic properties of it were investigated using density functional theory methods. The Al₄C unit formally transfers six electrons to the C₈₀ cage which induces stabilization of Al₄C@C₈₀. A favorable binding energy, relatively large HOMO-LUMO gap, electron affinities and ionization potentials suggested the Al₄C@C₈₀ is rather stable. The analysis of vertical ionization potential and vertical electron affinity indicate Al₄C@C₈₀ is a good electron acceptor.

Probing the specificity of gamma-glutamylamine cyclotransferase: an enzyme involved in the metabolism of transglutaminase-catalyzed protein crosslinks

γ-Glutamylamine cyclotransferase (gGACT) catalyzes the intramolecular cyclization of a variety of l-γ-glutamylamines producing 5-oxo-l-proline and free amines. Its substrate specificity implicates it in the downstream metabolism of transglutaminase products, and is distinct from that of γ-glutamyl cyclotransferase which acts on l-γ-glutamyl amino acids. To elucidate the mechanism by which gGACT distinguishes between l-γ-glutamylamine and amino acid substrates, the specificity of the rabbit kidney enzyme for the amide region of substrates was probed through the kinetic analysis of a series of l-γ-glutamylamines. The isodipeptide N ^?-(l-γ-glutamyl)-l-lysine 1 was used as a reference. The kinetic constants of the l-γ-glutamyl derivative of n-butylamine 7, were nearly identical to those of 1. Introduction of a methyl or carboxylate group on the carbon adjacent to the side-chain amide nitrogen in l-γ-glutamylamine substrates resulted in a dramatic decrease in substrate properties for gGACT thus providing an explanation of why gGACT does not act on l-γ-glutamyl amino acids except for l-γ-glutamylglycine. Placement of substituents on carbons further removed from the side-chain amide nitrogen in l-γ-glutamylamines restored activity for gGACT, and l-γ-glutamylneohexylamine 19 had a higher specificity constant (k _cat /K _m) than 1. gGACT did not exhibit any stereospecificity in the amide region of l-γ-glutamylamine substrates. In addition, analogues (26–30) with heteroatom substitutions for the γ methylene position of the l-γ-glutamyl moiety were examined. Several thiocarbamoyl derivatives of l-cysteine (28–30) were excellent substrates for gGACT.     

In pursuit of negative Fukui functions: examples where the highest occupied molecular orbital fails to dominate the chemical reactivity

In our quest to explore molecules with chemically significant regions where the Fukui function is negative, we explored reactions where the frontier orbital that indicates the sites for electrophilic attack is not the highest occupied molecular orbital. The highest occupied molecular orbital (HOMO) controls the location of the regions where the Fukui function is negative, supporting the postulate that negative values of the Fukui function are associated with orbital relaxation effects and nodal surfaces of the frontier orbitals. Significant negative values for the condensed Fukui function, however, were not observed.

Ab initio molecular dynamics simulation on the formation process of He@C60 synthesized by explosion

The applications of endohedral non-metallic fullerenes are limited by their low production rate. Recently, an explosive method developed in our group shows promise to prepare He@C₆₀ at fairly high yield, but the mechanism of He inserting into C₆₀ cage at explosive conditions was not clear. Here, ab initio molecular dynamics analysis has been used to simulate the collision between C₆₀ molecules at high-temperature and high-pressure induced by explosion. The results show that defects formed on the fullerene cage by collidsion can effectively decrease the reaction barrier for the insertion of He into C₆₀, and the self-healing capability of the defects was also observed.

The effect of C-vacancy on hydrogen storage and characterization of H2 modes on Ti functionalized C60 fullerene A first principles study

Density functional theory calculations were performed to examine the effect of a C vacancy on the physisorption of H₂ onto Ti-functionalized C₆₀ fullerene when H₂ is oriented along the x-, y-, and z-axes of the fullerene. The effect of the C vacancy on the physisorption modes of H₂ was investigated as a function of H₂ binding energy within the energy window (?0.2 to ?0.6 eV) targeted by the Department of Energy (DOE), and as functions of a variety of other physicochemical properties. The results indicate that the preferential orientations of H₂ in the defect-free (i.e., no C vacancy) C₆₀TiH₂ complex are along the x- and y-axes of C₆₀ (with adsorption energies of ?0.23 and ?0.21 eV, respectively), making these orientations the most suitable ones for hydrogen storage, in contrast to the results obtained for defect-containing fullerenes. The defect-containing (i.e., containing a C vacancy) C₅₉TiH₂ complex do not exhibit adsorption energies within the targeted energy range. Charge transfer occurs from Ti 3d to C 2p of the fullerene. The binding of H₂ is dominated by the pairwise support–metal interaction energy E(i)^Cn...Ti, and the role of the fullerene is not restricted to supporting the metal. The C vacancy enhances the adsorption energy of Ti, in contrast to that of H₂. A significant reduction in the energy gap of the pristine C₆₀ fullerene is observed when TiH₂ is adsorbed by it. While the C _n fullerene readily participates in nucleophilic processes, the adjacent TiH₂ fragment is available for electrophilic processes.

Molecular simulation study of PAMAM dendrimer composite membranes

Pure polysulfone (PSF) and its composites with chitosan (CST), hyaluronic acid (HA), conventional poly(amidoamine), and hydroxyl poly(amidoamine) dendrimers as the membranes for separation of the gases, methane, carbon dioxide, hydrogen sulfide, nitrogen, and oxygen have been studied by molecular dynamics (MD) and grand canonical Monte Carlo (GCMC) simulations. The transport properties (solubility, diffusivity, and permeability) of pure and gas mixtures in the membranes were calculated and the results of the simulations were compared with the available experimental data. The simulated structural properties of the pure and composite PSF membranes including occupied volume, free volume, surface area, fractional free volume (FFV), and radius of gyration (R _g ) were evaluated and their effects on the separability of the gases by the membranes were analyzed and interpreted by the obtained results.

Synthesis and biological evaluation of cationic fullerene quinazolinone conjugates and their binding mode with modeled Mycobacterium tuberculosis hypoxanthine-guanine phosphoribosyltransferase enzyme

The present work reports a series of novel cationic fullerene derivatives bearing a substituted-quinazolinone moiety as a side arm. Fullerene-quinazolinone conjugates synthesized using the 1,3-dipolar cycloaddition reaction of C₆₀ with azomethine ylides generated from the corresponding Schiff bases of substituted quinazolinone were characterized by elemental analysis, FT-IR, ¹H NMR, ¹³C NMR and ESI-MS and screened for their antibacterial activity against Mycobacterium tuberculosis (H ₃₇ Rv strain). All the compounds exhibited significant activity with the most effective having MIC in the range of 1.562–3.125 μg/mL. Compound 9f exhibited good biological activity compared to standard drugs. We developed a computational strategy based on the modeled M. tuberculosis hypoxanthine-guanine phosphoribosyltransferase (HGPRT) using homology modeling techniques and studied its binding pattern with synthesized fullerene derivatives. We then explored the surface geometry of the protein to place the cage adjacent to the active site while optimizing its quinazolinone side arm to establish H bonding with active site residues.

First-principles vdW-DF investigation on the interaction between the oxazepam molecule and C60 fullerene

The interaction between oxazepam and C₆₀ fullerene was explored using first-principles vdW-DF calculations. It was found that oxazepam binds weakly to the fullerene cage via its carbonyl group. The binding of oxazepam to C₆₀ is affected drastically by nonlocal dispersion interactions, while vdW forces affect the corresponding geometries only a little. Furthermore, aqueous solution affects the geometries of the oxazepam approaching to fullerene slightly, while oxazepam binds slightly farther away from the nanocage. The results presented provide evidence for the applicability of the vdW-DF method and serve as a practical benchmark for the investigation of host–guest interactions in biological systems.