Journal of Computational Chemistry最新文献

Development of organic nonlinear optical materials has become progressively more important due to their emerging applications in new-generation photonic devices. A novel series of chromophores based on innovative thiophene and furan-fused cyclopentadienyl bridge with various powerful donor and acceptor moieties were designed and theoretically investigated for applications in nonlinear optics. To unravel the structure-property relationship between this new push-pull conjugated systems and their nonlinear optical property, multiple methods, including density of states analysis, coupled perturbed Kohn-Sham (CPKS) method, sum-over-states (SOS) model, the two-level model (TSM), hyperpolarizability density analysis, and the (hyper)polarizability contribution decomposition, were performed to comprehensively investigated the nonlinear optical and electronic properties of this new π-system. Due to excellent charge transfer ability of new bridge and distinctive structure of donor and acceptor, the designed chromophores exhibit deep HOMO levels, low excitation energy, high dipole moment difference and large hyperpolarizability, indicating the appealing air-stable property and remarkable electrooptic performance of them. Importantly, THQ-CS-A3 and PA-CS-A3 shows outstanding NLO response properties with β_tot value of 6953.9 × 10^-30 and 5066.0 × 10^-30 esu in AN, respectively. The influence of the push-pull strength, the heterocycle and the π-conjugation of new bridge on the nonlinear optical properties of this novel powerful systems are clarified. This new series of chromophores exhibit remarkable electro-optical Pockels and optical rectification effect. More interestingly, PA-CS-A3 and THQ-CS-A2 also show appealing SHG effect. This study will help people understand the nature of nonlinear optical properties of innovative heteroarene-fused based cyclopentadienyl chromophores and offer guidance for the rational design of chromophores with outstanding electrooptic (EO) performance in the future.

Chlorine is an important chemical which has long been produced in chlor-alkali process using dimensionally stable anodes (DSA). However, some serious drawbacks of DSA inspire the development of alternative anodes for chlorine evolution reaction (CER). In this study, we focused on the graphene- and carbon nanotube-supported platinum tetra-phenyl porphyrins as electrocatalysts for CER, which have been theoretically investigated based on density functional theory. Our results reveal that the supported substrates possess potential CER electrocatalytic activity with very low thermodynamic overpotentials (0.012-0.028 V) via Cl* pathway instead of ClO*. The electronic structures analyses showed that electron transfer from the support to the adsorbed chlorine via the Pt center leads to strong Pt-Cl interactions. Furthermore, the supported electrocatalysts exhibited excellent selectivity toward CER because of high overpotentials and reaction barriers of oxygen evolution process. Therefore, our results may pave the way for designing CER electrocatalyst utilizing emerging carbon nanomaterials.

A recent work by Marks et al. on the formation of carbamic acid in NH$��{\hspace{0.17em}}_{�3�}�$-CO$��{\hspace{0.17em}}_{�2�}�$ interstellar ices pointed out its stability in the gas phase and the concomitant production of its dimer. Prompted by these results and the lack of information on these species, we have performed an accurate structural, energetic and spectroscopic investigation of carbamic acid and its dimer. For the former, the structural and spectroscopic characterization employed composite schemes based on coupled cluster (CC) calculations that account for the extrapolation to the complete basis set limit and core correlation effects. A first important outcome is the definitive confirmation of the nonplanarity of carbamic acid, then followed by an accurate estimate of its rotational and vibrational spectroscopy parameters. As far as the carbamic acid dimer is concerned, the investigation started from the identification of its most stable forms. For them, structure and vibrational properties have been evaluated using density functional theory, while a composite scheme rooted in CC theory has been employed for the energetic characterization. Our results allowed us to provide a better interpretation of the feature observed in the recent experiment mentioned above.

This paper is the first to look at the structural, electronic, mechanical, optical, and thermodynamic properties of the ANiX (ASc, Ti, Y, Zr, Hf; XBi, Sn) half-Heusler (HH) using DFT based first principles method. The lattice parameters that we have calculated are very similar to those obtained in prior investigations with theoretical and experimental data. The positive phonon dispersion curve confirm the dynamical stability of ANiX (ASc, Ti, Y, Zr, Hf; XBi, Sn). The electronic band structure and DOS confirmed that the studied materials ANiX (ASc, Ti, Y, Zr, Hf; XBi, Sn) are direct band gap semiconductors. The investigation also determined significant constants, including dielectric function, absorption, conductivity, reflectivity, refractive index, and loss function. These optical observations unveiled our compounds potential utilization in various electronic and optoelectronic device applications. The elastic constants were used to fulfill the Born criteria, confirming the mechanical stability and ductility of the solids ANiX (ASc, Ti, Y, Zr, Hf; XBi, Sn). The calculated elastic modulus revealed that our studied compounds are elastically anisotropic. Moreover, ANiX (ASc, Ti, Y, Zr, Hf; XBi, Sn) has a very low minimum thermal conductivity (K_min), and a low Debye temperature (θ_D), which indicating their appropriateness for utilization in thermal barrier coating (TBC) applications. The Helmholtz free energy (F), internal energy (E), entropy (S), and specific heat capacity (Cv) are determined by calculations derived from the phonon density of states.

Cocrystals are assemblies of more than one type of molecule stabilized through noncovalent interactions. They are promising materials for improved drug formulation in which the stability, solubility, or biocompatibility of the active pharmaceutical ingredient (API) is improved by including a coformer. In this work, a range of density functional theory (DFT) and density functional tight binding (DFTB) models are systematically compared for their ability to predict the lattice enthalpy of a broad range of existing pharmaceutically relevant cocrystals. These range from cocrystals containing model compounds 4,4′-bipyridine and oxalic acid to those with the well benchmarked APIs of aspirin and paracetamol, all tested with a large set of alternative coformers. For simple cocrystals, there is a general consensus in lattice enthalpy calculated by the different DFT models. For the cocrystals with API coformers the cocrystals, enthalpy predictions depend strongly on the DFT model. The significantly lighter DFTB models predict unrealistic values of lattice enthalpy even for simple cocrystals.

Intermolecular interactions with polycyclic aromatic hydrocarbons (PAHs) represent an important area of physisorption studies. These investigations are often hampered by a size of interacting PAHs, which makes the calculation prohibitively expensive. Therefore, methods designed to deal with large molecules could be helpful to reduce the computational costs of such studies. Recently we have introduced a new systematic approach for the molecular fragmentation of PAHs, denoted as AROFRAG, which decomposes a large PAH molecule into a set of predefined small PAHs with a benzene ring being the smallest unbreakable unit, and which in conjunction with the Molecules-in-Molecules (MIM) approach provides an accurate description of total molecular energies. In this contribution we propose an extension of the AROFRAG, which provides a description of intermolecular interactions for complexes composed of PAH molecules. The examination of interaction energy partitioning for various test cases shows that the AROFRAG3 model connected with the MIM approach accurately reproduces all important components of the interaction energy. An additional important finding in our study is that the computationally expensive long-range electron-correlation part of the interaction energy, that is, the dispersion component, is well described at lower AROFRAG levels even without MIM, which makes the latter models interesting alternatives to existing methods for an accurate description of the electron-correlated part of the interaction energy.

Conformational ensemble generation and the search for the global minimum conformation are important problems in computational chemistry. In this work, a variant on the conformer-rotamer ensemble sampling tool (CREST) iterative metadynamics (iMTD) algorithm designed for determining structural ensembles and energetics of noncovalent clusters of flexible molecules is presented. We term this new algorithm a low-energy diversity-enhanced variant on CREST, or LEDE-CREST. As with CREST, the energies are evaluated using the semiempirical GFN2-xTB extended tight binding approach. The utility of the algorithm is highlighted by generating ensembles for a variety of noncovalent clusters of flexible or rigid monomers using both CREST and LEDE-CREST.

This study focuses on the systematic exploration of the emodepside conformations bound to monovalent K⁺ ion using quantum mechanical density functional theory (DFT) calculations at the M06-2X/6-31+G(d,p) level of theory. Nine conformers of emodepside and their complexes with K⁺ ion were characterized as stationary points on the potential energy surface. The conformational isomers were examined for their 3D structures, bonding, energetics, and interactions with the cation. A cavitand-like structure (CC) is identified to be the energetically most stable arrangement. To arrive at a better understanding of the K⁺ ion binding, calculations were initially performed on complexes formed by the K⁺ and Na⁺ ions with model ligands (methyl ester and N,N-dimethyl acetamide). Both the natural bond orbital (NBO) method and the block-localized wavefunction (BLW) energy decomposition approach was employed to assess the bonding and energetic contributions stabilizing the ion-bound model complexes. Finally, the solvent effect was evaluated through complete geometry optimizations and energy minimizations for the model ion-ligand complexes and the emodepside-K⁺ bound complexes using an implicit solvent model mimicking water and DMSO.