Fitting complex potential energy surfaces to simple model potentials: application of the simplex-annealing method

A stochastic method of optimization, which combines simulated annealing with simplex, is implemented to fit the parameters of a simple model potential. The main characteristic of the method is that it explores the whole space of the parameters of the model potential, and therefore it is very efficient in locating the global minimum of the cost function, in addition to being independent of the initial guess of the parameters. The method is employed to fit the complex intermolecular potential energy surface of the dimer of water, using as a reference the spectroscopic quality anisotropic site-site potential of Feller et al. The simple model potential chosen for its reparameterization is the MCY model potential of Clementi et al. The quality of the fit is assessed by comparing the geometry of the minimum, the harmonic frequencies, and the second virial coefficients of the parameterized potential with the reference one. Finally, to prove more rigorously the robustness of this method, it is compared with standard nonstochastic methods of optimization.     

Absolute rate calculations: atom and proton transfers in hydrogen-bonded systems.

We calculate energy barriers of atom- and proton-transfer reactions in hydrogen-bonded complexes in the gas phase. Our calculations do not involve adjustable parameters and are based on bond-dissociation energies, ionization potentials, electron affinities, bond lengths, and vibration frequencies of the reactive bonds. The calculated barriers are in agreement with experimental data and high-level ab initio calculations. We relate the height of the barrier with the molecular properties of the reactants and complexes. The structure of complexes with strong hydrogen bonds approaches that of the transition state, and substantially reduces the barrier height. We calculate the hydrogen-abstraction rates in H-bonded systems using the transition-state theory with the semiclassical correction for tunneling, and show that they are in excellent agreement with the experimental data. H-bonding leads to an increase in tunneling corrections at room temperature.     

Is Cd2 truly a van der Waals molecule? Analysis of rotational profiles recorded at the , transitions

Rotational profiles of the ²²⁸Cd₂ isotopomer recorded in the (υ′, υ″) = (26, 0), (27, 0), (42, 0), (45, 0), (46, 0), (48, 0) vibrational bands of the transition were analysed. As a result, the , , , , and excited- as well as the ground-state rotational constants of the (¹¹⁴Cd)₂ were determined. The analysis allowed determining the absolute values for the and excited- and ground-state bond lengths, respectively. The obtained result – the – distinctly shorter than that obtained with assumption of pure ground-state van der Waals bonding, supports a theoretical prediction of a covalent admixture to the bonding. Analysis of the partially-resolved rotational profile recorded in the (υ′, υ″) = (38, 0) band of the same isotopomer recorded at the transition allowed estimating the rotational constant in the B1_u state.     

A comparative quantum-chemical analysis of electronic structures and spectroscopic parameters of cycloalkanes and cyclopolysilanes

A comparative quantum-chemical analysis of the electronic structures and spectroscopic parameters of the cycloalkanes C₃H₆, C₄H₈, C₅H₁₀, and C₆H₁₂ and their silicon analogs Si₃H₆, Si₄H₈, Si₅H₁₀ and Si₆H₁₂ was performed in the framework of the SCF MO LCAO method in the INDO approximation. Qualitative interpretation of “abnormal” ionization potentials and energies of electronic absorption spectra of cyclopolysilanes has been given. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1105–1108, June, 1997.     

pH,Definition and measurement at high temperatures

In this review paper, the NBS scale and its limitations are briefly discussed. The magnitude of liquid junction potentials and some calculated values are presented. The use of a molality scale for hydrogen electrode concentration cells at high temperatures is described, and results from measurements on ionization equilibria are summarized. Use of this scale is also recommended for certain circumstances with cells without liquid junction. As an alternative activity scale, use of the Pitzer ion-interaction treatment for ions is recommended for special cases. Finally, reference data are presented for _±HCl in HCl(aq) to 350°C and (HCl+NaCl)(aq) to 200°C that were derived by use of the Pitzer ion-interaction treatment.Presented at the Second International Symposium on Chemistry in High Temperature Water, Provo, UT, August 1991.     

Spin polarisation in dual catalysts for the oxygen evolution and reduction reactions

Strongly correlated catalysts can be understood from precise quantum approximations. Incorporating properly electronic correlations thus let’s define Spin rules in catalysis, opening a new door towards optimum compositions for the most important reactions for a sustainable future.     

A Photoionization Reflectron Time-of-Flight Mass Spectrometric Study on the Detection of Ethynamine (HCCNH2) and 2H-Azirine (c-H2CCHN)

Ices of acetylene (C₂H₂) and ammonia (NH₃) were irradiated with energetic electrons to simulate interstellar ices processed by galactic cosmic rays in order to investigate the formation of C₂H₃N isomers. Supported by quantum chemical calculations, experiments detected product molecules as they sublime from the ices using photoionization reflectron time-of-flight mass spectrometry (PI-ReTOF-MS). Isotopically-labeled ices confirmed the C₂H₃N assignments while photon energies of 8.81 eV, 9.80 eV, and 10.49 eV were utilized to discriminate isomers based on their known ionization energies. Results indicate the formation of ethynamine (HCCNH₂) and 2H-azirine (c-H₂CCHN) in the irradiated C₂H₂:NH₃ ices, and the energetics of their formation mechanisms are discussed. These findings suggest that these two isomers can form in interstellar ices and, upon sublimation during the hot core phase, could be detected using radio astronomy.     

Adaptive‐numerical‐bias metadynamics

A metadynamics scheme is presented in which the free energy surface is filled with progressively adding adaptive biasing potentials, obtained from the accumulated probability distribution of the collective variables. Instead of adding Gaussians with assigned height and width in conventional metadynamics method, here we add a more realistic adaptive biasing potential to the Hamiltonian of the system. The shape of the adaptive biasing potential is adjusted on the fly by sampling over the visited states. As the top of the barrier is approached, the biasing potentials become wider. This decreases the problem of trapping the system in the niches, introduced by the addition of Gaussians of fixed height in metadynamics. Our results for the free energy profiles of three test systems show that this method is more accurate and converges more quickly than the conventional metadynamics, and is quite comparable (in accuracy and convergence rate) with the well‐tempered metadynamics method. © 2017 Wiley Periodicals, Inc.     

Density functional theory and CCSD(T) evaluation of ionization potentials,redox potentials,and bond energies related to zirconocene polymerization catalysts

Quantum-mechanical-based computational design of molecular catalysts requires accurate and fast electronic structure calculations to determine and predict properties of transition-metal complexes. For Zr-based molecular complexes related to polyethylene catalysis, previous evaluation of density functional theory (DFT) and wavefunction methods only examined oxides and halides or select reaction barrier heights. In this work, we evaluate the performance of DFT against experimental redox potentials and bond dissociation enthalpies (BDEs) for zirconocene complexes directly relevant to ethylene polymerization catalysis. We also examined the ability of DFT to compute the fourth atomic ionization potential of zirconium and the effect the basis set selection has on the ionization potential computed with CCSD(T). Generally, the atomic ionization potential and redox potentials are very well reproduced by DFT, but we discovered relatively large deviations of DFT-calculated BDEs compared to experiment. However, evaluation of BDEs with CCSD(T) suggests that experimental values should be revisited, and our CCSD(T) values should be taken as most accurate.