Search

Catalytic-materials design requires predictive modeling of the interaction between catalyst and reactants. This is challenging due to the complexity and diversity of structure-property relationships across the chemical space. Here, we report a strategy for a rational design of catalytic materials using the artificial intelligence approach (AI) subgroup discovery. We identify catalyst genes (features) that correlate with mechanisms that trigger, facilitate, or hinder the activation of carbon dioxide (CO2) towards a chemical conversion. The AI model is trained on first-principles data for a broad family of oxides. We demonstrate that surfaces of experimentally identified good catalysts consistently exhibit combinations of genes resulting in a strong elongation of a C-O bond. The same combinations of genes also minimize the OCO-angle, the previously proposed indicator of activation, albeit under the constraint that the Sabatier principle is satisfied. Based on these findings, we propose a set of new promising catalyst materials for CO2 conversion. Here the authors demonstrate an artificial-intelligence based approach to identify catalytic materials features that correlate with mechanisms that trigger, facilitate, or hinder CO2 catalytic reactions. [ABSTRACT FROM AUTHOR]

Titanium dioxide (TiO2) nanoclusters (NCs) and nanoparticles (NPs) have been the focus of intense research in recent years since they play a prominent role in photocatalysis. In particular, the properties of their excited states determine the photocatalytic activity. Among the requirements for photocatalytic activity, low excitation energy and large separation of the charge carriers are crucial. While information regarding the first is straightforward from either experiment or theory, the information regarding the second is scarce or missing. In the present work we fill this gap through a topological analysis of the first singlet excited state of a series of TiO2 NCs, and anatase and rutile derived NPs containing up to 495 atoms. The excited states of all these systems in vacuo have been obtained from time-dependent density functional theory (TDDFT) calculations using hybrid functionals and the influence of water was taken into account through a continuum model. Three different topological descriptors based on the attachment/detachment one-electron charge density, are scrutinized: (i) charge transfer degree, (ii) charge density overlap, and (iii) distance between centroids of charge. The present analysis shows that the charge separation in the excited state strongly depends on the NP size and shape. The character of the electronic excitations, as arising from the analysis of the canonical Kohn–Sham molecular orbitals (MOs) or from natural transition orbitals (NTOs), is also investigated. The understanding and prediction of charge transfer and recombination in TiO2 nanostructures may have implications in the rational design of these systems to boost their photocatalytic potential. [ABSTRACT FROM AUTHOR]

A database containing 17 multiplicity-changing valence and Rydberg excitation energies of p-block elements is used to test the performance of density functional theory (DFT) with approximate density functionals for calculating relative energies of spin states. We consider only systems where both the low-spin and high-spin state are well described by a single Slater determinant, thereby avoiding complications due to broken-symmetry solutions. Because the excitations studied involve a spin change, they require a balanced treatment of exchange and correlation, thus providing a hard test for approximate density functionals. We test three formalisms for predicting the multiplicity-changing transition energies. First is the ΔSCF method; we also test time-dependent density functional theory (TDDFT), both in its conventional form starting from the low-spin state and in its collinear spin-flip form starting from the high-spin state. Very diffuse basis functions are needed to give a qualitatively correct description of the Rydberg excitations. The scalar relativistic effect needs to be considered when quantitative results are desired, and we include it in the comparisons. With the ΔSCF method, most of the tested functionals give mean unsigned errors (MUEs) larger than 6 kcal/mol for valence excitations and MUEs larger than 3 kcal/mol for Rydberg excitations, but the performance for the Rydberg states is much better than can be obtained with time-dependent DFT. It is surprising to see that the long-range corrected functionals, which have 100% Hartree-Fock exchange at large inter-electronic distance, do not improve the performance for Rydberg excitations. Among all tested density functionals, ΔSCF calculations with the O3LYP, M08-HX, and OLYP functionals give the best overall performance for both valence and Rydberg excitations, with MUEs of 2.1, 2.6, and 2.7 kcal/mol, respectively. This is very encouraging since the MUE of the CCSD(T) coupled cluster method with quintuple zeta basis sets is 2.0 kcal/mol; however, caution is advised since many popular density functionals give poor results, and there can be very significant differences between the ΔSCF predictions and those from TDDFT. [ABSTRACT FROM AUTHOR]

The adsorption of small molecules such as NO or CO on surfaces of magnetic oxides containing transition metals is difficult to model by current density functional approximations. Two such oxides are NiO(100) and Ni-doped MgO(100). Here we compare the results of a theoretical model of the Ni-doped MgO(100) surface with experimental results on NiO(100), which introduces some uncertainty into a quantitative theory-experiment comparison. In the present work, we tested seven meta-GGA and hybrid metafunctionals, in particular, three developed by the Minnesota group (M05, M06-L, and M06), and TPSS, TPSSh, TPSSKCIS, and B1B95; six GGA functionals, including BP86, PBE, and four other functionals that are modifications of PBE (PBEsol, SOGGA, revPBE, and RPBE); five hybrid GGA functionals (B3LYP, PBE0, B97–2, B97–3, and MPWLYP1M); and one unconventional functional of the generalized gradient type with scaled correlation called MOHLYP. The Minnesota meta-GGA functionals were found in the past to be very good choices when transition metal atoms were present; the other functionals chosen are a selection from the most currently used and most promising sets of functionals for bulk solids and surfaces and for transition metals. The difficulty is due to the charge transfer between open shells in the case of NO and to the weak character of the interaction in the case of CO. It is shown that the M06 hybrid meta functional applied to NO or CO on a model of the Ni-doped MgO(100) surface is able to provide a good description of both adsorbate geometries and binding energies. The M06 vibrational frequency shifts are more accurate than for other functionals, but there is still room for improvement. [ABSTRACT FROM AUTHOR]

As gold clusters increase in size, the preferred structure changes from planar to three-dimensional and, for anionic clusters, Aun-, the two-dimensional(2D)→three-dimensional (3D) transition is found experimentally to occur between n=11 and n=12. Most density functionals predict that planar structures are preferred up to higher n than is observed experimentally, an exception being the local spin density approximation. Here we test four relatively new functionals for this feature, in particular, M05, M06-L, M06, and SOGGA. We find that M06-L, M06, and SOGGA all predict the 2D→3D transition at the correct value of n. Since the M06-L and M06 functionals have previously been shown to be reasonably accurate for transition metal bond energies, main group atomization energies, barrier heights, and noncovalent interaction energies, and, since they are here shown to perform well for the s-d excitation energy and ionization potential of Au atoms and for the size of Aun- clusters at which the 2D→3D transition occurs, they are recommended for simulating processes catalyzed by gold clusters. [ABSTRACT FROM AUTHOR]

Using previously developed potential energy surfaces and their couplings, non-Born–Oppenheimer trajectory methods are used to study the state-selected photodissociation of ammonia, prepared with up to six quanta of vibrational excitation in the symmetric (ν1) or antisymmetric (ν3) stretching modes of NH3(Ã). The predicted dynamics is mainly electronically nonadiabatic (that is, it produces ground electronic state amino radicals). The small probability of forming the excited-state amino radical is found, for low excitations, to increase with total energy and to be independent of whether the symmetric or antisymmetric stretch is excited; however some selectivity with respect to exciting the antisymmetric stretch is found when more than one quantum of excitation is added to the stretches, and more than 50% of the amino radical are found to be electronically excited when six quanta are placed in the antisymmetric stretch. These results are in contrast to the mechanism inferred in recent experimental work, where excitation of the antisymmetric stretch by a single quantum was found to produce significant amounts of excited-state products via adiabatic dissociation at total energies of about 7.0 eV. Both theory and experiment predict a broad range of translational energies for the departing H atoms when the symmetric stretch is excited, but the present simulations do not reproduce the experimental translational energy profiles when the antisymmetric stretch is excited. The sensitivity of the predicted results to several aspects of the calculation is considered in detail, and the analysis leads to insight into the nature of the dynamics that is responsible for mode selectivity. [ABSTRACT FROM AUTHOR]

The adsorption of CO on Mg(001) constitutes a challenge for current density functional approximations because of its weak interaction character. In the present work we show that the M06-2X and M06-HF exchange-correlation functionals are the first ones to provide a simultaneously satisfactory description of adsorbate geometry, vibrational frequency shift, and adsorption energy of CO on MgO(001). For a sufficiently large embedded cluster model, the three functionals of the M06 family—which contain a nonzero percentage of Hartree–Fock exchange (M06, M06-2X, and M06-HF)—all predict positive C–O vibrational shifts, in agreement with the experimental findings, while the local M06-L functional gives large negative shifts. Moreover, the shifts computed with the M06-2X and M06-HF potentials are in good agreement with the experimental shift of +14 cm-1. The interaction energy (De) calculated with M06-2X and M06-HF is ∼6.0 kcal/mol, which agrees well with the De value (∼4 kcal/mol) deduced from the D0 obtained in thermal desorption measurements on single-crystal surfaces. [ABSTRACT FROM AUTHOR]

The performance of the M06 family of exchange-correlation potentials for describing the electronic structure and the Heisenberg magnetic coupling constant (J) is investigated using a set of representative open-shell systems involving two unpaired electrons. The set of molecular systems studied has well defined structures, and their magnetic coupling values are known experimentally. As a general trend, the M06 functional is about equally as accurate as B3LYP or PBE0. The performance of local functionals is important because of their economy and convenience for large-scale calculations; we find that M06-L local functional of the M06 family largely improves over the local spin density approximation and the generalized gradient approximation. [ABSTRACT FROM AUTHOR]

Bromoacetyl chloride photodissociation has been interpreted as a paradigmatic example of a process in which nonadiabatic effects play a major role. In molecular beam experiments by Butler and co-workers [J. Chem. Phys. 95, 3848 (1991); J. Chem. Phys. 97, 355 (1992)], BrCH2C(O)Cl was prepared in its ground electronic state (S0) and excited with a laser at 248 nm to its first excited singlet state (S1). The two main ensuing photoreactions are the ruptures of the C–Cl bond and of the C–Br bond. A nonadiabatic model was proposed in which the C–Br scission is strongly suppressed due to nonadiabatic recrossing at the barrier formed by the avoided crossing between the S1 and S2 states. Recent reduced-dimensional dynamical studies lend support to this model. However, another interpretation that has been given for the experimental results is that the reduced probability of C–Br scission is a consequence of incomplete intramolecular energy redistribution. To provide further insight into this problem, we have studied the energetically lowest six singlet electronic states of bromoacetyl chloride by using an ab initio multiconfigurational perturbative electronic structure method. Stationary points (minima and saddle points) and minimum energy paths have been characterized on the S0 and S1 potential energy surfaces. The fourfold way diabatization method has been applied to transform five adiabatic excited electronic states to a diabatic representation. The diabatic potential energy matrix of the first five excited singlet states has been constructed along several cuts of the potential energy hypersurfaces. The thermochemistry of the photodissociation reactions and a comparison with experimental translational energy distributions strongly suggest that nonadiabatic effects dominate the C–Br scission, but that the reaction proceeds along the energetically allowed diabatic pathway to excited-state products instead of being nonadiabatically suppressed. This conclusion is also supported by the low values of the diabatic couplings on the C–Br scission reaction path. The methodology established in the present study will be used for the construction of global potential energy surfaces suitable for multidimensional dynamics simulations to test these preliminary interpretations. [ABSTRACT FROM AUTHOR]

High probabilities of energy transfer from translation to molecular rotations are observed in the scattering of n-D2 from LiF(001) at an incident beam energy of 85.3 meV. For the <100> incidence direction, close-coupling calculations yield ratios of the rotationally inelastic (j=0→2) and (j=1→3) peaks to the rotationally elastic specular peaks (G=0) that are in reasonable agreement with experiment, as are the ratios of the rotationally elastic diffraction peak intensities to the specular peak intensities. The agreement between theory and experiment is also quite good for the rotationally inelastic diffractive (-1-1) transitions for (j=1→3), but rather poor for (j=0→2). The calculations show that the interaction between the electrostatic field of the surface ions and the quadrupole moment of the D2 molecule efficiently promotes the (j=0→2) and (j=1→3) transitions. If this electrostatic interaction is excluded from the potential model, the ratios of the (j=0→2) and (j=1→3) rotationally inelastic peaks to the corresponding specular peaks show a large discrepancy with experiment, underlining the importance of this interaction. The close-coupling calculations show a somewhat worse agreement with experiment for the <110> incidence direction. In particular, the sharp peaks observed experimentally in the ratios of the peak intensities of the rotationally inelastic G=0 (j=0→2) and (j=1→3) to the rotationally elastic G=0 transitions as a function of incident angle are not reproduced by the calculations. The theoretical ratios of the peak intensities of the rotationally elastic diffraction to G=0 transitions are shifted to lower incidence angles with respect to experiment. The rotationally inelastic diffractive (-10) transitions present an interesting resonance phenomenon for the (j=0→2) rotational transition. This resonance is predicted by both theory and experiment, although at rather different incident angles. [ABSTRACT FROM AUTHOR]

Full- [six-dimensional (6-D)] and reduced-dimensional [five-dimensional (5-D)] quantum wave packet calculations have been performed for the title reaction to obtain reaction probabilities deriving from the ground rovibrational states of OH and CO with total angular momentum J=0. Three potential energy surfaces (PES) are studied, namely, those of Bradley and Schatz (BS), Yu, Muckerman, and Sears (YMS), and Lakin, Troya, Schatz, and Harding (LTSH). 6-D calculations are performed only for the BS PES, while 5-D results are reported for all three PES’. The 6-D results obtained in the present work improve on those previously reported, since a larger vibrational basis and a better representation of the OH and CO bonds has been introduced. In particular, we now employ a generalized Lanczos–Morse discrete variable representation for both the OH and CO vibrations. In a further improvement, the generalized discrete variable representation of the CO vibration is based on different CO intramolecular potentials for the asymptotic and product grids employed in our projection formalism. This new treatment of the vibrational bases allows for a large reduction in computation time with respect to our previous implementation of the wave packet method, for a given level of accuracy. As a result, we have been able to extend the range of collision energies for which we can obtain converged 6-D results to a higher energy (0.8 eV) than was possible before (0.5 eV). The comparison of the new 6-D and previous 5-D results for the BS PES shows good agreement of the general trend in the reaction probabilities over all collision energies considered (0.1–0.8 eV), while our previous 6-D calculation showed reaction probabilities that differed from the 5-D results by up to 10% between 0.5 and 0.8 eV. The 5-D reaction probabilities reveal interesting trends for the different PES’. In particular, at low energies (<0.2 eV) the LTSH PES gives rise to much larger reactivity than the other PES’, while at high energies (>0.3 eV) its reaction probability decreases with respect to the BS and YMS PES’, being more than a factor of 2 smaller at 0.8 eV. A 5-D calculation on a modified version of the LTSH surface shows that the van der Waals interaction in the entrance channel, which is not correctly described in the other PES’ is largely responsible for its larger reactivity at low energies. The large difference between the 5-D reaction probabilities for the YMS and LTSH PES’ serves to emphasize the importance of the van der Waals interaction for the reactivity at low energies, because most of the stationary point energies on the YMS and LTSH PES are rather similar, being in line with high-level ab initio information. © 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

Variational transition state theory, within the canonical unified statistical model including multidimensional tunneling corrections, has been applied to derive thermal rate constants in the temperature range 80-2800 K for the title reaction and its deuterated counterpart in the low- and high-pressure limits. The analytical potential energy surfaces of Bradley and Schatz (BS) and of Yu, Muckerman, and Sears have been employed. In the low-pressure limit, the results are comparable though slightly better for the BS potential surface at the lowest temperatures. In the high-pressure limit, the BS potential surface is notably closer to the experimental data, which are only available at temperatures close to and higher than 300 K. The kinetic isotope effect (KIE) at room temperature is better predicted by the BS surface, although both surfaces fail to reproduce the experimental dependence of the KIE on temperature. [ABSTRACT FROM AUTHOR]

A report is presented on the study of the photodissociation dynamics of the prototypical HBr system using semiclassical trajectory methods. The Li[Br.sup.+] model system is used to show that accurate spin-coupled potential curves can be constructed for odd-electron system

An extension of the fourfold way to construct diabatic representations suitable for spin-coupled systems reported. A new method is formulated for the case of even-electron systems that yield pairs of fragments with doublet spin multiplicity.

The optical gap (Ogap) of a set of (TiO2)n nanoclusters and nanoparticles with n = 10–563 and different morphologies such as spherical, octahedral, lamellar, or tubular finite structures is investigated based on a relativistic all-electron description along with a numerical atomic centered orbital basis set. Two different functionals are used, PBE and PBEx, the former corresponds to a standard implementation of the generalized gradient approximation (GGA) and the latter to a hybrid functional with 12.5% of Fock exchange which reproduces the band gap of bulk TiO2 anatase and rutile. It is shown that the inclusion of exchange Fock in the PBE functional promotes a systematic energy gap opening of 1.25 eV relative to the PBE values. Remarkably, a linear correlation is found between PBEx and PBE Ogap calculated values with concomitant similar correlations for the HOMO and LUMO orbital energies. However, it appears that PBEx induces a larger downshift on the HOMO orbital than the upshift observed on the LUMO one. The fact that the PBEx hybrid functional was shown to reproduce the experimental energy gaps of stoichiometric and reduced TiO2 bulk phases leads to a suitable and practical way to successfully estimate Ogap of TiO2 nanoparticles containing up to thousands of atoms with PBEx precision from computationally affordable PBE calculations. [ABSTRACT FROM AUTHOR]