Effect of substituent and solvent on cation-π interactions in benzene and borazine: a computational study

Cation-lone Pair Interaction in Alkali/Alkaline Earth Metal Ion-Heavier Borazine Analogue Complexes

The present study is the first report on the formation of alkali/alkaline earth metal ion-heavier borazine analogue complexes via cation-lone pair interaction. Density functional calculations are performed in scrutinizing the complex formation between alkali (Li+, Na+, K+)/alkaline earth (Be2+, Mg2+, Ca2+) metal ions and heavier borazine analogues (HBA) viz. B3P3H6, Al3N3H6, Al3P3H6, Al3As3H6, and Ga3P3H6. The complexes are found to be stable in gas phase with stabilization energies within the range 26.40-324.74 kcal mol-1. The stability can be attributed to the polarizing power of the involved metal ions. Presence of solvent phase exerted notable impact on the stability of the complexes; stability is reduced significantly with the increase in solvent polarity. The process of complexation is exothermic and spontaneous. QTAIM analysis indicated the presence of both ionic and covalent interaction between HBAs and metal ions. HOMO energy, Wiberg bond index, NCI-isosurface and RDG plot analysis revealed the major role of cation-lone pair interaction in the complexation process.

Computational insight into the formation of cation-π/cation-lone pair complexes between 3d-metal (II) ions and furan

CONTEXT

Cation-π and cation-lone pair interactions between 3d-metal (II) ions [Fe(II), Co(II), Ni(II) and Cu(II)] and furan are explored in the formation of 1:1 and 1:2 type complexes. Both cation-π (IEgas = -192.27 to -312.65 kcal mol-1) and cation-lone pair (IEgas = -163.13 to -271.76 kcal mol-1) interactions are reasonably strong and lead to the formation of stable 1:1 and 1:2 type complexes in gas phase. The complexes are also stable in solvent phase, but their stability is reduced significantly in presence of solvent dielectrics, especially in ethanol, DMSO and water. Formation of the complexes is thermodynamically favourable (exothermic and spontaneous). Charge transfer (Δq = 0.62 to 1.92 e-), Laplacian of electron density (∇2ρ = 0.1435 to 0.6628 au) and total electron energy density (H(r) = -0.0019 to -0.0436 au) analysis have argued in favour of partial ionic and partial covalent character of the interactions.

METHODS

Density functional theory (DFT) is exclusively used for the study. Polarizable continuum model (PCM) is used to perform solvent phase study. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) and non-covalent interaction (NCI) analyses are performed for understanding other aspects of complex formation.

Computational Library Enables Pattern Recognition of Noncovalent Interactions and Application as a Modern Linear Free Energy Relationship

A quantitative and predictive understanding of how attractive noncovalent interactions (NCIs) influence functional outcomes is a long-standing goal in mechanistic chemistry. In that context, better comprehension of how substituent effects influence NCI strengths, and the origin of those effects, is still needed. We sought to build a resource capable of elucidating fundamental origins of substituent effects in NCIs and diagnosing NCIs in chemical systems. To accomplish this, a library of 893 NCI energies was calculated encompassing cation-π, anion-π, CH-π, and π-π interactions across 60 different arenes and heteroarenes. The interaction energies (IEs) were calculated using symmetry-adapted perturbation theory (SAPT), which identifies electrostatic, inductive, exchange-repulsive, and dispersive contributions to total IE. This descriptor library provides a comprehensive platform for evaluating substituent effect trends beyond traditional molecular descriptors such as Hammett values, frontier molecular orbital energies, and electrostatic potential, thereby expanding the tools available to analyze modern chemical processes that involve NCIs. To demonstrate the application of this library, three case studies in asymmetric catalysis and supramolecular chemistry are presented. These case studies informed the development of an automated NCI analysis tool, which employs statistical analyses to diagnose a particular NCI in a chemical system of interest.

A theoretical study on 1H-indole-2,3-dione complexes with lithium, sodium, and potassium cations

CONTEXT

A comparative study of the change in different properties of electronic and structural of the free 1H-indole-2,3-dione molecule and its complexes has been obtained. HOMA analysis was performed to investigate the effects of lithium sodium and potassium cations on the aromaticity of lithium sodium and potassium complexes of 1H-indole-2,3-dione.

METHODS

Several 1H-indole-2,3-dione complexes with lithium, sodium, and potassium cations were optimized at the B3LYP/6-311G(d,p) level. The cation and π interaction has been investigated from different aspects, including interaction energy calculations, charge transfer values, and changes in the aromaticity of the ring upon complexation. The charge transfer and natural population analysis for the complexes were performed with the natural bond orbital (NBO) analysis. The properties of bond critical points in complexes were studied by applying the quantum theory of atoms in molecules (QTAIM). Finally, the aromaticity change of phenyl induced upon complex formation was evaluated by applying the harmonic oscillator model of aromaticity (HOMA). [Li-INa]+ and [[Li-INb]+ were optimized with the wB97XD function using a version of Grimme's D2 dispersion model, and the absorption energy was compared with the calculation made with the B3LYP functional.

Ab Initio Molecular Dynamics Investigation of the Solvation States of Hydrated Ions in Confined Water

Ionic transport in nanoscale channels with a critical size comparable to that of ions and solutes exhibits exceptional performance in water desalination, ion separation, electrocatalysts, and supercapacitors. However, the solvation states (SSs), i.e., the hydration structures and probability distribution, of hydrated ions in nanochannels differ from those in the bulk and the perspective of continuum theory. In this work, we conduct ab initio enhanced-sampling atomistic simulations to investigate the ion-specific SSs of monovalent ions (including Li+, Na+, K+, F-, Cl-, and I-) in the graphene channel with a width of 1 nm. Our findings highlight that the SSs of those ions are primarily determined by ion-water hydration, where ion-wall interactions play a minor role. The distribution of ions in layered confined water is a result of ion-specific hydration, which arises from the synergy of entropy and enthalpy. The free energy barriers for transitions between SSs are on the order of 1kBT, allowing for modulation through applying external fields or modifying surface properties. As the ion-wall interaction strengthens, as observed in vermiculite and carbides and nitrides of transition metal channels, the probability of near-wall SSs increases. These results help to improve the performance of nanofluidic devices and provide crucial insights for developing accurate force fields of molecular simulations or advanced theoretical approaches for ion dynamics in confined channels.

Bi-terminal fusion of intrinsically-disordered mussel foot protein fragments boosts mechanical strength for protein fibers

Microbially-synthesized protein-based materials are attractive replacements for petroleum-derived synthetic polymers. However, the high molecular weight, high repetitiveness, and highly-biased amino acid composition of high-performance protein-based materials have restricted their production and widespread use. Here we present a general strategy for enhancing both strength and toughness of low-molecular-weight protein-based materials by fusing intrinsically-disordered mussel foot protein fragments to their termini, thereby promoting end-to-end protein-protein interactions. We demonstrate that fibers of a ~60 kDa bi-terminally fused amyloid-silk protein exhibit ultimate tensile strength up to 481 ± 31 MPa and toughness of 179 ± 39 MJ*m^-3, while achieving a high titer of 8.0 ± 0.70 g/L by bioreactor production. We show that bi-terminal fusion of Mfp5 fragments significantly enhances the alignment of β-nanocrystals, and intermolecular interactions are promoted by cation-π and π-π interactions between terminal fragments. Our approach highlights the advantage of self-interacting intrinsically-disordered proteins in enhancing material mechanical properties and can be applied to a wide range of protein-based materials.

CIPDB: A biological structure databank for studying cation and π interactions

As major forces for modulating protein folding and molecular recognition, cation and π interactions are extensively identified in protein structures. They are even more competitive than hydrogen bonds in molecular recognition, thus, are vital in numerous biological processes. In this review, we introduce the methods for the identification and quantification of cation and π interactions, provide insights into the characteristics of cation and π interactions in the natural state, and reveal their biological function together with our developed database (Cation and π Interaction in Protein Data Bank; CIPDB; http://chemyang.ccnu.edu.cn/ccb/database/CIPDB). This review lays the foundation for the in-depth study of cation and π interactions and will guide the use of molecular design for drug discovery.

Exploring alkali metal cation⋯hydrogen interaction in the formation half sandwich complexes with cycloalkanes: a DFT approach

Gas and solvent phase stability of half sandwich complexes between cycloalkanes viz. cyclopropane, cyclobutane, cyclopentane, cyclohexane, bicyclo[2.2.2]octane and adamantane with alkali metal cations (Li+, Na+ and K+) are analysed using density functional theory (DFT). M06-2X/6-31++G(d,p) level is primarily used for the study. The studied half sandwich complexes are stable in gas phase (stabilization energy upto 26.55 kcal mol−1). Presence of solvent phase irrespective of its dielectric, imparts negative impact on the stability of the chosen complexes. The formation of the complexes is exothermic in nature. The process of complexation is both enthalpy (ΔH) and free energy (ΔG) driven. Variation in HOMO (highest occupied molecular orbital) energy also indicates towards the chemical stability of complexes. The interaction is non-covalent with primary contribution from induction component. NBO analysis indicates that C–H bond is the donor and antibonding metal orbital is the acceptor site in the process of complexation. Stability of the complexes depends on the size of the interacting monomers.

Harnessing the cation-π interactions of metalated gold monolayer-protected clusters to detect aromatic volatile organic compounds

The strong, non-covalent interactions between π-systems and cations have been the focus of numerous studies on biomolecule structure and catalysis. These interactions, however, have yet to be explored as a sensing mechanism for detecting trace levels of volatile organic compounds (VOCs). In this article, we provide evidence that cation-π interactions can be used to elicit sensitive and selective chemiresistor responses to aromatic VOCs. The chemiresistors are fitted with carboxylate-linked alkali metals bound to the surface of gold monolayer-protected clusters formulated on microfabricated interdigitated electrodes. Sensor responses to aromatic and non-aromatic VOCs are consistent with a model for cation-π interactions arising from association of electron-rich aromatic π-systems to metal ions with the relative strength of attraction following the order K⁺ > Na⁺ > Li⁺. The results point toward cation-π interactions as a promising research avenue to explore for developing aromatic VOC-selective sensors.

Formation of sandwich and multidecker complexes between O2 and alkali/alkaline earth metals: A DFT study

Abstract: Feasibility of formation of sandwich and multidecker complexes between O2 molecules and alkali/alkaline earth metal has been analyzed in the light of density functional theory (DFT). High value of...

Anion-Anion Interactions in Aerogen-Bonded Complexes. Influence of Solvent Environment

Ab initio calculations are applied to the question as to whether a AeX₅^- anion (Ae = Kr, Xe) can engage in a stable complex with another anion: F^-, Cl^-, or CN^-. The latter approaches the central Ae atom from above the molecular plane, along its C₅ axis. While the electrostatic repulsion between the two anions prevents their association in the gas phase, immersion of the system in a polar medium allows dimerization to proceed. The aerogen bond is a weak one, with binding energies less than 2 kcal/mol, even in highly polar aqueous solvent. The complexes are metastable in the less polar solvents THF and DMF, with dissociation opposed by a small energy barrier.

Substituent effects in the so-called cationπ interaction of benzene and its boron-nitrogen doped analogues: overlooked role of σ-skeleton

Despite massive efforts to pinpoint the substituent effects in the so-called cationπ systems, no consensus has been yet reached on how substituents exercise their effects in the interaction of the aromatic molecule with the metal ion. The π-polarization (the Hunter model) and the direct local effect (the Wheeler-Houk model) are two lines of thought applied to this problem, but the justification of both approaches is based on insufficiently proven assumptions and approximations. In order to shed more light on this issue we propose a new approach which enables us to gauge directly the energetic trends resulting from the interaction of the ring with the cation. In our method we add one more partitioning level to the interacting quantum atoms (IQA) scheme and decompose the IQA interaction energies into contributions resulting from σ and π electron densities of the aromatic ring. The new approach, which is named partitioned-IQA, abbreviated as p-IQA, has been applied to complexes of derivatives of benzene or azaborine interacting with a sodium cation. The p-IQA approach reveals that in these systems both σ and π electronic moieties are polarized. Interestingly, for the majority of cases the σ-polarization outweighs the π one, contrary to the Hunter model. However, the Wheeler-Houk model is not precise, either, since the σ-polarization shows some degree of non-locality. In addition, the substituents are found to have a negligible influence on the ring orbital-overlapping capability, i.e. the covalency. Therefore, the substituent effect in the cationπ interaction is a nonlocal classical effect, indicating that neither Hunter model nor Wheeler-Houk model is able to fully describe all the aspects of the substituent effects. The p-IQA conclusions for the considered systems have been compared with the results from the functional-group SAPT (F-SAPT) method. We believe that the presented partitioning in the IQA framework will provide a deeper insight into the substituent effects in the cationπ interactions, which is beyond the σ-π atomic charge population separation.

Computational Approach To Reveal the Structural Stability and Electronic Properties of Lithiated M/CNT (M = Si, Ge) Nanocomposites as Anodes for Lithium-Ion Batteries

This work is motivated to explore the structural stability and electronic and electrochemical properties of nanocomposites of M₄Li _n (M = Si and Ge)-carbon nanotube (CNT) by employing first-principles density functional theory calculations. By analyzing the structural stability of various M₄Li _n (n = 0-10) clusters, it is revealed that a tetrahedron-shaped M₄Li₄ Zintl cluster is found to be highly stable. Our study on the interaction between the lithiated clusters and CNT illustrates that the charge transfer from the former to latter plays a pivotal role in stabilizing these nanocomposites. The structural stability of those nanocomposites arises as a consequence of bonding between lithiated clusters and CNT, which is mediated through the cation-π interaction. The strength of the interaction between them is well reflected in electronic structure calculations by shifting the energy levels with respect to the Fermi energy. Further, the electrochemical properties of these nanocomposites are explored by forming an assembly of the cluster-inserted CNT. The calculated average intercalation voltage of the systems is found to be low (maximum ∼1.0 V for M = Si and 1.05 V for M = Ge), which demonstrates their anodic behavior.

On the formation of sandwich and multidecker complexes via π⋯π interaction: a DFT study

Sandwich and multidecker complexes via organic π–inorganic π interaction.

Molecule Channels Directed by Cation-Decorated Graphene Oxide Nanosheets and Their Application as Membrane Reactors

Highly selective macromembranes, fabricated by cation-decorated graphene oxide, exhibit an excellent selectivity toward a wide range of solvents. Mixed solvents are successfully separated, based on which a membrane reactor is designed to promote a series of chemical reactions. The cations bonding to the graphene oxide nanosheets are found to be responsible for this selectivity by cation-π, electrostatic interactions, and hydrogen bonding.

B-Hb⋯π interactions in benzene–borazine sandwich and multidecker complexes: a DFT study

Multidecker complexes of benzene/borazine can be formed through B-H_b⋯π interactions.

Interlocked benzenes in triangular π-architectures: anchoring groups dictate ion binding and transmission

Graphynes and graphenylenes – potential materials for lithium ion batteries.

Sandwiches of N-doped diamondoids and benzenevialone pair–cation and cation–pi interaction: a DFT study

Cation–lone pair and cation–pi interactions in the complexes of N-doped dimondoids.

Association of phenyldiboronic acids with hydrogen bond acceptors to form hydrogen bonded DD·AA-type complexes: a DFT study

DD·AA-type complexes of phenyldiboronic acids.

A stress tensor and QTAIM perspective on the substituent effects of biphenyl subjected to torsion

The Quantum Theory of Atoms in Molecules (QTAIM) defines quantities in 3D space that can be easily obtained from routine quantum chemical calculations. The present investigation shows that local properties can be related quantitatively to measures traditionally connected to experimental data, such as Hammett constants. We consider the specific case of substituted biphenyl to quantify the effects of a torsion φ, 0.0° ≤ φ ≤ 180.0°, of the C-C bond linking the two phenyl rings for C12 H9 -x, where x = N(CH3 )2 , NH2 , CH3 , CHO, CN, NO2, on the entire molecule. QTAIM interpreted Hammett constants, aΔH(rb ) are introduced and constructed using the difference between the H(rb ) value of C12 H9 -x and the C12 H9 -H, biphenyl which is the reference molecule, with a constant of proportionality a. This investigation unexpectedly yields very good or good agreement for the x groups with the Hammett para-, meta-, and ortho-substituent constants and is checked against para-substituted benzene. We then proceed to present the interpreted substituent constants of seven new biphenyl substituent groups, where tabulated Hammett substituent constant values are not available; y = SiH3 , ZnCl, COOCH3 , SO2 NH2 , SO2 OH, COCl, CB3 . Consistency is found for the QTAIM interpreted biphenyl substituent constants of the seven new groups y independently using the stress tensor polarizability Pσ . In addition, a selection of future applications is discussed that highlight the usefulness of this approach. © 2016 Wiley Periodicals, Inc.

Chiral Nanostructures from Helical Copolymer-Metal Complexes: Tunable Cation-π Interactions and Sergeants and Soldiers Effect

Poly(phenylacetylene) (PPA) copolymers containing (R)- or (S)-MPA as minor chiral pendant can be forced to selectively adopt the right- o left-handed helix, in the presence of small amounts of Na(+) or Ag(+) ("Sergeants and Soldiers Effect") by addition of a donor cosolvent. The helical sense depends exclusively on the chiral monomer/donor cosolvent ratio, and this allows a perfect on/off tuning of the helicity of the copolymer. When the amount of the donor cosolvent is low, the metal ion complex is stabilized by a cation-π interaction, which is selectively cleaved when the amount of cosolvent is higher. Macroscopically chiral nanospheres and nanotubes composed by helical copolymers with P or M helical sense are also described. Our results demonstrate that it is possible to obtain the two enantiomeric helical structures (P and M helicities) and the corresponding nanospheres and nanotubes from a single helical copolymer, by controlled activation/deactivation of the Sergeant and Soldiers Effect with a donor cosolvent.

Behaviour of cation–pi interaction in presence of external electric field

External electric field effects cation–π interaction.

Cation–π interaction in cofacial molecular dyads: a DFT and TDDFT study

Cation–π interaction in molecular tweezer.

Effect of the cation size on the framework structures of magnesium tungstate, A4Mg(WO4)3 (A = Na, K), R2Mg2(WO4)3 (R = Rb, Cs)

A series of alkali metal magnesium tungstates, A₄Mg(WO₄)₃ (A = Na, K), R₂Mg₂(WO₄)₃ (R = Rb, Cs), were synthesized from a high temperature solution, and their structures were determined by single-crystal X-ray diffraction.

The leading role of cation–π interactions in polymer chemistry: the control of the helical sense in solution

Cation–π interactions determine the helical sense adopted by a polyphenylacetylene bearing (R)-α-methoxy-α-phenylacetamide as a pendant group (poly-1).

Observation of novel oxygen⋯oxygen interaction in supramolecular assembly of cobalt(iii) Schiff base complexes: a combined experimental and computational study

Structural features of two newly synthesized mononuclear cobalt(iii) complexes have been examined by DFT calculations.

Substituent and Solvent Effects on the Absorption Spectra of Cation-π Complexes of Benzene and Borazine: A Theoretical Study

Time-dependent density functional theory (TDDFT) has been used to predict the absorption spectra of cation-π complexes of benzene and borazine. Both polarized continuum model (PCM) and discrete solvation model (DSM) and a combined effect of PCM and DSM on the absorption spectra have been elucidated. With decrease in size of the cation, the π → π* transitions of benzene and borazine are found to undergo blue and red shift, respectively. A number of different substituents (both electron-withdrawing and electron-donating) and a range of solvents (nonpolar to polar) have been considered to understand the effect of substituent and solvents on the absorption spectra of the cation-π complexes of benzene and borazine. Red shift in the absorption spectra of benzene cation-π complexes are observed with both electron-donating groups (EDGs) and electron-withdrawing groups (EWGs). The same trend has not been observed in the case of substituted borazine cation-π complexes. The wavelength of the electronic transitions corresponding to cation-π complexes correlates well with the Hammet constants (σ_p and σ_m). This correlation indicates that the shifting of spectral lines of the cation-π complexes on substitution is due to both resonance and inductive effect. On incorporation of solvent phases, significant red or blue shifting in the absorption spectra of the complexes has been observed. Kamlet-Taft multiparametric equation has been used to explain the effect of solvent on the absorption spectra of complexes. Polarity and polarizability are observed to play an important role in the solvatochromism of the cation-π complexes.