Quantum chemical exploration of formaldehyde clusters (H2 CO)n (n = 2-4)

Structure, Carbonyl Vibrational Frequencies, and Local Energy Decomposition of Binding Energy in Formaldehyde Clusters, (HCHO)n=1-10

In this work, structures, vibrational frequencies, and binding energies of formaldehyde clusters, (HCHO)_n=1-10, are investigated by using global optimization genetic algorithm followed by density functional theory calculations and local energy decomposition (LED) analysis at the high-end DLPNO-CCSD(T) level of theory. With the use of genetic algorithm, different structures of all clusters are generated, which are further refined by the quantum chemical geometry optimization technique. From those structures, the conformer with the lowest energy is chosen and used for further analysis. The variation in carbonyl stretching frequency with change in cluster size is discussed and compared with available experimental data. Furthermore, by using standard and advanced LED analysis, different components of the binding energy are studied in all clusters and their variations with cluster size are also discussed.

Ionization potentials for the H2CO trimer

In this work, a computational study on the ionization potentials (IPs) of the formaldehyde trimer, (H₂CO)₃, is presented. Twelve lowest-lying vertical IPs were determined through the use of the coupled-cluster level of theory using correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. Specifically, the equation-of-motion ionization potential coupled-cluster with single and double excitations method with the aug-cc-pVnZ and aug-cc-pCVnZ (n = D and T) basis sets was used. The Feller-Peterson-Dixon (FPD) composite approach was employed to provide accurate IPs, and eight conformations of (H₂CO)₃ were considered. The FPD IPs determined for (H₂CO)₃ were found to be systematically lower than those computed for the dimer and monomer of H₂CO in the pattern IP_(monomer) > IP_(dimer) > IP_(trimer) for a given IP. In addition, the IPs calculated when considering only the more stable conformation (C₀) are in good agreement with those obtained using the eight conformations of the H₂CO trimer, and thus, the actual conformation played only a minor role in determining such properties in the present case. By providing first accurate IP results for the H₂CO trimer, we hope to motivate future experimental and computational investigations (e.g., studies involving photoionization) that rely on such quantities.

High performance global exploration of isomers and isomerization channels on quantum chemical potential energy surface of H5 C2 NO2

High performance global exploration of isomers and isomerization channels on the quantum chemical potential energy surface (PES) is performed for H₅ C₂ NO₂ by using the scaled hypersphere search-anharmonic downward distortion following (SHS-ADDF) method. A multi-node operation, NeoGRRM, has achieved high performance exploration calculations for the large system by submitting SHS-ADDF sub-jobs into many cores in parallel and unifying the results of sub-jobs into the total lists of the main-job. Global exploration of equilibrium (EQ) and transition-state structures at the level of B3LYP/6-31G(d) gave 3210 EQs and 23278 TSs. Nine compounds were found in the low energy regions of 0-100 kJ/mol; the lowest energy compound is N-methylcarbamic acid, the second is methyl carbamate, and the third is glycine (the most fundamental amino acid). Interconversion pathways between the conformers of each of the low energy compounds were surveyed. Isomerization channels around glycine were explored in detail. The lowest energy barriers around some of the EQs turned to be negative after zero-point energy corrections. This indicates that those structures cannot exist as independent structures because they spontaneously collapse into more stable structures. The global PES search showed various interesting dissociating channels which indicate synthon reaction pathways in the reverse directions.

Probing the ionization potentials of the formaldehyde dimer

In this work, we present a computational investigation on the ionization potentials (IPs) of the formaldehyde dimer, (H₂CO)₂. Twelve lowest lying IPs (corresponding to the entire valence orbitals) for both C_2h and C_s symmetry conformers have been computed at the coupled cluster level of theory using large correlation consistent basis sets with extrapolation to the complete basis set limit and consideration of core electron correlation effects. Specifically, the equation-of-motion ionization potential coupled-cluster with single and double (EOMIP-CCSD) excitations method with the aug-cc-pVXZ and aug-cc-pCVXZ (X = T, Q, and 5) basis sets combined with the Feller-Peterson-Dixon approach was employed, as well as CCSD with perturbative triples [CCSD(T)] with the aug-cc-pVTZ basis sets. In general, excellent agreement was observed from the comparison between the results obtained through the use of these approaches. In addition, the IPs for the formaldehyde monomer were also obtained using such methodologies and the results compared with existing experimental data; excellent agreement was also observed in this case. To the best of our knowledge, this work represents the first of its kind to determine the IPs for all these systems using a high level theory approach and is presented to motivate experimental investigations, e.g., studies involving photoionization, particularly for the formaldehyde dimer. The equilibrium binding energy of the C_2h dimer is calculated in this work at the CCSD(T)/aug-cc-pVTZ level of theory to be -4.71 kcal/mol. At this same level of theory, the equilibrium isomerization energy between C_2h and C_s conformers is 0.76 kcal/mol (C_s conformer being more stable).

Quantum Chemical Exploration of Intermolecular Reactions of Acetylene

Quantum chemical explorations of potential energy surfaces showed that acetylene produces various products starting from molecular arrays in short distances of 1.3-2.5 Å. Arrays of (C₂ H₂ )₂ gave cyclobutadiene, tetrahedrane, and acetylene dimers. Arrays of (C₂ H₂ )₃ gave benzene, prismane, benzvalene, Dewar benzene, and acetylene trimers. Arrays of (C₂ H₂ )₄ gave cubane, cyclooctatetranene, and acetylene tetramers. Different forms of initial arrays yielded different sets of products; a parallel array of two monomers gave cyclobutadiene, whereas a cross array gave tetrahedrane. Initial molecular arrays with unusually close contacts were estimated to require local forces of 1-9 × 10^-8  N. © 2019 Wiley Periodicals, Inc.