Impact of dielectric constant of solvent on the formation of transition metal-ammine complexes.

The DFT-level computational investigations into Gibbs free energies (ΔG) demonstrate that as the dielectric constant of the solvent increases, the stabilities of [M(NH3 )n ]2+/3+ (n = 4, 6; M = selected 3d transition metals) complexes decrease. However, there is no observed correlation between the stability of the complex and the solvent donor number. Analysis of the charge transfer and Wiberg bond indices indicates a dative-bond character in all the complexes. The solvent effect assessed through solvation energy is determined by the change in the solvent accessible surface area (SASA) and the change in the charge distribution that occurs during complex formation. It has been observed that the SASA and charge transfer are different in the different coordination numbers, resulting in a variation in the solvent effect on complex stability in different solvents. This ultimately leads to a change between the relative stability of complexes with different coordination numbers while increasing the solvent polarity for a few complexes. Moreover, the findings indicate a direct relationship between ΔΔG (∆Gsolvent -∆Ggas ) and ΔEsolv , which enables the computation of ΔG for the compounds in a particular solvent using only ΔGgas and ΔEsolv . This approach is less computationally expensive.

The Stability of Hydrogen-Bonded Ion-Pair Complex Unexpectedly Increases with Increasing Solvent Polarity.

The generally observed decrease of the electrostatic energy in the complex with increasing solvent polarity has led to the assumption that the stability of the complexes with ion-pair hydrogen bonds decreases with increasing solvent polarity. Besides, the smaller solvent-accessible surface area (SASA) of the complex in comparison with the isolated subsystems results in a smaller solvation energy of the latter, leading to a destabilization of the complex in the solvent compared to the gas phase. In our study, which combines Nuclear Magnetic Resonance, Infrared Spectroscopy experiments, quantum chemical calculations, and molecular dynamics (MD) simulations, we question the general validity of this statement. We demonstrate that the binding free energy of the ion-pair hydrogen-bonded complex between 2-fluoropropionic acid and n-butylamine (CH3CHFCOO-…NH3But+) increases with increased solvent polarity. This phenomenon is rationalized by a substantial charge transfer between the subsystems that constitute the ion-pair hydrogen-bonded complex. This unexpected finding introduces a new perspective to our understanding of solvation dynamics, emphasizing the interplay between solvent polarity and molecular stability within hydrogen-bonded systems.

The Impact of the Solvent Dielectric Constant on A←NH3 Dative Bond Depends on the Nature of the Lewis Electron-Pair Systems.

The present work aims to determine to what extent the value of the dielectric constant of the solvent can influence the dative bond in Lewis electron pair bonding systems. For this purpose, two different systems, namely H3 B←NH3 and {Zn←(NH3 )}2+ , were studied in selected solvents with significantly different dielectric constants. Based on the results from state-of-the-art computational methods using DFT, constrained DFT, energy decomposition analyses, solvent accessible surface area, and charge transfer calculations, we found that the stability of the neutral H3 B←NH3 system increases with increasing solvent polarity. In contrast, the opposite trend is observed for the positively charged {Zn←(NH3 )}2+ . The observed changes are attributed to different charge redistributions in neutral and charged complexes, which are reflected by a different response to the solvent and are quantified by changes in solvation energies.

Trends in the stability of covalent dative bonds with variable solvent polarity depend on the charge transfer in the Lewis electron-pair system.

In general, the stability of neutral complexes with dative bonds increases as the polarity of the solvent increases. This is based on the fact that the dipole moment of the complex increases as the charge transferred from the donor to the acceptor increases. As a result, the solvation energy of the complex becomes greater than that of subsystems, causing an increase in the stabilization energy with increasing solvent polarity. Our research confirms this assumption, but only when the charge transfer is sufficiently large. If it is below a certain threshold, the increase in the complex's dipole moment is insufficient to result in a higher solvation energy than subsystems. Thus, the magnitude of the charge transfer in the Lewis electron-pair system determines the stability trends of dative bonds with varying solvent polarity. We used molecular dynamics (MD) simulations based on an explicit solvent model, which is considered more reliable, to verify the results obtained with a continuous solvent model.

Addition Reaction between Piperidine and C60 to Form 1,4-Disubstituted C60 Proceeds through van der Waals and Dative Bond Complexes: Theoretical and Experimental Study.

A combined computational and experimental study reveals the character of the C60 complexes with piperidine formed under different reaction conditions. The IR and NMR experiments detect the dative bond complex, which according to NMR, is stable in the oxygen-free environment and transforms to the adduct complex in the presence of O2. Computational studies on the character of reaction channels rationalize the experimental observations. They show that the piperidine dimer rather than a single piperidine molecule is required for the complex formation. The calculations reveal significant differences in the dative bond and adduct complexes' character, suggesting a considerable versatility in their electronic properties modulated by the environment. This capability offers new application potential in several fields, such as in energy storage devices.

Structure-directed formation of the dative/covalent bonds in complexes with C70piperidine.

The combined experimental-computational study has been performed to investigate the complexes formed between C70 carbon allotrope and piperidine. The results of FT-IR, H-NMR, and C-NMR measurements, together with the calculations based on the DFT approach and molecular dynamics simulations, prove the existence of dative/covalent bonding in C70piperidine complexes. The dative bond forms not only at the region of five- and six-membered rings, observed previously with C60, but also at the region formed of six-membered rings. The structure, i.e., nonplanarity, explains the observed dative bond formation. New findings on the character of interaction of secondary amines with C70 bring new aspects for the rational design of modified fullerenes and their applications in electrocatalysis, spintronics, and energy storage.

Cyclo[n]carbons Form Strong N → C Dative/Covalent Bonds with Piperidine.

The newly synthesized C18 ring is demonstrated as the smallest all-carbon acceptor that exhibits strong electron acceptance. This study provides a quantum-chemical investigation of the electron-acceptance behavior of monocyclic carbon rings with a particular emphasis on C18 through the formation of a dative bond with piperidine. The results show that Cn rings form strong dative bonds with piperidine, whereas the respective van der Waals (vdW) complexes are higher in energy. The main driving force is the release of angle strain of cyclo[n]carbons caused by the change in hybridization from sp to sp2 associated with the formation of the dative bond. On the contrary, other sp allotropes, diynes, favorably form vdW complexes. Molecular dynamics (MD) simulations support the stability of the dative bond throughout a simulation of 20 ps. This opens up the possibility of stabilizing highly reactive C18 through dative/covalent functionalization.

An Isolated Molecule of Iron(II) Phthalocyanin Exhibits Quintet Ground-State: A Nexus between Theory and Experiment.

Iron(II) phthalocyanine (FePc) is an important member of the phthalocyanines family with potential applications in the fields of electrocatalysis, magnetic switching, electrochemical sensing, and phototheranostics. Despite the importance of electronic properties of FePc in these applications, a reliable determination of its ground-state is still challenging. Here we present combined state of the art computational methods and experimental approaches, that is, Mössbauer spectroscopy and Superconducting Quantum Interference Device (SQUID) magnetic measurements to identify the ground state of FePc. While the nature of the ground state obtained with density functional theory (DFT) depends on the functional, giving mostly the triplet state, multi-reference complete active space second-order perturbation theory (CASPT2) and density matrix renormalization group (DMRG) methods assign quintet as the FePc ground-state in gas-phase. This has been confirmed by the hyperfine parameters obtained from 57 Fe Mössbauer spectroscopy performed in frozen monochlorobenzene. The use of monochlorobenzene guarantees an isolated nature of the FePc as indicated by a zero Weiss temperature. The results open doors for exploring the ground state of other metal porphyrin molecules and their controlled spin transitions via external stimuli.

The X40×10 Halogen Bonding Benchmark Revisited: Surprising Importance of (n-1)d Subvalence Correlation.

We have re-evaluated the X40×10 benchmark for halogen bonding using conventional and explicitly correlated coupled cluster methods. For the aromatic dimers at small separation, improved CCSD(T)-MP2 "high-level corrections" (HLCs) cause substantial reductions in the dissociation energy. For the bromine and iodine species, (n-1)d subvalence correlation increases dissociation energies and turns out to be more important for noncovalent interactions than is generally realized; (n-1)sp subvalence correlation is much less important. The (n-1)d subvalence term is dominated by core-valence correlation; with the smaller cc-pVDZ-F12-PP and cc-pVTZ-F12-PP basis sets, basis set convergence for the core-core contribution becomes sufficiently erratic that it may compromise results overall. The two factors conspire to generate discrepancies of up to 0.9 kcal/mol (0.16 kcal/mol RMS) between the original X40×10 data and the present revision.

Prediction of neutral noble gas insertion compounds with heavier pnictides: FNgY (Ng = Kr and Xe; Y = As, Sb and Bi).

A novel class of interesting insertion compounds obtained through the insertion of a noble gas atom into the heavier pnictides have been explored by various ab initio quantum chemical techniques. Recently, the first neutral noble gas insertion compounds, FXeY (Y = P, N), were theoretically predicted to be stable; the triplet state was found to be the most stable state, with a high triplet-singlet energy gap, by our group. In this study, we investigated another noble gas inserted compound, FNgY (Ng = Kr and Xe; Y = As, Sb and Bi), with a triplet ground state. Density functional theory (DFT), second order Møller-Plesset perturbation theory (MP2), coupled-cluster theory (CCSD(T)) and multi-reference configuration interaction (MRCI) based techniques have been utilized to investigate the structures, stabilities, harmonic vibrational frequencies, charge distributions and topological properties of these compounds. These predicted species, FNgY (Ng = Kr and Xe; Y = As, Sb and Bi) are found to be energetically stable with respect to all the probable 2-body and 3-body dissociation pathways, except for the 2-body channel leading to the global minimum products (FY + Ng). Nevertheless, the finite barrier height corresponding to the saddle points of the compounds connected to their respective global minima products indicates that these compounds are kinetically stable. The structural parameters, energetics, and charge distribution results as well as atoms-in-molecules (AIM) analysis suggest that these predicted molecules can be best represented as F(-)[(3)NgY](+). Thus, all the aforementioned computed results clearly indicate that it may be possible to experimentally prepare the most stable triplet state of FNgY molecules under cryogenic conditions through a matrix isolation technique.

What Are the Ground State Structures of C20 and C24? An Explicitly Correlated Ab Initio Approach.

A new benchmark study has been performed for six isomers of C20 and four isomers of C24 using explicitly correlated methods, together with coupled cluster theory with large basis sets and DFT with advanced functionals. The relative energy trends obtained are extremely sensitive to the methods used. Combining our best CCSD(T)-MP2 difference with our best MP2 basis set limit, the dehydrocorannulene bowl is found to be the most stable for C20, followed by the cage at about 8 kcal/mol, and the ring at about 46 kcal/mol. For C24, the D3d cage is found to be the most stable isomer, followed at only a few kilocalories per mole by dehydrocoronene, and at larger separations by then octahedral cage and the ring, respectively. This makes C24 the smallest classical fullerene. The estimated residual basis set error of the estimated CCSD(T) basis set limit is conservatively expected to be ±1 kcal/mol. In general, DFT exhibits large errors for relative energies with RMSD values in the 8-34 kcal/mol range. However, among the DFT functionals, the DSD-PBEP86-D3BJ double hybrid comes close to our best ab initio results, while the ωB97X-V range-separated hybrid is in semiquantitative agreement.

Prediction of a neutral noble gas compound in the triplet state.

Discovery of the HArF molecule associated with H-Ar covalent bonding [Nature, 2000, 406, 874-876] has revolutionized the field of noble gas chemistry. In general, this class of noble gas compound involving conventional chemical bonds exists as closed-shell species in a singlet electronic state. For the first time, in a bid to predict neutral noble gas chemical compounds in their triplet electronic state, we have carried out a systematic investigation of xenon inserted FN and FP species by using quantum chemical calculations with density functional theory and various post-Hartree-Fock-based correlated methods, including the multireference configuration interaction technique. The FXeP and FXeN species are predicted to be stable by all the computational methods employed in the present work, such as density functional theory (DFT), second-order Møller-Plesset perturbation theory (MP2), coupled-cluster theory (CCSD(T)), and multireference configuration interaction (MRCI). For the purpose of comparison we have also included the Kr-inserted compounds of FN and FP species. Geometrical parameters, dissociation energies, transition-state barrier heights, atomic charge distributions, vibrational frequency data, and atoms-in-molecules properties clearly indicate that it is possible to experimentally realize the most stable state of FXeP and FXeN molecules, which is triplet in nature, through the matrix isolation technique under cryogenic conditions.

Noble-Gas-Inserted Fluoro(sulphido)boron (FNgBS, Ng = Ar, Kr, and Xe): A Theoretical Prediction.

The possibility of the existence of a new series of neutral noble gas compound, FNgBS (where Ng = Ar, Kr, Xe), is explored theoretically through the insertion of a Ng atom into the fluoroborosulfide molecule (FBS). Second-order Møller-Plesset perturbation theory, density functional theory, and coupled cluster theory based methods have been employed to predict the structure, stability, harmonic vibrational frequencies, and charge distribution of FNgBS molecules. Through energetics study, it has been found that the molecules could dissociate into global minima products (Ng + FBS) on the respective singlet potential energy surface via a unimolecular dissociation channel; however, the sufficiently large activation energy barriers provide enough kinetic stability to the predicted molecules, which, in turn, prevent them from dissociating into the global minima products. Moreover, the FNgBS species are thermodynamically stable, owing to very high positive energies with respect to other two two-body dissociation channels, leading to FNg + BS and F(-) + NgBS(+), and two three-body dissociation channels, corresponding to the dissociation into F + Ng + BS and F(-) + Ng + BS(+). Furthermore, the Mulliken and NBO charge analysis together with the AIM results reveal that the Ng-B bond is more of covalent in nature, whereas the F-Ng bond is predominantly ionic in character. Thus, these compounds can be better represented as F(-)[NgBS](+). This fact is also supported by the detail analysis of bond length, bond dissociation energy, and stretching force constant values. All of the calculated results reported in this work clearly indicate that it might be possible to prepare and characterize the FNgBS molecules in cryogenic environment through matrix isolation technique by using a mixture of OCS/BF3 in the presence of large quantity of noble gas under suitable experimental conditions.

Theoretical prediction of noble gas inserted thioformyl cations: HNgCS⁺ (Ng = He, Ne, Ar, Kr, and Xe).

The existence of new interesting insertion compounds, HNgCS(+) (Ng = He-Xe), have been predicted theoretically through insertion of a noble gas atom into the thioformyl cation, HCS(+). Second-order Møller-Plesset perturbation theory (MP2), density functional theory (DFT), and coupled-cluster theory (CCSD(T)) based techniques have been used to explore the structure, energetics, charge distribution, and harmonic vibrational frequencies of these compounds. These predicted ions are found to be energetically stable with respect to all the possible 2-body and 3-body dissociation pathways, except the 2-body channel leading to the global minimum products (HCS(+) + Ng). Nevertheless, all these ions are found to be kinetically stable with a finite barrier height corresponding to their transition states, which are connected to their respective global minima products. The results obtained from charge distribution as well as atoms in molecules (AIM) analysis suggest that these ions can be best described as [HNg](+)CS. Strong covalent character in the H-Ng bond is supported by the high positive energy value corresponding to the 3-body dissociation pathways. Thus, it might be possible to prepare the HNgCS(+) ions in a glow discharge containing H2S, CO, and noble gas under cryogenic conditions through matrix isolation technique.

N-Heterocyclic Carbene-Catalyzed Facile Synthesis of Phthalidyl Sulfonohydrazones: Density Functional Theory Mechanistic Insights and Docking Interactions.

N-heterocyclic carbene catalysis reaction protocol is disclosed for the synthesis of phthalidyl sulfonohydrazones. A broad range of N-tosyl hydrazones react effectively with phthalaldehyde derivatives under open-air conditions, enabling the formation of a new C-N bond via an oxidative path. The reaction proceeds under mild reaction conditions with broad substrate scopes, wide functional group tolerance, and good to excellent yields. The mechanistic pathway is studied successfully using control experiments, competitive reactions, ESI-MS spectral analyses of the reaction mixture, and computational study by density functional theory. The potential use of one of the phthalidyl sulfonohydrazone derivatives as the inhibitor of ß-ketoacyl acyl carrier protein synthase I of Escherichia coli is investigated using molecular docking.

Theoretical prediction of XRgCO(+) ions (X = F, Cl, and Rg = Ar, Kr, Xe).

In this work we have predicted novel rare gas containing cationic molecules, XRgCO(+) (X = F, Cl and Rg = Ar, Kr, Xe) using ab initio quantum chemical methods. Detail structural, stability, vibrational frequency, and charge distribution values are reported using density functional theory, second-order Møller-Plesset perturbation theory, and coupled-cluster theory based methods. These ions are found to be metastable in nature and exhibit a linear geometry with C∞v symmetry in their minima energy structures, and the nonlinear transition state geometries are associated with Cs symmetry. Except for the two-body dissociation channel (Rg + XCO(+)), these ions are stable with respect to all other dissociation channels. However, the connecting transition states between the above-mentioned two-body dissociation channel products and the predicted ions are associated with sufficient energy barriers, which restricts the metastable species to transform into the global minimum products. Thus, it may be possible to detect and characterize these metastable ions using an electron bombardment technique under cryogenic conditions.

Theoretical prediction of rare gas containing hydride cations: HRgBF+ (Rg = He, Ar, Kr, and Xe).

The existence of rare-gas-containing hydride ions of boron (HRgBF(+)) has been predicted by using ab initio quantum chemical methods. The HRgBF(+) ions are obtained by inserting a rare gas (Rg) atom in between the H and B atoms of a HBF(+) ion, and the geometries are optimized for minima as well as transition states using second-order Møller-Plesset perturbation theory (MP2), density functional theory (DFT), and coupled-cluster theory (CCSD(T)) based techniques. The predicted HRgBF(+) ions are found to be metastable, and they exhibit a linear structure at the minima and a nonlinear planar structure at the transition state, corresponding to C∞v and Cs symmetries, respectively. All of the predicted HRgBF(+) ions show negative binding energies with respect to the two-body dissociation channel, leading to global minima (HBF(+) + Rg) on the singlet potential energy surface. In contrast, the dissociation energies corresponding to another two-body dissociation channel leading to HRg(+) + BF and two three-body dissociation channels corresponding to the dissociation into H + Rg + BF(+) and H(+) + Rg + BF show very high positive energies. Apart from positive dissociation energies, the predicted ions show finite barrier heights corresponding to the transition states involving a H-Rg-B bending mode, leading to the global minima products (HBF(+) + Rg). The finite barrier heights in turn would prevent the metastable HRgBF(+) species from transforming to global minima products. Structure, harmonic vibrational frequencies, stability, and Mulliken and natural bonding orbital (NBO) charge distribution values for all of the species are reported using the MP2 and DFT methods. Furthermore, the intrinsic reaction coordinate analysis confirms that the metastable minimum-energy structure and the global minimum products are connected through the corresponding transition state for each of the species on the respective singlet potential energy surface. Atoms-in-molecules (AIM) analysis indicates that the HRgBF(+) ions are best described as HRg(+)BF and are analogous to the isoelectronic HRgCO(+) and HRgN2(+) ions. The energetic along with charge redistribution and spectroscopic data strongly support the possible existence of HRgBF(+) ions. Hence, it might be possible to generate HRgBF(+) ions in the DC discharge plasma of a BF3/H2/Rg mixture at low temperature, and the predicted ions may be characterized using the magnetic field modulated infrared laser spectroscopic technique, which has been used earlier to characterize HBF(+) ions.

Theoretical prediction of rare gas inserted hydronium ions: HRgOH2(+).

A possibility of existence of new species through insertion of a rare gas atom in hydronium ion resulting into HRgOH2(+) cation (Rg = He, Ar, Kr, and Xe) has been explored by using various ab initio quantum chemical techniques. Structure, harmonic vibrational frequencies, stability, and charge distribution of HRgOH2(+) species as obtained using density functional theory, second order Møller-Plesset perturbation theory, and coupled-cluster theory based methods are reported in this work. All the calculated results suggest that the HRgOH2(+) species are stable enough with respect to all the dissociation channels, except the 2-body dissociation path (H3O(+) + Rg). Nevertheless, this 2-body dissociation channel connected through the relevant transition state is associated with a finite barrier, which in turn would prevent the metastable species in transforming to global minimum products. The calculated values of topological properties within the framework of quantum theory of atoms-in-molecules are found to be consistent with the bond length values. Structural and energetic parameters clearly suggest that it might be possible to prepare and characterize the HRgOH2(+) species (except HHeOH2(+)) using electron bombardment matrix isolation technique in a way similar to that of the preparation of (Rg2H)(+) or mixed (RgHRg('))(+) cations.

Coexistence of Electrochromism and Bipolar Nonvolatile Memory in a Single Viologen.

Viologens are fascinating redox-active organic compounds that have been widely explored in electrochromic devices (ECDs). However, the combination of electrochromic and resistive random-access memory in a single viologen remains unexplored. We report the coexistence of bistate electrochromic and single-resistor (1R) memory functions in a novel viologen. A high-performance electrochromic function is achieved by combining viologen (BzV2+2PF6) with polythiophene (P3HT), enabling a "push-pull" electronic effect due to the efficient intermolecular charge transfer in response to an applied bias. The ECDs show high coloration efficiency (ca. 1150 ± 10 cm2 C-1), subsecond switching time, good cycle stability (>103 switching cycles), and low-bias operation (±1.5 V). The ECDs require low power for switching the color states (55 µW cm-2 for magenta and 141 µW cm-2 for blue color). The random-access memory devices (p+2-Si/BzV2+2PF6/Al) exhibit distinct low and high resistive states with an ON/OFF ratio of ∼103, bipolar and nonvolatile characteristics that manifest good performances, and "Write"-"Read"-"Erase" (WRE) functions. The charge conduction mechanism of the RRAM device is elucidated by the Poole-Frenkel model where SET and RESET states arise at a low transition voltage (VT = ±1.7 V). Device statistics and performance parameters for both electrochromic and memory devices are compared with the literature data. Our findings on electrochromism and nonvolatile memory originated in the same viologen could boost the development of multifunctional, smart, wearable, flexible, and low-cost optoelectronic devices.