Benchmarking Gaussian Basis Sets in Quantum-Chemical Calculations of Photoabsorption Spectra of Light Atomic Clusters
Ontology highlight
ABSTRACT:
SUBMITTER:
PROVIDER: S-EPMC9798528 | biostudies-literature | 2022 Dec
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Hydrogen sulfide is toxic and corrosive gas abundantly available in nature. The activation of hydrogen sulfide to produce hydrogen and elemental sulfur is of great significance for possible applications in toxic pollutant control and hydrogen energy regeneration. The activation of H2S by transition metal atoms (M = Cr, Mn, and Fe) has been studied by low-temperature matrix isolation infrared spectroscopy and quantum chemical calculations. Experimental and theoretical results indicate that the reaction between ground-state M atoms and H2S is inhibited by the repulsive interactions between the reactants. After being excited upon photolysis, the corresponding excited-state M atoms react with H2S molecules spontaneously. The produced insertion product HMSH further decomposed to metal sulfides upon full-arc mercury lamp irradiation by the splitting of hydrogen.
| S-EPMC9404174 | biostudies-literature
Project description:It is customary in molecular quantum chemistry to adopt basis set libraries in which the basis sets are classified according to either their size (triple-ζ, quadruple-ζ, ...) and the method/property they are optimal for (correlation-consistent, linear-response, ...) but not according to the chemistry of the system to be studied. In fact the vast majority of molecules is quite homogeneous in terms of density (i.e., atomic distances) and types of bond involved (covalent or dispersive). The situation is not the same for solids, in which the same chemical element can be found having metallic, ionic, covalent, or dispersively bound character in different crystalline forms or compounds, with different packings. This situation calls for a different approach to the choice of basis sets, namely a system-specific optimization of the basis set that requires a practical algorithm that could be used on a routine basis. In this work we develop a basis set optimization method based on an algorithm-similar to the direct inversion in the iterative subspace-that we name BDIIS. The total energy of the system is minimized together with the condition number of the overlap matrix as proposed by VandeVondele et al. [VandeVondele et al. J. Chem. Phys. 2007, 227, 114105]. The details of the method are here presented, and its performance in optimizing valence orbitals is shown. As demonstrative systems we consider simple prototypical solids such as diamond, graphene sodium chloride, and LiH, and we show how basis set optimizations have certain advantages also toward the use of large (quadruple-ζ) basis sets in solids, both at the DFT and Hartree-Fock level.
| S-EPMC7997400 | biostudies-literature
Project description:The luminescence line shape is an important feature of semiconductor quantum dots (QDs) and affects performance in various optical applications. Here, we report a first-principles method to predict the luminescence spectrum of thousands of atom QDs. In our approach, neural network potential calculations are combined with density functional theory calculations to describe exciton-phonon coupling (EPC). Using the calculated EPC, the luminescence spectrum is evaluated within the Franck-Condon approximation. Our approach results in the luminescence line shape for an InP/ZnSe core/shell QD (3406 atoms) that exhibits excellent agreement with the experiments. From a detailed analysis of EPC, we reveal that the coupling of both acoustic and optical phonons to an exciton are important in determining the spectral line shapes of core/shell QDs, which is in contrast with previous studies. On the basis of the present simulation results, we provide guidelines for designing high-performance core/shell QDs with ultrasharp emission spectra.
| S-EPMC9888616 | biostudies-literature
Project description:The prediction of isomer shifts in 57Fe Mossbauer spectra is typically achieved by building calibration lines using the values of the density at the nuclear position. Using Slater-type orbital basis or large and specific Gaussian-type orbital basis has been thus far mandatory to achieve accurate predictions with density functional theory methods. In this work, we show that replacing the value of the density at the nucleus by the density integrated in a sphere of radius 0.06 au centered on the Fe nuclei yields excellent calibration lines (r2 = 0.976) with a high predictive power (q2 = 0.975, MAE = 0.055 mm·s-1) while using the conventional def2-TZVP basis set and X-ray geometrical parameters. Our data set comprises 69 57Fe-containing compounds and 103 signals. We also find B3LYP performing significantly better than the PW91 functional.
| S-EPMC8675134 | biostudies-literature
Project description:This paper presents the first converged frequency-dependent HF polarizability results for general molecules, on a set of 89 closed-shell atoms and molecules. The solver employs multiresolution analysis (MRA) in a multiwavelet basis to compute both ground and response states to a guaranteed precision, which are validated against independent numerical grid calculations on atoms and linear molecules. The MRA ground-state energies and response properties are used to evaluate results in correlation-consistent basis sets up to 5Z augmented with either single or double diffuse functions and core-polarization functions. Systematic trends are revealed through consideration of chemical composition as well as the use of machine learning to cluster convergence trends, the latter suggesting the possibility of learning and correcting basis-set error.
| S-EPMC11209939 | biostudies-literature
Project description:Multicomponent methods are a conceptually simple way to include nuclear quantum effects into quantum chemistry calculations. In multicomponent methods, the electronic molecular orbitals are described using the linear combination of atomic orbitals approximation. This requires the selection of a one-particle electronic basis set which, in practice, is commonly a correlation-consistent basis set. In multicomponent method studies, it has been demonstrated that large electronic basis sets are required for quantum hydrogen nuclei to accurately describe electron-nuclear correlation. However, as we show in this study, much of the need for large electronic basis sets is due to the correlation-consistent electronic basis sets not being optimized to describe nuclear properties and electron-nuclear correlation. Herein, we introduce a series of correlation-consistent electronic basis sets for hydrogen atoms called cc-pVnZ-mc with additional basis functions optimized to reproduce multicomponent density functional theory protonic densities. These new electronic basis sets are shown to yield better protonic densities with fewer electronic basis functions than the standard correlation-consistent basis sets and reproduce other protonic properties such as proton affinities and protonic excitation energies, even though they were not optimized for these purposes. The cc-pVnZ-mc basis sets should enable multicomponent many-body calculations on larger systems due to the improved computational efficiency they provide for a given level of accuracy.
| S-EPMC9910068 | biostudies-literature
Project description:We have calculated redox potentials of the two metal clusters in Mo-nitrogenase with quantum mechanical (QM) calculations. We employ an approach calibrated for iron-sulfur clusters with 1-4 Fe ions, involving QM-cluster calculations in continuum solvent and large QM systems (400-500 atoms), based on structures from combined QM and molecular mechanics (QM/MM) geometry optimisations. Calculations on the P-cluster show that we can reproduce the experimental redox potentials within 0.33 V. This is similar to the accuracy obtained for the smaller clusters, although two of the redox reactions involve also proton transfer. The calculated P1+/PN redox potential is nearly the same independently of whether P1+ is protonated or deprotonated, explaining why redox titrations do not show any pH dependence. For the FeMo cluster, the calculations clearly show that the formal oxidation state of the cluster in the resting E0 state is MoIIIFe3IIFe4III , in agreement with previous experimental studies and QM calculations. Moreover, the redox potentials of the first five E0-E4 states are nearly constant, as is expected if the electrons are delivered by the same site (the P-cluster). However, the redox potentials are insensitive to the formal oxidation states of the Fe ion (i.e., whether the added protons bind to sulfide or Fe ions). Finally, we show that the later (E4-E8) states of the reaction mechanism have redox potential that are more positive (i.e., more exothermic) than that of the E0/E1 couple.
| S-EPMC9822455 | biostudies-literature
Project description:We report absolute photoabsorption cross sections for gas-phase 2- and 5-bromopyrimidine in the 3.7-10.8 eV energy range, in a joint theoretical and experimental study. The measurements were carried out using high-resolution vacuum ultraviolet synchrotron radiation, with quantum chemical calculations performed through the nuclear ensemble approach in combination with time-dependent density functional theory, along with additional Franck-Condon Herzberg-Teller calculations for the first absorption band (3.7-4.6 eV). The cross sections of both bromopyrimidines are very similar below 7.3 eV, deviating more substantially from each other at higher energies. In the 7.3-9.0 eV range where the maximum cross-section is found, a single and broad band is observed for 5-bromopyrimidine, while more discernible features appear in the case of 2-bromopyrimidine. Several π* ← π transitions account for the most intense bands, while weaker ones are assigned to transitions involving the nitrogen and bromine lone pairs, the antibonding σ*Br orbital, and the lower-lying Rydberg states. A detailed comparison with the available photo-absorption data of bromobenzene is also reported. We have found significant differences regarding the main absorption band, which is more peaked in bromobenzene, becoming broader and shifting to higher energies in both bromopyrimidines. In addition, there is a significant suppression of vibrational structures and of Rydberg states in the pair of isomers, most noticeably for 2-bromopyrimidine.
| S-EPMC8235550 | biostudies-literature
Project description:A set of atomic polarizability parameters for a new polarizable Gaussian model (pGM) has been developed with the goal to accurately reproduce the polarizability anisotropy, taking advantage of its ability to attenuate all short-range electrostatic interactions, by fitting the ab initio molecular polarizability tensors ( A pq) calculated at the B3LYP/aug-cc-pVTZ level. For comparison, we also rederived the parameters for three Thole models in which the 1-2 (bonded), 1-3 (separated by two bonds), and 1-4 (separated by three bonds) interactions are fully included. The average percent errors (APEs) of molecular polarizability tensors for 4842 molecules or dimers are 2.98, 3.76, 3.28, and 3.82% for the pGM, Thole linear, Thole exponential, and Thole Amoeba models, respectively, with atom-type independent, universal screening factors (USF). The APEs are reduced further to 2.30, 2.69, 2.25, and 2.48% for the four corresponding polarizable models with atom-type dependent, variable screening factors (VSF). It is encouraging that the pGM with variable screening factors achieved APEs of 1.83 for 1155 amino acid analogs, dipeptides, and tetrapeptides, 1.39 for 28 nucleic acid bases, 0.708 for 1464 water clusters, and 1.99 for 85 dimers of water and biological building blocks. Compared to the new set of models, the APEs of the old Thole models that were fitted to isotropic molecular polarizabilities are 8.7% for set A (without the 1-2 and 1-3 interactions) and 6.3% for set D (with the 1-2 and 1-3 interactions) models, respectively. MPAD, a metric of molecular polarization anisotropy difference based on the diagonal terms of molecular polarizability tensors was defined and applied to assess the polarizable models in reproducing the ab initio molecular polarization anisotropy. The MPADs are 3.71, 4.70, 4.11, and 4.77% for the pGM, Thole linear, Thole exponential, and Thole Amoeba USF models, respectively. The APEs are reduced further to 2.85, 3.58, 2.90, and 3.15% for the four corresponding VSF models. Thus, the new pGM and Thole models notably improve molecular polarization anisotropy. Since pGM attenuates all short-range electrostatic interactions, its application is expected to improve stability in charge fitting, energy, and force calculations and the accuracy of multibody polarization.
| S-EPMC7197406 | biostudies-literature
Project description:Ligand enrichment among top-ranking hits is a key metric of molecular docking. To avoid bias, decoys should resemble ligands physically, so that enrichment is not simply a separation of gross features, yet be chemically distinct from them, so that they are unlikely to be binders. We have assembled a directory of useful decoys (DUD), with 2950 ligands for 40 different targets. Every ligand has 36 decoy molecules that are physically similar but topologically distinct, leading to a database of 98,266 compounds. For most targets, enrichment was at least half a log better with uncorrected databases such as the MDDR than with DUD, evidence of bias in the former. These calculations also allowed 40x40 cross-docking, where the enrichments of each ligand set could be compared for all 40 targets, enabling a specificity metric for the docking screens. DUD is freely available online as a benchmarking set for docking at http://blaster.docking.org/dud/.
| S-EPMC3383317 | biostudies-literature