Journal of Computational Chemistry最新文献

The previously introduced workflow to achieve an energetically and structurally optimized description of frontier bonds in quantum mechanical/molecular mechanics (QM/MM)-type applications was extended into the regime of computational material sciences at the example of a layered carbon model systems. Optimized QM/MM link bond parameters at HSEsol/6-311G(d,p) and self-consistent density functional tight binding (SCC-DFTB) were derived for graphitic systems, enabling detailed investigation of specific structure motifs occurring in graphene-derived structures $��v�i�a�$ quantum-chemical calculations. Exemplary molecular dynamics (MD) simulations in the isochoric-isothermic (NVT) ensemble were carried out to study the intercalation of lithium and the properties of the Stone–Thrower–Wales defect. The diffusivity of lithium as well as hydrogen and proton adsorption on a defective graphene surface served as additional example. The results of the QM/MM MD simulations provide detailed insight into the applicability of the employed link-bond strategy when studying intercalation and adsorption properties of graphitic materials.

DFT functionals are of paramount importance for an accurate electronic and structural description of transition metal systems. In this work, a systematic analysis using some well-known and commonly used DFT functionals is performed. A comparison of the structural and energetic parameters calculated with the available experimental data is made in order to find the adequate functional for an accurate description of the TiO₂ bulk and surface of both anatase and rutile structures. In the absence of experimental data on the surface energy, the theoretical predictions obtained using the high-accuracy HSE06 functional were used as a reference to compare against the surface energy values calculated with the other DFT functionals. A clear improvement in the electronic description of both anatase and rutile was observed by introducing the Hubbard U correction term to PBE, PW91, and OptPBE functionals. The OptPBE-U4 functional was found to offer a good compromise between accurately describing the structural and electronic properties of titania.

A symbolic C++ library—Tinned—has been developed for symbolic differentiation and manipulation in response theory. By recognizing different key building blocks in the density matrix-based (Thorvaldsen et al., J. Chem. Phys. 2008, 129, 214108) and coupled-cluster response theories, we have implemented their corresponding C++ symbolic classes, including but not limited to one- and two-electron operators, exchange-correlation energy and potential, and coupled-cluster operator. Formulas of response theory can be well expressed in terms of the symbolic classes in the library Tinned. Their high-order perturbation-strength derivatives can be straightforwardly computed and extracted afterwards for numerical evaluation. The library Tinned will greatly facilitate the development work of response theory and may lead to a unified framework for response theory at different levels of electronic structure theory.

The regioselective radical CH trifluoromethylation of aromatic compounds have been shown to proceed in good yield and high regioselectivity when cyclodextrin (CD) is present. Yet, the reaction mechanism and the role of CD during the reaction have remained obscure. To this end, here we performed density functional theory (DFT) calculations to the conformations obtained by semiempirical quantum mechanical molecular dynamics calculations to reveal the reaction mechanism and the role of CD in controlling regioselectivity. The results show that metal salt increases the yield but do not affect the regioselectivity, which we further confirmed by an experiment. In contrast, multiple CD-substrate complex conformations and reaction pathways were obtained, and CD was shown to contribute to improving the regioselectivity by stabilizing the intermediate state via encapsulation. The present study indicates that CDs can increase the regioselectivity by stabilizing the intermediate and product states while only marginally affecting the transition state.

Theoretical modeling of the solid-state photocatalysis is one of the important issues as various useful photocatalysts have been developed to date. In this work, we investigated the mechanism of the alcohol photooxidation on niobium oxide (Nb₂O₅) which was experimentally developed, using the density functional theory (DFT)/time-dependent (TD)DFT calculations based on the cluster model. The alcohol adsorption and the first hydrogen transfer from hydroxy group to surface occur in the ground state, while the second hydrogen transfer from CH proceeds in the excited states during the photoirradiation of UV or visible light. The spin crossing was identified and the low-lying triplet states were solved for the reaction pathway. The photoabsorption in the visible light region was characterized as the charge transfer transition from O 2p of alcohol to Nb 4d of the Nb₂O₅ surface. The spin density and the natural population analysis indicated the generation of spin density in the moiety of carbonyl compound and its dissipation to the interface of the surface, which partly explains the electron paramagnetic resonance measurement. It was confirmed that the rate determining step is the desorption of carbonyl compound and water molecule in agreement with the experimental rate equation analysis. The present findings with the theoretical modeling will provide useful information for the further studies of the solid-state photocatalysis.

This computational study explores the copper (I) chloride catalyzed synthesis of (E)-1-(2,2-dichloro-1-phenylvinyl)-2-phenyldiazene (2Cl-VD) from readily available hydrazone derivative and carbon tetrachloride (CCl₄). 2Cl-VD has been extensively utilized to synthesize variety of heterocyclic organic compounds in mild conditions. The present computational investigations primarily focus on understanding the role of copper (I) and N¹,N¹,N²,N²-tetramethylethane-1,2-diamine (TMEDA) in this reaction, TMEDA often being considered a proton scavenger by experimentalists. Considering TMEDA as a ligand significantly alters the energy barrier. In fact, it is only 8.3 kcal/mol higher compared to the ligand-free (LF) route for the removal of a chlorine atom to form the radical ·CCl₃ but the following steps are almost barrierless. This intermediate then participates in attacking the electrophilic carbon in the hydrazone. Crucially, the study reveals that the overall potential energy surface is thermodynamically favorable, and the theoretical turnover frequency (TOF) value is higher in the case of Cu(I)-TMEDA complex catalyzed pathway.

Bandgap is a key property that determines electrical and optical properties in materials. Modulating the bandgap thus is critical in developing novel materials particularly semiconductors with improved features. This study examines the bandgap, highest occupied molecular orbital (HOMO), and lowest unoccupied molecular orbital (LUMO) energy level trends in a metal organic framework, metal–organic framework 5 (MOF-5), as a function of Hammett substituent effect (with the constant σ_m in the meta-position of the benzene ring) and solvent dielectric effect (with the constant ε). Specifically, experimental design and response surface methodologies helped to assess the significance of trends and correlations between these molecular properties with σ_m and ε. While the HOMO and LUMO decrease with increasing σ_m, the LUMO exhibits greater sensitivity to the substituent's electron withdrawing capability. The relative difference in these trends helps to explain why the bandgap tends to decrease with increasing σ_m.

In this article, we explore the construction of Hamiltonians with long-range interactions and their corrections using the short-range behavior of the wave function. A key aspect of our investigation is the examination of the one-particle potential, kept constant in our previous work, and the effects of its optimization on the adiabatic connection. Our methodology involves the use of a parameter-dependent potential dependent on a single parameter to facilitate practical computations. We analyze the energy errors and densities in a two-electron system (harmonium) under various conditions, employing different confinement potentials and interaction parameters. The study reveals that while the mean-field potential improves the expectation value of the physical Hamiltonian, it does not necessarily improve the energy of the system within the bounds of chemical accuracy. We also delve into the impact of density variations in adiabatic connections, challenging the common assumption that a mean field improves results. Our findings indicate that as long as energy errors remain within chemical accuracy, the mean field does not significantly outperform a bare potential. This observation is attributed to the effectiveness of corrections based on the short-range behavior of the wave function, a universal characteristic that diminishes the distinction between using a mean field or not.

We show herein that 1,10-dicyano substitution restricts the paragon fluxionality of bullvalene to just 14 isomers which isomerize along a single cycle. The restricted fluxionality of 1,10-dicyanobullvalene (DCB) is investigated by means of: (i) Bonding analyses of the isomer structures using the adaptive natural density partitioning (AdNDP). (ii) Quantum dynamical simulations of the isomerizations along the cyclic intrinsic reaction coordinate of the potential energy surface (PES). The PES possesses 14 equivalent potential wells supporting 14 isomers which are separated by 14 equivalent potential barriers supporting 14 transition states. Accordingly, at low temperatures, DCB appears as a hindered molecular rotor, without any delocalization of the wavefunction in the 14 potential wells, without any nuclear spin isomers, and with completely negligible tunneling. These results are compared and found to differ from those for molecular boron rotors. (iii) Born-Oppenheimer molecular dynamics (BOMD) simulations of thermally activated isomerizations. (iv) Calculations of the rate constants in the frame of transition state theory (TST) with reasonable agreement achieved with the BOMD results. (v) Simulations of the equilibration dynamics using rate equations for the isomerizations with TST rate coefficients. Accordingly, in the long-time limit, isomerizations of the 14 isomers, each with C_s symmetry, approach the “14 C_s → C_7v” thermally averaged structure. This is a superposition of the 14 equally populated isomer structures with an overall C_7v symmetry. By extrapolation, the results for DCB yield working hypotheses for so far un-explored properties e.g. for the equilibration dynamics of C₁₀H₁₀.

Graphene is the newest form of elemental carbon and it is becoming rapidly a potential candidate in the framework of nano-bio research. Many reports confirm the successful use of graphene-based materials as carriers of anticancer drugs having relatively high loading capacities compared with other nanocarriers. Here, the outcomes of a systematic study of the adsorption behavior of FDA approved Pt^II drugs cisplatin, oxaliplatin, and carboplatin on surface models of pristine, holey, and nitrogen-doped holey graphene are reported. DFT investigations in water solvent have been carried out considering several initial orientations of the drugs with respect to the surfaces. Adsorption free energies, calculated including basis set superposition error (BSSE) corrections, result to be significantly negative for many of the drug@carrier adducts indicating that tested layers could be used as potential carriers for the delivery of anticancer Pt^II drugs. The reduced density gradient (RDG) analysis allows to show that many kinds of non-covalent interactions, including canonical H-bond, are responsible for the stabilization of the formed adducts.