RESUMO

An accurate potential energy surface (PES) for the HSiS system based on MRCI+Q calculations extrapolated to the complete basis set limit is presented. Modeled with the double many-body expansion (DMBE) method, the PES provides an accurate description of the long-range interactions, including electrostatic and dispersion terms decaying as R-4, R-5, R-6, R-8, R-10 that are predicted from dipole moments, quadrupole moments, and dipolar polarizabilities, which are also calculated at the MRCI+Q level. The novel PES is then used in quasiclassical trajectory calculations to predict the rate coefficients of the Si + SH → SiS + H reaction, which has been shown to be a major source of the SiS in certain regions of the interstellar medium. An account of the zero-point energy leakage based on various nonactive models is also given. It is shown that the reaction is dominated by long-range forces, with the mechanism Si + SH → SiSH → SSiH → SiS + H being the most important one for all temperatures studied. Although SSiH corresponds to the global minimum of the PES, the contribution from the direct reaction Si + SH → SSiH → SiS + H is less than 0.5% for temperatures higher than 500 K. The rovibrational distributions of the products are calculated using the momentum Gaussian binning method and show that as the temperature is increased the average vibrational quantum number decreases while the rotational distribution spreads up to larger values.

RESUMO

A realistic double many-body expansion potential energy surface (PES) is developed for the 2A″ state of the carbon-nitrogen-oxygen (CNO) system based on MRCI-F12/cc-pVQZ-F12 ab initio energies. The new PES reproduces the fitted points with chemical accuracy (root mean square deviation up to 0.043 eV) and explicitly incorporates long range energy terms that can accurately describe the electrostatic and dispersion interactions. Thermal rate coefficients were computed for the C(3P) + NO(2Π) reaction for temperatures ranging from 15 K to 10 000 K, and the values are compared to previously reported results. The differences are rationalized, and the major importance of long range forces in predicting the rate coefficients for barrierless reactions is emphasized.

RESUMO

The problem of obtaining the geometrical configuration of a molecule that minimizes its potential energy is a very complicated one for a series of applications, ranging from determining the structure of biological macromolecules to nanoclusters of atoms. Global optimization tools are available for this task, and many of them are based in performing successive local optimizations, where the starting geometries for these steps are determined by an intelligent algorithm. Here we develop a method to save computing time in the optimization of nanoclusters by predicting if a given minimum has been previously visited during local optimization steps. Our application to Cu-Al nanoalloys indicates that it is possible to save a substantial amount of computational cost. The application also reveals new promising AlxCuy clusters and explain their stabilities in terms of the jellium model.

RESUMO

A global potential energy surface for the lowest quartet state of HN2 is reported for the first time from accurate multireference ab initio calculations extrapolated to the complete basis set limit using the double many-body expansion method. All its stationary points are characterized, with the lowest quartet of HN2 predicted to have a bent global minimum 36 kcal mol-1 below the N(4S) + NH(X̃3Σ-) asymptote, from which it is barrierlessly achievable. The entire set of calculated ab initio points has been fitted for energies up to 1000 kcal mol-1 above the global minimum with an RMSD of 0.89 kcal mol-1, a gap comprising all identified stationary points. Special care is taken in modeling the involved long-range forces and cusps caused by crossing seams. The novel PES prompts for the calculation of rate constants for several unexplored reactions that are relevant for combustion, plasma, and atmospheric chemistry.

RESUMO

The predicted rate constants for C + NO and O + CN collisions in three potential energy surfaces (PESs) for the 2A' state of the CNO molecule are compared using quasiclassical trajectories. Different temperature dependencies are obtained for the C + NO reaction, which are explained in terms of the long-range properties of the PESs. Recommended values and mechanistic details are also reported. For O + CN collisions, a better agreement between the theoretical results is found, except for temperatures below 100 K.

RESUMO

The DFT potential energy hypersurfaces of closed-shell nitrogen clusters up to ten atoms are explored via a genetic algorithm (GA). An atom-atom distance threshold parameter, controlled by the user, and an "operator manager" were added to the standard evolutionary procedure. Both B3LYP and PBE exchange-correlation functionals with 6-31G basis set were explored using the GA. Further evaluation of the structures generated were performed through reoptimization and vibrational analysis within MP2 and CCSD(T) levels employing larger correlation consistent basis set. The binding energies of all stable structures found are calculated and compared, as well as their energies relative to the dissociation into N2, [Formula: see text] and [Formula: see text] molecules. With the present approach, we confirmed some previously reported polynitrogen structures and predicted the stability of new ones. We can also conclude that the energy surface profile clearly depends on the calculation method employed.

RESUMO

We report a new global double many-body expansion potential energy surface for the ground state of the CNO(2A') manifold, calculated by the explicit correlation multireference configuration interaction method. The functional form was accurately fitted to 3701 ab initio points with a root mean squared deviation of 0.99 kcal mol-1. All stationary points reported in previous forms are systematically described and improved, in addition to three new ones and a characterization of an isomerization transition state between the CNO and NCO minima. The novel proposed form gives a realistic description of both short-range and long-range interactions and hence is commended for dynamics studies.

RESUMO

The transitions between quartet and doublet states of the NH2 molecule are studied for the first time, allowing the evaluation of the N((4)S) + H2 reactive channel. High level ab initio calculations of the spin-orbit coupling are performed over the whole configurational space of the NH2 molecule and fitted to a proposed analytic form. Quasiclassical trajectories coupled with the surface hopping method are employed to calculate reaction cross section and rate constants. The reaction is largely affected by the initial rovibrational states of H2, while the formation of long-lived complexes enhances the reaction probability.

RESUMO

A method for modeling cusps on adiabatic potential energy surfaces without the need for any adiabatic-to-diabatic transformation is presented and shown to be successfully applied to the (2)A″ state of NO2. The more complicated case of a system with permutationally equivalent crossing seams is also examined and illustrated by considering the two first (2)A' states of the nitrogen trimer.

RESUMO

The potential energy hypersurface associated with sodium-potassium alloy clusters is explored via an enhanced genetic algorithm, where two different operators are added to the standard evolutionary procedure. Based on the recent result that the empirical Gupta many-body potential yields reasonable results for clusters with more than seven atoms, we have employed this function in the evaluation of the energies. Agglomerates from seven to the well-established 55-atom structure are studied, and their second-order energy difference and excess energies are calculated. It is found that the most stable alloys (compared to the homonuclear counterparts) are found with the proportion of sodium atoms in the range of 30 to 40%. The experimental propensity of core-shell segregation is successfully predicted by the current approach.

Assuntos

RESUMO

The electronic quenching reaction N((2)D) + N2 → N((4)S) + N2 is studied using the trajectory surface hopping method and employing two doublet and one quartet accurate potential energy surfaces. State-specific properties are analyzed, such as the dependence of the cross section on the initial quantum state of the reactants, vibrational energy transfer, and rovibrational distribution of the product N2 molecule in thermalized conditions. It is found that rotational energy on the reactant N2 molecule is effective in promoting the reaction, whereas vibrational excitation tends to reduce the reaction probability. For initial states and collision energy thermalized in an initial bath, it is found that the products are "hotter", both vibration and rotation wise.