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.

Similarities and Differences of Hydridic and Protonic Hydrogen Bonding.

Ab initio calculations were employed to investigate the interactions between selected electron-donating groups, characterized by M-H bonds (where M represents a transition metal and H denotes a hydridic hydrogen), and electron-accepting groups featuring both σ- and π-holes. The study utilized the ωB97X-D3BJ/def2-TZVPPD level of theory. Hydridic hydrogen complexes were found in all complexes with σ- and π-holes. A comparative analysis was conducted on the properties hydridic H-bond complexes, presented here and those studied previously, alongside an extended set of protonic H-bonds complexes. While the stabilization energies changes in M-H bond lengths, vibrational frequencies, intensities of the spectral bands, and charge transfer for these complexes are comparable, the nature of hydridic and protonic H-bonds fundamentally differ. In protonic H-bond complexes, the main stabilization forces arise from electrostatic contributions, while in hydridic H-bond complexes, dispersion energy, is the primary stabilization factor due to the excess of electrons and thus larger polarizability at hydridic H. The finding represents an important characteristic that distinguishes hydridic H-bonding from protonic H-bonds.

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.

Controlling On-Surface Photoactivity: The Impact of π-Conjugation in Anhydride-Functionalized Molecules on a Semiconductor Surface.

On-surface synthesis has become a prominent method for growing low-dimensional carbon-based nanomaterials on metal surfaces. However, the necessity of decoupling organic nanostructures from metal substrates to exploit their properties requires either transfer methods or new strategies to perform reactions directly on inert surfaces. The use of on-surface light-induced reactions directly on semiconductor/insulating surfaces represents an alternative approach to address these challenges. Here, exploring the photochemical activity of different organic molecules on a SnSe semiconductor surface under ultra-high vacuum, we present a novel on-surface light-induced reaction. The selective photodissociation of the anhydride group is observed, releasing CO and CO2. Moreover, we rationalize the relationship between the photochemical activity and the π-conjugation of the molecular core. The different experimental behaviour of two model anhydrides was elucidated by theoretical calculations, showing how the molecular structure influences the distribution of the excited states. Our findings open new pathways for on-surface synthesis directly on technologically relevant substrates.

The unusual stability of H-bonded complexes in solvent caused by greater solvation energy of complex compared to those of isolated fragments.

Here, the effect of solvent on the stability of non-covalent complexes, was studied. These complexes were from previously published S22, S66, and X40 datasets, which include hydrogen-, halogen- and dispersion-bonded complexes. It was shown that the charge transfer in the complex determines whether the complex is stabilized or destabilized in solvent.

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.

Stability and Reactivity of the Phenalene and Olympicene Isomers.

The phenalene (triangulene) and olympicene molecules belong to the polycyclic aromatic hydrocarbon class, which have attracted substantial technological interest due to their unique electronic properties. Electronic structure calculations serve as a valuable tool in investigating the stability and reactivity of these molecular systems. In the present work, the multireference calculations, namely, the complete active space second-order perturbation theory and multireference averaged quadratic coupled cluster (MR-AQCC), were employed to study the reactivity and stability of phenalene and olympicene isomers, as well as their modified structures where the sp3-carbon at the borders were removed. The harmonic oscillator model of aromaticity (HOMA) and the nucleus-independent chemical shift as geometric and magnetic indexes calculated with density functional theory were utilized to assess the aromaticity of the studied molecules. These indexes were compared with properties such as the excitation energy and natural orbitals occupation. The reactivity analyzed using the HOMA index combined with MR-AQCC revealed the radical character of certain structures as well as the weakening of their aromaticity. Moreover, the results suggest that the removal of sp3-carbon atoms and the addition of hydrogen atoms did not alter the π network and the excitation energies of the phenalene molecules.

Pathways to fluorescence via restriction of intramolecular motion in substituted tetraphenylethylenes.

The design of materials with enhanced luminescence properties is a fast-developing field due to the potential applicability of these materials as light-emitting diodes or for bioimaging. A transparent way to enhance the emission properties of interesting molecular candidates is blocking competing and unproductive non-radiative relaxation pathways by the restriction of intramolecular motions. Rationalized functionalization is an important possibility to achieve such restrictions. Using time-dependent density functional theory (TD-DFT) based on the ωB97XD functional and the semiempirical tight-binding method including long-range corrections (TD-LC-DFTB), this work investigates the effect of functionalization of the paradigmatic tetraphenylethylene (TPE) on achieving restricted access to conical intersections (RACI). Photodynamical surface hopping simulations have been performed on a larger set of compounds including TPE and ten functionalized TPE compounds. Functionalization has been achieved by means of electron-withdrawing groups, bulky groups which block the relaxation channels via steric hindrance and groups capable of forming strong hydrogen bonds, which restrict the motion via the formation of hydrogen bond channels. Most of the investigated functionalized TPE candidates show ultrafast deactivation to the ground state due to their still existing structural flexibility, but two examples, one containing -CN and -CF3 groups and a second characterized by a network of hydrogen bonds, have been identified as interesting candidates for creating efficient luminescence properties in solution.

Unexpected Strengthening of the H-Bond Complexes in a Polar Solvent Due to a More Efficient Solvation of the Complex Compared to Isolated Monomers.

It is generally assumed that hydrogen-bonded complexes are less stable in solvents than in the gas phase and that their stability decreases with increasing solvent polarity. This assumption is based on the size of the area available to the solvent, which is always smaller in the complex compared to the subsystems, thereby reducing the solvation energy. This reduction prevails over the amplification of the electrostatic hydrogen bond by the polar solvent. In this work, we show, using experimental IR spectroscopy and DFT calculations, that there are hydrogen-bonded complexes whose stability becomes greater with increasing solvent polarity. The explanation for this surprising stabilization is based on the analysis of the charge redistribution in the complex leading to increase of its dipole moment and solvation energy. Constrained DFT calculations have shown a dominant role of charge transfer over polarization effects for dipole moment and solvation energy.

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.

On-Surface Strain-Driven Synthesis of Nonalternant Non-Benzenoid Aromatic Compounds Containing Four- to Eight-Membered Rings.

The synthesis of polycyclic aromatic hydrocarbons containing various non-benzenoid rings remains a big challenge facing contemporary organic chemistry despite a considerable effort made over the last decades. Herein, we present a novel route, employing on-surface chemistry, to synthesize nonalternant polycyclic aromatic hydrocarbons containing up to four distinct kinds of non-benzenoid rings. We show that the surface-induced mechanical constraints imposed on strained helical reactants play a decisive role leading to the formation of products, energetically unfavorable in solution, with a peculiar ring current stabilizing the aromatic character of the π-conjugated system. Determination of the chemical and electronic structures of the most frequent product reveals its closed-shell character and low band gap. The present study renders a new route for the synthesis of novel nonalternant polycyclic aromatic hydrocarbons or other hydrocarbons driven by internal stress imposed by the surface not available by traditional approaches of organic chemistry in solution.

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.

The Existence of a N→C Dative Bond in the C60 -Piperidine Complex.

The complexes formed between carbon allotropes (C20 , C60 fullerenes, graphene, and single-wall carbon nanotubes) and piperidine have been investigated by means of computational quantum chemical and experimental IR and NMR techniques. Alongside hydrogen bonds, the C⋅⋅⋅N tetrel bond, and lone-pair⋅⋅⋅π interactions, the unexpected N→C dative/covalent bond has been detected solely in complexes of fullerenes with piperidine. Non-planarity and five-member rings of carbon allotropes represent the key structural prerequisites for the unique formation of a dative N→C bond. The results of thermodynamics calculations, molecular dynamics simulations, and NMR and FTIR spectroscopy explain the specific interactions between C60 and piperidine. The differences in behavior of individual carbon allotropes in terms of dative bonding formation brings a new insight into their controllable organic functionalization.

Mechanisms of Orthogonal Photodecarbonylation Reactions of 3-Hydroxyflavone-Based Acid-Base Forms.

Carbon monoxide is a naturally occurring gasotransmitter combining inherent toxicity with a remarkable therapeutic potential and arduous administration. Photoactivatable carbon monoxide-releasing molecules (photoCORMs) are chemical agents that allow for precise spatial and temporal control over the CO release. In this work, we present a comprehensive mechanistic study of the photochemical CO release from 3-hydroxy-2-phenyl-4H-chromen-4-one, a π-extended 3-hydroxyflavone photoCORM, in methanol using steady-state and transient absorption spectroscopies and quantum chemical calculations. The multiplicity of the productive excited states and the role of oxygen (O2) in the CO production are emphasized, revealing a photoreaction dichotomy of the 3-hydroxyflavone acid and base forms. The utilization of three major orthogonal mechanistic pathways, all of which lead to the CO release, can fuel future endeavors to improve the CO release efficacy of 3-hydroxyflavone-based derivatives and refine their potential medical applications as photoCORMs.

Multireference Approaches for Excited States of Molecules.

Understanding the properties of electronically excited states is a challenging task that becomes increasingly important for numerous applications in chemistry, molecular physics, molecular biology, and materials science. A substantial impact is exerted by the fascinating progress in time-resolved spectroscopy, which leads to a strongly growing demand for theoretical methods to describe the characteristic features of excited states accurately. Whereas for electronic ground state problems of stable molecules the quantum chemical methodology is now so well developed that informed nonexperts can use it efficiently, the situation is entirely different concerning the investigation of excited states. This review is devoted to a specific class of approaches, usually denoted as multireference (MR) methods, the generality of which is needed for solving many spectroscopic or photodynamical problems. However, the understanding and proper application of these MR methods is often found to be difficult due to their complexity and their computational cost. The purpose of this review is to provide an overview of the most important facts about the different theoretical approaches available and to present by means of a collection of characteristic examples useful information, which can guide the reader in performing their own applications.

Tuning the UV spectrum of PAHs by means of different N-doping types taking pyrene as paradigmatic example: categorization via valence bond theory and high-level computational approaches.

Tuning of the electronic spectra of carbon dots by means of inserting heteroatoms into the π-conjugated polycyclic aromatic hydrocarbon (PAH) system is a popular tool to achieve a broad range of absorption and emission frequencies. Especially nitrogen atoms have been used successfully for that purpose. Despite the significant progress achieved with these procedures, the prediction of specific shifts in the UV-vis spectra and the understanding of the electronic transitions is still a challenging task. In this work, high-level quantum chemical methods based on multireference (MR) and single-reference (SR) methods have been used to predict the effect of different nitrogen doping patterns inserted into the prototypical PAH pyrene on its absorption spectrum. Furthermore, a simple classification scheme based on valence bond (VB) theory and the Clar sextet rule in combination with the harmonic oscillator measure of aromaticity (HOMA) index was applied to arrange the different doping structures into groups and rationalize their electronic properties. The results show a wide variety of mostly redshifts in the spectra as compared to the pristine pyrene case. The most interesting doping structures with the largest red shifts leading to absorption energies below one eV could be readily explained by the occurrence of diradical VB structures in combination with Clar sextets. Moreover, analysis of the electronic transitions computed with MR methods showed that several of the low-lying excited states possess double-excitation character, which cannot be realized by the popular SR methods and, thus, are simply absent in the calculated spectra.

Ground state of the Fe(ii)-porphyrin model system corresponds to quintet: a DFT and DMRG-based tailored CC study.

Fe(ii)-porphyrins play an important role in many reactions relevant to material science and biological processes, due to their closely lying spin states. Although the prevalent opinion is that these systems posses the triplet ground state, the recent experiment on Fe(ii)-phthalocyanine under conditions matching those of an isolated molecule points toward the quintet ground state. We present a thorough DFT and DMRG-based tailored CC study of Fe(ii)-porphyrin model, in which we address all previously discussed correlation effects. We examine the importance of geometrical parameters, the Fe-N distances in particular, and conclude that the system possesses the quintet ground state.

Doping Capabilities of Fluorine on the UV Absorption and Emission Spectra of Pyrene-Based Graphene Quantum Dots.

Functionalization of quantum carbon dots (QCDs) and graphene quantum dots (GQDs) is a popular way to tune their optical spectra increasing their potential applicability in material science and biorelated disciplines. Based on the experimental observation, functionalization by fluorine atoms induces substantial shifts in absorption and emission spectra and an intensity increase. Understanding of the effects due to fluorine functionalization at the atomic scale level is still challenging due to the complex structure of fluorinated QCDs. In this work, the effect of covalent edge-fluorination and fluorine anion doping on absorption and emission spectra of prototypical polycyclic aromatic hydrocarbons pyrene and circum-pyrene has been investigated. The ways to achieve efficient red-shifts in the UV spectra and obtaining reasonable intensities stood in the focus of the work. High-level quantum chemical methods based on density functional theory/multireference configuration interaction (DFT/MRCI) and single-reference second-order algebraic diagrammatic construction (ADC(2)) and density functional theory (DFT) using the CAM-B3LYP functional have been used for this purpose. The calculations show that doping with the fluoride anion can have significant effects on the electronic spectrum. However, the effect of the fluoride ion is strongly dependent on its position with respect to the QCD. The localization above the GQDs causes large red-shifts to both the absorption and emission of spectra of GQDs, while in-plane localization leads to only negligible shifts and a tendency to dissociation after electronic excitation. Thus, large red-shifts, observed in complexes with F-, are obtained due to the introduction of new excited states with large CT character not yet been considered previously in this context, although they have the potential to significantly influence the photophysics of quantum dots. Doping by edge fluorination redshifts the spectra only slightly. This study provides insights on fluorine-doped GQDs, which is conducive to promoting its rational design and controllable synthesis.

The generality of the GUGA MRCI approach in COLUMBUS for treating complex quantum chemistry.

The core part of the program system COLUMBUS allows highly efficient calculations using variational multireference (MR) methods in the framework of configuration interaction with single and double excitations (MR-CISD) and averaged quadratic coupled-cluster calculations (MR-AQCC), based on uncontracted sets of configurations and the graphical unitary group approach (GUGA). The availability of analytic MR-CISD and MR-AQCC energy gradients and analytic nonadiabatic couplings for MR-CISD enables exciting applications including, e.g., investigations of π-conjugated biradicaloid compounds, calculations of multitudes of excited states, development of diabatization procedures, and furnishing the electronic structure information for on-the-fly surface nonadiabatic dynamics. With fully variational uncontracted spin-orbit MRCI, COLUMBUS provides a unique possibility of performing high-level calculations on compounds containing heavy atoms up to lanthanides and actinides. Crucial for carrying out all of these calculations effectively is the availability of an efficient parallel code for the CI step. Configuration spaces of several billion in size now can be treated quite routinely on standard parallel computer clusters. Emerging developments in COLUMBUS, including the all configuration mean energy multiconfiguration self-consistent field method and the graphically contracted function method, promise to allow practically unlimited configuration space dimensions. Spin density based on the GUGA approach, analytic spin-orbit energy gradients, possibilities for local electron correlation MR calculations, development of general interfaces for nonadiabatic dynamics, and MRCI linear vibronic coupling models conclude this overview.