A Hessian free method to prevent zero-point energy leakage in classical trajectory simulation

The problem associated with the zero-point energy (ZPE) leak in classical trajectory calculations is well known. Since ZPE is a manifestation of the quantum uncertainty principle, there are no restrictions on energy during the classical propagation of nuclei. This phenomenon can lead to unphysical results, such as forming products without the ZPE in the internal vibrational degrees of freedom (DOFs). The ZPE leakage also permits reactions below the quantum threshold for the reaction. We have developed a new Hessian-free method, inspired by the Lowe-Andersen thermostat model, to prevent energy dipping below a threshold in the local-pair (LP) vibrational DOFs. The idea is to pump the leaked energy to the corresponding local vibrational mode, taken from the other vibrational DOFs. We have applied the new correction protocol on the ab initio ground-state molecular dynamics simulation of the water dimer (H20)2, which dissociates due to unphysical ZPE spilling from the high-frequency OH modes. The LP-ZPE method has been able to prevent the ZPE spilling of the OH stretching modes by pumping back the leaked energy into the corresponding modes while this energy is taken from the other modes of the dimer itself, keeping the system as a microcanonical ensemble.

Vibrational Analysis of Benziodoxoles and Benziodazolotetrazoles

Tetrazoles are well known for their high positive enthalpy of formation which makes them attractive as propellants, explosives, and energetic materials. As a step towards a deeper understanding of the stability of benziodazolotetrazole (BIAT)-based materials compared to their benziodoxole (BIO) counterparts, we investigated in this work electronic structure features and bonding properties of two monovalent iodine precursors: 2-iodobenzoic acid and 5-(2-iodophenyl)tetrazole and eight hypervalent iodine (III) compounds: I-hydroxybenzidoxolone, I-methoxybenziodoxolone, I-ethoxybenziodoxolone, I-iso-propoxybenziodoxolone and the corresponding I-hydroxyben ziodazolotetrazole, I-methoxybenziodazolotetrazole, I-ethoxybenziodazolotetrazole and I-iso- propoxybenziodazolotetrazole. As an efficient tool for the interpretation of the experimental IR spectra and for the quantitative assessment of the I−C, I−N, and I−O bond strengths in these compounds reflecting substituent effects, we used the local vibrational mode analysis, originally introduced by Konkoli and Cremer, complemented by electron density and natural bond orbital analyses. Based on the hypothesis that stronger bonds correlate with increased stability, we predict that, for both series, i.e., substituted benziodoxoles and benziodazolotetrazoles, the stability increases as follows: I-iso-propoxy < I-ethoxy < I-methoxy < I-hydroxy. In particular, the I−N bonds in the benziodazolotetrazoles could be identified as the so-called trigger bonds being responsible for the initiation of explosive decomposition in benziodazolotetrazoles. The new insight gained by this work will allow for the design of new benziodazolotetrazole materials with controlled performance or stability based on the modulation of the iodine bonds with its three ligands. The local mode analysis can serve as an effective tool to monitor the bond strengths, in particular to identify potential trigger bonds. We hope that this article will foster future collaboration between the experimental and computational community being engaged in vibrational spectroscopy.

Assessing the Intrinsic Strengths of Ion–Solvent and Solvent–Solvent Interactions for Hydrated Mg2+ Clusters

Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.

Exceptionally Long Covalent CC Bonds—A Local Vibrational Mode Study

For decades one has strived to synthesize a compound with the longest covalent C−C bond applying predominantly steric hindrance and/or strain to achieve this goal. On the other hand electronic effects have been added to the repertoire, such as realized in the electron deficient ethane radical cation in its D3d form. Recently, negative hyperconjugation effects occurring in diamino-o-carborane analogs such as di-N,N-dimethylamino-o-carborane have been held responsible for their long C−C bonds. In this work we systematically analyzed CC bonding in a diverse set of 53 molecules including clamped bonds, highly sterically strained complexes such as diamondoid dimers, electron deficient species, and di-N,N-dimethylamino-o-carborane to cover the whole spectrum of possibilities for elongating a covalent C−C bond to the limit. As a quantitative intrinsic bond strength measure, we utilized local vibrational CC stretching force constants ka(CC) and related bond strength orders BSO n(CC), computed at the ωB97X-D/aug-cc-pVTZ level of theory. Our systematic study quantifies for the first time that whereas steric hindrance and/or strain definitely elongate a C−C bond, electronic effects can lead to even longer and weaker C−C bonds. Within our set of molecules the electron deficient ethane radical cation, in D3d symmetry, acquires the longest C−C bond with a length of 1.935 Å followed by di-N,N-dimethylamino-o-carborane with a bond length of 1.930 Å. However, the C−C bond in di-N,N-dimethylamino-o-carborane is the weakest with a BSO n value of 0.209 compared to 0.286 for the ethane radical cation; another example that the longer bond is not always the weaker bond. Based on our findings we provide new guidelines for the general characterization of CC bonds based on local vibrational CC stretching force constants and for future design of compounds with long C−C bonds.

The Interstitial Carbon–Dioxygen Center in Irradiated Silicon

We investigated, experimentally as well as theoretically, defect structures in electron irradiated Czochralski-grown silicon (Cz-Si) containing carbon. Infrared spectroscopy (IR) studies observed a band at 1020 cm−1 arisen in the spectra around 300 °C. Its growth occurs concomitantly with the decay out of the well-known vacancy-oxygen (VO) defect, with a Local Vibrational Mode (LVM) at 830 cm−1 and carbon interstitial-oxygen interstitial (CiOi) defect with a LVM at 862 cm−1, in silicon (Si). The main purpose of this work is to establish the origin of the 1020 cm−1 band. One potential candidate is the carbon interstitial-dioxygen (CiO2i) defect since it is expected to form upon annealing out of the CiOi pair. To this end, systematic density functional theory (DFT) calculations were used to predict the lowest energy structure of the (CiO2i) defect in Si. Thereafter, we employed the dipole–dipole interaction method to calculate the vibrational frequencies of the structure. We found that CiO2i defect has an LVM at ~1006 cm−1, a value very close to our experimental one. The analysis and study of the results lead us to tentatively correlate the 1020 cm−1 band with the CiO2i defect.

Equilibrium Geometries, Adiabatic Excitation Energies and Intrinsic C=C/C–H Bond Strengths of Ethylene in Lowest Singlet Excited States Described by TDDFT

Seventeen singlet excited states of ethylene have been calculated via time-dependent density functional theory (TDDFT) with the CAM-B3LYP functional and the geometries of 11 excited states were optimized successfully. The local vibrational mode theory was employed to examine the intrinsic C=C/C–H bond strengths and their change upon excitation. The natural transition orbital (NTO) analysis was used to further analyze the C=C/C–H bond strength change in excited states versus the ground state. For the first time, three excited states including πy′ → 3s, πy′ → 3py and πy′ → 3pz were identified with stronger C=C ethylene double bonds than in the ground state.