Data-driven stabilization of NimPdn-m nanoalloys: a study using density functional theory and data mining approaches.

Green hydrogen, generated through the electrolysis of water, is a viable alternative to fossil fuels, although its adoption is hindered by the high costs associated with the catalysts. Among a wide variety of potential materials, binary nickel-palladium (NiPd) systems have garnered significant attention, particularly at the nanoscale, for their efficacious roles in catalyzing hydrogen and oxygen evolution reactions. However, our atom-level understanding of the descriptors that drive their energetic stability at the nanoscale remains largely incomplete. Here, we investigate by density functional theory calculations the descriptors that drives the stability of the NimPdn-m clusters for different sizes (n = 13, 27, 41) and compositions. To achieve our goals, a large number of trial configurations were generated and selected using data mining algorithms (k-means, t-SNE) and genetic algorithms, while the most important physical-chemical descriptors were identified using Spearman correlation analysis. We have found that core-shell formation, with the smaller Ni atoms lying in the center of the particle, plays a major role in the stabilization of the nanoalloys, and this effect causes the alloys to assume a icosahedral-fragment configuration (as the unary nickel cluster) instead of a fcc fragment (as the unary palladium cluster). However, the core-shell formation in this alloy is unique in that Pd poor compositions exhibit scattered Pd atoms on the surface. As the palladium content increases, this gives rise to the complete Pd shell. This stabilization mechanism is quantitatively supported by the different correlations observed in the number of Ni-Ni and Pd-Pd bonds with energy, in which the latter tends to decrease alloy stability. Furthermore, a notable trend is the correlation between the coordination number of Ni atoms with alloy stabilization, while the coordination of Pd atoms shows an inverse correlation.

Correction: Quantum and semiclassical studies of nonadiabatic electronic transitions between N(4S) and N(2D) by collisions with N2.

Correction for 'Quantum and semiclassical studies of nonadiabatic electronic transitions between N(4S) and N(2D) by collisions with N2' by Dandan Lu et al., Phys. Chem. Chem. Phys., 2023, 25, 15656-15665, https://doi.org/10.1039/D3CP01429K.

Computational screening of silver-based single-atom alloys catalysts for CO2 reduction.

Electrocatalytically reducing CO2 into value-added products is a challenging but promising process. Catalysts have been proposed to reduce the potential necessary for the reaction to occur, among which single-atom alloys (SAAs) are particularly promising. Here, we employ density functional theory calculations and the computational electrode model to predict whether silver-based SAAs have the potential to be effective electrocatalysts to convert CO2 into C1 products. We take into account surface defects by using the Ag(211) surface as a model. We also verify whether the proposed materials are prone to OH poisoning or enhance the competing hydrogen evolution reaction. Our calculations predict that these materials show weak mixing between the host and the dopant, characterized by a sharp peak in the density of states near the Fermi energy, except when copper (also a coinage metal) is used as the dopant. This affects the adsorption energy of the different intermediate molecules, yielding different reaction profiles for each substrate. As non-doped silver, copper-doped SAA tends to spontaneously desorb carbon monoxide (CO) instead of proceeding with its reduction. Other elements of the fourth period (Fe, Co, and Ni) tend to bind to the CO molecule but do not favor more reduced products. These metals also tend to enhance the hydrogen evolution reaction. On the contrary, we show that the Ir and Rh dopants have significant potential as electrocatalysts, which favors the reduction of CO over its desorption while also suppressing the hydrogen evolution reaction at potentials lower than those required by copper. They have also been shown to not be prone to poisoning by OH radicals.

Theoretical study of the O(3P) + CN(X2Σ+) → CO(X1Σ+) + N(2D)/N(4S) reactions.

The barrierless exothermic reactions between atomic oxygen and the cyano radical, O(3P) + CN(X2Σ+) → CO(X1Σ+) + N(2D)/N(4S), play a significant role in combustion, astrochemistry, and hypersonic environments. In this work, their dynamics and kinetics are investigated using both wave packet (WP) and quasi-classical trajectory (QCT) methods on recently developed potential energy surfaces of the 12A', 12A,″ and 14A″ states. The product state distributions in the doublet pathway obtained with the WP method for a few partial waves show extensive internal excitation in the CO product. This observation, combined with highly oscillatory reaction probabilities, signals a complex-forming mechanism. The statistical nature of the reaction is confirmed by comparing the WP results with those from phase space theory. The calculated rate coefficients using the WP (with a J-shifting approximation) and QCT methods exhibit agreement with each other near room temperature, 1.77 × 10-10 and 1.31 × 10-10 cm3 molecule-1 s-1, but both are higher than the existing experimental results. The contribution of the quartet pathway is small at room temperature due to a small entrance channel bottleneck. The QCT rate coefficients are further compared with experimental results above 3000 K, and the agreement is excellent.

High-accuracy DMBE potential energy surface for CNO(A''4) and the rate coefficients for the C + NO reaction in the A'2, A''2, and A''4 states.

An accurate potential energy surface (PES) for the lowest lying A''4 state of the CNO system is presented based on explicitly correlated multi-reference configuration interaction calculations with quadruple zeta basis set (MRCI-F12/cc-pVQZ-F12). The ab initio energies are fitted using the double many-body expansion method, thus incorporating long-range energy terms that can accurately describe the electrostatic and dispersion interactions with physically motivated decaying functions. Together with the previously fitted lowest A'2 and A''2 states using the same theoretical framework, this constitutes a new set of PESs that are suitable to predict rate coefficients for all atom-diatom reactions of the CNO system. We use this set of PESs to calculate thermal rate coefficients for the C(P3) + NO(Π2) reaction and compare the temperature dependence and product branching ratios with experimental results. The comparison between theory and experiment is shown to be improved over previous theoretical studies. We highlight the importance of the long-range interactions for low-temperature rate coefficients.

Underlying mechanisms of gold nanoalloys stabilization.

Gold nanoclusters have attracted significant attention due to their unique physical-chemical properties, which can be tuned by alloying with elements such as Cu, Pd, Ag, and Pt to design materials for various applications. Although Au-nanoalloys have promising applications, our atomistic understanding of the descriptors that drive their stability is far from satisfactory. To address this problem, we considered 55-atom model nanoalloys that have been synthesized by experimental techniques. Here, we combined data mining techniques for creating a large sample of representative configurations, density functional theory for performing total energy optimizations, and Spearman correlation analyses to identify the most important descriptors. Among our results, we have identified trends in core-shell formation in the AuCu and AuPd systems and an onion-like design in the AuAg system, characterized by the aggregation of gold atoms on nanocluster surfaces. These features are explained by Au's surface energy, packing efficiency, and charge transfer mechanisms, which are enhanced by the alloys' preference for adopting the structure of the alloying metal rather than the low-symmetry one presented by Au55. These generalizations provide insights into the interplay between electronic and structural properties in gold nanoalloys, contributing to the understanding of their stabilization mechanisms and potential applications in various fields.

Quantum and semiclassical studies of nonadiabatic electronic transitions between N(4S) and N(2D) by collisions with N2.

The dynamics and kinetics of spin-forbidden transitions between N(2D) and N(4S) via collisions with N2 molecules are investigated using a quantum wave packet (WP) method and the semi-classical coherent switches with decay of mixing (CSDM) method. These electronic transition processes are competing with exchange reaction channels on both the doublet and quartet potential energy surfaces. The WP and CSDM quenching rate coefficients are found in reasonable agreement with each other, and both reproduce the previous theoretical results. For the excitation process, the agreement between the two approaches is dependent on the treatment of the zero-point energy (ZPE) in the product, because the high endoergicity of this process leads to severe violation of the vibrational ZPE. The Gaussian-binning (GB) method is found to improve the agreement with the quantum result. The excitation rate coefficients are found to be two orders of magnitude smaller than that of the adiabatic exchange reaction, underscoring the inefficient intersystem crossing due to the weak spin-orbit coupling between the two spin manifolds of the N3 system.

Study of the Counter Cation Effects on the Supramolecular Structure and Electronic Properties of a Dianionic Oxamate-Based {NiII2} Helicate.

Herein, we describe the synthesis, crystal structure, and electronic properties of {[K2(dmso)(H2O)5][Ni2(H2mpba)3]·dmso·2H2O}n (1) and [Ni(H2O)6][Ni2(H2mpba)3]·3CH3OH·4H2O (2) [dmso = dimethyl sulfoxide; CH3OH = methanol; and H4mpba = 1,3-phenylenebis(oxamic acid)] bearing the [Ni2(H2mpba)3]2- helicate, hereafter referred to as {NiII2}. SHAPE software calculations indicate that the coordination geometry of all the NiII atoms in 1 and 2 is a distorted octahedron (Oh) whereas the coordination environments for K1 and K2 atoms in 1 are Snub disphenoid J84 (D2d) and distorted octahedron (Oh), respectively. The {NiII2} helicate in 1 is connected by K+ counter cations yielding a 2D coordination network with sql topology. In contrast to 1, the electroneutrality of the triple-stranded [Ni2(H2mpba)3] 2- dinuclear motif in 2 is achieved by a [Ni(H2O)6]2+ complex cation, where the three neighboring {NiII2} units interact in a supramolecular fashion through four R22(10) homosynthons yielding a 2D array. Voltammetric measurements reveal that both compounds are redox active (with the NiII/NiI pair being mediated by OH- ions) but with differences in formal potentials that reflect changes in the energy levels of molecular orbitals. The NiII ions from the helicate and the counter-ion (complex cation) in 2 can be reversibly reduced, resulting in the highest faradaic current intensities. The redox reactions in 1 also occur in an alkaline medium but at higher formal potentials. The connection of the helicate with the K+ counter cation has an impact on the energy levels of the molecular orbitals; this experimental behavior was further supported by X-ray absorption near-edge spectroscopy (XANES) experiments and computational calculations.

The role of single-atom Rh-dopants in the adsorption properties of OH and CO on stepped Ag(211) surfaces.

Several chemical reactions with commercial and environmental importance can benefit from the development of more active or selective heterogeneous catalysts. Particularly, those catalyzed by metallic surfaces are usually impacted by the presence of defects such as kinks and dopants. Here, we employed density functional theory calculations within van der Waals correction to investigate the effects of single-atom Rh-dopants in the adsorption properties of OH and CO on stepped Ag(211) surfaces. From our calculations and analyses, we found that the dopant is more energetically stable when replacing more coordinated (and less exposed to the vacuum) sites of the surface. However, in the presence of both molecules, this trend is inverted, and the dopant is more stable in the least coordinated site (step). While OH presents high adsorption energies on both doped and non-doped silver surfaces, CO binds weakly to the noble metal, and strongly on doped sites. The results are relevant for understanding single-atom catalysts on noble-metal surfaces, where the difference in selectivity and activity between the host metal and dopants is exploited. The charge redistribution caused by the dopant, and the appearance of a sharp peak in the density of states of the surface are used to rationalize the results and provide insights into the interactions involved in the adsorption of both molecules.

Formation of phosphorus monoxide through the [Formula: see text] reaction.

Phosphorus is a key and vital element for a diverse set of important biological molecules, being indispensable for life as we know. A deeper comprehension of its role in astrochemistry and atmospheric chemistry may aid in finding answers to how this element became available on Earth. The PO molecule is one of the main reservoirs of phosphorus in the interstellar medium (ISM), and a better understanding of the mechanisms and rate coefficients for its formation in the ISM is important for modelling its abundances. In this work, we perform multireference configuration interaction calculations on the formation of PO via the [Formula: see text] reaction, analyzing its potential energy surface and rate coefficients for the global reaction on both doublet and quartet states. We also perform DFT (M06-2X) and CCSD(T) calculations, in order to compare the results. We found that the OPO system possesses a high multiconfigurational character, making DFT and CCSD methodologies not suitable for its potential energy landscape calculation. The rate coefficients have been calculated using the master equation system solver (MESS) package, and the results compared to recent experimental data. It is shown that the quartet state contributes for temperatures higher than 700K. The computed rate coefficient can be described by a modified Arrhenius equation [[Formula: see text]] with [Formula: see text], [Formula: see text] and [Formula: see text] K.

Gas-Phase Preparation of Silyl Cyanide (SiH3CN) via a Radical Substitution Mechanism.

The silyl cyanide (SiH3CN) molecule, the simplest representative of a fully saturated silacyanide, was prepared in the gas phase under single-collision conditions via a radical substitution mechanism. The chemical dynamics were direct and revealed a pronounced backward scattering as a consequence of a transition state with a pentacoordinated silicon atom and almost colinear geometry of the attacking cyano radical and leaving hydrogen. Compared to the isovalent cyano (CN)-methane (CH4) system, the CN-SiH4 system dramatically reduces the energy of the transition state to silyl cyanide by nearly 100 kJ mol-1, which reveals a profound effect on the chemical bonding and reaction mechanism. In extreme high-temperature environments including circumstellar envelopes of IRC +10216, this versatile radical substitution mechanism may synthesize organosilicon molecules via reactions of silane with doublet radicals. Overall, this study provides rare insights into the exotic reaction mechanisms of main-group XIV elements in extreme environments and affords deeper insights into fundamental molecular mass growth processes involving silicon in our universe.

Directed gas-phase preparation of the elusive phosphinosilylidyne (SiPH2, X2A'') and cis/trans phosphinidenesilyl (HSiPH; X2A') radicals under single-collision conditions.

The reaction of the D1-silylidyne radical (SiD; X2Π) with phosphine (PH3; X1A1) was conducted in a crossed molecular beams machine under single collision conditions. Merging of the experimental results with ab initio electronic structure and statistical Rice-Ramsperger-Kassel-Marcus (RRKM) calculations indicates that the reaction is initiated by the barrierless formation of a van der Waals complex (i0) as well as intermediate (i1) formed via the barrierless addition of the SiD radical with its silicon atom to the non-bonding electron pair of phosphorus of the phosphine. Hydrogen shifts from the phosphorous atom to the adjacent silicon atom yield intermediates i2a, i2b, i3; unimolecular decomposition of these intermediates leads eventually to the formation of trans/cis-phosphinidenesilyl (HSiPH, p2/p4) and phosphinosilylidyne (SiPH2, p3) via hydrogen deuteride (HD) loss (experiment: 80 ± 11%, RRKM: 68.7%) and d-trans/cis-phosphinidenesilyl (DSiPH, p2'/p4') plus molecular hydrogen (H2) (experiment: 20 ± 7%, RRKM: 31.3%) through indirect scattering dynamics via tight exit transition states. Overall, the study reveals branching ratios of p2/p4/p2'/p4' (trans/cis HSiPH/DSiPH) to p3 (SiPH2) of close to 4 : 1. The present study sheds light on the complex reaction dynamics of the silicon and phosphorous systems involving multiple atomic hydrogen migrations and tight exit transition states, thus opening up a versatile path to access the previously elusive phosphinidenesilyl and phosphinosilylidyne doublet radicals, which represent potential targets of future astronomical searches toward cold molecular clouds (TMC-1), star forming regions (Sgr(B2)), and circumstellar envelopes of carbon rich stars (IRC + 10216).

Nonadiabatic reaction dynamics to silicon monosulfide (SiS): A key molecular building block to sulfur-rich interstellar grains.

Sulfur- and silicon-containing molecules are omnipresent in interstellar and circumstellar environments, but their elementary formation mechanisms have been obscure. These routes are of vital significance in starting a chain of chemical reactions ultimately forming (organo) sulfur molecules-among them precursors to sulfur-bearing amino acids and grains. Here, we expose via laboratory experiments, computations, and astrochemical modeling that the silicon-sulfur chemistry can be initiated through the gas-phase reaction of atomic silicon with hydrogen sulfide leading to silicon monosulfide (SiS) via nonadiabatic reaction dynamics. The facile pathway to the simplest silicon and sulfur diatomic provides compelling evidence for the origin of silicon monosulfide in star-forming regions and aids our understanding of the nonadiabatic reaction dynamics, which control the outcome of the gas-phase formation in deep space, thus expanding our view about the life cycle of sulfur in the galaxy.

A Crossed Molecular Beams and Computational Study of the Formation of the Astronomically Elusive Thiosilaformyl Radical (HSiS, X2A').

The formation pathways to silicon- and sulfur-containing molecules are crucial to the understanding of silicon-sulfur chemistry in interstellar and circumstellar environments. While multiple silicon- and sulfur-containing species have been observed in deep space, their fundamental formation mechanisms are largely unknown. The crossed molecular beams technique combined with electronic structure and Rice-Ramsperger-Kassel-Marcus (RRKM) calculations was utilized to study the bimolecular reaction of atomic silicon (Si(3Pj)) with thiomethanol (CH3SH, X1A') leading to the thiosilaformyl radical (HSiS, X2A') via an exclusive methyl radical (CH3, X2A2″) loss via indirect scattering dynamics which involves barrierless addition and hydrogen migration in an overall exoergic reaction, indicating the possibility that HSiS can form in cold molecular clouds. The astronomically elusive thiosilaformyl radical may act as a tracer of an exotic silicon-sulfur chemistry to be deciphered toward, for example, the star-forming region SgrB2, thus leading to a better understanding of the formation of silicon-sulfur bonds in deep space.

Chemical dynamics study on the gas-phase reaction of the D1-silylidyne radical (SiD; X2Π) with deuterium sulfide (D2S) and hydrogen sulfide (H2S).

The reactions of the D1-silylidyne radical (SiD; X2Π) with deuterium sulfide (D2S; X1A1) and hydrogen sulfide (H2S; X1A1) were conducted utilizing a crossed molecular beams machine under single collision conditions. The experimental work was carried out in conjunction with electronic structure calculations. The elementary reaction commences with a barrierless addition of the D1-silylidyne radical to one of the non-bonding electron pairs of the sulfur atom of hydrogen (deuterium) sulfide followed by possible bond rotation isomerization and multiple atomic hydrogen (deuterium) migrations. Unimolecular decomposition of the reaction intermediates lead eventually to the D1-thiosilaformyl radical (DSiS) (p1) and D2-silanethione (D2SiS) (p3) via molecular and atomic deuterium loss channels (SiD-D2S system) along with the D1-thiosilaformyl radical (DSiS) (p1) and D1-silanethione (HDSiS) (p3) through molecular and atomic hydrogen ejection (SiD-H2S system) via indirect scattering dynamics in barrierless and overall exoergic reactions. Our study provides a look into the complex dynamics of the silicon and sulfur chemistries involving multiple deuterium/hydrogen shifts and tight exit transition states, as well as insight into silicon- and sulfur-containing molecule formation pathways in deep space. Although neither of the non-deuterated species - the thiosilaformyl radical (HSiS) and silanethione (H2SiS) - have been observed in the interstellar medium (ISM) thus far, astrochemical models presented here predict relative abundances in the Orion Kleinmann-Low nebula to be sufficiently high enough for detection.

Non-Adiabatic Reaction Dynamics in the Gas-Phase Formation of Phosphinidenesilylene, the Isovalent Counterpart of Hydrogen Isocyanide, under Single-Collision Conditions.

The phosphinidenesilylene (HPSi; X1A') molecule was prepared via a directed gas-phase synthesis in the bimolecular reaction of ground-state atomic silicon (Si; 3P) with phosphine (PH3; X1A1) under single-collision conditions. The chemical dynamics are initiated on the triplet surface via addition of a silicon atom to the non-bonding electron pair of phosphine, followed by non-adiabatic dynamics and surface hopping to the singlet manifold, accompanied by isomerization via atomic hydrogen shift and decomposition to phosphinidenesilylene (HPSi, X1A') along with molecular hydrogen. Statistical calculations predict that silylidynephosphine (HSiP, X1Σ+) is also formed, albeit with lower yields. The barrier-less route to phosphinidenesilylene opens up a multipurpose mechanism to access the hitherto obscure class of phosphasilenylidenes through silicon-phosphorus coupling via reactions of atomic silicon with alkylphosphines under single-collision conditions in the absence of successive reactions of the reaction products, which are not feasible to prepare by traditional synthetic routes.

Reliability of semiempirical and DFTB methods for the global optimization of the structures of nanoclusters.

In this work, we explore the possibility of using computationally inexpensive electronic structure methods, such as semiempirical and DFTB calculations, for the search of the global minimum (GM) structure of chemical systems. The basic prerequisite that these inexpensive methods will need to fulfill is that their lowest energy structures can be used as starting point for a subsequent local optimization at a benchmark level that will yield its GM. If this is possible, one could bypass the global optimization at the expensive method, which is currently impossible except for very small molecules. Specifically, we test our methods with clusters of second row elements including systems of several bonding types, such as alkali, metal, and covalent clusters. The results reveal that the DFTB3 method yields reasonable results and is a potential candidate for this type of applications. Even though the DFTB2 approach using standard parameters is proven to yield poor results, we show that a re-parametrization of only its repulsive part is enough to achieve excellent results, even when applied to larger systems outside the training set.

An ab initio investigation of the adsorption properties of water on binary AlSi clusters.

The potential of doped aluminium clusters as catalysts for the water splitting reaction has attracted considerable scientific effort, however, the water-cluster interactions, which are a key step in the overall mechanism, are not fully understood. Here, we report an ab initio investigation of water adsorption on AlSi clusters at the MP2 level to elucidate the bonding and structural properties employing unary and binary 8- and 13-atom clusters, namely, Si8, Al2Si6, Al4Si4, Al8, Si13, Al2Si11, Al12Si, and Al13, which were selected by their relevance and energetic stability. We found that H2O binds via the O atom near to the on-top sites of the Si or Al atoms; in particular, there is a strong preference for the Al sites on the binary AlSi clusters, which is supported by the strong adsorption energy. Furthermore, we found a large enhancement of the adsorption energy on the Al2Si6 and Al2Si11 clusters, which can be explained by the cationic character of the Al site, which increases the Coulomb contribution to the Al+-O- interaction.

What Electronic Structure Method Can Be Used in the Global Optimization of Nanoclusters?

The problem of obtaining the spatial structure of nanoclusters is known to be very difficult due to the large number of local minima associated with their potential energy surfaces (isomers). In global optimization approaches, such as basin hopping and genetic algorithms, the problem is normally tackled by first using a low-level and affordable method to evaluate the energy. Afterward, the putative global minimum (and often a few others) is refined with calculations using higher level methods and larger basis sets. There is no guarantee, however, that the structure obtained at the lower level method will be the global minimum at the refined one. In this work, we have performed benchmark coupled cluster calculations at the complete basis set limit for a large number of different isomers of representative clusters of third row elements. Such calculations are then employed to check the hypothesis that lower level methods can be used in the global optimization with reliable results. For this, we have developed a methodology that allows us to compare a large number of minima obtained at different calculation levels. The results indicate that, if the global optimization is capable of reaching not only the global minimum but also a reduced number of low lying structures, most of the tested density functional theory (DFT) functionals are good choices, with emphasis on TPSSh. Besides giving a more solid ground to this commonly used approach, this work helps guiding such global optimizations. The use of the MP2 method and several scaled variants is also assessed, from where it is concluded that the scaled variants yield better results than standard MP2 or DFT approaches, except for one system where a large number of van der Waals structures exist.

Emerging contaminants removal by granular activated carbon obtained from residual Macauba biomass.

The removal of emergent contaminants via adsorption on granular activated carbon, prepared from Macauba palm, has been studied, contributing to the recovery of the residual biomass, endocarp, obtained in the Macauba palm oil extraction process. The material was characterized by different techniques, such as Raman spectroscopy, thermal analysis, adsorption/desorption of N2, zeta potential, and scanning electron microscopy. The N2 adsorption studies showed that the material presents wide micropores and narrow mesopores, and has a surface area of 907.0 m2 g-1. Its maximum adsorption capacity towards the three main emerging contaminants (bisphenol A, ethinylestradiol, and amoxicillin) is much higher than that obtained with benchmark adsorbents (0.148, 0.104, and 0.072 mmol g-1, respectively). The influence of temperature and pH on the adsorption was also analyzed, allowing an improved description of the adsorption mechanism and showing very promising results.