First-principles vdW-DF investigation on the interaction between the oxazepam molecule and C60 fullerene

The interaction between oxazepam and C₆₀ fullerene was explored using first-principles vdW-DF calculations. It was found that oxazepam binds weakly to the fullerene cage via its carbonyl group. The binding of oxazepam to C₆₀ is affected drastically by nonlocal dispersion interactions, while vdW forces affect the corresponding geometries only a little. Furthermore, aqueous solution affects the geometries of the oxazepam approaching to fullerene slightly, while oxazepam binds slightly farther away from the nanocage. The results presented provide evidence for the applicability of the vdW-DF method and serve as a practical benchmark for the investigation of host–guest interactions in biological systems.

Simple benzene derivatives adsorption on defective single-walled carbon nanotubes: a first-principles van der Waals density functional study

We have investigated the interaction between open-ended zig-zag single-walled carbon nanotube (SWCNT) and a few benzene derivatives using the first-principles van der Waals density functional (vdW-DF) method, involving full geometry optimization. Such sp ²-like materials are typically investigated using conventional DFT methods, which significantly underestimate non-local dispersion forces (vdW interactions), therefore affecting interactions between respected molecules. Here, we considered the vdW forces for the interacting molecules that originate from the interacting π electrons of the two systems. The ?0.54 eV adsorption energy reveals that the interaction of benzene with the side wall of the SWCNT is typical of the strong physisorption and comparable with the experimental value for benzene adsorption onto the graphene sheet. It was found that aromatics are physisorbed on the sidewall of perfect SWCNTs, as well as at the edge site of the defective nanotube. Analysis of the electronic structures shows that no orbital hybridization between aromatics and nanotubes occurs in the adsorption process. The results are relevant in order to identify the potential applications of noncovalent functionalized systems.

Structure and spectral characteristics of diquat-cucurbituril complexes from density functional theory

Electronic structure, ¹H NMR and infrared spectra of diquat (6,7-dihydrodipyrido[1,2-b:1′,2′-e] pyrazine-5,8-diium or DQ²⁺) encapsulated by cucurbit[n]uril (n?=?7,8) hosts are obtained using the density functional theory. Theoretical calculations have shown that both CB[7] or CB[8] host possesses strong affinity toward DQ²⁺ compared to its reduced cation or neutral species. Calculated ¹H NMR spectra reveal that H_α protons on bi-pyridinium rings of DQ²⁺@CB[8] complex are de-shielded owing to C=O?H interactions. On the other hand aromatic (H_β and H_δ) of DQ²⁺ within the CB[8] cavity exhibit significant shielding. The complexation of CB[8] with DQ²⁺ splits the carbonyl stretching vibration (1788 cm^?1) into two distinct vibrations which correspond to 1765 cm^?1 arising from hydrogen bonded carbonyls and the 1792 cm^?1 band from non-interacting ones. Further, the CN stretching vibration in DQ²⁺ exhibits a frequency blue-shift of 6 cm^?1 on its encapsulation within the CB[8] cavity. The direction of frequency shift has been explained on the basis of natural bond orbital analyses.

Quantitative characterization of protein tertiary motifs

Dimerization of 2-naphthalenecarbonitrile (2-NpCN) mediated by cucurbit[8]uril (CB[8]) has been investigated employing the density functional theory. Different structures of 2-NpCN dimers were generated by combining monomers in anti-head-to-head (A), anti-head-to-tail (B) and syn-head-to-tail (C) fashion. All these dimeric structures possess rigid cube-like architecture. On confinement within the CB[8] dimer A turns out to be the lowest energy structure. Calculated ¹H NMR spectra revealed that the 2-NpCN dimer exhibits large shielding for aromatic protons consistent with the experiment. The protons attached to cubane moiety on the other hand, led to down-field signals. Dimerization mediated with CB[8] cavitand is further accompanied by the frequency up-shift (blue shift) of methylene stretching vibration in its infrared spectra.

A comparative DFT study on aquation and nucleobase binding of ruthenium (II) and osmium (II) arene complexes

The potential energy surfaces of the reactions of organometallic arene complexes of the type [(η ⁶-arene)M^II(pic)Cl] (where pic = 2-picolinic acid, M = Ru or Os) were examined by a DFT computational study. Among the seven density functional methods, hybrid exchange functional B3LYP outperforms the others to explain the aquation of the complexes. The reactions and binding energies of Ru^II and Os^II arene complexes with both 9EtG and 9EtA were studied to gain insight into the reactivity of these types of organometallic complexes with DNA. The obtained data rationalize experimental observation, contributing to partly understanding the potential biological and medical applications of organometallic complexes.

Structural elucidation of supramolecular alpha-cyclodextrin dimer/aliphatic monofunctional molecules complexes

The structural elucidation of 2α-cyclodextrin/1-octanethiol, 2α-cyclodextrin/1-octylamine and 2α-cyclodextrin/1-nonanoic acid inclusion complexes by nuclear magnetic resonance (NMR) spectroscopy and molecular modeling has been achieved. The detailed spatial configurations are proposed for the three inclusion complexes based on 2D NMR method. ROESY experiments confirm the inclusion of guest molecules inside the α-cyclodextrin (α-CD) cavity. On the other hand, the host-guest ratio observed was 2:1 for three complexes. The detailed spatial configuration proposed based on 2D NMR methods were further interpreted using molecular modeling studies. The theoretical calculations are in good agreement with the experimental data.

Theoretical investigation on switchable second-order nonlinear optical (NLO) properties of novel cyclopentadienylcobalt linear [4]phenylene complexes

As a kind of novel organometallic complexes, the cyclopentadienylcobalt (CpCo) linear [4]phenylene complexes (4 = number of benzene rings) display efficient switchable nonlinear optical (NLO) response when CpCo reversibly migrates along the linear [4]phenylene triggered by heating or lighting. In this paper, the second-order NLO properties for CpCo linear [4]phenylene complexes were calculated by using the density functional theory (DFT) methods with four functionals. All of the functionals yield the same order of β _tot values: 1<2<4<3. The effect of solvent on second-order NLO properties has been studied using polarized continuum model (PCM) in the tetrahydrofuran (THF) solution. The solvent leads to a slight enhancement of the NLO responses for the studied complexes relevant to their NLO responses in vacuo. The electronic absorption spectra were investigated by the TDDFT methods. The TDDFT calculations indicate that the maximum absorption peaks of complexes 2–4 in the near-infrared spectrum area show the bathochromic shift together with a decreasing intensity compared to complex 1. We have also found that the cobalt (Co) atom acts as a donor in all the organometallic complexes and the d → π* and π → π* charge transfer (CT) transitions contribute to the enhancement of second-order NLO response. Furthermore, two experimentally existing complexes 1 and 3 are found to have a large difference in β _tot values. It is our expectation that this difference may stimulate the search for a new type of switchable NLO material based on CpCo linear [4]phenylene complexes.

Modeling the effect of H-bonding interactions and molecular packing on the molecular structure of [Ag(ethylnicotinate)2]NO3 complex

The gas phase molecular structure of a single isolated molecule of [Ag(Etnic)₂NO₃];1 where Etnic = Ethylnicotinate was calculated using B3LYP method. The H-bonding interaction between 1 with one (complex 2) and two (complex 3) water molecules together with the dimeric formula [Ag(Etnic)₂NO₃]₂;4 and the tetrameric formula [Ag(Etnic)₂NO₃]₄;5 were calculated using the same level of theory to model the effect of intermolecular interactions and molecular packing on the molecular structure of the titled complex. The H-bond dissociation energies of complexes 2 and 3 were calculated to be in the range of 12.220–14.253 and 30.106–31.055 kcal?mol^?1, respectively, indicating the formation of relatively strong H-bonds between 1 and water molecules. The calculations predict bidentate nitrate ligand in the case of 1 and 2, leading to distorted tetrahedral geometry around the silver ion with longer Ag–O distances in case of 2 compared to 1, while 3 has a unidentate nitrate ligand leading to a distorted trigonal planar geometry. The packing of two [Ag(Etnic)₂NO₃] complex units; 4 does not affect the molecular geometry around Ag(I) ion compared to 1. In the case of 5, the two asymmetric units of the formula [Ag(Etnic)₂NO₃] differ in the bonding mode of the nitrate group, where the geometry around the silver ion is distorted tetrahedral in one unit and trigonal planar in the other. The calculations predicted almost no change in the charge densities at the different atomic sites except at the sites involved in the C–H?O interactions as well as at the coordinated nitrogen of the pyridine ring.

A molecular dynamics study on sI hydrogen hydrate

A molecular dynamics simulation is carried out to explore the possibility of using sI clathrate hydrate as hydrogen storage material. Metastable hydrogen hydrate structures are generated using the LAMMPS software. Different binding energies and radial distribution functions provide important insights into the behavior of the various types of hydrogen and oxygen atoms present in the system. Clathrate hydrate cages become more stable in the presence of guest molecules like hydrogen.

Stereoselectivity of chalcone isomerase with chalcone derivatives: a computational study

Chalcone isomerase (CHI) catalyzes the intramolecular cyclization of chalcones into flavonoids. The activity of CHI is essential for the biosynthesis of flavonoids precursors of floral pigments and phenylpropanoid plant defense compounds. In the present study, we explored the detailed binding structures and binding free energies for two different active site conformations of CHI with s-cis/s-trans conformers of three chalcone compounds by performing molecular dynamics (MD) simulations and binding free energy calculations. The computational results indicate that s-cis/s-trans conformers of chalcone compounds are orientated in the similar binding position in the active site of CHI and stabilized by the different first hydrogen bond network and the same second hydrogen bond network. The first hydrogen bond network results in much lower binding affinity of s-trans conformer of chalcone compound with CHI than that of s-cis conformer. The conformational change of the active site residue T48 from indirectly interacting with the substrate via the second hydrogen bond network to directly forming the hydrogen bond with the substrates cannot affect the binding mode of both conformers of chalcone compounds, but remarkably improves the binding affinity. These results show that CHI has a strong stereoselectivity. The calculated binding free energies for three chalcone compounds with CHI are consistent with the experimental activity data. In addition, several valuable insights are suggested for future rational design and discovery of high-efficiency mutants of CHI.

Nucleus-independent chemical shift criterion for aromaticity in π-extended tetraoxa[8]circulenes

Recently synthesized π-extended symmetrical tetraoxa[8]circulenes that exhibit electroluminescent properties were calculated at the density functional theory (DFT) level using the quantum theory of atoms in molecules (QTAIM) approach to electron density distribution analysis. Nucleus-independent chemical shift (NICS) indices were used to characterize the aromaticity of the studied molecules. The tetraoxa[8]circulene molecules were found to consist of two antiaromatic perimeters (according to the Hückel “4n” antiaromaticity rule) that include 8 and 24 π-electrons. Conversely, NICS calculations demonstrated the existence of a common π-extended system (distributed like a flat ribbon) in the studied tetraoxa[8]circulene molecules. Thus, these symmetrical tetraoxa[8]circulene molecules provide examples of diatropic systems characterized by the presence of induced diatropic ring currents.

Theoretical investigation on the structure and thermodynamic properties of the 2,4-dinitroimidazole complex with methanol

The structure and thermodynamic properties of the 2, 4-dinitroimidazole complex with methanol were investigated using the B3LYP and MP2(full) methods with the 6-31++G(2d,p) and 6-311++G(3df,2p) basis sets. Four types of hydrogen bonds [N–H?O, C–H?O, O–H?O (nitro oxygen) and O–H?π] were found. The hydrogen-bonded complex having the highest binding energy had a N–H?O hydrogen bond. Analyses of natural bond orbital (NBO) and atoms-in-molecules (AIM) revealed the nature of the intermolecular hydrogen-binding interaction. The changes in thermodynamic properties from monomers to complexes with temperatures ranging from 200.0 to 800.0 K were investigated using the statistical thermodynamic method. Hydrogen-bonded complexes of 2,4-dinitroimidazole with methanol are fostered by low temperatures.

Single crystal architecture and absorption spectra of octathio[8]circulene and sym-tetraselenatetrathio[8]circulene: QTAIM and TD-DFT approach

The single crystal architecture of the high-symmetry octathio[8]circulene and sym-tetraselenatetrathio[8]circulene is studied at the density functional theory (DFT) level with the quantum theory of atoms in molecules (QTAIMs) approach to the electron density distribution analysis. The presence of stabilizing intermolecular C---C, C---S and C---Se contacts in the longitudinal and transversal projections of the single crystals is postulated on the grounds of the previous high-resolution X-ray data for octathio[8]circulene; it is supported by the present QTAIM calculations and also predicted in some new details for both circulenes. We suggest that the appearance of the observed red color for the monocrystalline octathio[8]circulene is caused by strong intermolecular interactions between the molecules in the single crystal. However, the intermolecular interactions for the sym-tetraselenatetrathio[8]circulene crystal fragment are weaker and molecular layers are more friable in comparison to octathio[8]circulene crystal structure. These lead to the absence of visible absorption for the sym-tetraselenatetrathio[8]circulene crystal.

Microsolvation of Mg2+, Ca2+: strong influence of formal charges in hydrogen bond networks

A stochastic exploration of the quantum conformational spaces in the microsolvation of divalent cations with explicit consideration of up to six solvent molecules [Mg (H ₂ O) _n )]²⁺, (n?=?3, 4, 5, 6) at the B3LYP, MP2, CCSD(T) levels is presented. We find several cases in which the formal charge in Mg²⁺ causes dissociation of water molecules in the first solvation shell, leaving a hydroxide ion available to interact with the central cation, the released proton being transferred to outer solvation shells in a Grotthus type mechanism; this particular finding sheds light on the capacity of Mg²⁺ to promote formation of hydroxide anions, a process necessary to regulate proton transfer in enzymes with exonuclease activity. Two distinct types of hydrogen bonds, scattered over a wide range of distances (1.35–2.15 Å) were identified. We find that in inner solvation shells, where hydrogen bond networks are severely disturbed, most of the interaction energies come from electrostatic and polarization+charge transfer, while in outer solvation shells the situation approximates that of pure water clusters.

Relationship between structure and entropy contributions in an anthraquinone mercapto derivative

The structural and thermodynamic properties of an anthraquinone derivative were studied by means of quantum-chemical calculations. Conformational analysis using ab initio and density functional theory methods revealed 14 low-energy conformers. In order to discuss similarities and differences in entropy of the conformers, the rotational and vibrational contributions to entropy were correlated with changes in conformer structure. The component of the moment of inertia perpendicular to the molecular plane gives significant input to ΔS _rot , whereas the largest contributions to the ΔS _vib have vibrations associated with the τ _S1C20 coordinate.

Theoretical study of photophysical properties of 1,4-dihydropyrrolo[3,2-b]pyrrole-cored branched molecules with thienylenevinylene arms toward broad absorption spectra for solar cells

A series of oligo(thienylenevinylene) derivatives with 1,4-dihydropyrrolo[3,2-b]pyrrole as core has been investigated at the PBE0/6-31G(d) and the TD-PBE0/6-31+G(d,p) levels to design materials with high performances such as broad absorption spectra and higher balance transfer property. The results show that position and amount of arm affect the electronic density contours of frontier molecular orbitals significantly. The molecule with four arms owns the narrowest energy gap and the largest maximum absorption wavelength, and the molecule with two arms in positions a and c has the broadest absorption region among the designed molecules. Calculated reorganization energies of the designed molecules indicate that the molecules with two arms can be good potential ambipolar transport materials under proper operating conditions.

DFT studies of the conversion of four mesylate esters during reaction with ammonia

The energetics of the Menshutkin-like reaction between four mesylate derivatives and ammonia have been computed using B3LYP functional with the 6-31+G** basis set. Additionally, MPW1K/6-31+G** level calculations were carried out to estimate activation barrier heights in the gas phase. Solvent effect corrections were computed using PCM/B3LYP/6-31+G** level. The conversion of the reactant complexes into ion pairs is accompanied by a strong energy decrease in the gas phase and in all solvents. The ion pairs are stabilized with two strong hydrogen bonds in the gas phase. The bifurcation at C2 causes a significant activation barrier increase. Also, bifurcation at C5 leads to noticeable barrier height differentiation. Both B3LYP/6-31+G** and MPW1K/6-31+G** activation barriers suggest the reaction 2 (2a?+?NH₃) to be the fastest in the gas phase. The reaction 4 is the slowest one in all environments.