Search

A previously developed spectrometer for broadband electron paramagnetic resonance (EPR) spectroscopy of dilute randomly oriented systems has been considerably modified to extend the frequency reach down to the hundred MHz range and to boost concentration sensitivity by 1 to 2 orders of magnitude. The instrument is now suitable for the study of biological systems in particular metalloproteins. As a proof of concept, examples from the class of low-spin ferric hemoproteins are studied in terms of frequency-dependent changes in their EPR spectra. Mono-heme cytochrome c EPR is determined by g-strain over a wide frequency range, whereas a combination of unresolved ligand hyperfine interaction and concentration-dependent intermolecular dipolar interaction becomes dominant at very low frequencies. In the four heme containing cytochrome c3, g-strain combines with intramolecular dipolar interaction over the full-studied frequency range of 0.23-12.0 GHz. It is concluded that the point-dipole approach is inappropriate to describe magnetic interactions between low-spin ferric heme systems and that a body of literature on redox interactions in multi-heme proteins will be affected by this conclusion.

Directional internal rotation in molecular systems, generally controlled by chirality, is known to occur in natural and artificial systems driven by light or fueled chemically, but spontaneous directional molecular rotation is believed to be forbidden. We have designed a molecular rotor, whereby ferrocene and triptycene linked by a methylene bridge provide two rotational degrees of freedom. On the basis of experimental observations, in conjunction with computational data, we show that the two different modes of rotation are strongly coupled and the spatial orientation of the bistable ferrocene moiety controls the barrier to its own rotation about the triptycene axis. It is proposed that the barrier to clockwise 120° rotation across each individual triptycene blade is lower in the M-enantiomer and for counterclockwise 120° rotation, it is lower in its P-counterpart. These findings demonstrate the possibility of locally preferred thermal directional intramolecular rotation for each dynamically interconverting enantiomer.

We present a theoretical study of a two-photon absorption (2PA) process in dipolar and quadrupolar systems containing two BF2 units. For this purpose, we considered 13 systems studied by Ponce-Vargas et al. [J. Phys. Chem. B2017, 121, 10850−1085829136383] and performed linear and quadratic response theory calculations based on the RI-CC2 method to obtain the 2PA parameters. Furthermore, using the recently developed generalized few-state model, we provided an in-depth view of the changes in 2PA properties in the molecules considered. Our results clearly indicate that suitable electron-donating group substitution to the core BF2 units results in a large red-shift of the two-photon absorption wavelength, thereby entering into the desired biological window. Furthermore, the corresponding 2PA strength also increases significantly (up to 30-fold). This makes the substituted systems a potential candidate for biological imaging.

Moving beyond the six-atomic benchmark systems, we discuss the new age and future of first-principles reaction dynamics, which investigates complex, multichannel chemical reactions. We describe the methodology starting from the benchmark ab initio characterization of the stationary points, followed by full-dimensional potential energy surface (PES) developments and reaction dynamics computations. We highlight our composite ab initio approach providing benchmark stationary-point properties with subchemical accuracy, the Robosurfer program system enabling automatic PES development, and applications for the Cl + C2H6, F + C2H6, and OH- + CH3I post-six-atom reactions focusing on ab initio issues and their solutions as well as showing the excellent agreement between theory and experiment.

The properties of a material depend on how its atoms are arranged, and predicting these arrangements from first principles is a longstanding challenge. Orbital-free density functional theory provides a quantum-mechanical model based solely on the electron density, not individual wave functions. The resulting speedups make it attractive for random structure searching, whereby random configurations of atoms are relaxed to local minima in the energy landscape. We use this strategy to map the low-energy crystal structures of Li, Na, Mg, and Al at zero pressure. For Li and Na, our searching finds numerous close-packed polytypes of almost-equal energy, consistent with previous efforts to understand their low-temperature forms. For Mg and Al, the searching identifies the expected ground state structures unambiguously, in addition to revealing other low-energy structures. This new role for orbital-free density functional theory-particularly as continued advances make it accurate for more of the periodic table-will expedite crystal structure prediction over wide ranges of compositions and pressures.

A dispersion function Das in the form of a damped atom–atom asymptotic expansion fitted to ab initio dispersion energies from symmetry-adapted perturbation theory was improved and extended to systems containing heavier halogen atoms. To illustrate its performance, the revised Das function was implemented in the multipole first-order electrostatic and second-order dispersion (MED) scoring model. The extension has allowed applications to a much larger set of biocomplexes than it was possible with the original Das. A reasonable correlation between MED and experimentally determined inhibitory activities was achieved in a number of test cases, including structures featuring nonphysically shortened intermonomer distances, which constitute a particular challenge for binding strength predictions. Since the MED model is also computationally efficient, it can be used for reliable and rapid assessment of the ligand affinity or multidimensional scanning of amino acid side-chain conformations in the process of rational design of novel drugs or biocatalysts.

Independently of the preparation method, for cluster cations of aliphatic amino acids, the protonated form MnH+ is always the dominant species. This is a surprising fact considering that in the gas phase, they dissociate primarily by the loss of 45 Da, i.e., the loss of the carboxylic group. In the present study, we explore the dissociation dynamics of small valine cluster cations Mn+ and their protonated counterparts MnH+ via collision-induced dissociation experiments and ab initio calculations with the aim to elucidate the formation of MnH+-type cations from amino acid clusters. For the first time, we report the preparation of valine cluster cations Mn+ in laboratory conditions, using a technique of cluster ion assembly inside He droplets. We show that the Mn+ cations cooled down to He droplet temperature can dissociate to form both Mn-1H+ and [Mn-COOH]+ ions. With increasing internal energy, the Mn-1H+ formation channel becomes dominant. Mn-1H+ ions then fragment nearly exclusively by monomer loss, describing the high abundance of protonated clusters in the mass spectra of amino acid clusters.

Absorption and magnetic circular dichroism (MCD) spectra have been measured and theoretically simulated for a series of palladium octaethylporphyrins substituted at the meso positions with phenyl groups (n = 0–4). Analysis of the spectra included the perimeter model and time-dependent density functional theory (TDDFT) calculations. With the increasing number of phenyl substituents, the molecule is transformed from a positive hard (ΔHOMO > ΔLUMO) to a soft (ΔHOMO ≈ ΔLUMO) chromophore. This is manifested by a drastic decrease of the absorption intensity in the 0–0 region of the Q-band and by the strongly altered ratio of MCD intensities in the Q and Soret regions. Such behavior can be readily predicted using perimeter model, by analyzing frontier orbital shifts caused by various perturbations: alkyl and aryl substitution, insertion of a metal, and deviations from planarity. TDDFT calculations confirm the trends predicted by the perimeter model, but they fail in cases of less symmetrical derivatives to properly reproduce the MCD spectra in the Soret region. Our results confirm the power of the perimeter model in predicting absorption and MCD spectra of large organic molecules, porphyrins in particular. We also postulate, contrary to previous works, that the isolated porphyrin dianion is not a soft chromophore, but rather a strongly positive-hard one.

Recent electrochemical experiments have shown that the reduction of O2 can be driven backward to water oxidation, which is the first case that has been successfully demonstrated for an enzyme. To understand this ability of the enzyme, both the forward reduction and backward oxidation have been studied here. For the forward reaction, a mechanism similar to earlier studies was obtained. All steps of the full catalytic cycle were obtained for the first time, and it was shown that the explicit reduction steps contribute significantly to the rate-limiting step of the O–O bond cleavage. For the backward oxidation reaction, it was found that the mechanism of the O–O bond formation is not just the reverse of the reduction step where the O–O bond is cleaved for a protonated peroxide. The formation of two fully deprotonated oxo groups was found to be important, which leads to a large radical character for one of the oxo groups. For this possibility, it is important that the pKa of the water bound to the cofactor is quite high.

The Waddington mechanism, or the Waddington-type reaction pathway, is crucial for low-temperature oxidation of both alkenes and alcohols. In this study, the Waddington mechanism in the oxidation chemistry of butene and butanol isomers was systematically investigated. Fundamental quantum chemical calculations were conducted for the rate constants and thermodynamic properties of the reactions and species in this mechanism. Calculations were performed using two different ab initio solvers: Gaussian 09 and Orca 4.0.0, and two different kinetic solvers: PAPR and MultiWell, comprehensively. Temperature- and pressure-dependent rate constants were performed based on the transition state theory, associated with the Rice Ramsperger Kassel Marcus and master equation theories. Temperature-dependent thermochemistry (enthalpies of formation, entropy, and heat capacity) of all major species was also conducted, based on the statistical thermodynamics. Of the two types of reaction, dissociation reactions were significantly faster than isomerization reactions, while the rate constants of both reactions converged toward higher temperatures. In comparison, between two ab initio solvers, the barrier height difference among all isomerization and dissociation reactions was about 2 and 0.5 kcal/mol, respectively, resulting in less than 50%, and a factor of 2-10 differences for the predicted rate coefficients of the two reaction types, respectively. Comparing the two kinetic solvers, the rate constants of the isomerization reactions showed less than a 32% difference, while the rate of one dissociation reaction (P1 ↔ WDT12) exhibited 1-2 orders of magnitude discrepancy. Compared with results from the literature, both reaction rate coefficients (R4 and R5 reaction systems) and species' thermochemistry (all closed shell molecules and open shell radicals R4 and R5) showed good agreement with the corresponding values obtained from the literature. All calculated results can be directly used for the chemical kinetic model development of butene and butanol isomer oxidation.

Recent advances in the calculation of activation energies are shedding new light on the dynamical time scales of liquid water. In this Perspective, we examine how activation energies elucidate the central, but not singular, role of the exchange of hydrogen-bond (H-bond) partners that rearrange the H-bond network of water. The contributions of other motions to dynamical time scales and their associated activation energies are discussed along with one case, vibrational spectral diffusion, where H-bond exchanges are not mechanistically significant. Nascent progress on outstanding challenges, including descriptions of non-Arrhenius effects and activation volumes, are detailed along with some directions for future investigations.

As a key diagnostic property of benzenoids and other polycyclic hydrocarbons, induced ring current has inspired diverse approaches for calculation, modeling, and interpretation. Grid-based methods include the ipsocentric ab initio calculation of current maps, and its surrogate, the pseudo-π model. Graph-based models include a family of conjugated-circuit (CC) models and the molecular-orbital Huckel-London (HL) model. To assess competing claims for physical relevance of derived current maps for benzenoids, a protocol for graph-reduction and comparison was devised. Graph reduction of pseudo-π grid maps highlights their overall similarity to HL maps, but also reveals systematic differences. These are ascribed to unavoidable pseudo-π proximity limitations for benzenoids with short nonbonded distances, and to poor continuity of pseudo-π current for classes of benzenoids with fixed bonds, where single-reference methods can be unreliable. Comparison between graph-based approaches shows that the published CC models all shadow HL maps reasonably well for most benzenoids (as judged by L1-, L2-, and L∞-error norms on scaled bond currents), though all exhibit physically implausible currents for systems with fixed bonds. These comparisons inspire a new combinatorial model (Model W) based on cycle decomposition of current, taking into account the two terms of lowest order that occur in the characteristic polynomial. This improves on all pure-CC models within their range of applicability, giving excellent adherence to HL maps for all Kekulean benzenoids, including those with fixed bonds (halving the rms discrepancy against scaled HL bond currents, from 11% in the best CC model, to 5% for the set of 18 360 Kekulean benzenoids on up to 10 hexagonal rings). Model W also has excellent performance for open-shell systems, where currents cannot be described at all by pure CC models (4% rms discrepancy against scaled HL bond currents for the 20112 non-Kekulean benzenoids on up to 10 hexagonal rings). Consideration of largest and next-to-largest matchings is a useful strategy for modeling and interpretation of currents in Kekulean and non-Kekulean benzenoids (nanographenes).

Quasi-classical direct dynamics simulations, performed with the B3LYP-D3/cc-pVDZ electronic structure theory, are reported for vibrational relaxation of the three NH stretches of the -NH3+ group of protonated tryptophan (TrpH+), excited to the n = 1 local mode states. The intramolecular vibrational energy relaxation (IVR) rates determined for these states, from the simulations, are in good agreement with the experiment. In accordance with the experiment, IVR for the free NH stretch is slowest, with faster IVR for the remaining two NH stretches which have intermolecular couplings with an O atom and a benzenoid ring. For the free NH and the NH coupled to the benzenoid ring, there are beats (i.e., recurrences) in their relaxations versus time. For the free NH stretch, 50% of the population remained in n = 1 when the trajectories were terminated at 0.4 ps. IVR for the free NH stretch is substantially slower than for the CH stretch in benzene. The agreement found in this study between quasi-classical direct dynamics simulations and experiments indicates the possible applicability of this simulation method to larger biological molecules. Because IVR can drive or inhibit reactions, calculations of IVR time scales are of interest, for example, in unimolecular reactions, mode-specific chemistry, and many photochemical processes.

Accurate predictions of hyperfine structure (HFS) constants are important in many areas of chemistry and physics, from the determination of nuclear electric and magnetic moments to benchmarking of new theoretical methods. We present a detailed investigation of the performance of the relativistic coupled cluster method for calculating HFS constants within the finite-field scheme. The two selected test systems are 133Cs and 137BaF. Special attention has been paid to construct a theoretical uncertainty estimate based on investigations on basis set, electron correlation and relativistic effects. The largest contribution to the uncertainty estimate comes from higher order correlation contributions. Our conservative uncertainty estimate for the calculated HFS constants is ∼5.5%, while the actual deviation of our results from experimental values is

We present the first experimental study of the frontier orbitals in an ultrathin film of the novel hexa-carbene photosensitizer [Fe(btz)3]3+, where btz is 3,3′-dimethyl-1,1′-bis(p-tolyl)-4,4′-bis(1,2,3-triazol-5-ylidene). Resonant photoelectron spectroscopy (RPES) was used to probe the electronic structure of films where the molecular and oxidative integrities had been confirmed with optical and X-ray spectroscopies. In combination with density functional theory calculations, RPES measurements provided direct and site-selective information about localization and interactions of occupied and unoccupied molecular orbitals. Fe 2p, N 1s, and C 1s measurements selectively probed the metal, carbene, and side-group contributions revealing strong metal–ligand orbital mixing of the frontier orbitals. This helps explain the remarkable photophysical properties of iron-carbenes in terms of unconventional electronic structure properties and favorable metal–ligand bonding interactions—important for the continued development of these type of complexes toward light-harvesting and light-emitting applications.

The rotational spectrum of 2,5-dichlorothiophene (DCT) was measured for the first time using Fourier transform microwave spectroscopy from 5.5-19 GHz. Dense hyperfine splitting patterns due to the two quadrupolar chlorine nuclei (I = 3/2) were resolved and assigned for the 35Cl-35Cl, 37Cl-35Cl, and 37Cl-37Cl isotopologues and for the two 13C and one 34S analogues with two 35Cl atoms, allowing derivation of their respective nuclear quadrupole coupling tensors. The rotational constants obtained from fitting the spectra of the six isotopic species allowed derivation of the experimental geometry of DCT for comparison with the equilibrium structure computed at the MP2/aug-cc-pVTZ level. This revealed that the electron-withdrawing effect of chlorine causes small distortions in the ring geometry relative to thiophene, including a 1.1° increase in the two S-C-C angles and a 0.012 A increase in the two S-C bond lengths.

Quantum mechanics/molecular mechanics (QM/MM) calculations are widely used embedding techniques to computationally investigate enzyme reactions. In most QM/MM computations, the quantum mechanical region is treated through density functional theory (DFT), which offers the best compromise between chemical accuracy and computational cost. Nevertheless, to obtain more accurate results, one should resort to wave function-based methods, which however lead to a much larger computational cost already for relatively small QM subsystems. To overcome this drawback, we propose the coupling of our QM/ELMO (quantum mechanics/extremely localized molecular orbital) approach with molecular mechanics, thus introducing the three-layer QM/ELMO/MM technique. The QM/ELMO strategy is an embedding method in which the chemically relevant part of the system is treated at the quantum mechanical level, while the rest is described through frozen ELMOs. Since the QM/ELMO method reproduces results of fully QM computations within chemical accuracy and with a much lower computational effort, it can be considered a suitable strategy to extend the range of applicability and accuracy of the QM/MM scheme. In this paper, other than briefly presenting the theoretical bases of the QM/ELMO/MM technique, we will also discuss its validation on the well-tested deprotonation of acetyl coenzyme A by aspartate in citrate synthase.

Superalkalis and superhalogens are atomic clusters that mimic the chemistry of alkali and halogen atoms, respectively; the ionization energies of the superalkalis are less than those of alkali atoms, while the electron affinities of superhalogens are larger than those of the halogen atoms. These superions can serve as the building blocks of a new class of supersalts with applications in solar cells, metal-ion batteries, multiferroic materials, and so on. While considerable progress has been made in the design and synthesis of superhalogens, a similar understanding of superalkalis is lacking. Using density functional theory with hybrid exchange-correlation functional and Gaussian basis sets, we have systematically studied the role of size and composition on the properties of two different classes of clusters whose stabilities are governed by the Wade-Mingos polyhedral skeletal electron pair theory. One class belongs to the closo-borane family LimBnXn (m = 1, 2, 3; n = 6, 12; X = H, F, CN), while the other to the Zintl ions Lim[Be@Ge9]. We show that Li3BnXn and Li3[Be@Ge9] clusters are superalkalis with ionization energies as low as 2.84 eV in Li3B6H6. However, contrary to expectation, the ionization energies do not decrease with increasing cluster volume. Instead, ionization energies are linked to the X ligands' electron affinities; the larger the electron affinity, the higher is the ionization energy. The understanding gained here will help in the discovery of superalkalis and, hence, enrich the library of supersalts.

Fluid-fluid interfacial free energy can be measured accurately and can also be calculated from molecular simulations. However, it is challenging to measure solid-fluid interfacial free energy directly. Accurate computation has not yet been advanced by molecular simulations. In this study, we derive working expressions for estimating solid-fluid interfacial free energy based on the free-energy perturbation method with consideration of solid deformation. A Lennard-Jones solid-fluid system is simulated. Our derivations indicate that the effect of solid deformation is pronounced on solid-fluid interfacial free energy, and the results may be significantly different from the conventional test area method. Our results reveal that the contribution of the solid deformation highly depends on the stress conditions in the solid, which can be either positive or negative. Adsorption of fluids onto the solid surface has a significant effect on interfacial free energy. In weak adsorption, the interfacial free energy is close to the solid-vacuum surface free energy. Strong adsorption results in a significant reduction in interfacial free energy.

The reactions of 1-tetradecene and 1-pentadecene, the C14 and C15 linear 1-alkenes, with OH radicals in the presence of NOx were investigated in a series of environmental chamber experiments. Particle-phase β-hydroxynitrates, dihydroxynitrates, dihydroxycarbonyls, and 1,4-hydroxynitrates and gas-phase aldehydes were sampled and then identified and quantified using a suite of offline analytical techniques that included derivatization, gas and liquid chromatography, and multiple types of mass spectrometry. Measured molar yields of products formed by OH radical addition to the C═C double bond, including β-hydroxynitrates, dihydroxynitrates, dihydroxycarbonyls (which have not been previously directly quantified with high accuracy), and aldehydes were 0.125 ± 0.01, 0.048 ± 0.005, 0.240 ± 0.04, and 0.268 ± 0.03 (0.264 ± 0.02 and 0.271 ± 0.04 for the formaldehyde and tridecanal/tetradecanal co-products of β-hydroxyalkoxy radical decomposition), respectively. These values give a total molar yield of 0.681 ± 0.05, which agrees very well with the results of kinetics measurements that indicate that the fraction of reaction that occurs by OH radical addition is 0.70. The yields were used to calculate branching ratios for all OH radical addition pathways, including a value of 0.18 for the formation of dihydroxynitrates from the reaction of dihydroxyperoxy radicals with NO and values of 0.47 and 0.53 for β-hydroxyalkoxy radical decomposition and isomerization. The results were used with literature data on the yields of aldehydes measured for similar reactions of smaller alkenes, a model for the effect of carbon number on branching ratios for organic nitrate formation, and a mechanism for H atom abstraction derived from studies of linear alkanes to achieve a complete, quantitative gas-phase reaction mechanism for 1-alkenes. The results should also be useful for constructing mechanisms for more complex reactions of volatile organic compounds.

Diffusion Monte Carlo (DMC) provides a powerful method for understanding the vibrational landscape of molecules that are not well-described by conventional methods. The most computationally demanding step of these calculations is the evaluation of the potential energy. In this work, a general approach is developed in which a neural network potential energy surface is trained by using data generated from a small-scale DMC calculation. Once trained, the neural network can be evaluated by using highly parallelizable calls to a graphics processing unit (GPU). The power of this approach is demonstrated for DMC simulations on H2O, CH5+, and (H2O)2. The need to include permutation symmetry in the neural network potentials is explored and incorporated into the molecular descriptors of CH5+ and (H2O)2. It is shown that the zero-point energies and wave functions obtained by using the neural network potentials are nearly identical to the results obtained when using the potential energy surfaces that were used to train the neural networks at a substantial savings in the computational requirements of the simulations.

Light-induced blinking, an inherent feature of many forms of super-resolution microscopy, has been linked to transient reduction of the fluorescent cyanine dye used as an imaging agent. There is, however, only scant literature information related to one-electron reduced cyanine dyes, especially in an aqueous environment. Here, we examine a small series of cyanine dyes, possessing disparate π-conjugation lengths, under selective reducing or oxidizing conditions. The experiment allows recording of both differential absorption spectra and decay kinetics of the resultant one-electron reduced or oxidized transient species in water. Relative to the ground state, absorption transitions for the various radicals are weak and somewhat broadened but do allow correlation with the π-conjugation length. In all cases, absorption maxima lie to the blue of the main ground-state transition. Under anaerobic conditions, the transient species decay on the microsecond to millisecond time scale, with the mean lifetime depending on molecular structure, radiation dose, and dye concentration. The experimental absorption spectra recorded for the one-electron reduced radicals and the presumed dimer cation radical compare well to spectra obtained from time-dependent density functional theory calculations. The results allow conclusions to be drawn regarding the plausibility of the reduced species being responsible for light-induced blinking in direct stochastic optical reconstruction microscopy.

A numerical approach is developed to capture the effect of rotation-vibration coupling in a practically affordable way. In this approach only a limited number of adjacent rotational components are considered to be coupled, while the couplings to other rotational components are neglected. This partially coupled (PC) approach permits to reduce the size of Hamiltonian matrix significantly, which enables the calculations of ro-vibrational states above dissociation threshold (scattering resonances) for large values of total angular momentum. This method is employed here to reveal the role of the Coriolis effect in the ozone formation reaction at room temperature, dominated by large values of total angular momentum states, on the order of J = 24 and 28. We found that, overall, the effect of ro-vibrational coupling is not minor for large J. Compared to the results of symmetric top rotor approximation, where the ro-vibrational coupling is neglected, we found that the widths of scattering resonances, responsible for the lifetimes of metastable ozone states, remain nearly the same (on average), but the number of these states increases by as much as 20%. We also found that these changes are nearly the same in symmetric and asymmetric ozone isotopomers 16O18O16O and 16O16O18O. Therefore, based on the results of these calculations, the Coriolis coupling does not seem to favor the formation of asymmetric ozone molecules and thus cannot be responsible for symmetry-driven mass-independent fractionation of oxygen isotopes.

The gas-phase rotational spectrum of (cyanomethylene)cyclopropane, (CH2)2C═CHCN, generated by a Wittig reaction between the hemiketal of cyclopropanone and (cyanomethylene)triphenylphosphorane, is presented for the first time. This small, highly polar nitrile is a cyclopropyl-containing structural isomer of pyridine. The rotational spectra of the ground state and two vibrationally excited states were observed, analyzed, and least-squares fit from 130 to 360 GHz. Over 3900 R-, P-, and Q-branch, ground-state rotational transitions were fit to low-error, partial octic, A- and S-reduced Hamiltonians, providing precise determinations of the spectroscopic constants. The two lowest-energy vibrationally excited states, ν17 and ν27, form a Coriolis-coupled dyad displaying small a- and b-type resonances. Transitions for these two states were measured and least-squares fit to a two-state, partial octic, A-reduced Hamiltonian in the Ir representation with nine Coriolis-coupling terms (Ga, GaJ, GaK, GaJJ, Fbc, FbcJ, FbcK, Gb, and GbJ). The observation of many resonant transitions and nine nominal interstate transitions enabled a very accurate and precise energy difference between ν17 and ν27 to be determined: ΔE17,27 = 29.8975453 (33) cm-1. The spectroscopic constants presented herein provide the foundation for future astronomical searches for (cyanomethylene)cyclopropane.

Ceria (CeO2) is a well-known catalytic oxide with many environmental, energy production, and industrial applications, most of them involving water as a reactant, byproduct, solvent, or simple spectator. In this work, we parameterized a Ce/O/H ReaxFF for the study of ceria and ceria/water interfaces. The parameters were fitted to an ab initio training set obtained at the DFT/PBE0 level, including the structures, cohesive energies, and elastic properties of the crystalline phases Ce, CeO2, and Ce2O3; the O-defective structures and energies of vacancy formation on CeO2 bulk and CeO2 (111) surface, as well as the absorption and reaction energies of H2 and H2O molecules on CeO2 (111). The new potential reproduced reasonably well all the fitted properties as well as the relative stabilities of the different ceria surfaces, the oxygen vacancies formation, and the energies and structures of associative and dissociative water molecules on them. Molecular dynamics simulations of the liquid water on the CeO2 (111) and CeO2 (100) surfaces were carried out to study the coverage and the mechanism of water dissociation. After equilibration, on average, 35% of surface sites of CeO2 (111) are hydroxylated whereas 15% of them are saturated with molecular water associatively adsorbed. As for the CeO2 (100) surface, we observed that water preferentially dissociates covering 90% of the available surface sites in excellent agreement with recent experimental findings.

The calculation of acid dissociation constants (pKa) is an important task in computational chemistry and chemoinformatics. Theoretically and with minimal empiricism, this is possible from computed acid dissociation free energies via so-called linear free-energy relationships. In this study some modifications are introduced to the latter, providing a straightforward, broadly applicable protocol with an adjustable degree of sophistication for quantum chemistry-based calculations of pKa in water. It targets a wide pKa range (∼70 units) and medium-sized, flexible molecules. Herein, a focus is set on the recently published r2SCAN-3c and related efficient composite density functionals and the semiempirical GFN2-xTB method, including a newly introduced energy correction for heterolytic dissociation, both in combination with implicit solvation models. The performance is evaluated in comparison with experimental data, showing mean errors often smaller than a targeted 1 pKa unit accuracy. Larger deviations are observed only upon inclusion of challenging highly negative ( 15) pKa values. Among all those tested, it is found that B97-3c is the best performing functional, although rather independently of the density functional theory (DFT) method used; low root-mean-square errors of 0.8-1.0 pKa units for typical drugs are obtained. For optimal performance, it is recommended to employ DFT functional specific free-energy relationship parameters. Additionally, a significant conformational dependence of the pKa values is revealed and quantified for some nonrigid drug molecules.

The dispersive optical activity of two saturated cyclic amines, (R)-2-methylpyrrolidine (R-2MPY) and (S)-2-methylpiperidine (S-2MPI), has been interrogated under isolated and solvated conditions to elucidate the roles of large-amplitude motion associated with nitrogen-center inversion and ring-puckering dynamics. Experimental optical rotatory dispersion profiles were almost mirror images of one another and displayed parallel solvent dependencies. Quantum-chemical analyses built on density-functional and coupled-cluster methods revealed four low-lying conformers for each molecule, which are distinguished by axial/equatorial orientations of their amino hydrogens and methyl substituents. Chiroptical signatures predicted for these species were combined through an independent-conformer ansatz to simulate the ensemble-averaged response, with a polarizable continuum model (PCM) being used to treat implicit solute-solvent interactions. The intrinsic behavior observed for isolated (gaseous) R-2MPY and S-2MPI was reproduced best by merging coupled-cluster (CCSD) estimates of rotatory powers with thermal population fractions deduced from complete basis set (CBS-APNO) free-energy calculations. Although prior claims of sizable chiroptical contributions arising from helically twisted (chiral) heterocyclic frameworks could be discounted, less satisfactory agreement between experiment and theory was realized for solution phases. Response properties sustained modest isomer-dependent changes in the presence of PCM solvation, but the corresponding energy metrics showed systematic trends, whereby structures having larger electric-dipole moments were stabilized preferentially in media of high polarity. Despite the fact that R-2MPY conformations were predicted to undergo a progressive reordering of their relative energies across the six solvents of interest, S-2MPI was found to exhibit more pronounced solvent-induced perturbations at long wavelengths (viz., in regions far removed from electronic resonances). Experimental results are discussed in terms of the distinct ring-puckering mechanisms for R-2MPY and S-2MPI, which are expected to be dominated by hindered pseudorotation among envelope/twist motifs and semi-inversion between chairlike antipodes, respectively.

The Coulomb explosion dynamics following strong field ionization of chlorocarbonylsulfenyl chloride was studied using multimass coincidence detection and covariance imaging analysis, supported by density functional theory calculations. These results show evidence of multiple dissociation channels from various charge states. Double ionization to low-lying electronic states leads to a dominant C-S cleavage channel, while higher states can alternatively correlate to the loss of Cl+. Triple ionization leads to a double dissociation channel, the observation of which is confirmed via three-body covariance analysis, while further ionization leads primarily to atomic or diatomic fragments whose relative momenta depend strongly on the starting structure of the molecule.

Carbon clusters exhibit a broad diversity of topologies and shapes, encompassing fullerene-like cages, graphene-like flakes, and more disordered pretzel-like and branched structures. Here, we examine computationally their infrared spectra in relation with these structures from a statistical perspective. Individual spectra for broad samples of isomers were determined by means of the self-consistent charge density functional-based tight-binding method, and an interpolation scheme is designed to reproduce the spectral features by regression on a much smaller subset of the sample. This interpolation proceeds by encoding the structures using appropriate descriptors and selecting them through principal component analysis, Gaussian regression or inverse distance weighting providing the nonlinear weighting functions. Metric learning is employed to reduce the global error on a preselected testing set. The interpolated spectra satisfactorily reproduce the specific spectral features and their dependence on the size and shape, enabling quantitative prediction away from the testing set. Finally, the classification of structures within the four proposed families is critically discussed through a statistical analysis of the sample based on iterative label spreading.

The thermal unimolecular decomposition of 2-methyltetrahydrofuran (2-MTHF) was studied behind reflected shock waves in a single-pulse shock tube over the temperature range of 1179-1361 K and pressure range of 9-17 atm. Methane, ethylene, ethane, 1,3-butadiene, propylene, acetaldehyde, and acetylene were identified as products in the decomposition of 2-MTHF. A reaction scheme was proposed to explain the mechanism for the observed products. The experimentally determined rate coefficients were best fit to an Arrhenius expression for the overall decomposition and is represented as ktotalexp(1179-1361 K) = (3.23 ± 0.59) × 1011 s-1 exp(-51.3 ± 1.4 kcal mol-1/RT). Quantum chemistry methods were used to calculate the energetics and kinetics of various possible unimolecular dissociation pathways involved in the thermal decomposition of 2-MTHF. The initial decomposition of 2-MTHF occurs predominantly via ring-methyl (C-CH3) single bond fission, leading to the formation of tetrahydrofuran (C4H7O) radical, and methyl radical was found to be the major reaction compared to all the possible initial bond fission, ring opening, and molecular elimination channels. The temperature-dependent rate coefficients for the unimolecular dissociation of 2-MTHF were calculated using the RRKM (Rice-Ramsperger-Kassel-Marcus) theory in combination with the CCSD(T)/cc-pVTZ//B3LYP/cc-pVTZ level of electronic structure calculations over the temperature range of 800-1500 K. The computed high-pressure limiting rate coefficients for the initial decomposition of 2-MTHF through C-CH3 single bond fission channel were found to be ∼2 times higher in the temperatures between 800 and 900 K, and above this temperature, they agree well with the values reported in the literature.

The ultraviolet (UV) photodissociation dynamics of the jet-cooled cyclohexyl (c-C6H11) radical is studied using the high-n Rydberg atom time-of-flight (HRTOF) technique. The cyclohexyl radical is produced by the 193 nm photodissociation of chlorocyclohexane and bromocyclohexane and is examined in the photolysis wavelength region of 232-262 nm. The H-atom photofragment yield (PFY) spectrum contains a broad peak centered at 250 nm, which is in good agreement with the UV absorption spectrum of the cyclohexyl radical and assigned to the 3p Rydberg states. The translational energy distributions of the H-atom loss product channel, P(ET)'s, are bimodal, with a slow (low ET) component peaking at ∼6 to 7 kcal/mol and a fast (high ET) component peaking at ∼44-48 kcal/mol. The fraction of the average translational energy in the total excess energy, ⟨fT⟩, is in the range of 0.16-0.25 in the photolysis wavelength region of 232-262 nm. The H-atom product angular distribution of the slow component is isotropic, while that of the fast component is anisotropic with an anisotropy parameter of β ≈ 0.5-0.7. The bimodal product translational energy and angular distributions indicate two dissociation pathways to the H + C6H10 products in cyclohexyl. The high-ET anisotropic component is from a repulsive, prompt dissociation on a repulsive potential energy surface coupling with the Rydberg excited states to produce H + cyclohexene. The low-ET isotropic component is consistent with the unimolecular dissociation of hot radical on the ground electronic state after internal conversion from the Rydberg states. The similarity of the photodissociation dynamics of the cyclohexyl radical to the previously studied small linear and branched alkyls expands on the understanding of the dissociation dynamics of alkyl radicals to include larger cyclic alkyl radicals.

Pyrrole-2-carbaldehyde (Pa) forms one of the unnatural nucleic acid bases, and as a base pair with 7-(2-thienyl)imidazo[4,5-b]pyridine (Ds), it has been known to be stable in DNA. The Ds-Pa pair is stabilized in DNA via van der Waals' interaction and shape fitting. There are some studies on the origin of its stability and reactivity in the ground state. However, for a successful unnatural base pair, it needs to be stable not only in the ground state but also upon irradiation with UV-visible light. To understand the photoinduced reactivity, we investigate the excited-state properties of the Pa base and understand the energetically feasible photoprocesses that can occur upon excitation in the UV region. Two distinct pathways are obtained. One of the pathways involves an out-of-plane mode and has some similarities with the deactivation channels in the natural pyrimidine bases. On the other hand, the second pathway involves an excited-state proton transfer.

The energetics of small cationic tantalum clusters and their gas-phase adsorption and dehydrogenation reaction pathways with methane are investigated with ion-trap experiments and spin-density-functional-theory calculations. Tan+ clusters are exposed to methane under multicollision conditions in a cryogenic ring electrode ion-trap. The cluster size affects the reaction efficiency and the number of consecutively dehydrogenated methane molecules. Small clusters (n = 1-4) dehydrogenate CH4 and concurrently eliminate H2, while larger clusters (n > 4) demonstrate only molecular adsorption of methane. Unique behavior is found for the Ta+ cation, which dehydrogenates consecutively up to four CH4 molecules and is predicted theoretically to promote formation of a [Ta(CH2-CH2-CH2)(CH2)]+ product, exhibiting C-C coupled groups. Underlying mechanisms, including reaction-enhancing couplings between potential energy surfaces of different spin-multiplicities, are uncovered.

A time-dependent (TD) formulation of equation-of-motion coupled-cluster (EOM-CC) theory can provide excited-state information over an arbitrarily wide energy window with a reduced memory footprint relative to conventional, frequency-domain EOM-CC theory. However, the floating-point costs of the time-integration required by TD-EOM-CC are generally far larger than those of the frequency-domain form of the approach. This work considers the potential of the short iterative Lanczos (SIL) integration scheme [J. Chem. Phys.1986, 85, 5870-5876] to reduce the floating-point costs of TD-EOM-CC simulations. Low-energy and K-edge absorption features for small molecules are evaluated using TD-EOM-CC with single and double excitations, with the time-integrations carried out via SIL and fourth-order Runge-Kutta (RK4) schemes. Spectra derived from SIL- and RK4-driven simulations are nearly indistinguishable, and with an appropriately chosen subspace dimension, the SIL requires far fewer floating-point operations than are required by RK4. For K-edge spectra, SIL is the more efficient scheme by an average factor of 7.2.

An iterative approach is introduced, which allows the efficient solution of the hierarchical equations of motion (HEOM) for the steady-state of open quantum systems. The approach combines the method of matrix equations with an efficient preconditioning technique to reduce the numerical effort of solving the HEOM. Illustrative applications to simulate nonequilibrium charge transport in single-molecule junctions demonstrate the performance of the method.

Secondary organic aerosol formation in the atmospheric aqueous/particulate phase by photosensitized reactions is currently subject to uncertainties. To understand the impact of photosensitized reactions, photophysical and -chemical properties of photosensitizers, kinetic data, and reaction mechanisms of these processes are required. The photophysical properties of acetophenones, benzaldehydes, benzophenones, and naphthalenes were investigated in aqueous solution using laser flash excitation. Quantum yields of excited photosensitizers were determined giving values between 0.06-0.80 at 298 K and pH = 5. Molar absorption coefficients (emax(3PS*) = (0.8-13) × 104 L mol-1 cm-1), decay rate constants in water (k1st = (9.4 ± 0.5) × 102 to (2.2 ± 0.1) × 105 s-1), and quenching rate constants with oxygen (kq(O2) = (1.7 ± 0.1-4.4 ± 0.4) × 109 L mol-1 s-1) of the excited triplet states were determined at 298 K and pH = 5. Photosensitized reactions of carboxylic acids and alkenes show second-order rate constants in the range of (37 ± 7.0-0.55 ± 0.1) × 104 and (27 ± 5.0-0.04 ± 0.01) × 108 L mol-1 s-1. The results show that different compound classes act differently as a photosensitizer and can be a sink for certain organic compounds in the atmospheric aqueous phase.

We provide the exact analytical solution of the ground-state hydrogenic problem with soft Coulomb potential in 1-, 2- and 3-D. We show that the wave function is an analytical function of the inverse of the soft Coulomb potential and identify a power term, an exponentially decaying term and a mildly varying modulator function on the exponential. In approaching the bare Coulomb limit, only the exponentially decaying term survives in 2D and 3D and converges to the well-known result. This is in contrast with the 1D case, where the wave function shrinks to a delta function with a total energy of minus infinity. The asymptotic behavior of the energy in such limit has been analyzed. Moreover, by analyzing the solution in different dimensions, we find that the total energy increases with dimension and scales linearly rather than quadratically with the nuclear charge Z in the large Z limit.

The complex formed from trimethylamine ((CH3)3N) and trifluoromethanesulfonic acid (triflic acid, CF3SO3H) has been observed by Fourier transform microwave spectroscopy in a supersonic jet. Spectroscopic data, most notably 14N nuclear quadrupole coupling constants, are combined with computational results at several levels of theory to unambiguously demonstrate complete or near-complete proton transfer from the triflic acid to the trimethylamine upon complexation. Thus, the system is best regarded as a trimethylammonium triflate ion pair in the gas phase. The formation of an isolated ion pair in a 1:1 complex of a Bronsted acid and base is unusual and likely arises due to the strong acidity of triflic acid. Simple energetic arguments based on proton affinities and the Coulomb interaction energy can be used to rationalize this result.

Bond dissociation energy (BDE) has been calculated for a series of compounds that contain N-O bonds. These structures encompass model N,N,O-trisubstituted hydroxylamines that include O-methoxy, O-acyl, and O-phenyl hydroxylamines. The calculations used three accurate composite methods, CBS-QB3, CBS-APNO, and G4 methods and the computationally more affordable M06-2X/6-311+G(3df,2p) density functional theory (DFT) functional. The calculated N-O single-bond BDEs are 5-15 kcal/mol higher than a generic N-O BDE of 48 kcal/mol quoted in the literature and in textbooks. The M06-2X DFT functional provides BDEs that are in excellent agreement with the higher-level composite methods. We also provided a comparison of the N-O BDE for pyridine-N-oxide to simple trialkylamine oxides. Based on an experimental BDE of 63.3 ± 0.5 kcal/mol for pyridine-N-oxide, our best estimate gives 56.7 ± 0.9 kcal/mol N-O BDE for trimethylamine-N-oxide and 59.0 ± 0.8 kcal/mol for triethylamine-N-oxide.

Formates are a class of organic molecules emitted into the atmosphere from fuel additives and industrial solvents. Formate-derived esters can undergo a vast range of chemical reactions in the atmosphere, most of which are initiated by oxidation by hydroxyl radicals. One potential reaction upon their interaction with atmospheric water is proton transfer to form protonated formates. The goal of the present work is to explore the dissociation of these protonated species and thus their possible atmospheric fate. Tandem mass spectrometry was employed to study the unimolecular dissociation of the protonated forms of methyl- (1), ethyl- (2), isopropyl- (3), tert-butyl- (4), and phenylformate (5). 1 and 2 lose CO as a common fragmentation product, forming a protonated alcohol, and 2 also generates neutral ethanol (forming protonated CO). 3 and 4 readily lose the stable isopropyl and tert-butyl radicals as well as neutral alkenes propene and isobutene. Methanol loss is also observed from both ions. 5 exhibits both phenyl radical loss (similar to 3 and 4) and CO loss (like 1 and 2). Density functional theory was used to explore the observed minimum energy reaction pathways for each ion, and CBS-QB3 single-point energy calculations provided reliable energetics.

International audience; This work investigates the potential of the particle swarm algorithm for the optimization of detailed kinetic mechanisms. To that end, empirical analysis has beenconducted to evaluate the efficiency of this algorithm in comparison with the genetic algorithm. Both algorithms are built on evolutionary processes according to which a randomly defined population will evolve, over the iterations, towards an optimal solution. The genetic algorithm is driven by crossover and mutation operators and by a selection method. The PSO approach is based on the experience of each individual and on the group experience to control the direction of its evolution. The success of the application of an algorithm can be sensitive to the choice of operators and the relative importance attributed to them. Therefore, to make the comparison as rigorous as possible, about a dozen strategies were proposed for each algorithm and the performances were evaluated. A degraded version of the GRI-Mech 3.0 mechanism (i.e. containing some of the kinetic constants randomly modified) was generated and then optimized by the two evolutionary algorithms to recover the predictive character of the original mechanism. The results show that, for the majority of the proposed strategies, PSO is more efficient than the GA, whereas the latter is generally much more used for the optimization of detailed kinetic mechanisms.

Ab initio calculations were carried out to understand the structural, electronic, and energetic properties of molybdenum oxide clusters, MomOn (m = 1-6; n = 1-3m), to understand the relationships between size, composition, and reactivity. In clusters with a low oxygen-to-molybdenum ratio, there are bridge-bonded and linearly bonded oxygen atoms on a molybdenum core, while at higher ratios, Mo atoms are separated from each other and oxygen atoms located between the molybdenum atoms. The energy gap between the highest occupied molecular orbitals (HOMOs) and lowest unoccupied molecular orbitals (LUMOs) widens with n, i.e., at a high oxygen-to-molybdenum ratio. Stoichiometric MomO3m clusters (m > 1) have a HOMO-LUMO gap that ranges from 2.6 to 3.4 eV in neutral conditions and less than 0.6 eV in ionic states. The ionization potential of MomO3m clusters is higher than 10 eV. MomOn clusters qualitatively and quantitatively exhibit a similar electronic structure to the bulk. The energy of the reduction reaction, MomOn → MomOn-1 + 1/2O2, is on average lower in clusters with high oxygen content; for example, the reduction energies of Mo6O18 and Mo6O9 are 2.23 and 5.19 eV, respectively. In the fragmentation of MomOn clusters, the general trend for clusters with a low oxygen-to-molybdenum ratio is the separation of a Mo atom or a Mo2 dimer from the cluster, while clusters with higher oxygen content mostly form stoichiometric MoO3, Mo2O6, and Mo3O9 clusters.

It is well known that the presence of an electron-withdrawing substituent (EWS) placed near the halogen (X) atom on a Lewis acid molecule amplifies the ability of this unit to engage in a halogen bond with a base. Quantum calculations are applied to examine how quickly these effects fade as the EWS is moved further and further from the X atom. Conjugated alkene and alkyne chains of varying lengths with a terminal C-I first facilitate analysis as to how the number of these multiple bonds affects the strength of CI··N XB to NH3. Then, electron-withdrawing F and C≡N substituents are placed on the opposite end of the chain, and their effects on the XB properties are monitored as a function of their distance from I. These same EWSs are added to the ortho, meta, and para positions of aromatic iodobenzene. It is found that the XB grows in strength as more triple bonds are added to the alkyne, but there is little change caused by elongating an alkene. The cyano group has a much stronger effect than does F. While F strengthens the XB, its effects are quickly attenuated as it is moved further from I. The consequences of C≡N substitution are stronger and extend over a longer distance. Placement of an EWS on the phenyl ring diminishes with distance: o > m > p, and the effects of disubstitution are nearly additive. These trends apply not only to energetics but also to geometries, properties of the wave function, σ-hole depth, and NMR shielding.

The collision diameter σ for a large set of molecular species is related to the static electronic polarizability αel. A remarkable correlation between these quantities conceptually similar to the analogous one previously identified for atoms is revealed. Our recommended model is the function σ(αel) = p1 + p2αel1/3, where p1 = 0.768 A and p2 = 2.168 are the fitting parameters providing the best overall match to the reference data for collision diameter (181 data points). The obtained correlation allows one to easily find the collision diameter of molecules from the known polarizability and vice versa. These findings can be useful for many applications, where there is a need for inexpensive assessments of the collision diameters or electronic polarizabilities, for example, when developing the transport property databases for modeling of chemically reacting flows.

We demonstrate that 2,4-bis(4,5-diphenyl-1H-imidazol-2-yl)phenol (2,4-bImP) undergoes photoinduced conversion into the so-called "π-conjugated zwitterion" after causing an excited-state intramolecular proton transfer (ESIPT) reaction. The powder sample of 2,4-bImP exhibits largely Stokes-shifted fluorescence characteristics to ESIPT fluorophores. On the other hand, its originally colorless solutions become colored when exposed to UV light for several minutes, whose color depends on the type of solvent. In particular, the CHCl3 solution rapidly turns dark green with the absorption maximum around 700 nm, and the colored solution is nearly restored to original by alternating addition of acid and base. To explain such drastic and reversible color changes, we hypothesized that the occurrence of ESIPT (i.e., deprotonation of the phenol and protonation of the imidazolyl group at its 2-position) triggered the charge-separated structure between the negatively charged phenolate and the positively charged imidazoliumyl group at its 4-position, which allowed resonance with the neutral p-quinoid structure. The formation of this π-conjugated zwitterion was strongly supported by the results of 1H and 15N NMR and Raman measurements.

Polycyclic hydrocarbons are often used to understand the electronic structure of nanographene systems. Among them, indeno[1,2b]fluorene and indeno[1,2c]fluorene isomers present a central p-quinodimethane unit leading to unique optical properties. In this work, we characterized the absorption spectra of indeno[1,2b]fluorene and [2,1-c]diindeno[n]thiophene derivatives with (spin-flip) simplified time-dependent density functional theory [(SF-)sTD-DFT] methods. Note that the SF-sTD-DFT level of theory allows a computationally efficient treatment for large diradicals. To interpret spectra, we implemented natural transition orbitals (NTOs) at both SF-sTD-DFT and sTD-DFT levels. This compact and method-independent representation of the electronic excitation provides a simple interpretation for the low-lying excited states of this set of molecules in terms of three different types of NTOs: "quinoid", "aromatic", and "π-bonded". When comparing with experiment, we found that only one molecule of this set is actually a high-spin triplet diradical. Others are almost closed-shell molecules with a very small contribution from a doubly excited configuration that only the spin-flip method could capture. The small amount of static correlation recovered by the spin-flip active space provides a linear relation between the first visible theoretical and experimental excitation energies among this set.

Photodissociation of vibrationally and electronically excited sulfur dimer molecules (S2) has been studied in a combined experimental and computational quantum chemistry study in order to characterize bound-continuum transitions. Ab initio quantum chemistry calculations are carried out to predict the potential energy curves, spin-orbit coupling, transition moments, and bound-continuum spectra of S2 for comparison with the experimental data. The experiment uses velocity map imaging to measure S-atom production following S2 photoexcitation in the ultraviolet region (320-205 nm). A pulsed electric discharge in H2S produces ground-state S2 X3ςg -(v = 0-15) as well as electronically excited singlet sulfur and b1ςg +(v = 0, 1), and evidence is presented for the production and photodissociation of S2 a1Δg. In a previous paper, we reported threshold photodissociation of S2X3ςg -(v = 0) in the 282-266 nm region. In the present study, S(3PJ) fine structure branching and angular distributions for photodissociation of S2 (X3ςg -(v = 0), a1Δg and b1ςg +) via the B″3IIu, B3ςu - and 11IIu excited states are reported. In addition, photodissociation of the X3ςg -(v = 0) state of S2 to the second dissociation limit producing S(3P2) + S(1D) is characterized. The present results on S2 photodynamics are compared to those of the well-studied electronically isovalent O2 molecule.

Small iridium nanoclusters are prominent subnanometric systems for catalysis-related applications, mainly because of a large surface-to-volume ratio, noncoalescence feature, and tunable properties, which are completely influenced by the number of atoms, geometry, and molecular interaction with the chemical environment. Herein, we investigate the interaction between Irn nanoclusters (n = 2-7) and polluting molecules, CO, NO, and SO, using van der Waals D3 corrected density functional theory calculations. Starting from a representative structural set, we determine the growth pattern of the lowest energy unprotected Irn nanoclusters, which is based on open structural motifs, and from the adsorption of a XO (X = C, N, and S) molecule, the preferred high-symmetric adsorption sites were determined, dominated by the onefold top site. For protected systems, 4XO/Ir4 and 6XO/Ir6, we found a reduction in the total magnetic moment, while the equilibrium bonds of the nanoclusters expanded (contracted) due to mCO and mNO (mSO) adsorption, with exceptions for systems with large structural distortions (4SO/Ir4 and 6NO/Ir6). Meanwhile, the C-O and N-O (S-O) bond strength decreases (increases) following an increase (decrease) in the C-O and N-O (S-O) distances upon adsorption. We show, through energetic analysis, that for the different chemical environments, relative stability changes occur from the most stable unprotected nanoclusters, planar square (Ir4), and prism (Ir6) to higher energy isomers. The change in the stability order between the two competing protected systems is feasible if the balance between the interaction energy (additive term) and distortion energies (nonadditive terms) compensates for the relative total energies of the unprotected configurations. For all systems, the interaction energy is the main reason responsible for stability alterations, except for 4SO/Ir4, where the main contribution is from a small penalty due to Ir4 distortions upon adsorption, and for 4NO/Ir4, where the energetic effects from the adsorption do not overcome the difference between the binding energies of the unprotected nanoclusters. Finally, from energy decomposition and Hirshfeld charge analysis, we find a predominant covalent nature of the physical contributions in mOX···Irn interactions with a cationic core (Irn) and an anionic shell (XO coverage).

The bonding structures of the ground state and the lowest five excited states of rhodium monoboride are identified by determining the quasi-atomic orbitals in full valence space MCSCF wave functions and the interactions between these orbitals. A quadruple bond, namely two π-bonds and two σ-bonds, is identified and characterized for the X1Σ+ ground state, in agreement with a previous report (Cheung J. Phys. Chem. Lett. 2020, 11, 659-663). However, in all excited states, the bonding is predicted to be weaker because, in these states, one of the σ-bonding interactions has a small magnitude. In the a3Δ and A1Δ states, the bond order is between a triple and quadruple bond. In the b3Σ+ state, the Rh-B linkage is a triple bond. In the c3Π and B1Π states, the atoms are linked by a double bond due to an additional weakening of the two π-bonds. The decreases in the predicted bond strengths are reflected in the decreases of the predicted binding energies and in the increases of the predicted bond lengths from the X1Σ+ ground state to the c3Π and the B1Π excited states. Notably, the 5pσ orbital of rhodium, which is vacant in the ground state of the atom, plays a significant role in the molecule.

Density functional theory (DFT) calculations were conducted to gain insight into the reaction mechanism of the Bronsted acid-catalyzed unsymmetrical 1,2,4,5-tetrazine synthesis. Various possible reaction pathways were considered, and the most favorable one can be characterized via sequential six steps, including addition of DCM to hydrazine 1 giving complex IM4, N-H bond activation in IM4 mediated by sulfur, AcOH-assisted substitution of 3 with sulfur-activated hydrazine 2, HNO2-assisted addition of nitrile to intermediate 8, cyclization, and intramolecular elimination leading to the final product 7. Among the six steps, sulfur activation of IM4 N-H bond is found to be the rate-determining step (RDS). The mechanism rationalizes the experimental observation that 2 equiv of sulfur leads to the best yield of product. Furthermore, we disclosed that the Bronsted acid additives (i.e., acetic acid and nitrous acid) served triple roles as catalyst, proton shuttle, and hydrogen bond donor and acceptor in the whole catalysis.