No Arabic abstract
We developed an interface program between a program suite for an automated search of chemical reaction pathways, GRRM, and a program package of semiempirical methods, MOPAC. A two-step structural search is proposed as an application of this interface program. A screening test is first performed by semiempirical calculations. Subsequently, a reoptimization procedure is done by ab initio or density functional calculations. We apply this approach to ion adsorption on cellulose. The computational efficiency is also shown for a GRRM search. The interface program is suitable for the structural search of large molecular systems for which semiempirical methods are applicable.
High-throughput computational materials design promises to greatly accelerate the process of discovering new materials and compounds, and of optimizing their properties. The large databases of structures and properties that result from computational searches, as well as the agglomeration of data of heterogeneous provenance leads to considerable challenges when it comes to navigating the database, representing its structure at a glance, understanding structure-property relations, eliminating duplicates and identifying inconsistencies. Here we present a case study, based on a data set of conformers of amino acids and dipeptides, of how machine-learning techniques can help addressing these issues. We will exploit a recently developed strategy to define a metric between structures, and use it as the basis of both clustering and dimensionality reduction techniques showing how these can help reveal structure-property relations, identify outliers and inconsistent structures, and rationalise how perturbations (e.g. binding of ions to the molecule) affect the stability of different conformers.
A study of the precision of the semiempirical methods used in the determination of the chemical abundances in gas-rich galaxies is carried out. In order to do this the oxygen abundances of a total of 438 galaxies were determined using the electronic temperature, the $R_{23}$ and the P methods. The new calibration of the P method gives the smaller dispersion for the low and high metallicity regions, while the best numbers in the turnaround region are given by the $R_{23}$ method. We also found that the dispersion correlates with the metallicity. Finally, it can be said that all the semiempirical methods studied here are quite insensitive to metallicity with a value of $8.0pm0.2$ dex for more than 50% of the total sample. keywords{ISM: abundances; (ISM): H {sc ii} regions}
Using the simple (symmetric) Hubbard dimer, we analyze some important features of the $GW$ approximation. We show that the problem of the existence of multiple quasiparticle solutions in the (perturbative) one-shot $GW$ method and its partially self-consistent version is solved by full self-consistency. We also analyze the neutral excitation spectrum using the Bethe-Salpeter equation (BSE) formalism within the standard $GW$ approximation and find, in particular, that i) some neutral excitation energies become complex when the electron-electron interaction $U$ increases, which can be traced back to the approximate nature of the $GW$ quasiparticle energies; ii) the BSE formalism yields accurate correlation energies over a wide range of $U$ when the trace (or plasmon) formula is employed; iii) the trace formula is sensitive to the occurrence of complex excitation energies (especially singlet), while the expression obtained from the adiabatic-connection fluctuation-dissipation theorem (ACFDT) is more stable (yet less accurate); iv) the trace formula has the correct behavior for weak (ie, small $U$) interaction, unlike the ACFDT expression.
Eumelanin is regarded to be an attractive candidate material for biomedical applications. Despite many theoretical studies exploring the structure of eumelanin, an exact mapping of the energetic landscape of the very large phase space of eumelanin is still elusive. In this work, we implement a piecewise Ising Model to predict formation enthalpies of Eumelanin single and double tetramers, and demonstrate its superior predictive and generalizable capabilities. We believe this model will prove very useful in theoretically characterizing the many unique properties attributed to its disorder. The modular nature of the predictive Ising model built up in this work is well-suited for analysis and characterization of a larger phase space of eumelanin polymers, including hexamers and octomers, as well as larger stacked structures, such as potential triple and quadruple eumelanin tetramers. Absorbance data can be incorporated with population-wide predictions of polymer abundance to produce weighted-average predictions of broadband absorbance of bulk eumelanin.
The concept of chemical bonding can ultimately be seen as a rationalization of the recurring structural patterns observed in molecules and solids. Chemical intuition is nothing but the ability to recognize and predict such patterns, and how they transform into one another. Here we discuss how to use a computer to identify atomic patterns automatically, so as to provide an algorithmic definition of a bond based solely on structural information. We concentrate in particular on hydrogen bonding -- a central concept to our understanding of the physical chemistry of water, biological systems and many technologically important materials. Since the hydrogen bond is a somewhat fuzzy entity that covers a broad range of energies and distances, many different criteria have been proposed and used over the years, based either on sophisticate electronic structure calculations followed by an energy decomposition analysis, or on somewhat arbitrary choices of a range of structural parameters that is deemed to correspond to a hydrogen-bonded configuration. We introduce here a definition that is univocal, unbiased, and adaptive, based on our machine-learning analysis of an atomistic simulation. The strategy we propose could be easily adapted to similar scenarios, where one has to recognize or classify structural patterns in a material or chemical compound.