Hydrodesulfurization of Dibenzothiophene: A Machine Learning Approach.

The hydrodesulfurization (HDS) process is widely used in the industry to eliminate sulfur compounds from fuels. However, removing dibenzothiophene (DBT) and its derivatives is a challenge. Here, the key aspects that affect the efficiency of catalysts in the HDS of DBT were investigated using machine learning (ML) algorithms. The conversion of DBT and selectivity was estimated by applying Lasso, Ridge, and Random Forest regression techniques. For the estimation of conversion of DBT, Random Forest and Lasso offer adequate predictions. At the same time, regularized regressions have similar outcomes, which are suitable for selectivity estimations. According to the regression coefficient, the structural parameters are essential predictors for selectivity, highlighting the pore size, and slab length. These properties can connect with aspects like the availability of active sites. The insights gained through ML techniques about the HDS catalysts agree with the interpretations of previous experimental reports.

Conformational preference of dipeptide zwitterions in aqueous solvents.

Proper description of solvent effects is challenging for theoretical methods, particularly if the solute is a zwitterion. Here, a series of theoretical procedures are used to determine the preferred solvated conformations of twelve hydrophobic dipeptides (Leu-Leu, Leu-Phe, Phe-Leu, Ile-Leu, Phe-Phe, Ala-Val, Val-Ala, Ala-Ile, Ile-Ala, Ile-Val, Val-Ile and Val-Val) in the zwitterionic state. First, the accuracy of density functional theory (DFT), combined with different implicit solvent models, for describing zwitterions in aqueous solvent is assessed by comparing the predicted against the experimental glycine tautomerization energy, i.e., the energetic difference between canonical and zwitterionic glycine in aqueous solvents. It is found that among the tested solvation schemes, the charge-asymmetric nonlocally determined local-electric solvation model (CANDLE) predicts an energetic difference in excellent agreement with the experimental value. Next, DFT-CANDLE is used to determine the most favorable solvated conformation for each of the investigated dipeptide zwitterions. The CANDLE-solvated structures are obtained by exploring the conformational space of each dipeptide zwitterion concatenating DFT calculations, in vacuum, with classical molecular dynamics simulations, in explicit solvents, and DFT calculations including explicit water molecules. It is found that the energetically most favorable conformations are similar to those of the dipeptide zwitterions in their respective crystal structures. Such structural agreement is indicative of the DFT-CANDLE accomplishment of the description of solvated zwitterions, and suggests that these biomolecules self-assemble as quasi-rigid objects.

Activated layered double hydroxides: assessing the surface anion basicity and its connection with the catalytic activity in the cyanoethylation of alcohols.

Layered double hydroxides (LDHs) act as catalysts in several reactions like in the cyanoethylation of alcohols with acrylonitrile to produce alkoxypropionitriles. Here we report an experimental and theoretical study in which it is shown that the experimental catalytic activity of LDHs in the cyanoethylation of 2-propanol and methanol correlates with the predicted strength of the basicity of the adsorbed surface species. First, it is shown that using activated LDHs containing Mg2+ and Al3+ (MgAl-LDH), Mg2+ and Ga3+ (MgGa-LDH), and Mg2+, Al3+ and Ga3+ (MgAlGa-LDH) great conversions to alkoxypropionitriles in high yields are obtained. Next, the basicity of these LDHs is estimated by means of the local softness, a local reactivity index calculated using density functional theory and appropriate surface models. For that, the adsorption of hydroxide and methoxide anions at the (001) surface of MgAl and MgGa-LDHs is investigated. We include LDHs containing Zn2+ and Al3+ (ZnAl-LDH) and Zn2+ and Ga3+ (ZnGa-LDH) in this part of the study to account for the effect of changing the divalent and trivalent metal composition on the basicity. It is found that hydroxide anions adsorbed on the MgGa-LDH surface and methoxide anions adsorbed on the MgAl-LDH surface are the most basic ones. This basicity trend correlates with our experimental findings about the catalytic activity of the activated LDHs. Further analyzing the connection between the LDH composition and the anion basicity, it is argued that the key steps dictating the LDH catalytic activity are the alcohol deprotonation in the cyanoethylation of 2-propanol, as it has been previously suggested, and the methoxide anion attack to the acrylonitrile double bond in the methanol cyanoethylation reaction.

Electronic structure behavior of PbO2, IrO2, and SnO2 metal oxide surfaces (110) with dissociatively adsorbed water molecules as a function of the chemical potential.

A detailed analysis of the electronic structure of three different electrochemical interfaces as a function of the chemical potential (µ) is performed using the grand canonical density functional theory in the joint density functional theory formulation. Changes in the average number of electrons and the density of states are also described. The evaluation of the global softness, which measures the tendency of the system to gain or lose electrons, is straightforward under this formalism. The observed behavior of these quantities depends on the electronic nature of the electrochemical interfaces.

Establishing the Relationship between Quantum Capacitance and Softness of N-Doped Graphene/Electrolyte Interfaces within the Density Functional Theory Grand Canonical Kohn-Sham Formalism.

The joint density functional theory (JDFT) is applied in the context of the grand canonical Kohn-Sham theory to calculate the global and local softness of pristine and N-substituted graphene structures. A comparison is established between the different theoretical approaches to evaluate total capacitance, revealing that the JDFT approach presents the closest result of this property with experimental data. A model of series capacitors is used to determine the quantum and nonquantum contributions of total capacitance, which enables us to determine the limitations of the rigid band approximation for the studied systems. It is found that global chemical softness is proportional to the total capacitance measured in the experiments, when the geometry relaxation is neglected. In this context, it is possible to obtain quantum and total capacitance (and consequently softness) from an average number of electrons vs applied potential plots and the model of series capacitors. Likewise, the relation of capacitance and softness gives rise to a new definition of local capacitance within the JDFT formalism. The evaluation of global and local softness paves the way to analyze electrochemical surface reactivity as a function of applied potential for a solid-electrolyte interface.

Using local softness to reveal oxygen participation in redox processes in cathode materials.

In this paper, the use of chemical local softness s(r) is proposed as an alternative way of analyzing the initial redox processes that occur in cathode materials used for lithium-ion batteries. It is shown that the chemical local softness is a quantity able to capture the same effects as the standard analysis based on the projected density of states. Because of its own nature, the local softness reveals the atomic sites involved in charge-transfer events and allows a quantitative comparative analysis between different materials. As pointed out by Johannes et al. (Solid State Ion 286:83-89, 2016), this analysis can be used as an indicator of stability of cathode materials for Li-ion batteries.

Nature of Hydrogen Bonds and S···S Interactions in the l-Cystine Crystal.

The intermolecular interactions that govern the stability of the l-cystine crystal were studied. This task is accomplished by using density-functional theory (DFT) with the generalized-gradient approximation (GGA) and including many-body dispersion (MBD) interactions. The strengths of the different interactions within the molecular crystal were obtained by a decomposition of the total interaction energy in two-, three-, and four-body contributions. It was determined that most of the hydrogen bonds formed within the crystal are strong (13, 15, and 19 kcal/mol) and the van der Waals nature of the S···S interaction is fully confirmed. Also, the presence of strong repulsive three-body contributions is determined. The results obtained support the idea of designing crystal growth inhibitors for this system in such a way that, when inserted in the crystal, they maintain the disulfide bridge environment but its capacity of generate hydrogen-bond networks is removed.

Sensing polarization effects through the analysis of the effective C6 dispersion coefficients in NaCl solutions.

Density functional theory based ab initio molecular dynamics is used to obtain microscopic details of the interactions in sodium chloride solutions. By following the changes in the atomic C6 coefficients under the Tkatchenko-Scheffler's scheme, we were able to identify two different coordination situations for the Cl(-) ion with significant different capabilities to perform dispersion interactions. This capability is enhanced when the ion-ion distance corresponds to the contact ion-pair situation. Also, the oxygen and hydrogen atoms of the water molecules change their aptitudes to interact through van der Waals like terms when they are close to the cation region of the ion-pair. These results have interesting implications on the design of force fields to model electrolyte solutions.

Peptide sr11a from Conus spurius is a novel peptide blocker for Kv1 potassium channels.

More than a hundred conotoxins are known today and from them, only seven conopeptides have been identified to target voltage-gated potassium channels (Kv). Conotoxin sr11a belongs to the I(2)-superfamily which is characterized by four disulfide bridges and provokes muscle stiffness when injected intracranially in mice. The aim of this work was to test the biological activity of sr11a on recombinant voltage-gated Kv1 potassium channels expressed in Xenopus laevis oocytes. Peptide sr11a was purified by high-performance liquid chromatography from the venom of the vermivorous Conus spurius. We found that peptide sr11a inhibits the delayed rectifiers Kv1.2 and Kv1.6 but had not effect on the slowly inactivating Kv1.3 channel. The functional dyad composed of a basic Lys and a hydrophobic amino acid residue is a crucial structural element, regarding the binding properties and blocking activities of more than a hundred K(+) channel toxins. Peptide sr11a does not contain Lys residues and then, it lacks the functional dyad. Molecular modeling of peptide sr11a reveals the presence of exposed basic residues of Arg and suggests that Arg17 and Arg29 are important on its biological activity.

Theoretical investigations of oxygen-17 NMR chemical shifts to discriminate among helical forms.

17O, 15N, 13C, and 1H NMR chemical shieldings are calculated using density functional theory to differentiate among the three primarily helical forms, 310, alpha, and pi in polyalanine peptides under periodic boundary conditions. This study suggests 17O as the best observable, as it has been demonstrated to be sensitive to hydrogen bonding and highly affected by small changes in the polypeptide in helix conformations. This theoretical study seeks to characterize the subtle conformational differences of helical structures by NMR chemical shift observables which may lead to important questions in experimental structure determination on the basis of using chemical shifts to identify protein secondary structures.

Dimerization of thymol blue in solution: Theoretical evidence.

The possibility of dimerization of thymol blue was addressed by ab initio and force field calculations. In agreement with experimental information, a dimer forming symmetrical chemical environments for hydrogen bond formation was determined. This dimer is stable in vacuum and aqueous media and corresponds to the same protonated state proposed by the experiment. A comparison of the CVFF and MM3 force fields and ab initio results shows the suitability of CVFF to qualitatively describe this system.

Potential energy landscape of monolayer-surface systems governed by repulsive lateral interactions: the case of (3 x 3)-I-Pt(111).

Combined density functional theory (DFT) and Monte Carlo (MC) approach is applied to study the potential energy landscape of four iodine atoms adsorbed on the Pt(111) surface in a (3 x 3) unit cell. Three critical points were identified: (3 x 3)-sym and (3 x 3)-asym, corresponding to structures well known from experimental studies, while the third one (3 x 3)-zigzag is a new structure not reported before. An interaction model fitted to DFT calculations allows us to explain the difference between arrangements of iodine monolayer in vacuum, air, and solution environments as a result of different repulsion regimes.

Effect of double bonds on the conducting properties of ciguatoxin 3C and tetrahydropyrane-based polymers: a theoretical study.

The electronic structure of the ciguatoxin 3C is analyzed through the Kohn-Sham model by using two different kinds of basis sets: localized basis set (Gaussian functions) and nonlocalized basis set (plane wave functions). With the localized basis functions, two approximations are used for the exchange-correlation functional: the local density approximation and the generalized gradient approximation. With the nonlocalized basis set, just the local density approximation is used. The energy gap, obtained from the frontier molecular orbitals, for this molecule predicts that this system is a semiconductor, even when the number of double bonds is increased inside the structure. However, as large molecules built with the basic unit--the tetrahydropyrane--of the ciguatoxin 3C are found in nature, it suggests studying the gap in polymeric systems built with the basic unit of this molecule. It is demonstrated that the presence of double bonds reduces considerably the gap, indicating the possibility of forming conducting materials by introducing double bonds in this kind of molecular systems. Thus, molecules strongly linked with biological systems can be used as precursor to build electric conducting systems.

Structural transitions in the polyalanine alpha-helix under uniaxial strain.

We analyzed the response to strain of an infinite polyalanine chain in the alpha-helical conformation using density functional theory. Under compressive strain the alpha-helix is found to undergo structural transitions to a pi-helix when the length of the helix is reduced by more than 10%. Under tensile strain the structure changes into a 3(10)-helix when the length is stretched by more than 10%. Our analysis of these transitions shows that they proceed essentially in two steps: At first there is mainly a length change, and only with some delay the helix twist adjusts.

Substituent effects.

Inductive and resonance effects of different substituents (CH3, Cl, NH3, CN, NO2) on the reactivity of the nitrogen and oxygen centers of several saturated and unsaturated hydroxylamine derivates toward hard and soft electrophiles are analyzed in terms of various conceptual DFT based descriptors calculated using B3LYP/6-311G** method. In most cases, the +/-I and +/-M effects are correctly described by these descriptors. The way the substituent effect dies down as the distance between the substituent and the active center increases is also analyzed. It is observed that more than one effect is to be considered to properly understand the effect of a substituent on reactivity.

Basis set effects on frontier molecular orbital energies and energy gaps: a comparative study between plane waves and localized basis functions in molecular systems.

In order to study the Kohn-Sham frontier molecular orbital energies in the complete basis limit, a comparative study between localized functions and plane waves, obtained with the local density approximation exchange-correlation functional is made. The analyzed systems are ethylene and butadiene, since they are theoretical and experimentally well characterized. The localized basis sets used are those developed by Dunning. For the plane-waves method, the pseudopotential approximation is employed. The results obtained by the localized basis sets suggest that it is possible to get an estimation of the orbital energies in the limit of the complete basis set, when the basis set size is large. It is shown that the frontier molecular orbital energies and the energy gaps obtained with plane waves are similar to those obtained with a large localized basis set, when the size of the supercell and the plane-wave expansion have been appropriately calibrated.