Hardness and Polarizability Profiles for Intramolecular Proton Transfer in Water Dimer Radical Cation

Ab Initio Study of Ionized Water Radical Cation (H2O)8+ in Combination with the Particle Swarm Optimization Method

The structures of cationic water clusters (H₂O)₈⁺ have been globally explored by the particle swarm optimization method in combination with quantum chemical calculations. Geometry optimization and vibrational analysis for the 15 most interesting clusters were computed at the MP2/aug-cc-pVDZ level and infrared spectrum calculation at MPW1K/6-311++G** level. Special attention was paid to the relationships between their configurations and energies. Both MP2 and B3LYP-D3 calculations revealed that the cage-like structure is the most stable, which is different from a five-membered ring lowest energy structure but agrees well with a cage-like structure in the literature. Furthermore, our obtained cage-like structure is more stable by 0.87 and 1.23 kcal/mol than the previously reported structures at MP2 and B3LYP-D3 levels, respectively. Interestingly, on the basis of their relative Gibbs free energies and the temperature dependence of populations, the cage-like structure predominates only at very low temperatures, and the most dominating species transforms into a newfound four-membered ring structure from 100 to 400 K, which can contribute greatly to the experimental infrared spectrum. By topological analysis and reduced density gradient analysis, we also investigated the structural characteristics and bonding strengths of these water cluster radical cations.

Ab initio investigation of structure, stability, thermal behavior, bonding, and infrared spectra of ionized water cluster (H2O)6. J Chem Phys 2016;145:154307. [PMID: 27782468 DOI: 10.1063/1.4964860]

The low-lying isomers of cationic water cluster (H₂O)₆⁺ have been globally explored by using particle swarm optimization algorithm in conjunction with quantum chemical calculations. Compared with previous results, our searching method covers a wide range of structural isomers of (H₂O)₆⁺ and therefore turns out to be more effective. With these local minima, geometry optimization and vibrational analysis are performed for the most interesting clusters at second-order Møller-Plesset (MP2)/aug-cc-pVDZ level, and their energies are further refined at MP2/aug-cc-pVTZ and coupled-cluster theory with single, double, and perturbative triple excitations/aug-cc-pVDZ level. The interaction energies using the complete basis set limits at MP2 level are also reported. The relationships between their structure arrangement and their energies are discussed. Based on the results of thermal simulation, structural change from a four-numbered ring to a tree-like structure occurs at T ≈ 45 K, and the relative population of six lowest-free-energy isomers is found to exceed 4% at some point within the studied temperature range. Studies reveal that, among these six isomers, two new-found isomers constitute 10% of isomer population at 180 K, and the experimental spectra can be better explained with inclusions of the two isomers. The molecular orbitals for six representative cationic water clusters are also studied. Through topological and reduced density gradient analysis, we investigated the structural characteristics and the bonding strengths of these water cluster radical cations.

Symmetry breaking in the axial symmetrical configurations of enolic propanedial, propanedithial, and propanediselenal: pseudo Jahn–Teller effect versus the resonance-assisted hydrogen bond theory

The origin of the symmetry breaking in the axial symmetrical configurations of enolic propanedial (1), propanedithial (2), and propanediselenal (3) have been investigated by means of time-dependence density functional theory and natural bond orbital interpretations. The results obtained at the quantum chemistry composite (G2MP2, CBS-QB3), ab initio molecular orbital (MP2/6-311++G**), and hybrid density functional theory (B3LYP/6-311++G**) levels of theory showed that the hydrogen-centered synchronous axial symmetrical (C2v) configurations of compounds 1–3 possessing the maximum π-electron delocalization within the M1=C2–C3=C4–M5–H6 keto-enol groups are less stable than their corresponding plane symmetrical (Cs) forms. Importantly, the symmetry breaking in the C2v configurations of the enol forms of compounds 1–3 to their corresponding plane symmetrical Cs configurations is due to the pseudo Jahn–Teller effect (PJTE) by mixing the ground A1 and excited B2 electronic states resulting in a PJT (A1 + B2) ⊗ b2 problem. We may expect that by the decrease of the energy gaps between reference states in the C2v forms that are involved in the PJTE decrease from compound 1 to compound 3, the PJT stabilization energy (PJTSE) may increase but the results obtained showed that the corresponding PJTSEs decrease. This fact can be justified by the increase of the electron delocalizations from the nonbonding orbitals of the C=M moieties to the antibonding orbitals of the H–M bonds, which leads to an increase of the π-electron delocalization within the M1=C2–C3=C4–M5–H6 keto-enol groups. In confrontation between the impacts of the resonance-assisted hydrogen bond and PJTE in the structural and configurational properties of compounds 1–3, PJTE has an overwhelming contribution and causes the symmetry breaking of the C2v configurations to their corresponding Cs forms. The correlations between the structural parameters, synchronicity indices, natural charges, PJTSEs, electron delocalizations, and the hardness of compounds 1–3 have been investigated.

Imaginary Vibrational Modes in Polycyclic Aromatic Hydrocarbons: A Challenging Test for the Hardness Profiles

In a very recent article (J. Am. Chem. Soc. 2006, 128, 9342-9343) Moran et al. found that electron-correlated methodologies using popular Pople basis sets lead to spurious nonplanar polycyclic aromatic hydrocarbon (PAH) equilibrium structures. Furthermore, some of the present authors have shown that the hardness profiles along a reaction path can be a useful tool to find spurious stationary points in the potential energy surface. Herein, we test the performance of the hardness profiles to detect shortcomings in energy profiles for the challenging case of nonplanar PAHs. The results obtained show that in 41 of the 42 imaginary vibrational modes studied, the hardness profiles indicate the wrong number and type of the potential energy surface stationary points.

Structural Investigation of Asymmetrical Dimer Radical Cation System (H2O−H2S)+: Proton-Transferred or Hemi-Bonded?

Ab initio molecular orbital and hybrid density functional methods have been employed to characterize the structure and bonding of (H2O-H2S)+, an asymmetrical dimer radical cation system. A comparison has been made between the two-center three-electron (2c-3e) hemi-bonded system and the proton-transferred hydrogen-bonded systems of (H2O-H2S)+. Geometry optimization of these systems was carried out using unrestricted Hartree Fock (HF), density functional theory with different functionals, and second-order Møller-Plesset perturbation (MP2) methods with 6-311++G(d,p) basis set. Hessian calculations have been done at the same level to check the nature of the equilibrium geometry. Energy data were further improved by calculating basis set superposition error for the structures optimized through MP2/6-311++G(d,p) calculations. The calculated results show that the dimer radical cation structure with H2O as proton acceptor is more stable than those structures in which H2O acts as a proton donor or the 2c-3e hemi-bonded (H2O thereforeSH2)+ system. This stability trend has been further confirmed by more accurate G3, G3B3, and CCSD(T) methods. On the basis of the present calculated results, the structure of H4OS+ can best be described as a hydrogen-bonded complex of H3O+ and SH with H2O as a proton acceptor. It is in contrast to the structure of neutral (H2O...H2S) dimer where H2O acts as a proton donor. The present work has been able to resolve the ambiguity in the nature of bonding between H2O and H2S in (H2O-H2S)+ asymmetrical dimer radical cation.

Stiffness and Raman Intensity: a Conceptual and Computational DFT Study

A DFT-based reactivity descriptor, the nuclear stiffness, is related to the Raman scattering intensity, which is experimentally accessible. The application of this new relationship obtained within certain approximations has been checked in two different sets of molecules. First, we study a favorable case, where the contribution of the anisotropy to the Raman intensity is zero (symmetric stretching mode in 15 tetrahedral molecules). Second, we consider a "worst" case scenario, where the anisotropy contribution can be expected to be important (stretching mode in 32 diatomic molecules). The numerical results clearly show a relationship between stiffness and Raman intensity reflecting the expected anisotropy influence.

The hardness profile as a tool to detect spurious stationary points in the potential energy surface

In the present work, we have computed the energy and hardness profiles for a series of inter and intramolecular conformational changes at several levels of calculation. All processes studied have in common the fact that the choice of a weak methodology or a poor basis set results in the presence of spurious stationary points in the energy profile. At variance with the energy profiles, the hardness profiles calculated as the difference between the vertical ionization potential and electron affinity always show the correct number of stationary points independently of the basis set and methodology used. For this reason, we have concluded that hardness profiles can be used to check the reliability of the energy profiles for those chemical systems that, because of their size, cannot be treated with high level ab initio methods.