Journal of Computational Chemistry最新文献

Herein, we present a density functional theory with dispersion correction (DFT-D) calculations that focus on the intercalation of ionic liquids (ILs) electrolytes into the two-dimensional (2D) Ti₃C₂T_x MXenes. These ILs include the cation 1-ethyl-3-methylimidazolium (Emim⁺), accompanied by three distinct anions: bis(trifluoromethylsulfonyl)imide (TFSA⁻), (fluorosulfonyl)imide (FSA⁻) and fluorosulfonyl(trifluoromethanesulfonyl)imide (FTFSA⁻). By altering the surface termination elements, we explore the intricate geometries of IL intercalation in neutral, negative, and positive pore systems. Accurate estimation of charge transfer is achieved through five population analysis models, such as Hirshfeld, Hirshfeld-I, DDEC6 (density derived electrostatic and chemical), Bader, and VDD (voronoi deformation density) charges. In this work, we recommend the DDEC6 and Hirshfeld-I charge models, as they offer moderate values and exhibit reasonable trends. The investigation, aimed at visualizing non-covalent interactions, elucidates the role of cation-MXene and anion-MXene interactions in governing the intercalation phenomenon of ionic liquids within MXenes. The magnitude of this role depends on two factors: the specific arrangement of the cation, and the nature of the anionic species involved in the process.

Mechanistic investigations at the density functional theory level of organic and organometallic reactions in solution are now broadly accessible and routinely implemented to complement experimental investigations. The selection of an appropriate functional among the plethora of developed ones is the first challenge on the way to reliable energy barrier calculations. To provide guidelines for the choice of an initial and reliable computational level, the performances of commonly used non-empirical (PBE, PBE0, PBE0-DH) and empirical density functionals (BLYP, B3LYP, B2PLYP) were evaluated relative to experimental activation enthalpies. Most reactivity databases to assess density functional performances are primarily based on high level calculations, here a set of experimental activation enthalpies of organic and organometallic reactions performed in solution were selected from the literature. As a general trend, the non-empirical functionals outperform the empirical ones. The most accurate energy barriers are obtained with hybrid PBE0 and double-hybrid PBE0-DH density functionals, both providing similar performance. Regardless of the functional under consideration, the addition of the GD3-BJ empirical dispersion correction does not enhance the accuracy of computed energy barriers.

In this proof-of-concept paper, we show how exchange-correlation effects can be simply recovered for interatomic energies within the interacting quantum atoms decomposition when local, gradient generalized, or meta-gradient generalized approximations are used in density functional theory (DFT) calculations. We also demonstrate how inhomogeneity and non-local effects can be introduced even from a pure local scheme, without resorting to any orbital information. Finally, we provide numerical evidence on a database of selected energetic molecules that this decomposition scheme can be efficiently used to build accurate models for the prediction of molecular energies from an initial “cheap” DFT calculation.

This study investigates the formation of partial sigma (σ) covalent bonds in experimentally synthesizable biradicals formed from hydrogenated and fluorinated C₈, C₂₀, and C₆₀ cage structures, by assessing their stability, geometry, and bonding character in singlet and triplet states using restricted B3LYP-D3/6–31+G(d,p) theory, natural bond orbital (NBO) analysis, and complete active space self-consistent field (CASSCF) method. The results show that these partial σ_CC bonds have Wiberg bond orders of 0.38 to 0.48 and bond lengths ranging from 2.62 Å to 5.93 Å. Cage size influences the characteristics of the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), with electrons favoring more antibonding orbitals in smaller cages where electrons reside more on the exterior of the cage and favoring bonding orbitals in larger ones where electrons are more in the interior. Fluorination enhances electron density on bonding orbitals. The analysis further clarified that the differentiation between antibonding and bonding features of HOMOs and LUMOs extends beyond merely electron transfer from s- to p-atomic orbitals, also noting possible interactions of the same symmetry repel. The study also introduces hyperconjugation from α-position CH bonds as a factor in stabilizing partial σ-bond formation. The results also caution against the use of broken symmetry methodology in unrestricted SCF wavefunctions for biradicals, such as those in this study as it may cause large spin contamination and thus errors in the calculated electronic properties results.

Elucidating the speciation (in terms of concentration versus pH) and understanding the formation mechanisms of polyoxometalates remains a significant challenge, both in experimental and computational domains. POMSimulator is a new methodology that tackles this problem from a purely computational perspective. The methodology uses results from quantum mechanics based methods to automatically set up the chemical reaction network, and to build speciation models. As a result, it becomes possible to predict speciation and phase diagrams, as well as to derive new insights into the formation mechanisms of large molecular clusters. In this work we present the main features of the first open-source version of the software. Since the first report [Chem. Sci. 2020, 11, 8448-8456], POMSimulator has undergone several improvements to keep up with the growing challenges that were tackled. After four years of research, we recognize that the source code is sufficiently stable to share a polished and user-friendly version. The Python code, manual, examples, and install instructions can be found at https://github.com/petrusen/pomsimulator.

Oxazine dyes act as reporters of their near environment by the response of their fluorescence spectra. At the same time, their fluorescence spectra exhibit a pronounced vibrational progression. In this work, we computationally investigate the impact of near-environment models consisting of aggregated water as well as betaine molecules on the vibrational profile of fluorescence spectra of different oxazine derivatives. For aggregated betaine and a water molecule located above the plane of the dyes, we observe a distinct modification of the vibrational profile, which is more pronounced than the effect of a continuum description of a solvent environment. Our analysis shows that this effect cannot be explained by a pure change in the electronic structure, but that also vibrational degrees of freedom of the environment can be decisive for the vibrational profile and should, hence, not generally be neglected.

Several theoretical studies at different levels of theory have attempted to calculate the absolute position of the SnO₂ conduction band, whose knowledge is key for its effective application in optoelectronic devices such us, for example, perovskite solar cells. However, the predicted band edges fall outside the experimentally measured range. In this work, we introduce a computational scheme designed to calculate the conduction band minimum values of SnO₂, yielding results aligned with experiments. Our analysis points out the fundamental role of encompassing surface oxygen vacancies to properly describe the electronic profile of this material. We explore the impact of both bridge and in-plane oxygen vacancy defects on the structural and electronic properties of SnO₂, explaining from an atomistic perspective the experimental observables. The results underscore the importance of simulating both types of defects to accurately predict SnO₂ features and provide new fundamental insights that can guide future studies concerning design and optimization of SnO₂-based materials and functional interfaces.

The fragment molecular orbital (FMO) scheme is one of the popular fragmentation-based methods and has the potential advantage of making the circuit shallow for quantum chemical calculations on quantum computers. In this study, we used a GPU-accelerated quantum simulator (cuQuantum) to perform the electron correlation part of the FMO calculation as unitary coupled-cluster singles and doubles (UCCSD) with the variational quantum eigensolver (VQE) for hydrogen-bonded (FH)$��{\hspace{0.17em}}_{�3�}�$ and (FH)$��{\hspace{0.17em}}_{�2�}�$-H$��{\hspace{0.17em}}_{�2�}�$O systems with the STO-3G basis set. VQE-UCCSD calculations were performed using both canonical and localized MO sets, and the results were examined from the point of view of size-consistency and orbital-invariance affected by the Trotter error. It was found that the use of localized MO leads to better results, especially for (FH)$��{\hspace{0.17em}}_{�2�}�$-H$��{\hspace{0.17em}}_{�2�}�$O. The GPU acceleration was substantial for the simulations with larger numbers of qubits, and was about a factor of 6.7–7.7 for 18 qubit systems.

Density functional theory (DFT) calculations explore the stability of a single platinum atom on various flat, stepped, and defective ceria surfaces, in the context of single-atom catalysts (SACs) for the water–gas shift (WGS) reaction. The adsorption properties and diffusion kinetics of the metal strongly depend on the support termination with large stability on metastable and stepped CeO₂(100) and (210) surfaces where the diffusion of the platinum atom is hindered. At the opposite, the more stable CeO₂(111) and (110) terminations weakly bind the platinum atom and can promote the growth of metallic clusters thanks to fast diffusion kinetics. The adsorption of carbon monoxide on the single platinum atom supported on the various ceria terminations is also sensitive to the surface structure. Carbon monoxide weakly binds to the single platinum atom supported on reduced CeO₂(111) and (211) terminations. The desorption of the CO₂ formed during the WGS reaction is thus facilitated on the latter terminations. A vibrational analysis underlines the significant changes in the calculated scaled anharmonic CO stretching frequency on these catalysts.

According to the Hofmeister series, thiocyanate is the strongest “salting in” anion. In fact, it has a strong denaturant activity against the native state of globular proteins. A molecular level rationalization of the Hofmeister series is still missing, and therefore the denaturant activity of thiocyanate also awaits a robust explanation. In the last years, different types of experimental studies have shown that thiocyanate is capable to directly interact with both polar and nonpolar groups of polypeptide chains. This finding has been scrutinized via a careful computational procedure based on density functional theory approaches. The results indicate that thiocyanate is able to make H-bonds via both the nitrogen and sulfur atom, and to make strong van der Waals interactions with almost all the groups of polypeptide chains, regardless of their polarity.