A multi-descriptor analysis of substituent effects on the structure and aromaticity of benzene derivatives: π-Conjugation versus charge effects.

This work provides a detailed multi-component analysis of aromaticity in monosubstituted (X = CH3, C H 2 - , C H 2 + , NH2, NH-, NH+, OH, O-, and O+) and para-homodisubstituted (X = CH3, CH2, NH2, NH, OH, and O) benzene derivatives. We investigate the effects of substituents using single-reference (B3LYP/DFT) and multireference (CASSCF/MRCI) methods, focusing on structural (HOMA), vibrational (AI(vib)), topological (ELFπ), electronic (MCI), magnetic (NICS), and stability (S0-T1 splitting) properties. The findings reveal that appropriate π-electron-donating and π-electron-accepting substituents with suitable size and symmetry can interact with the π-system of the ring, significantly influencing π-electron delocalization. While the charge factor has a minimal impact on π-electron delocalization, the presence of a pz orbital capable of interacting with the π-electron delocalization is the primary factor leading to a deviation from the typical aromaticity characteristics observed in benzene.

Stability and Reactivity of the Phenalene and Olympicene Isomers.

The phenalene (triangulene) and olympicene molecules belong to the polycyclic aromatic hydrocarbon class, which have attracted substantial technological interest due to their unique electronic properties. Electronic structure calculations serve as a valuable tool in investigating the stability and reactivity of these molecular systems. In the present work, the multireference calculations, namely, the complete active space second-order perturbation theory and multireference averaged quadratic coupled cluster (MR-AQCC), were employed to study the reactivity and stability of phenalene and olympicene isomers, as well as their modified structures where the sp3-carbon at the borders were removed. The harmonic oscillator model of aromaticity (HOMA) and the nucleus-independent chemical shift as geometric and magnetic indexes calculated with density functional theory were utilized to assess the aromaticity of the studied molecules. These indexes were compared with properties such as the excitation energy and natural orbitals occupation. The reactivity analyzed using the HOMA index combined with MR-AQCC revealed the radical character of certain structures as well as the weakening of their aromaticity. Moreover, the results suggest that the removal of sp3-carbon atoms and the addition of hydrogen atoms did not alter the π network and the excitation energies of the phenalene molecules.

Electronic Structure of Small Isolated and Supported Manganese Oxide Clusters.

In the present work, possible molecular models of the isolated manganese oxides and supported Mn3Ox/Al2O3 structures were built based on small clusters of passivated MnOx. The support was represented as a simplified model of the alumina tetramer cluster based on small fragments of AlOxHy. Combinations of MnOxHy and AlOxHy clusters were made to form both the isolated and supported manganese oxides clusters. The electronic structure of these systems was characterized by ab initio methods (DFT and CASPT2). It was observed that the vertical excitation energy of the isolated and supported Mn3OxHy clusters is significantly lower than that of the alumina cluster model, while both the isolated and supported Mn3OxHy wave function characters are qualitatively similar with respect to the ground state and electronic transition processes, suggesting that the alumina cluster behaves as an inert support, since there is little contribution of this component in the description of the low-lying electronic states. The present study also reports for the first time the spectroscopic parameters of several clusters containing the manganese transition metal atom.

Structural stability and the low-lying singlet and triplet states of BN-n-acenes, n = 1-7.

The chemical stability and the low-lying singlet and triplet excited states of BN-n-acenes (n = 1-7) were studied using single reference and multireference methodologies. From the calculations, descriptors such as the singlet-triplet splitting, the natural orbital (NO) occupations and aromaticity indexes are used to provide structural and energetic analysis. The boron and nitrogen atoms form an isoelectronic pair of two carbon atoms, which was used for the complete substitution of these units in the acene series. The structural analysis confirms the effects originated from the insertion of a uniform pattern of electronegativity difference within the molecular systems. The covalent bonds tend to be strongly polarized which does not happen in the case of a carbon-only framework. This effect leads to a charge transfer between neighbor atoms resulting in a more strengthened structure, keeping the aromaticity roughly constant along the chain. The singlet-triplet splitting also agrees with this stability trend, maintaining a consistent gap value for all molecules. The BN-n-acenes molecules possess a ground state with monoconfigurational character indicating their electronic stability. The low-lying singlet excited states have charge transfer character, which proceeds from nitrogen to boron.

Electronic structure and physicochemical properties of the metal and semimetal oxide nanoclusters.

Clusters are physical entities composed of a few to thousands of atoms with capabilities to develop novel materials, like cluster-assembled materials. In this sense, knowing the electronic structure and physicochemical properties of the isolated clusters can be useful to understand how they interact with other chemical species by intermolecular forces, as free, embedded, and saturated clusters, and by intramolecular forces, acting as support clusters. In this way, in the present work, the electronic structure and physicochemical properties of metal oxide nanoclusters (MgO, Al2O3, SiO2, and TiO2) were studied by highly correlated molecular quantum chemistry methods. Through the electronic state's characterization, a semiconductor aspect was found for the titania oxide nanocluster (Te < 0.8 eV) while the other agglomerates showed a characteristic of insulating material (Te > 3.3 eV). From the stability index, the following stability order can be characterized: (SiO2)4 > (Al2O3)4 > (MgO)4 > (TiO2)3. Initial information of intermolecular and intramolecular forces caused by the studied clusters was calculated through the relative electrophilicity index, which classified the (MgO)4 and (TiO2)3 clusters as the more reactive ones, in which the (MgO)4 cluster was identified as a nucleophilic species, while the (TiO2)3 cluster as an electrophilic molecule.

Multireference study of ionic/covalent electronic states of MF (M = Be, Mg and Ca).

Ultracold environments composed by atoms or molecules offer an opportunity to study chemical reactions at the quantum-state level, for simulation of solid-state systems, as qubits in quantum computing, and for test fundamental symmetries. Those ultracold conditions formed by molecules can be obtained from cryogenic buffer gas, via supersonic expansion, followed by deceleration or from the laser cooling process. Diatomic alkaline earth monofluoride molecules have been shown as great candidates for the laser cooling process. In this sense, the present work focuses on the characterization of the low-lying doublet electronic states correlated to the first dissociation channel of the alkaline earth monofluorides diatomic molecules MF (M = Be, Mg and Ca). The developed state-of-the-art methodology was based on a qualitative analysis of the diatomic electronic structure, employing a hypothetical potential energy curve or by a simple molecular orbital diagram combined with bond order analysis. The potential energy curves, excitation and dissociation energies, and various sets of spectroscopic parameters were calculated by the MRCI/cc-pV5Z methodology. Transition probabilities for emission and radiative lifetimes among the characterized electronic states were also calculated for the (A)2Π âŸ¶ (X)2Σ+ electronic transition. Comparing the spectroscopy properties, we were able to indicate the CaF molecule as the best candidate molecule for laser cooling devices among the studied molecules.

Identification of Magic Numbers in Homonuclear Clusters: The ε3 Stability Ranking Function.

With the rise of cluster-assembled materials, an index that is able to rank and identify stable clusters or molecules is of great interest in materials sciences and engineering. In the present work, we applied a stability ranking function (ε3) in nanoclusters formed by simple metals (Na, Mg), main group elements (Al), or transition metals (Ti, Cu). The ε3 function parameters are molecular properties derived from the wave function. These parameters can be divided into kinetic and thermodynamic descriptors, in which the kinetic descriptors are the ionization potential and electronic excitation energy, while the atomization free Gibbs energy is the thermodynamic one. This simple ε3 function was able to identify the possible magic numbers of the studied clusters across the periodic table in a good agreement with previous experimental and theoretical works.

Stability and Reactivity of Silicon Magic Numbers Doped with Aluminum and Phosphorus Atoms.

The progressive scaling down of the silicon-based electronics has allowed to develop devices at nanometer scale, requiring new engineering techniques guided by fundamental chemical and physical concepts. Particularly, the nanostructured cluster systems are promising materials since their physical-chemical properties are sensitive to its shape, size, and chemical components, such that completely different materials can be produced by the simple addition or removal of a single atom. These size-tunable properties can open a new area in materials science and engineering. In the present work, quantum chemical methods were used to study the chemical substitution effects caused by subvalent (aluminum) and supervalent (phosphorus) atoms in the physical-chemical properties of some small silicon clusters, which demonstrate high stability, called magic numbers. The changes in the electronic structure and chemical acceptance to the dopants were evaluated with respect to ionization potential, electronic excitation energy, stability, and reactivity parameters. Taken together, these results enable to identify the most stable silicon-doped clusters. Regarding electrophilic reactions, Si10P is the most favorable system, while for nucleophilic reactions, none of the doped clusters resulted in higher stability.

A quantitative tool to establish magic number clusters, ε3, applied in small silicon clusters, Si2-11.

The present work focuses on establishing a function to rank the stability of small silicon clusters to characterize their magic numbers. This function is composed by a thermodynamic descriptor, the atomization Gibbs free energy, and indirect kinetic descriptors, the highest occupied molecular orbital energy and the lowest excitation energy of each system. The silicon clusters geometries were optimized using density functional theory within a hybrid meta-GGA approximation (M06), while the electronic energy was corrected by single-point calculation using CASPT2 level of theory to obtain the molecular properties. Both methodologies were combined with polarized diffused triple zeta, 6-311++G(3df,3pd), basis set for all atoms. Some molecular properties and their combinations were considered to create the aforementioned function to represent the clusters chemical stability and their magic numbers. The chosen stability ranking function, called ε3, presents results in agreement with the previous mass spectrometry experimental data identifying 4, 6, 7 and 10 as magic numbers for small silicon clusters. We believe this stability ranking function can be useful to study other intramolecular atomic and molecular clusters. Graphical abstract Stability ranking function, Îµ31, applied on Sin (n = 2 - 11) clusters showing Fukui's functions for the Sin (n = 2 - 11) obtained by the electronic density difference through CASPT2//M06/6-311++G(3df,3pd) with an isosurface value equal to 0.003.

On the importance of non-covalent interactions for porous membranes: unraveling the role of pore size.

Membrane-based gas separation technology is of crucial importance in the current economy and nanoporous graphene, given its single-atomic layer, is an essential building-block material to achieve efficiency towards permeability and selectivity for such processes. Classically, pore size is the main feature that governs the diffusion energy barrier. Its nature, nevertheless, is also affected by other non-negligible physical mechanisms not yet discussed. Here we propose a theoretical study on the role of non-covalent interactions towards H2 diffusion through two graphene-based membranes. Symmetry-Adapted Perturbation Theory (SAPT) was used to investigate the total interaction energy and its physically meaningful components (electrostatics, exchange, induction and dispersion). The study reveals the importance of quantum effects such as polarization and electron delocalization in order to counterbalance the abiding idea of pore size being the dominant factor accounting for the energy barrier. These results have important implications for the rational design of efficient nanoporous devices for separation applications.

How to efficiently tune the biradicaloid nature of acenes by chemical doping with boron and nitrogen.

Acenes are fascinating polyaromatic compounds that combine impressive semiconductor properties with an open-shell character by varying their molecular sizes. However, the increasing chemical instabilities related to their biradicaloid structures pose a great challenge for synthetic chemistry. Modifying the π-bond topology through chemical doping allows modulation of the electronic properties of graphene-related materials. In spite of the practical importance of these techniques, remarkably little is known about the basic question - the extent of the radical character created or quenched thereby. In this work, we report a high-level computational study on two acene oligomers doubly-doped with boron and nitrogen, respectively. These calculations demonstrate precisely which the chemical route is in order to either quench or enhance the radical character. Moving the dopants from the terminal rings to the central ones leads to a remarkable variation in the biradicaloid character (and thereby also in the chemical stability). This effect is related to a π-charge transfer involving the dopants and the radical carbon centers at the zigzag edges. This study also provides specific guidelines for a rational design of large polyaromatic compounds with enhanced chemical stability.

Singlet La and Lb Bands for N-Acenes (N = 2-7): A CASSCF/CASPT2 Study.

In this work CASPT2 calculations of polyacenes (from naphthalene to heptacene) were performed to find a methodology suitable for calculations of the absorption spectra, in particular of the La (B2u state) and Lb (B3u state) bands, of more extended systems. The effect of the extension of the active space and of freezing σ orbitals was investigated. The MCSCF excitation energy of the B2u state is not sensitive to the size of the active space used. However, the CASPT2 results depend strongly on the amount of σ orbitals frozen reflecting the ionic character of the B2u state. On the other hand, the excitation energies of the B3u state are much more sensitive to the size of the active space used in the calculations reflecting its multiconfigurational character. We found a good agreement with experimental data for both bands by including 14 electrons in 14 π orbitals in the active space followed by the CASPT2(14,14) perturbation scheme in which both σ and π orbitals are included.

Thermochemical and Kinetics of CH3SH + H Reactions: The Sensitivity of Coupling the Low and High-Level Methodologies.

The reaction system formed by the methanethiol molecule (CH3SH) and a hydrogen atom was studied via three elementary reactions, two hydrogen abstractions and the C-S bond cleavage (CH3SH + H → CH3S + H2 (R1); → CH2SH + H2 (R2); → CH3 + H2S (R3)). The stable structures were optimized with various methodologies of the density functional theory and the MP2 method. Two minimum energy paths for each elementary reaction were built using the BB1K and MP2 methodologies, and the electronic properties on the reactants, products, and saddle points were improved with coupled cluster theory with single, double, and connected triple excitations (CCSD(T)) calculations. The sensitivity of coupling the low and high-level methods to calculate the thermochemical and rate constants were analyzed. The thermal rate constants were obtained by means of the improved canonical variational theory (ICVT) and the tunneling corrections were included with the small curvature tunneling (SCT) approach. Our results are in agreement with the previous experimental measurements and the calculated branching ratio for R1:R2:R3 is equal to 0.96:0:0.04, with kR1 = 9.64 × 10-13 cm3 molecule-1 s-1 at 298 K.

Thermochemical and Kinetics of Hydrazine Dehydrogenation by an Oxygen Atom in Hydrazine-Rich Systems: A Dimer Model.

The kinetics of the reaction of N2H4 with oxygen depends sensitively on the initial conditions used. In oxygen-rich systems, the rate constant shows a conventional positive temperature dependence, while in hydrazine-rich setups the dependence is negative in certain temperature ranges. In this study, a theoretical model is presented that adequately reproduces the experimental results trend and values for hydrazine-rich environment, consisting of the hydrogen abstraction from the hydrazine (N2H4) dimer by an oxygen atom. The thermochemical properties of the reaction were computed using two quantum chemical approaches, the coupled cluster theory with single, double, and noniterative triple excitations (CCSD(T)) and the M06-2X DFT approach with the aug-cc-pVTZ and the maug-cc-pVTZ basis sets, respectively. The kinetic data were calculated with the improved canonical variational theory (ICVT) using a dual-level methodology to build the reaction path. The tunneling effects were considered by means of the small curvature tunneling (SCT) approximation. Potential wells on both sides of the reaction ((N2H4)2 + O → N2H4···N2H3 + OH) were determined. A reaction path with a negative activation energy was found leading, in the temperature range of 250-423 K, to a negative dependence of the rate constant on the temperature, which is in good agreement with the experimental measurements. Therefore, the consideration of the hydrazine dimer model provides significantly improved agreement with the experimental data and should be included in the mechanism of the global N2H4 combustion process, as it can be particularly important in hydrazine-rich systems.

Comparative study of small boron, silicon and germanium clusters: B(m)Si(n) and B(m)Ge(n) (m + n = 2-4).

Chemically speaking, atomic clusters are very rich, allowing their application in a broad range of technological areas such as developing functional materials, heterogeneous catalysis, and building optical devices. In this work, high level computational chemistry methods were used in a systematic manner to improve the characterization of small clusters formed by boron, silicon, germanium, mixed boron/silicon, and mixed boron/germanium. Calculations were carried out with both ab initio [MP2 and CCSD(T)] and density functional (B3LYP) methods with extended basis sets. The CCSD(T) results were then extrapolated to the complete basis set (CBS) limit. Finally, geometrical parameters, vibrational frequencies, and relative energies were then obtained and compared to data presented in the literature. Graphical Abstract Small boron, silicon and germanium clusters: BmSin and BmGen (m + n = 2-4).

Hydrogen abstraction from the hydrazine molecule by an oxygen atom.

Thermochemical and kinetics properties of the hydrogen abstraction from the hydrazine molecule (N2H4) by an oxygen atom were computed using high-level ab initio methods and the M06-2X DFT functional with aug-cc-pVXZ (X = T, Q) and maug-cc-pVTZ basis sets, respectively. The properties along the reaction path were obtained using the dual-level methodology to build the minimum energy path with the potential energy surface obtained with the M06-2X method and thermochemical properties corrected with the CCSD(T)/CBS//M06-2X/maug-cc-pVTZ results. The thermal rate constants were calculated in the framework of variational transition-state theory. Wells on both sides of the reaction (reactants and products) were found and considered in the chemical kinetics calculations. Additionally, the product yields were investigated by means of a study of the triplet and singlet surfaces of the N2H4 + O → N2H2 + H2O reaction.

Thermochemical and kinetics studies of the CH3SH+S (3P) hydrogen abstraction and insertion reactions.

Sulfur-containing molecules have a significant impact on atmosphere and biosphere. In this work we studied, from the point of view of electronic structure and chemical kinetics methods, the elementary reactions between a methanethiol molecule and a sulfur atom leading to hydrogen abstraction C-S bond cleavage (CH(3)SH+S; R1:→ CH(3)S+SH; R2: → CH(2)SH+SH; R3:→ CH(3)+HS(2)). The geometrical structures of the reactants, products, and saddle points for the three reaction paths were optimized using the BB1K method with the aug-cc-pV(T+d)Z basis set. The thermochemical properties were improved using single point coupled-cluster (CCSD(T)) calculations on the BB1K geometries followed by extrapolation to the complete basis set (CBS) limit. This methodology was previously applied and has given accurate values of thermochemical and kinetics properties when compared to benchmark calculations and experimental data. For each reaction, the thermal rate constants were calculated using the improved canonical variational theory (ICVT) including the zero-curvature (ICVT/ZCT) and small-curvature (ICVT/SCT) tunneling corrections. For comparison, the overall ICVT/SCT reaction rate constant at 300 K obtained with single-point CCSD(T)/CBS calculations for the CH(3)SH+S reaction is approximately 1400 times lower than the isovalent CH(3)SH+O reaction, obtained with CVT/SCT. The reaction path involving the hydrogen abstraction from the thiol group is the most important reactive path in all temperatures.

A multireference configuration interaction study of CuB and CuAl molecular constants and photoionization spectra.

Accurate potential energy curves and molecular constants for the low-lying electronic states of CuX(y) (X = B, Al; y = 0, +1) were investigated using the complete active space self-consistent field/multireference configuration interaction (MRCI) methodology with aug-cc-pV5Z basis set. The photoionization spectra of CuX were computed, showing electron detachment in the region of far ultraviolet. The results complement the previous theoretical characterizations and the few experimental studies. A comparative analysis was carried out concerning the different choices of reference configuration state functions in the MRCI calculations with and without the contribution of scalar relativistic effects. The results obtained with a small reference set adequately constructed are competitive to those using a much larger number of configuration state functions, and also the scalar relativistic effects improve significantly the molecular constants in this kind of system containing a 3d metal atom.

A product branching ratio controlled by vibrational adiabaticity and variational effects: kinetics of the H + trans-N2H2 reactions.

The abstraction and addition reactions of H with trans-N(2)H(2) are studied by high-level ab initio methods and density functional theory. Rate constants were calculated for these two reactions by multistructural variational transition state theory with multidimensional tunneling and including torsional anharmonicity by the multistructural torsion method. Rate constants of the abstraction reaction show large variational effects, that is, the variational transition state yields a smaller rate constant than the conventional transition state; this results from the fact that the variational transition state has a higher zero-point vibrational energy than the conventional transition state. The addition reaction has a classical barrier height that is about 1 kcal∕mol lower than that of the abstraction reaction, but the addition rates are lower than the abstraction rates due to vibrational adiabaticity. The calculated branching ratio of abstraction to addition is 3.5 at 200 K and decreases to 1.2 at 1000 K and 1.06 at 1500 K.