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.

High-Level Multireference Investigations on the Electronic States in Single-Vacancy (SV) Graphene Defects Using a Pyrene-SV Model.

The nonplanar character of graphene with a single carbon vacancy (SV) defect is investigated utilizing a pyrene-SV model system by way of complete-active-space self-consistent field theory (CASSCF) and multireference configuration interaction singles and doubles (MR-CISD) calculations. Planar structures were optimized with both methods, showing the 3B1 state to be the ground state with three energetically close states within an energy range of 1 eV. These planar structures constitute saddle points. However, following the out-of-plane imaginary frequency yields more stable (by 0.22 to 0.53 eV) but nonplanar structures of Cs symmetry. Of these, the 1A' structure is the lowest in energy and is strongly deformed into an L shape. Following a further out-of-plane imaginary frequency in the nonplanar structures leads to the most stable but most deformed singlet structure of C1 symmetry. In this structure, a bond is formed between the carbon atom with the dangling bond and a carbon of the cyclopentadienyl ring. This bond stabilizes the structure by more than 3 eV compared to the planar 3B1 structure. Higher excited states were calculated at the MR-CISD level, showing a grouping of four states low in energy and higher states starting around 3 eV.

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.

Accurate rate constants for elementary reactions of molecular hydrogen and carbon monoxide mixtures and the role of the H2 rich environment.

This investigation provides accurate rate constant values for a set of elementary reactions relevant to mixtures between molecular hydrogen (H2) and carbon monoxide (CO) such as syngas. We considered intermediates and products including formaldehyde (H2CO), hydroxymethylene (c-HCOH and t-HCOH) and methanol (CH3OH). The calculations were performed employing the improved canonical variational transition state theory with small-curvature tunneling corrections based on high-level electronic structure results. This study demonstrates for the first time that H2 can act as an effective catalyst to the reaction from t-HCOH to H2CO. In this case, the adiabatic barrier height for the reaction decreases from 30.6 kcal⋅mol- 1 to 18.1 kcal⋅mol- 1 in the presence of H2. The results obtained here can improve the comprehension regarding processes such as the combustion of hydrogen-rich syngas.

The influence of the environment in chemical reactivity: the HCOOH formation from the H2O + CO reaction.

The reaction between carbon monoxide and water was studied occurring in an aerosol medium rich in methanol. This environment is plausible for the primitive and prebiotic Earth atmosphere. The chemical environment is expressed in terms of dielectric constant (ε) and the chemical system was modeled employing the polarizable continuum model (PCM). The main results were acquired from calculations employing the M06-2X density functional for the electronic structure calculations and the canonical variational theory with small curvature tunneling for the chemical kinetic calculations. The rise of ε affects both the thermochemistry and the kinetics of the reaction, increasing the barrier height and decreasing the rate constant for the reaction occurring at room temperature. For example, the rate constant at 300 K is 5-10× 10- 53 cm3 ⋅molecule- 1 ⋅s- 1 for low dielectric constant (ε < 3) and around 2-4× 10- 53 cm3 ⋅molecule- 1 ⋅s- 1 for ε between 7 and 40. Our results indicate that the ε variation allows a fine tuning to the rate of the reaction.

Relating Bond Strength and Nature to the Thermodynamic Stability of Hypervalent Togni-Type Iodine Compounds.

The bond strength and nature of a set of 32 Togni-like reagents have been investigated at the M062X/def2-TZVP(D) level of theory in acetonitrile described with the SMD continuum solvent model, to rationalize the main factors responsible for their thermodynamic stability in different conformations, and trifluoromethylation capabilities. For the assessment of bond strength, we utilized local stretching force constants and associated bond strength orders, complemented with local features of the electron density to access the nature of the bonds. Bond dissociation energies varied from 31.6 to 79.9 kcal/mol depending on the polarizing power of the ligand trans to CF3 . Based on the analysis of the Laplacian of the density, we propose that the charge-shift bond character plays an important role in the stability of the molecules studied, especially for those containing I-O bonds. New insights on the trans influence and on possible ways to fine-tune the stability of these reagents are provided.

The relativistic effects on the methane activation by gold(I) cations.

The reactivity of gold has been investigated for a long time. Here, we performed an in-depth analysis of relativistic effects over the chemical kinetic properties of elementary reactions associated with methane activation by gold(I) cations, CH4 + Au+ ↔ AuCH2 + + H2. The global reaction is modeled as a two-step process, CH4 + Au+ ↔ HAuCH3 + ↔ AuCH2 + + H2. Moreover, the barrierless dissociation of the initial adduct between reactants, AuCH4 +, is discussed as well. Higher-order relativistic treatments are used to provide corrections beyond the commonly considered scalar effects of relativistic effective core potentials (RECPs). Although the scalar relativistic contributions predominate, lowering the forward barrier heights by 48.4 and 36.1 kcal mol-1, the spin-orbit coupling effect can still provide additional reductions of these forward barrier heights by as much as 9% (1.0 and 2.2 kcal mol-1). The global reaction proceeds rapidly at low temperatures to the intermediate attained after the first hydrogen transfer, HAuCH3 +. The relativistic corrections beyond the ones from RECPs are still able to double the rate constant of the CH4 + Au+ → HAuCH3 + process at 300 K, while the reverse reaction becomes five times slower. The formation of global products from this intermediate only becomes significant at much higher temperatures (∼1500 K upward). The scalar relativistic contributions decrease the dissociation energy of the initial adduct, AuCH4 +, into the global products by 105.8 kcal mol-1, while the spin-orbit effect provides an extra lowering of 1.8 kcal mol-1.

Tunneling Enhancement of the Gas-Phase CH + CO2 Reaction at Low Temperature.

The rates of numerous activated reactions between neutral species increase at low temperatures through quantum mechanical tunneling of light hydrogen atoms. Although tunneling processes involving molecules or heavy atoms are well known in the condensed phase, analogous gas-phase processes have never been demonstrated experimentally. Here, we studied the activated CH + CO2 → HCO + CO reaction in a supersonic flow reactor, measuring rate constants that increase rapidly below 100 K. Mechanistically, tunneling is shown to occur by CH insertion into the C-O bond, with rate calculations accurately reproducing the experimental values. To exclude the possibility of H-atom tunneling, CD was used in additional experiments and calculations. Surprisingly, the equivalent CD + CO2 reaction accelerates at low temperature as zero-point energy effects remove the barrier to product formation. In conclusion, heavy-particle tunneling effects might be responsible for the observed reactivity increase at lower temperatures for the CH + CO2 reaction, while the equivalent effect for the CD + CO2 reaction results instead from a submerged barrier with respect to reactants.

Accurate Rate Constants for the Forward and Reverse H + CO ↔ HCO Reactions at the High-Pressure Limit.

The forward and reverse H + CO ↔ HCO reactions are important for combustion chemistry and have been studied from a wide variety of theoretical and experimental techniques. However, most of the chemical kinetic investigations concerning these processes are focused on low pressures or fall-off regions. Hence, a high-level electronic structure treatment was employed here in order to provide accurate rate constant values for these reactions at the high-pressure limit along temperatures from 50 to 4000 K. In relative terms, the variational effects on rate constants are shown to be almost as important at high temperatures as quantum tunneling corrections at the lowest temperatures investigated. The activation energies fitted by using modified and traditional Arrhenius' equations for the forward rate constants from 298 to 2000 K are, respectively, equal to 2.64 and 3.89 kcal mol-1, while similar fittings provided, respectively, 1.96 and 3.22 kcal mol-1, considering now forward rate constants from a temperature range of 298-373 K. This last activation energy determination (3.22 kcal mol-1) is in better agreement with the commonly referenced experimental value of 2.0 ± 0.4 kcal mol-1, also obtained from traditional fittings in the range 298-373 K, than the value attained from a broader temperature range fitting (3.89 kcal mol-1). However, some additional care must be considered along these comparisons once the experimental reaction rate measurements have been done for the trimolecular H + CO + M → HCO + M reaction instead. Anyway, the usage of appropriate temperature ranges is fundamental for reliable activation energy comparisons, although the remaining deviation between theory and experiment is still large and is possibly caused by the different pressure regimes assessed in each case. Finally, we roughly estimated that the high-pressure limit for the HCO decomposition into H and CO can be achieved along temperatures ranging from ∼246 and ∼255 K downward, respectively, at pressures of 1.1 and 9.6 atm, although further experiments should be carried out in order to improve these estimates. On the other hand, pressures larger than 9.8 × 104 atm are required for the aforementioned dissociation reaction to attain the high-pressure limit at 700 K. Therefore, the rate constants determined here are probably applicable in combustion studies at lower temperatures.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.

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.

Metal-Halogen Bonding Seen through the Eyes of Vibrational Spectroscopy.

Incorporation of a metal center into halogen-bonded materials can efficiently fine-tune the strength of the halogen bonds and introduce new electronic functionalities. The metal atom can adopt two possible roles: serving as halogen acceptor or polarizing the halogen donor and acceptor groups. We investigated both scenarios for 23 metal-halogen dimers trans-M(Y2)(NC5H4X-3)2 with M = Pd(II), Pt(II); Y = F, Cl, Br; X = Cl, Br, I; and NC5H4X-3 = 3-halopyridine. As a new tool for the quantitative assessment of metal-halogen bonding, we introduced our local vibrational mode analysis, complemented by energy and electron density analyses and electrostatic potential studies at the density functional theory (DFT) and coupled-cluster single, double, and perturbative triple excitations (CCSD(T)) levels of theory. We could for the first time quantify the various attractive contacts and their contribution to the dimer stability and clarify the special role of halogen bonding in these systems. The largest contribution to the stability of the dimers is either due to halogen bonding or nonspecific interactions. Hydrogen bonding plays only a secondary role. The metal can only act as halogen acceptor when the monomer adopts a (quasi-)planar geometry. The best strategy to accomplish this is to substitute the halo-pyridine ring with a halo-diazole ring, which considerably strengthens halogen bonding. Our findings based on the local mode analysis provide a solid platform for fine-tuning of existing and for design of new metal-halogen-bonded materials.

A Proposal for the Mechanism of the CH + CO2 Reaction.

Few experimental studies on the CH + CO2 global reaction propose H, CO, and HCO as major products. However, the reaction mechanisms behind this process have not yet been elucidated. Moreover, some intriguing kinetic particularities were noticed in these previous investigations. The advanced theoretical study performed here shows that a CH insertion mechanism is capable of explaining all the experimental data available. Hence, the strong deviations from a traditional Arrhenius behavior ascribed to the rate-determining elementary reaction (the CH insertion step) account for the kinetic particularities observed experimentally. A change in the preferred product channel as temperatures increase (from HCO + CO to H + 2CO) is also predicted to occur due to the HCO decomposition, although the CH depletion rates in typical conditions are not affected by this additional step.

Pushing 3c-4e Bonds to the Limit: A Coupled Cluster Study of Stepwise Fluorination of First-Row Atoms.

To better understand why hypervalent F, O, N, C, and B compounds are rarely stable, we carried out a systematic study of 28 systems, including anionic, cationic, and neutral molecules, held together by covalent, hypervalent, and noncovalent bonds. Molecular geometries, frequencies, atomic charges, electrostatic potentials, energy and electron densities, Mayer bond orders, local stretching force constants, and bond strength orders (BSOs) were derived from high accuracy CCSD(T) calculations and utilized to compare the strength and nature of hypervalent bonds with other types of bonds. All hypervalent molecules studied in this work were found to be either first-order transition states or unstable to dissociation, with F3- and OF3- as the only exceptions. For several systems, we found that a weak noncovalent bonded complex is more stable than a hypervalent one, due to the high energetic cost to accommodate an extra ligand, which can surpass the stability gained by 3c-4e bonding.

Emission Energies and Stokes Shifts for Single Polycyclic Aromatic Hydrocarbon Sheets in Comparison to the Effect of Excimer Formation.

Emission spectra of paradigmatic single-sheet polycyclic aromatic hydrocarbons (PAHs), pyrene, circum-1-pyrene, coronene, circum-1-coronene, and circum-2-coronene and Stokes shifts were computed and compared with previously calculated comparable data for relaxed excimer structures using the SOS-ADC(2), TD-B3LYP, and TD-CAM-B3LYP methods with multireference DFT/MRCI data as the benchmark. Vertical emission transitions and Stokes shifts were extrapolated to infinite PAH size. Comparison of Stokes shifts computed from theoretical monomer and dimer data confirms assumptions that relaxed excimers are responsible for the unusually large Stokes shifts in carbon dots observed in experimental investigations.

Potential Energy Curves for Formation of the CH2O2 Criegee Intermediate on the 3CH2 + 3O2 Singlet and Triplet Potential Energy Surfaces.

The potential energy curves (PECs) for the interaction of 3CH2 with 3O2 in singlet and triplet potential energy surfaces (PESs) leading to singlet and triplet Criegee intermediates (CH2OO) are studied using electronic structure calculations. The bonding mechanism is interpreted by analyzing the ground state multireference configuration interaction (MRCI) wave function of the reacting species and at all points along the PES. The interaction of 3CH2 with 3O2 on the singlet surface leads to a flat long-range attractive PEC lacking any maxima or minima along the curve. The triplet surface stems into a maximum along the curve resulting in a transition state with an energy barrier of 5.3 kcal/mol at CASSCF(4,4)/cc-pVTZ level. The resulting 3CH2OO is less stable than the 1CH2OO. In this study, the biradical character (ß) is used as a measure to understand the difference in the topology of the singlet and triplet PECs and the relation of the biradical nature of the species with their structures. The 3CH2OO has a larger biradical character than 1CH2OO, and because of the larger bond order of 1CH2OO, the C-O covalent bond becomes harder to break, thereby stabilizing 1CH2OO. Thus, this study provides insights into the shape of the PEC obtained from the reaction between 3CH2 and 3O2 in terms of their bonding nature and from the shape of the curves, the temperature dependence or independence of the rate of the reaction is discussed.