Structure to Property: Chemical Element Embeddings for Predicting Electronic Properties of Crystals.

We present a new general-purpose machine learning model that is able to predict a variety of crystal properties, including Fermi level energy and band gap, as well as spectral ones such as electronic densities of states. The model is based on atomic representations that enable it to effectively capture complex information about each atom and its surrounding environment in a crystal. The accuracy achieved for band gaps exceeds results previously published. By design, our model is not restricted to the electronic properties discussed here but can be extended to fit diverse chemical descriptors. Its advantages are (a) its low computational requirements, making it an efficient tool for high-throughput screening of materials; and (b) the simplicity and flexibility of its architecture, facilitating implementation and interpretation, especially for researchers in the field of computational chemistry.

A polarizable valence electron density based force field for high-energy interactions between atoms and molecules.

High-accuracy molecular force field models suited for hot gases and plasmas are not as abundant as those geared toward ambient pressure and temperature conditions. Here, we present an improved version of our previous electron-density based force field model that can now account for polarization effects by adjusting the atomic valence electron contributions to match ab initio calculated Mulliken partial charges. Using a slightly modified version of the Hohenberg-Kohn theorem, we also include an improved theoretical formulation of our model when applied to systems with degenerate ground states. We present two variants of our polarizable model, fitted from ab initio reference data calculated at CCSD(T)/cc-pVTZ and CCSD(T)/CEP-31G levels of theory, that both accurately model water dimer interaction energies. Further improvements include the additional interaction components with fictitious non-spherically symmetric, yet atom-centered, electron densities and fitting the exchange and correlation coefficients against analytical expressions. The latter removes all unphysical oscillations that are observed in the previous non-polarizable variant of our force field.

Theoretical study of fructose adsorption and conversion to trioses on metal-organic frameworks.

The conversion of chemically modified biomass into more valuable chemicals has recently gained significant attention from industry. In this study, we investigate the adsorption of fructose and its conversion into two trioses, glyceraldehyde (GLA) and dihydroxyacetone (DHA), on metal-organic frameworks using density functional theory calculations. The reaction mechanism proceeds through two main steps: first, the opening of the fructose ring; second, the retro-aldol fragmentation, which is favored over intramolecular hydrogen shifts. The substitution of a tetravalent metal in the metal-organic framework leads to different adsorption strengths in the order Hf-NU-1000 > Zr-NU-1000 > Ti-NU-1000. The catalytic activities of Hf-NU-1000 and Zr-NU-1000 are found to be similar. Both are more active than Ti-NU1000, corresponding to their relative Lewis acidity. It was found that functionalization of the organic linkers of the Hf-NU-1000 MOF does not improve its catalytic activity. The catalytic activity follows the order Hf-MOF-808 > Hf-NU-1000 > Hf-UIO-66 when based on either the overall activation energy or the turnover frequency (TOF).

Investigation of the Suzuki-Miyaura cross-coupling reaction on a palladium H-beta zeolite with DFT calculations.

Metal or metal cluster-doped zeolites catalyse a wide variety of reactions. In this work, a coupling reaction between bromobenzene and phenylboronic acid to yield biphenyl with the Pd-H-Beta zeolite catalyst was investigated with density functional theory (DFT) calculations. Utilizing a model system with tetrahedral Pd4 clusters within the H-Beta zeolite, it was demonstrated that the catalyst exhibited notable reactivity by effectively reducing the activation energy barrier for the reaction. Our investigation revealed that the zeolite framework facilitated electron transfer to the Pd cluster, thereby increasing the reaction activity. The coupling reaction was shown to be exothermic and comprise three main steps: oxidative addition of bromobenzene (C6H5Br), transmetallation with phenylboronic acid (C6H5B(OH)2), and reductive elimination of biphenyl (C12H10). Specifically, in the transmetallation step, which was the rate-determining step, the C-B bond breaking in phenylboronic acid (C6H5B(OH)2) and the phenylboronate anion (C6H5B(OH)3-) were compared under neutral and basic conditions, respectively. This comprehensive study clarifies the mechanism for the reaction with the modified Pd zeolite catalyst and highlights the essential role of the zeolite framework.

Combined Experimental and Theoretical Study of the Synthesis of 5,7-Dihydroxy-4-methylcoumarin via a Pechmann Condensation in the Presence of UiO-66-SO3H Catalysts.

An efficient synthesis of 5,7-dihydroxy-4-methylcoumarin from phloroglucinol with ethyl acetoacetate in the UiO-66-SO3H metal-organic framework is reported. The potential of UiO-66-SO3H as a solid catalyst was determined through optimized-condition experiments and quantum molecular calculations. The optimal conditions for the synthesis of 5,7-dihydroxy-4-methylcoumarin with UiO-66-SO3H were as follows: phloroglucinol/ethyl acetoacetate molar ratio = 1:1.6, reaction time = 4 h, and temperature = 140 °C, for which the reaction yield reached 66.0%. The reusability of UiO-66-SO3H catalysts for Pechmann condensation was examined. The activation energy of the reaction occurring on a sulfonic group of the UiO-66-SO3H catalyst was 12.6 kcal/mol, which was significantly lower than 22.6 kcal/mol of the same reaction on the UiO-66 catalyst. To comprehend the reaction mechanism, density functional theory with the ONIOM approach was applied for the synthesis of coumarin on the UiO-66-SO3H and UiO-66 clusters. A possible reaction mechanism was proposed involving three steps: a trans-esterification step, an intramolecular hydroxyalkylation step, and a dehydration step. The rate-determining step was suggested to be the first step which acquired an activation energy of 15.7 and 29.5 kcal/mol, respectively. Information from this study can be used as guidelines to develop more efficient catalytic metal-organic frameworks for various organic syntheses.

Coordination polymer-forming liquid Cu(2-isopropylimidazolate).

The structure of the melt state of one-dimensional (1D) coordination polymer crystal Cu(isopropylimidazolate) (melting temperature T m = 143 °C) was characterized by DSC, variable temperature PXRD, solid-state NMR (SSNMR), viscoelastic measurements, XAS, and DFT-AIMD calculations. These analyses suggested "coordination polymer-forming liquid" formation with preserved coordination bonds above T m. Variable chain configurations and moderate cohesive interaction in adjacent chains are the keys to the rarely observed polymer-forming liquid. The melt structure is reminiscent of the common 1D organic polymer melts such as entanglement or random coil structures.

Microcracking of Ni-rich layered oxide does not occur at single crystal primary particles even abused at 4.7 V.

There is a controversial issue based on the particle cracking of the Ni-rich layered oxide cathode materials whether it occurs at the primary particles or the grain boundary. Herein, we found that the microcracking of NMC811 does not occur at single crystalline primary particles even abused at a severe upper cell voltage of 4.7 V having a lot of gas evolution since the single-crystal NMC811 has superior mechanical stability. The capacity retentions determined at 1C rate and a 100% state of charge (SOC) are 80% and 50% after 1000 cycles for single crystal and polycrystal NMC811, respectively.

Sensing the ortho Positions in C6Cl6 and C6H4Cl2 from Cl2- Formation upon Molecular Reduction.

The geometrical effect of chlorine atom positions in polyatomic molecules after capturing a low-energy electron is shown to be a prevalent mechanism yielding Cl2-. In this work, we investigated hexachlorobenzene reduction in electron transfer experiments to determine the role of chlorine atom positions around the aromatic ring, and compared our results with those using ortho-, meta- and para-dichlorobenzene molecules. This was achieved by combining gas-phase experiments to determine the reaction threshold by means of mass spectrometry together with quantum chemical calculations. We also observed that Cl2- formation can only occur in 1,2-C6H4Cl2, where the two closest C-Cl bonds are cleaved while the chlorine atoms are brought together within the ring framework due to excess energy dissipation. These results show that a strong coupling between electronic and C-Cl bending motion is responsible for a positional isomeric effect, where molecular recognition is a determining factor in chlorine anion formation.

Understanding the interactions between lithium polysulfides and anchoring materials in advanced lithium-sulfur batteries using density functional theory.

Lithium-sulfur batteries (LSBs) are promising energy storage devices because of their high theoretical capacity and energy density. However, the "shuttle" effect in lithium polysulfides (LiPSs) is an unresolved issue that can hinder their practical commercial application. Research on LSBs has focused on finding appropriate materials that suppress this effect by efficiently anchoring the LiPSs intermediates. Quantum chemical computations are a useful tool for understanding the mechanistic details of chemical interaction involving LiPSs, and they can also offer strategies for the rational design of LiPSs anchoring materials. In this perspective, we highlight computational and theoretical work performed on this topic. This includes elucidating and characterizing the adsorption mechanisms, and the dominant types of interactions, and summarizing the binding energies of LiPSs on anchoring materials. We also give examples and discuss the potential of descriptors and machine learning approaches to predict the adsorption strength and reactivity of materials. We believe that both approaches will become indispensable in modelling future LSBs.

Metal triflate formation of C12-C22 phenolic compounds by the simultaneous C-O breaking and C-C coupling of benzyl phenyl ether.

Catalytic pathways to produce high carbon number compounds from benzyl phenyl ether require multiple steps to break the aryl etheric carbon-oxygen bonds; these steps are followed by energy-intensive processes to remove oxygen atoms and/or carbon-carbon coupling. Here, we show an upgrading strategy to transform benzyl phenyl ether into large phenolic (C12-C22) compounds by a one-step C-O breaking and C-C coupling catalyzed by metal triflates under a mild condition (100 °C and 1 bar). Hafnium triflate was the most selective for the desired products. In addition, we measured the effect of solvent polarity on the catalytic performance. Solvents with a polarity index of less than 3.4 promoted the catalytic activity and selectivity to C12-C22 phenolic products. These C12-C22 phenolic compounds have potential applications for phenol-formaldehyde polymers, diesel/jet fuels, and liquid organic hydrogen carriers.

Adsorption and dehydration of ethanol on isomorphously B, Al, and Ga substituted H-ZSM-5 zeolite: an embedded ONIOM study.

Dehydration reactions are important in the petroleum and petrochemical industries, especially for the feedstock production. In this work, the catalytic activity of zeolites with different acidities for the dehydration of ethanol to ethylene and diethylether is investigated by density functional calculations on cluster models of three isomorphous B, Al, and Ga substituted H-ZSM-5 zeolites. Both unimolecular and bimolecular mechanisms are investigated. Detailed reaction profiles for the dehydration reaction, assuming either a stepwise or a concerted mechanism, were calculated by using the ONIOM(MP2:M06-2X) + SCREEP method. The adsorption energies of ethanol are -21.6, -28.1, and -27.7 kcal mol-1 on H-[B]-ZSM-5, H-[Al]-ZSM-5, and H-[Ga]-ZSM-5 zeolites, respectively. The activation energies for the rate-determining step of the unimolecular concerted mechanism for the ethylene formation are 48.5, 42.6, and 43.6 kcal mol-1 on H-[B]-ZSM-5, H-[Al]-ZSM-5, and H-[Ga]-ZSM-5 zeolites, respectively. The activation energies for the ethoxy formation as the rate-determining step for the bimolecular formation of diethylether are 42.3, 40.0, and 41.1 kcal mol-1 on H-[B]-ZSM-5, H-[Al]-ZSM-5, and H-[Ga]-ZSM-5 zeolites, respectively. The results indicate that the catalytic activities for the dehydration of ethanol decrease in the order H-[Al]-ZSM-5 ~ H-[Ga]-ZSM-5 > H-[B]-ZSM-5. Besides the acid strength, the zeolite framework affects the reaction by stabilizing the reaction intermediates, leading to more stable adsorption complexes and lower activation barriers.

Electronic structure and reactivity of tirapazamine as a radiosensitizer.

Tirapazamine (TP) has been shown to enhance the cytotoxic effects of ionizing radiation in hypoxic cells, thus making it a candidate for a radiosensitizer. This selective behavior is often directly linked to the abundance of O2. In this paper, we study the electronic properties of TP in vacuum, micro-hydrated from one up to three molecules of water and embedded in a continuum of water. We discuss electron affinities, charge distribution, and bond dissociation energies of TP, and find that these properties do not change significantly upon hydration. In agreement with its large electron affinity, and bond breaking triggered by electron attachment requires energies higher than 2.5 eV, ruling out the direct formation of bioactive TP radicals. Our results suggest, therefore, that the selective behavior of TP cannot be explained by a one-electron reduction from a neighboring O2 molecule. Alternatively, we propose that TP's hypoxic selectivity could be a consequence of O2 scavenging hydrogen radicals.

Density Functional Investigation of the Conversion of Furfural to Furfuryl Alcohol by Reaction with i-Propanol over UiO-66 Metal-Organic Framework.

Carbonyl C═O bond reduction via catalytic transfer hydrogenation (CTH) is one of the essential processes for biomass conversion to valuable chemicals and fuels. Here, we investigate the CTH of furfural to furfuryl alcohol with i-propanol on UiO-66 metal-organic frameworks using density functional theory calculations and linear scaling relations. Initially, the reaction over two defect sites presented on Zr-UiO-66, namely, dehydrated and hydrated sites, have been compared. The hydrated active site is favored over that on the dehydrated active site since the activation free energy of the rate-determining reaction step occurring on the hydrated active site is lower than that occurring on the dehydrated active site (14.9 vs 17.9 kcal/mol). The catalytic effect of exchanged tetravalent metals (Hf and Ti) on Zr-UiO-66 is also considered. We found that Hf-UiO-66 (13.5 kcal/mol) provides a lower activation energy than Zr-UiO-66 (14.9 kcal/mol) and Ti-UiO-66 (19.4 kcal/mol), which corresponds to it having a higher Lewis acidity. The organic linkers of UiO-66 MOFs play a role in stabilizing all of the species on potential energy surfaces. The linear scaling relationship also reveals the significant role of the UiO-66 active site in activating the carbonyl C═O of furfural, and strong relationships are observed between the activation free energy, the charge of the metal at the MOF active sites, and the complexation energies in reaction coordinates.

Effects of single and double active sites of Cu oxide clusters over the MFI zeolite for direct conversion of methane to methanol: DFT calculations.

In this work, we investigate the effect of various species of Cu oxide clusters including single and double active sites incorporated in the MFI zeolite framework for the direct conversion of methane to methanol. An M06-2X density functional calculation is employed to fine-tune the suitable number and species of active sites and to provide insights into the effect of the active sites on the reaction mechanism of methane to methanol. Two models, single and double active sites of Cu oxide clusters, have been chosen, in which the single active site of Cu oxide clusters, (mono(µ-oxo)dicopper(ii)), is located at the Al1'-Al12' pair ([Cu(µ-O)Cu]2+@Al1'-Al12'/MFI) or at the Al6-Al7 pair ([Cu(µ-O)Cu]2+@Al6-Al7/MFI) in the MFI framework. For the double active sites of Cu oxide clusters, two species of double active sites of Cu oxide are considered. The first one is the double active site of mono(µ-oxo)dicopper(ii) containingtwo Al-Al pairs (Al1'-Al12' and Al6-Al7 pairs) in the MFI framework (2[Cu(µ-O)Cu]2+/MFI) and the other is the double active site of trans-µ-1,2-peroxo dicopper(ii), which occupies two Al-Al pairs (Al1'-Al12' and Al6-Al7 pairs) in the MFI framework (2[Cu(µ-1,2-peroxo)Cu]2+/MFI). Furthermore, the activation energy for C-H bond dissociation in direct methane conversion to methanol is considered. Compared with the single active site of [Cu(µ-O)Cu]2+/MFI, the double active sites, in particular (2[Cu(µ-O)Cu]2+/MFI), exhibited the lowest activation energy, approximately 12.5 kcal mol-1. The high charge transfer between activated methane and Cu oxide active sites and also the high negative partial charge at the bridging oxygen of Cu oxide active sites, which directly interact with the methane molecule and abstracts its H atom, are considered as the important factors which affect the catalytic activity of Cu oxide clusters for direct methane conversion to methanol. These findings strongly support that the number and species of Cu oxide active sites incorporated in the MFI framework can highly affect the reaction mechanism of methane to methanol.

Effect of Intercalants inside Birnessite-Type Manganese Oxide Nanosheets for Sensor Applications.

Hydrazine is a common reducing agent widely used in many industrial and chemical applications; however, its high toxicity causes severe human diseases even at low concentrations. To detect traces of hydrazine released into the environment, a robust sensor with high sensitivity and accuracy is required. An electrochemical sensor is favored for hydrazine detection owing to its ability to detect a small amount of hydrazine without derivatization. Here, we have investigated the electrocatalytic activity of layered birnessite manganese oxides (MnO2) with different intercalants (Li+, Na+, and K+) as the sensor for hydrazine detection. The birnessite MnO2 with Li+ as an intercalant (Li-Bir) displays a lower oxidation peak potential, indicating a catalytic activity higher than the activities of others. The standard heterogeneous electron transfer rate constant of hydrazine oxidation at the Li-Bir electrode is 1.09- and 1.17-fold faster than those at the Na-Bir and K-Bir electrodes, respectively. In addition, the number of electron transfers increases in the following order: K-Bir (0.11 mol) < Na-Bir (0.17 mol) < Li-Bir (0.55 mol). On the basis of the density functional theory calculation, the Li-Bir sensor can strongly stabilize the hydrazine molecule with a large adsorption energy (-0.92 eV), leading to high electrocatalytic activity. Li-Bir also shows the best hydrazine detection performance with the lowest limit of detection of 129 nM at a signal-to-noise ratio of ∼3 and a linear range of 0.007-10 mM at a finely tuned rotation speed of 2000 rpm. Additionally, the Li-Bir sensor exhibits excellent sensitivity, which can be used to detect traces of hydrazine without any effect of interference at high concentrations and in real aqueous-based samples, demonstrating its practical sensing applications.

Chiral Macroporous MOF Surfaces for Electroassisted Enantioselective Adsorption and Separation.

The development of surfaces with chiral features is a fascinating challenge for modern materials science, especially when they are used for chiral separation technologies. In this contribution, the design of hierarchically structured chiral macroporous zeolitic imidazolate framework-8 (ZIF-8) electrodes is presented. They are elaborated by an electrochemical deposition-dissolution technique based on the electrodeposition of metal through a colloidal crystal template, followed by controlled electrooxidation. This generates locally metal cations, which can interact with a chiral ligand present in the solution to form metal-organic frameworks (MOFs). The macroporous structure facilitates the access of the chiral recognition sites, located in the mesoporous MOF, and thus helps to overcome mass transport limitations. The efficiency of the designed functional materials for chiral adsorption and separation can be fine-tuned by applying an adjustable electric potential to the electrode surfaces. This hierarchical chiral ZIF-8 structure was deposited at the walls of a microfluidic device and used as a stationary phase for enantioselective separation. The potential-controlled interaction between the stationary phase and the chiral analytes allows baseline separation of two enantiomers. This opens up interesting perspectives for using hierarchically structured chiral MOFs as an efficient material for the selective adsorption and separation of chiral compounds.

Dehydrogenation of ethanol to acetaldehyde with nitrous oxide over the metal-organic framework NU-1000: a density functional theory study.

The conversion of ethanol to more valuable hydrocarbon compounds receives great attention in chemical industries because it could diminish the dependency on petroleum as raw material. We investigate the catalytic performance of Fe-supported MOF NU-1000 for the dehydrogenation of ethanol to acetaldehyde with nitrous oxide (N2O) by deriving the relevant reaction profiles with density functional theory calculations. In the proposed mechanism, the activation barrier of the rate-determining step is almost four times lower in the presence of N2O than without it. The supported NU-1000 framework plays also important role since it facilitates electron transfers and stabilizes all species along the reaction coordinate. When considering the catalytic activity of tetravalent metal centers (Zr, Hf and Ti) substituted into NU-1000 it is found that their activity decreases in the order Hf ≥ Zr > Ti, based on activation energies and turnover frequencies (TOF). Concerning MOF linkers, we show that the catalytic activity is not further improved by functionalizing NU-1000 with either electron-donating or electron-withdrawing organic groups.

Theoretical study of methane adsorption and C─H bond activation over Fe-embedded graphene: Effect of external electric field.

The effect of an external electric field (EF) on the methane adsorption and its activation on iron-embedded graphene (Fe-GPs) are investigated by using the M06-L density functional method. The EF is applied in the perpendicular direction to the graphene in the range of -0.015 to +0.015 a.u. with the interval of 0.005 a.u. The effects of EF on the adsorption, transition state and product complexes of the methane activation reaction are revealed. The binding energies of methane on Fe site in Fe-GPs are increased from -12.9 to -15.3, -18.1 and -21.5 kcal/mol for the negative EF of -0.005, -0.010 and -0.015, respectively. By applying positive EF, the activation barriers for methane activation are reduced in range of 3-8 kcal/mol (around 12-31%) and the reaction energies are more exothermic. The positive EF kinetically favors the reaction compared to the system without EF. The adsorption and activation of methane on Fe-GPs can be easily tuned by adjusting the external electric field for various applications. © 2019 Wiley Periodicals, Inc.

Nanocavity effects of various zeolite frameworks on n-pentane cracking to light olefins: combination studies of DFT calculations and experiments.

Better control of the product selectivity of light olefins (e.g., ethylene and propylene) obtained from the n-pentane catalytic cracking process has attracted considerable attention from both scientific and petrochemical industrial points of view. In this context, we report insights into the effects of the nanocavities of various zeolite frameworks, including H-FER, H-ZSM-5, and H-FAU, representing small, medium, and large cavities, on the reaction mechanism of n-pentane cracking to light olefins by using M06-2X/6-31G(d,p) density functional calculations, eventually leading to fine-tuning the product distribution of light olefins. The reaction mechanism consists of the following two main steps: (i) the protolytic cracking of n-pentane to form a pentonium intermediate; and (ii) the subsequent dissociation of the intermediate to either ethane-propylene or ethylene-propane. The key reaction pathways controlling the product distribution of light olefins relate to the dissociation of the pentonium intermediate, which can produce selectively either propylene (P) or ethylene (E), resulting in a controllable P/E ratio. The differences in the activation energies for ethylene production compared with those of propylene production over H-FER, H-ZSM-5, and H-FAU are 6.7, 5.0, and 0.5 kcal mol-1, respectively. Compared with H-ZSM-5 and H-FAU, the higher difference in the activation energy of these two pathways over H-FER implies that the preferable production of ethane-propylene compared with ethylene-propane is more pronounced. It is therefore reasonable to conclude that a smaller pore zeolite such as H-FER eventually leads to a high ratio of production of propylene to ethylene, in accordance with experimental observations.

Phenol Tautomerization Catalyzed by Acid-Base Pairs in Lewis Acidic Beta Zeolites: A Computational Study.

We investigate the tautomerization of phenol catalyzed by acid-base active pair sites in Lewis acidic Beta zeolites by means of density functional calculations using the M06-L functional. An analysis of the catalytic mechanism shows that hafnium on the Beta zeolite causes the strongest absorption of phenol compared to zirconium, tin, and germanium. This can be rationalized by the highest delocalization of electron density between the Lewis site and the oxygen of phenol which can in turn be quantified by the perturbative E(2) stabilization energy. The reaction is assumed to proceed in two steps, the phenol O-H bond dissociation and the protonation of the intermediate to form the cyclohexa-2,4-dien-1-one product. The rate determining step is the first one with a free activation energy of 26.3, 25.0, 22.1 and 22.7 kcal mol-1 for Ge-Beta, Sn-Beta, Zr-Beta, and Hf-Beta zeolites, respectively. The turnover frequencies follow these reaction barriers. Hence, the intrinsic catalytic activity of the Lewis acidic Beta zeolites studied here is in the order of Hf-Beta≈Zr-Beta>Sn-Beta> Ge-Beta for the tautomerization of phenol.