Effect of Terminal Fluorination on Chiroptical Properties of Carbo[5-8]helicenes: A Systematic Computational Study at the RI-ADC(2) Level

In this article, effects of di-, tetra-, and octafluorination on the structural and chiroptical properties of carbo[5-8]helicenes are reported. Three fluorinated derivatives are designed from each parent carbohelicene by substituting either one, two, or four hydrogens at each terminal ring with fluorine atoms. Excited states properties such as UV-vis and CD spectra of all the six fluorinated carbohelicenes are computed at the ADC(2)/def2-TZVP level, and the results are compared against the results of their respective parent carbohelicene. In addition, CPL properties are also computed at the same level of theory. In the case of carbo[5]helicene (5H), g_CPL decreases with an increase in the degree of fluorination. A similar observation is made in carbo[6]helicene (6H) too, although the value for tetrafluorinated 6H is slightly larger than for difluorinated 6H. Di- and tetrafluorination in carbo[7]helicene (7H), and all types of considered fluorination in carbo[8]helicene (8H) produce improved g_CPL results. Results for fluorescence rate constants are also shown. Results are analyzed in terms of the transition dipole moment vectors and the angles between those.

Classical Electrostatics Remains the Driving Force for Interanion Hydrogen and Halogen Bonding

Interanion hydrogen bonding (IAHB) and halogen bonding (IAXB) have emerged as a counterintuitive linker in a range of fascinating applications. Despite the overall repulsive (positive) binding energy, anions are trapped in a local minimum with its corresponding transition state (TS) preventing dissociation. In other words, the adduct of anions is metastable. Seemingly, the electrostatic paradigm and force field description of hydrogen/halogen bonding (HB/XB) are challenged, because of the preconceived Coulombic repulsion. Aiming at an insightful understanding of these interanion phenomena, we employed the energy decomposition approach based on the block-localized wavefunction method (BLW-ED) to investigate a series of exemplary interanion complexes. As expected, the key distinction from the conventional HB/XB lies in the electrostatic interaction, which is not increasingly repulsive as anions gradually approach to each other. Rather, there is a Coulombic barrier at a certain point. After this point, the electrostatic repulsion diminishes with the decreasing distance between anions. Differently, other energy components vary monotonically just like in conventional cases. The nonmonotonic characteristic of the electrostatic interaction in interanion complexes was reproduced using the multipole expansion in AMOEBA polarizable force field in which the state-specified atomic multipoles were adopted. This suggests that the nonmonotonicity can be well interpreted by classical electrostatic theory and there is no conceptual difference between conventional HB/XB and IAHB/IAXB. The stability of IAHB/IAXB depends on the competition between the local attractive HB/XB and the global Coulombic repulsion of net charges, though there is cooperativity between these two contrasting forces. This concise model was supported by the attractive IAHB/IAXB in modified molecular capsules, which exhibit strong quadruple HB/XBs and a considerable distance between charged substituents.

Resonance-assisted/impaired anion-π interaction: towards the design of novel anion receptors

Substituents alter the electron density distribution in benzene in various ways, depending on their electron withdrawing and donating capabilities, as summarized by the empirical Hammett equation. The change of the π electron density distribution subsequently impacts the interaction of substituted benzenes or other cyclic conjugated rings with anions. Currently the design and synthesis of conjugated cyclic receptors capable of binding anions is an active field due to their applications in the sensing and removal of environmental contaminants and molecular recognition. By using the block-localized wavefunction (BLW) method, which is a variant of ab initio valence bond (VB) theory and can derive the reference resonance-free state self-consistently, we quantified the resonance-assisted (RA) or resonance-impaired (RI) phenomena in anion–π interactions from both structural and energetic perspectives. The frozen interaction, in which the electrostatic attraction is involved, has been shown to be the governing factor for the RA or RI interactions with anions. Energy analyses based on the empirical point charge (EPC) model indicated that the anion–π interactions can be simplified as the attraction between a negative point charge (anion) and a group of local dipoles, affected by the enriched or diminished π-cloud due to the resonance between the substituents and the conjugated ring. Hence, two strategies for the design of novel anion receptors can be envisioned. One is the enhancement of the magnitudes and/or numbers of local dipoles (polarized σ bonds), and the other is the reduction of π electron density in conjugated rings. For cases with the RI characteristics, “curved” aromatic molecules are preferred to be anion receptors. Indeed, extremely strong binding was found in complexes formed with fluorinated corannulene (F-CDD) and fluorinated [5]cycloparaphenylene (F-[5]CPP). Inspired by the RA phenomenon, complexes of p-, o- and m-benzoquinones with halides were revisited.

Substituents alter the electron density distribution in benzene in various ways, depending on their electron withdrawing and donating capabilities, as summarized by the empirical Hammett equation.

Origin of the Diastereoselectivity of the Heterogeneous Hydrogenation of a Substituted Indolizine

In this work, the stereoselective heterogeneous hydrogenation of a tetrasubstituted indolizine was studied. Partial hydrogenation products were obtained in three steps from a substituted pyridine-2-carboxaldehyde prepared from commercial pyridoxine hydrochloride. The hydrogenation of the indolizine ring was shown to be diastereoselective, forming trans-6b and cis-9. Theoretical calculations (ab initio and DFT) were used to rationalize the unusual trans stereoselectivity for 6b, and a keto–enol tautomerism under kinetic control has been proposed as the source of diastereoselectivity.

DL-3-Aminoisobutyric acid: vibrational, NBO and AIM analysis of N-H⋯O bonded-zwitterionic dimer model

A zwitterionic dimer model constructed of inter-molecular -N-H⋯O bonding has been proposed for the solid sample of DL-3-Aminoisobutyric acid consistent with IR absorption and Raman spectral features measured in the 3500-400/50 cm^-1. This zwitterionic dimer model in water as solvent has been computed at B3LYP/6-311++G(d,p) and B3LYP-D3/6-311++G(d,p) levels including Grimme's dispersion correction associated with the -N-H⋯O interaction and SCRF-SMD method. Of the several possible monomer and dimer conformational structures, the most stable dimer constructed of two zwitterion monomer units has produced vibrational modes due to the -NH₃ ⁺ cation and -CO₂‾ anion involved in the -N-H⋯O bonding in fair agreement with the observed broad but composite IR modal features near the 3500-2000 cm^-1. Except for the frequency of asymmetric stretching mode of the -NH₃ ⁺ cation, its symmetric and bending modes agree with the observed values. As for the -CO₂‾ anion, the frequencies of all of its modes are in good agreement with the experiment. Natural bond orbital (NBO), molecular electrostatic potential (MEP), atoms-in-molecules (AIM) and non-covalent interaction (NCI) analyses have been used to understand electronic characterization of the -N-H⋯O bonding.

Interplay among Electrostatic, Dispersion, and Steric Interactions: Spectroscopy and Quantum Chemical Calculations of π-Hydrogen Bonded Complexes

π-Hydrogen bonding interactions are ubiquitous in both materials and biology. Despite their relatively weak nature, great progress has been made in their investigation by experimental and theoretical methods, but this becomes significantly more complicated when secondary intermolecular interactions are present. In this study, the effect of successive methyl substitution on the supramolecular structure and interaction energy of indole⋅⋅⋅methylated benzene (ind⋅⋅⋅n-mb, n=1-6) complexes is probed through a combination of supersonic jet experiments and benchmark-quality quantum chemical calculations. It is demonstrated that additional secondary interactions introduce a subtle interplay among electrostatic and dispersion forces, as well as steric repulsion, which fine-tunes the overall structural motif. Resonant two-photon ionization and IR-UV double-resonance spectroscopy techniques are used to probe jet-cooled ind⋅⋅⋅n-mb (n=2, 3, 6) complexes, with redshifting of the N-H IR stretching frequency showing that increasing the degree of methyl substitution increases the strength of the primary N-H⋅⋅⋅π interaction. Ab initio harmonic frequency and binding energy calculations confirm this trend for all six complexes. Electronic spectra of the three dimers are broad and structureless, with quantum chemical calculations revealing that this is likely to be due to multiple tilted conformations of each dimer possessing similar stabilization energies.

Stabilizing and Modulating Color by Copigmentation: Insights from Theory and Experiment

Natural anthocyanin pigments/dyes and phenolic copigments/co-dyes form noncovalent complexes, which stabilize and modulate (in particular blue, violet, and red) colors in flowers, berries, and food products derived from them (including wines, jams, purees, and syrups). This noncovalent association and their electronic and optical implications constitute the copigmentation phenomenon. Over the past decade, experimental and theoretical studies have enabled a molecular understanding of copigmentation. This review revisits this phenomenon to provide a comprehensive description of the nature of binding (the dispersion and electrostatic components of π-π stacking, the hydrophobic effect, and possible hydrogen-bonding between pigment and copigment) and of spectral modifications occurring in copigmentation complexes, in which charge transfer plays an important role. Particular attention is paid to applications of copigmentation in food chemistry.

Direct Spectroscopic Evidence for an n→π* Interaction

The n→π* interaction is an extremely weak but very important noncovalent interaction. Although this interaction is widely present in biomolecules and materials, its existence is counterintuitive and so has been debated extensively. Herein, we have reported direct spectroscopic evidence for an n→π* interaction for the first time by probing the carbonyl stretching frequency in phenyl formate using isolated gas-phase IR spectroscopy. This result also demonstrates that the conformational preference for the cis conformer of phenyl formate compared to the trans conformer arises due to the presence of the n→π* interaction in the former. The direct proof reported herein for this controversial but important noncovalent interaction should stimulate further experimental and theoretical investigation on this intriguing research topic.

Soft hydrogen bonds to alkenes: the methanol-ethene prototype under experimental and theoretical scrutiny

An FTIR spectroscopic study of the elusive hydrogen-bonded methanol-ethene complex, the most elementary example for weak intermolecular alcohol hydrogen bonding to a π cloud, is presented. By isolating the complex in a supersonic jet, the rigorous comparability to high-level quantum chemical calculations is ensured. In stark contrast to classical hydrogen bonds, experimental overtone analysis reveals the harmonic oscillator approximation for the OH red shift to be accurate. Harmonic calculations up to explicitly correlated local coupled-cluster level are thus found to agree very well with experiment. The experimental OH values for the red shift (45 cm^-1), the small change in diagonal anharmonicity (-3 cm^-1) and the overtone intensity attenuation (2 × 10²-fold) together with theoretical predictions for the preferred structural arrangement and the zero-point-corrected dissociation energy (8 kJ mol^-1) may thus be regarded as definitive reference values for related systems and for more approximate computational methods. In particular, MP2 calculations are shown to fail for this kind of weak intermolecular interaction.

Role of size and shape selectivity in interaction between gold nanoclusters and imidazole: a theoretical study

We present a theoretical study on the structure, stability, spectra and electronic properties of imidazole (Im) adsorbed on gold nanoclusters (Aun, n = 2, 4, 6, 8, 10, and 20). These computations were performed using various density functional theories with and without inclusion of Grimme's (D3) dispersion correction. For small clusters, we also carried out wavefunction-based ab initio (MP2 and SCS-MP2) computations for comparison. Vibrational, atoms in molecules (AIM) and natural bond orbital (NBO) analyses clearly reveal the occurrence of charge transfer (CT) through covalent (N1-Au) and noncovalent interactions that play important roles in the stability of the Im@Aun complexes with anchor assisted H-bonds (Cα-H · Au). Therefore, gold clusters can act as H-bond acceptors with biomolecules for development of new materials and applications. Our study establishes also the ability and reliability of PBE0 and M05-2X functionals compared to B3LYP and PBE for an accurate description of covalent and noncovalent interactions between Im and gold clusters since they lead to close agreement with MP2. Finally, we show that the Au8 cluster may be viewed as large enough to mimic the 3D gold surface.

Structure and intermolecular vibrations of perylene·trans-1,2-dichloroethene, a weak charge-transfer complex

The vibronic spectra of strong charge-transfer complexes are often congested or diffuse and therefore difficult to analyze. We present the spectra of the π-stacked complex perylene trans-1,2-dichloroethene, which is in the limit of weak charge transfer, the electronic excitation remaining largely confined to the perylene moiety. The complex is formed in a supersonic jet, and its S0 ↔ S1 spectra are investigated by two-color resonant two-photon ionization (2C-R2PI) and fluorescence spectroscopies. Under optimized conditions, vibrationally cold (T(vib) ≈ 9 K) and well resolved spectra are obtained. These are dominated by vibrational progressions in the “hindered-rotation” Rc intermolecular vibration with very low frequencies of 11 (S0) and 13 cm(–1) (S1). The intermolecular Tz stretch and the Ra and Rb bend vibrations are also observed. The normally symmetry-forbidden intramolecular 1a(u) “twisting” vibration of perylene also appears, showing that the π- stacking interaction deforms the perylene moiety, lowering its local symmetry from D2h to D2. We calculate the structure and vibrations of this complex using six different density functional theory (DFT) methods (CAM-B3LYP, BH&HLYP, B97-D3, ωB97X-D, M06, and M06-2X) and compare the results to those calculated by correlated wave function methods (SCS-MP2 and SCS-CC2). The structures and vibrational frequencies predicted with the CAM-B3LYP and BH&HLYP methods disagree with the other calculations and with experiment. The other four DFT and the ab initio methods all predict a π-stacked “centered” structure with nearly coplanar perylene and dichloroethene moieties and intermolecular binding energies of D(e) = −20.8 to −26.1 kJ/mol. The 000 band of the S0 → S1 transition is red-shifted by δν = −301 cm(–1) relative to that of perylene, implying that the D(e) increases by 3.6 kJ/mol or 15% upon electronic excitation. The intermolecular vibrational frequencies are assigned to the calculated Rc, Tz, Ra, and Rb vibrations by comparing to the observed/calculated frequencies and S0 ↔ S1 Franck–Condon factors. Of the three TD-DFT methods tested, the hybrid-meta-GGA functional M06-2X shows the best agreement with the experimental electronic transition energies, spectral shifts, and vibronic spectra, closely followed by the ωB97X-D functional, while the M06 functional gives inferior results.

Ring opening of 2-aza-3-borabicyclo[2.2.0]hex-5-ene, the Dewar form of 1,2-dihydro-1,2-azaborine: stepwise versus concerted mechanisms

The ring opening of the Dewar form of 1,2-dihydro-1,2-azaborine, 2-aza-3-borabicyclo[2.2.0]hex-5-ene (3) is investigated by theoretical methods by using multiconfiguration SCF (CASSCF) and coupled cluster theory [CCSD(T)] with basis sets up to polarised quadruple-zeta quality. The title compound was previously reported to form photochemically in cryogenic noble gas matrices from 1,2-dihydro-1,2-azaborine (4). Four reaction paths for the thermal ring opening of 3 to 4 could be identified. These are the conventional disrotatory and conrotatory electrocyclic ring-opening pathways where the BN unit is only a bystander. Two more favourable paths are stepwise and involve 1,3-boron-carbon interactions. The lowest energy barrier for the isomerisation reaction, 22 kcal mol(-1), should be high enough for an experimental observation in solution. However, in solution the dimerisation of 3 is computed to have a very low barrier (3 kcal mol(-1)), and thus 3 is expected to be a short-lived reactive intermediate.

Insight on the formation of chitosan nanoparticles through ionotropic gelation with tripolyphosphate

This work reports details pertaining to the formation of chitosan nanoparticles that we prepare by the ionic gelation method. The molecular interactions of the ionic cross-linking of chitosan with tripolyphosphate have been investigated and elucidated by means of all-electron density functional theory. Solvent effects have been taken into account using implicit models. We have identified primary-interaction ionic cross-linking configurations that we define as H-link, T-link, and M-link, and we have quantified the corresponding interaction energies. H-links, which display high interaction energies and are also spatially broadly accessible, are the most probable cross-linking configurations. At close range, proton transfer has been identified, with maximum interaction energies ranging from 12.3 up to 68.3 kcal/mol depending on the protonation of the tripolyphosphate polyanion and the relative coordination of chitosan with tripolyphosphate. On the basis of our results for the linking types (interaction energies and torsion bias), we propose a simple mechanism for their impact on the chitosan/TPP nanoparticle formation process. We introduce the β ratio, which is derived from the commonly used α ratio but is more fundamental since it additionally takes into account structural details of the oligomers.

A reinvestigation of the dimer of para-benzoquinone and pyrimidine with MP2, CCSD(T), and DFT using functionals including those designed to describe dispersion

We reevaluate the interaction of pyridine and p-benzoquinone using functionals designed to treat dispersion. We compare the relative energies of four different structures: stacked, T-shaped (identified for the first time), and two planar H-bonded geometries using these functionals (B97-D, ωB97x-D, M05, M05-2X, M06, M06L, and M06-2X), other functionals (PBE1PBE, B3LYP, X3LYP), MP2, and CCSD(T) using basis sets as large as cc-pVTZ. The functionals designed to treat dispersion behave erratically as the predictions of the most stable structure vary considerably. MP2 predicts the experimentally observed structure (H-bonded) to be the least stable, while single-point CCSD(T) at the MP2 optimized geometry correctly predicts the observed structure to be the most stable. We have confirmed the assignment of the experimental structure using new calculations of the vibrational frequency shifts previously used to identify the structure. The MP2/cc-pVTZ vibrational calculations are in excellent agreement with the observations. All methods used to calculate the energies provide vibrational shifts that agree with the observed structure even though most do not predict this structure to be most stable. The implications for evaluating possible π-stacking in biologically important systems are discussed.

The MP2-F12 method in the Turbomole program package

A detailed description of the explicitly correlated second-order Møller-Plesset perturbation theory (MP2-F12) method, as implemented in the TURBOMOLE program package, is presented. The TURBOMOLE implementation makes use of density fitting, which greatly reduces the prefactor for integral evaluation. Methods are available for the treatment of ground states of open- and closed-shell species, using unrestricted as well as restricted (open-shell) Hartree-Fock reference determinants. Various methodological choices and approximations are discussed. The performance of the TURBOMOLE implementation is illustrated by example calculations of the molecules leflunomide, prednisone, methotrexate, ethylenedioxytetrafulvalene, and a cluster model for the adsorption of methanol on the zeolite H-ZSM-5. Various basis sets are used, including the correlation-consistent basis sets specially optimized for explicitly correlated calculations (cc-pVXZ-F12).

Estimating the hydrogen bond energy

First, different approaches to detect hydrogen bonds and to evaluate their energies are introduced newly or are extended. Supermolecular interaction energies of 256 dimers, each containing one single hydrogen bond, were correlated to various descriptors by a fit function depending both on the donor and acceptor atoms of the hydrogen bond. On the one hand, descriptors were orbital-based parameters as the two-center or three-center shared electron number, products of ionization potentials and shared electron numbers, and the natural bond orbital interaction energy. On the other hand, integral descriptors examined were the acceptor-proton distance, the hydrogen bond angle, and the IR frequency shift of the donor-proton stretching vibration. Whereas an interaction energy dependence on 1/r(3.8) was established, no correlation was found for the angle. Second, the fit functions are applied to hydrogen bonds in polypeptides, amino acid dimers, and water cluster, thus their reliability is demonstrated. Employing the fit functions to assign intramolecular hydrogen bond energies in polypeptides, several side chain CH...O and CH...N hydrogen bonds were detected on the fly. Also, the fit functions describe rather well intermolecular hydrogen bonds in amino acid dimers and the cooperativity of hydrogen bond energies in water clusters.