Theoretical investigation of solvent effects on tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol

The density functional theory has been used to study the tautomeric equilibrium of 2‐diazo‐4,6‐dinitrophenol(DDNP) in the gas phase and in 14 solvents at the B3LYP/6‐31G* level. The solvent effects on the tautomeric equilibria were investigated by the self‐consistent reaction field theory (SCRF) based on conductor polarized continuum model (CPCM) in apolar and polar solvents and by the hybrid continuum‐discrete model in protic solvent, respectively. Solvent effects on the computed molecular properties, such as molecular geometries, dipole moments, E_LUMO, E_HOMO, total energies for DDNP tautomers and transition state, tautomerization energies and solvation energies have been found to be evident. The tautomeric equilibrium of DDNP is solvent‐dependent to a certain extent. The tautomer I (cyclic azoxy form) is preferred in the gas phase, while in nonpolar solvents tautomer I and II (quinold form) exist in comparable amounts, and in highly polar solvents, the tautomeric equilibrium is shifted in favor of the more polar tautomer II . © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Three‐dimensional reference interaction site model self‐consistent field analysis of solvent and substituent effects on the absorption spectra of Brooker's merocyanine

Solvent and substituent effects on the absorption spectra of Brooker's merocyanine (BM) are investigated using the three‐dimensional reference interaction site model self‐consistent field method and time‐dependent density functional theory. The π–π* excitation energies are computed for BM and its derivative 2,6‐di‐tert‐butyl (di‐t‐Bu) BM. The behaviors of the computed excitation energies with increasing solvent polarity are in good agreement with those of the corresponding experimental measurements. In addition, analysis of the solute–solvent interaction energies and spatial distribution functions reveals that the effects of the solvent on the absorption spectra are reduced by the steric hindrance of the t‐Bu groups. Furthermore, from the difference in the solute–solvent interaction energies of BM and di‐t‐Bu BM, it is shown that the effect of the t‐Bu substituents on the absorption spectrum is greater in high‐polarity solvents. © 2015 Wiley Periodicals, Inc.     

Study of the formamide–methanol dimer with ab initio and density functional theory methods

For the first time, the structures and energies for the hydrogen bonding of a 1:1 complex formed between formamide and methanol molecules have been computed with various pure and hybrid density functional theory (DFT) and ab initio methods at varied basis set levels from 6‐31g to 6‐31+g(d,p). Five reasonable geometries on the potential energy surface of methanol and formamide system are considered and their relative stability is discussed. The infrared (IR) spectrum frequencies, IR intensities, and vibrational frequency shifts are reported. From the systematic studies, it is found that all the DFT methods selected here correctly compute the dimerization energies and geometries, with the B3P86 method predicting the hydrogen bond lengths relatively shorter and BPW91 yielding the interaction energies relatively lower. Finally, the solvent effects on the geometries of the formamide–methanol complexes have also been investigated using self‐consistent reaction field (SCRF) calculations with five different DFT methods at the 6‐31+g(d,p) basis set level. The results indicate that the polarity of the solvent has played an important role on the structures and relative stabilities of different isomers. Moreover, the basis set superposition error correction is critical to the interaction energies in the polar solvents. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

Theoretical study of the syn- and anti-conformers of 2,2ʹ-bifuran derivatives: rotational barrier and solvent effect

The optimised molecular structures, vibrational wavenumbers and corresponding vibrational assignments of the syn- and anti-conformers of 2,2?-bifuran and its nitro, fluoro, methyl and hydroxyl derivatives were obtained using density functional theory. The starting structures with C₂ symmetry of all the ground state structures were considered and the transition state arising from syn-anti isomerisation was also modelled. All structures were fully optimised, and the geometries, dipole moments, charge, thermodynamic properties, and energies are reported. The calculated vibrational wavenumbers were assigned to the various fundamental modes of vibrations. The integral Equation Formalism Polarisation Continuum Model (IEF-PCM) was used to calculate the optimised geometry and the vibrational wavenumbers for all the compounds in different solvents. The results indicate that in the gas phase, the syn-conformer is more stable while in solution phases the conformational preference depends on the polarity of the solvent.     

Tautomeric equilibria of 2‐ and 4‐thiouracil in gas phase and in solvent: A density functional study

The relative stabilities of the five favored tautomers of 2‐ and 4‐thiouracil in gas phase and in water solution were determined by density functional theory employing the Becke, Lee, Yang, and Parr (B3LYP) exchange–correlation potential and the three 6‐31G(d,p), 6‐311++G(d,p), and triple‐zeta valence (TZVP) basis sets. Zero‐point vibrational corrections were also computed. Bulk solvent effects were studied in the framework of the self‐consistent reaction field approach by the polarizable continuum model. All calculations indicate that the most stable tautomer for both species, in the gas phase as well as in solution, has the oxo‐thione form, in full agreement with the previous ab initio and experimental studies. The tautomeric stability orders obtained in the aqueous solution are sensibly different from that in the gas phase. At B3LYP/6‐311++G(d,p) level in the gas phase, the following orders of stability for 2‐ and 4‐thiouracil tautomers were observed, respectively: S2U1>S2U2>S2U4>S2U5>S2U3 and S4U1>S4U2>S4U3>S4U4>S4U5. The corresponding trends in the aqueous phase are S2U1>S2U3>S2U2>S2U5>S2U4 and S4U1>S4U2>S4U3>S4U5>S4U4. On the basis of the computed energy differences we can hypothesize that only the oxo‐thione forms of 2‐ and 4‐thiouracil should exist in the gas phase and in water solution. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 82: 44–52, 2001     

Solvent‐Polarity‐Tuned Morphology and Inversion of Supramolecular Chirality in a Self‐Assembled Pyridylpyrazole‐Linked Glutamide Derivative: Nanofibers,Nanotwists, Nanotubes,and Microtubes

The self‐assembly of a low‐molecular‐weight organogelator into various hierarchical structures has been achieved for a pyridylpyrazole linked L ‐glutamide amphiphile in different solvents. Upon gel formation, supramolecular chirality was observed, which exhibited an obvious dependence on the polarity of the solvent. Positive supramolecular chirality was obtained in nonpolar solvents, whereas it was inverted into negative supramolecular chirality in polar solvents. Moreover, the gelator molecules self‐assembled into a diverse array of nanostructures over a wide scale range, from nanofibers to nanotubes and microtubes, depending on the solvent polarity. Such morphological changes could even occur for the xerogels in the solvent vapors. We found that the interactions between the pyridylpyrazole headgroups and the solvents could subtly change the stacking of the molecules and, hence, their self‐assembled nanostructures. This work exemplifies that organic solvents can significantly involve the gelation, as well as tune the structure and properties, of a gel.     

Conventional strain energies of azetidine and phosphetane: Can density functional theory yield reliable results?

The conventional strain energies for azetidine and phosphetane are determined within the isodesmic, homodesmotic, and hyperhomodesmotic models. Optimum equilibrium geometries, harmonic vibrational frequencies, and corresponding electronic energies and zero‐point vibrational energies are computed for all pertinent molecular systems using self‐consistent field theory, second‐order perturbation theory, and density functional theory and using the correlation consistent basis sets cc‐pVDZ, cc‐pVTZ, and cc‐pVQZ. Single point fourth‐order perturbation theory, CCSD, and CCSD(T) calculations using the cc‐pVTZ and the cc‐pVQZ basis sets are computed using the MP2/cc‐pVTZ and MP2/cc‐pVQZ optimized geometries, respectively, to ascertain the contribution of higher order correlation effects and to determine if the quadruple‐zeta valence basis set is needed when higher order correlation is included. In the density functional theory study, eight different functionals are used including B3LYP, wB97XD, and M06‐2X to determine if any functional can yield results similar to those obtained at the CCSD(T) level. © 2012 Wiley Periodicals, Inc.     

On the vibrational linear and nonlinear optical properties of compounds involving noble gas atoms: HXeOXeH,HXeOXeF, and FXeOXeF

The vibrational (hyper)polarizabilities of some selected Xe derivatives are studied in the context of Bishop–Kirtman perturbation theory (BKPT) and numerical finite field methodology. It was found that for this set of rare gas compounds, the static vibrational properties are quite large, in comparison to the corresponding electronic ones, especially those of the second hyperpolarizability. This also holds for the dc‐Pockels β(?ω;ω,0), Kerr γ(?ω;ω,0,0) and electric field second harmonic generation γ (?2ω;ω,ω,0) effects, although the computed nuclear relaxation (nr) vibrational contributions are smaller in magnitude than the static ones. HXeOXeH was used to study the effects of electron correlation, basis set, and geometry. Geometry effects were found to lead to noticeable changes of the vibrational and electronic second hyperpolarizability. A limited study of the effect of Xe insertion to the nr vibrational properties is also reported. Assessment of the results revealed that Xe insertion has a remarkable effect on the nr (hyper)polarizabilities. In terms of the BKPT, this is associated with a remarkable increase of the electrical and mechanical anharmonicity terms. The latter is consistent with the anharmonic character of several vibrational modes reported for rare gas compounds. © 2013 Wiley Periodicals, Inc.     

Computer simulation of the linear and nonlinear optical susceptibilities of p-nitroaniline in cyclohexane, 1,4-dioxane, and tetrahydrofuran in quadrupolar approximation. I. Molecular polarizabilities and hyperpolarizabilities

This is the first part of a study of the local field effects on (non)linear optical susceptibilities of solutions of para-nitroaniline (pNA) in three different solvents, cyclohexane (CH), 1,4-dioxane (DI), and tetrahydrofuran (THF), using a discrete molecular representation of the condensed phase. To account for dipolar and quadrupolar effects, the latter of which are especially important for DI solution, all the electric properties necessary to compute the local fields and local field gradients in quadrupolar approximation as well as the dipolar hyperpolarizabilities for the four molecules are computed, including frequency dispersion and vibrational contributions to the dipolar properties. The convergence of the perturbation treatment for the pure vibrational (PV) contributions is examined by comparison of the values obtained at the lowest order with those of partially computed second order in mechanical and electrical anharmonicity. For pNA, for which previous computations of the hyperpolarizabilities have generally found poor agreement with experimental results, a thorough investigation of the effects of solvent-induced geometry changes, dynamic and static correlation, frequency dispersion, and classical thermal averaging over the torsional modes of the substituent groups and the inversion mode of the amino group on the dipolar properties is carried out. Computations using self-consistent continuum reaction field models show that the amino group is substantially less pyramidalized in polar solvents than in the gas phase. With all the effects taken into account, reasonable agreement with the experimental electric-field induced second harmonic generation (EFISH) result on pNA vapor of Kaatz, Donley, and Shelton is obtained.     

Testing computational models of hyperpolarizability in a merocyanine dye using spectroscopic and DFT methods

The structural and electronic properties of a highly solvatochromic merocyanine dye, 2-(3-cyano-5,5-dimethyl-4-(3-(1-octadecylpyridin-4(1H)-ylidene)prop-1-enyl)furan-2(5H)-ylidene)malononitrile (pyr3pi), have been investigated using UV-vis, NMR, hyper-Rayleigh scattering, and Raman spectroscopies and further interpreted using computational chemistry. Spectroscopic data indicate that pyr3pi exists in its zwitterionic form even in low polarity solvents with electronic absorption spectra showing a hypsochromic shift with an increase in solvent polarity and NMR experiments indicating an increasingly zwitterionic structure in chloroform as the temperature is lowered. Raman spectra in increasingly polar solvents show small variations of the structure that are consistent with a change toward a structure with more zwitterionic character. However, comparison of the calculated and experimental vibrational energies and intensities and comparison of NMR coupling constants with calculated bond order indicate that calculations underestimate the amount of charge separation seen in low polarity solvents. Although for this system density functional theory (DFT) calculations and the two-state model qualitatively reproduce negative solvatochromism, they fail to reproduce the trends in hyperpolarizability seen experimentally. This is attributed to solvent field DFT calculations underestimating the degree of charge separation in reaction fields representing low polarity solvents.     

Theoretical study on the solvent effect of the nitrosyl cation (NO+) generating reaction

Nitrosyl cation (NO⁺) generating reaction HONO + H⁺ → NO⁺ + H₂O has been theoretically investigated by B3LYP and high‐electron‐correlation QCISD methods with 6‐31G (d,p) basis set. The solvent effects on the geometries, reaction path properties, energies, thermodynamic, and kinetic characters in four solvents (benzene, tetrahydrofuran, acetonitrile, and water) have been calculated using self‐consistent reaction field (SCRF) approach with the polarizable continuum model (PCM). The results show that the activation energy barriers and the relative energies of the products are decreased with increase of the polarities of the solvents, and the reaction is favored in polar solvents thermodynamically and kinetically. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Solvent effects on the thioamide rotational barrier: an experimental and theoretical study.

The solvent effect on the C-N rotational barriers of N,N-dimethylthioformamide (DMTF) and N,N-dimethylthioacetamide (DMTA) has been investigated using ab initio theory and NMR spectroscopy. Selective inversion recovery NMR experiments were used to measure rotational barriers in a series of solvents. These data are compared to ab initio results at the G2(MP2) theoretical level. The latter are corrected for large amplitude vibrational motions to give differences in free energy. The calculated gas phase barriers are in very good agreement with the experimental values. Solvation effects were calculated using reaction field theory. This approach has been found to give barriers that are in good agreement with experiment for many aprotic, nonaromatic solvents that do not engage in specific interactions with the solute molecules. The calculated solution-phase barriers for the thioamides using the above solvents are also in good agreement with the observed barriers. The solvent effect on the thioamide rotational barrier is larger than that for the amides because the thioamides have a larger ground-state dipole moment, and there is a larger change in dipole moment with increasing solvent polarity. The transition-state dipole moments for the amides and thioamides are relatively similar. The origin of the C-N rotational barrier and its relation to the concept of amide "resonance" is examined.     

Anharmonic vibrational probe for N-methylacetamide in different micro-environments

In the present study, anharmonic vibrational properties of the amide modes in N-methylacetamide (NMA), a model molecule for peptide vibrational spectroscopy, are examined by DFT calculations. The 3N-6 normal mode frequencies, diagonal and off-diagonal anharmonicities are evaluated by means of the second order vibrational perturbation theory (VPT2). Good performance of B3LYP/6-31+G** is found for predicting vibrational frequencies in comparison with gas phase experimental data. The amide vibrational modes are assigned through potential energy distribution analysis (PED). The solvation effect on the amide vibrational modes is modeled within the PCM method. From gas phase to polar solvents, red shifts are observed for both harmonic and anharmonic vibrational frequency of amide I mode while the CO bond length increases upon the solvent polarity. Cubic and quartic force constants are further calculated to evaluate the origin of the anharmonicity for the amide I mode of NMA in different micro-environments.     

First Principles Investigation of Noncovalent Complexation: A [2.2.2]‐Cryptand Ion‐Binding Selectivity Study

A first principles methodology, aimed at understanding the roles of molecular conformation and energetics in host–guest binding interactions, is developed and tested on a system that pushes the practical limits of ab initio methods. The binding behavior between the [2.2.2]‐cryptand host (4,7,13,16,21,24‐hexaoxa‐1,10‐diaza‐bicyclo[8.8.8]hexacosane) and alkali metal cations (Li⁺, Na⁺, and K⁺) in gas, water, methanol, and acetonitrile is characterized. Hartree–Fock and density functional theory methods are used in concert with crystallographic information to identify gas phase, energy‐minimized conformations. Gas phase free energies of binding, with vibrational contributions, are compared to solution‐state binding constants through relative binding selectivity analysis. Calculated relative binding free energies qualitatively correlated with solution state experiments only after gas phase metal desolvation is considered. The B3LYP exchange–correlation functional improves theoretical correlations with experimental relative binding free energies. The relevance of gas phase calculations towards understanding binding in condensed phases is discussed. Natural bond orbital methods highlights previously unidentified intramolecular and intermolecular M⁺(222) chemistries, such as an intramolecular n′_O→σ*_CH hydrogen bond.     

A First‐Principles Study on the Electronic,Vibrational, and Thermodynamic Properties of Jadeite and its Tentative Low‐Density Polymorph

The high‐pressure clinopyroxene mineral jadeite (chemical composition NaAlSi₂O₆) was studied by density‐functional theory with respect to its electronic, vibrational, and thermodynamic properties, correctly reproducing the available experimental data. At a larger‐than‐normal volume, however, a low‐density alumosilicate phase with tetrahedral instead of octahedral Al–O coordination was identified. This low‐density phase was investigated theoretically, too, and the results were compared to jadeite and experimentally observed properties of what has been dubbed jadeite glass. It turned out that the theoretically obtained properties of this hypothetical polymorph, such as the bulk modulus, the molar volume and the vibrational frequencies, agree well with the corresponding properties of the glass phase. Hence, from an inverse structure‐property relationship we propose to model jadeite glass with the aforementioned low‐density alumosilicate phase and tetrahedral Al–O motifs, as suggested from first principles and corroborated from experiment. The possibilities as well as the limitations of electronic, phonon and thermodynamic properties calculations applied to such a polymorph are also discussed.     

Solvent‐Induced Structural Transition of Self‐Assembled Dipeptide: From Organogels to Microcrystals

Organogels that are self‐assembled from simple peptide molecules are an interesting class of nano‐ and mesoscale soft matter with simplicity and functionality. Investigating the precise roles of the organic solvents and their effects on stabilization of the formed organogel is an important topic for the development of low‐molecular‐weight gelators. We report the structural transition of an organogel self‐assembled from a single dipeptide building block, diphenylalanine (L ‐Phe‐L ‐Phe, FF), in toluene into a flower‐like microcrystal merely by introducing ethanol as a co‐solvent; this provides deeper insights into the phase transition between mesostable gels and thermodynamically stable microcrystals. Multiple characterization techniques were used to reveal the transitions. The results indicate that there are different molecular‐packing modes formed in the gels and in the microcrystals. Further studies show that the co‐solvent, ethanol, which has a higher polarity than toluene, might be involved in the formation of hydrogen bonds during molecular self‐assembly of the dipeptide in mixed solvents, thus leading to the transition of organogels into microcrystals. The structural transformation modulated by the co‐solvent might have a potential implication in controllable molecular self‐assembly.     

Detection of Significant Aprotic Solvent Effects on the Conformational Distribution of Methyl 4‐Nitrophenyl Sulfoxide: From Gas‐Phase Rotational to Liquid‐Crystal NMR Spectroscopy

The conformational equilibrium of methyl 4‐nitrophenyl sulfoxide (MNPSO) was experimentally investigated in the gas phase by using microwave spectroscopy and in isotropic and nematic liquid‐crystal solutions, in which the solvents are nonaqueous and aprotic, by using NMR spectroscopy; moreover, it was theoretically studied in vacuo and in solution at different levels of theory. The overall set of results indicates a significant dependence of the solute conformational distribution on the solvent dielectric permittivity constant: when dissolved in low‐polarity media, the most stable conformation of MNPSO proved to be strongly twisted with respect to that in more polar solvents, in which the conformational distribution maximum essentially coincides with that obtained in the gas phase. We discuss a possible explanation of this behavior, which rests on electrostatic solute–solvent interactions and is supported by calculations of the solute electric dipole moment as a function of the torsional angle. This function shows that the least polar conformation of MNPSO is located at a twist angle close to that of the conformational distribution maximum found in less‐polar solvents. This fact, associated with a relatively flat torsional potential, can justify the stabilization of the twisted conformation by the less‐polar solvents.     

Self‐Assembled Aggregates and Molecular Bilayer Films of a Double‐Chain Fullerene Lipid: Structure and Electrochemistry

The self‐aggregation behavior of C₆₀ fullerenes that bear two octadecyl chains (lipid 1 ) as well as the structures and electrochemical properties of cast films of 1 are described. We also examined the self‐aggregation behavior in organic solvents of three previously reported compounds: C₆₀ with three each of hexadecyl (lipid 2 ), tetradecyl (lipid 3 ), or dodecyl (lipid 4 ) chains. The fullerene lipids in alcohols spontaneously formed spherical aggregates, whose diameters are related to the alkyl‐chain lengths, concentrations of the fullerene lipids, and the solvent polarity. The morphologies of the aggregates showed temperature dependence. Cast films of 1 formed multimolecular bilayer structures that undergo a phase transition typical of lipid bilayer membranes. The electrochemistry of cast films of 1 on an electrode in aqueous medium exhibits temperature dependence.     

QM/MM calculation of solvent effects on absorption spectra of guanine

Electronic spectra of guanine in the gas phase and in water were studied by quantum mechanical/molecular mechanical (QM/MM) methods. Geometries for the excited‐state calculations were extracted from ground‐state molecular dynamics (MD) simulations using the self‐consistent‐charge density functional tight binding (SCC‐DFTB) method for the QM region and the TIP3P force field for the water environment. Theoretical absorption spectra were generated from excitation energies and oscillator strengths calculated for 50 to 500 MD snapshots of guanine in the gas phase (QM) and in solution (QM/MM). The excited‐state calculations used time‐dependent density functional theory (TDDFT) and the DFT‐based multireference configuration interaction (DFT/MRCI) method of Grimme and Waletzke, in combination with two basis sets. Our investigation covered keto‐N7H and keto‐N9H guanine, with particular focus on solvent effects in the low‐energy spectrum of the keto‐N9H tautomer. When compared with the vertical excitation energies of gas‐phase guanine at the optimized DFT (B3LYP/TZVP) geometry, the maxima in the computed solution spectra are shifted by several tenths of an eV. Three effects contribute: the use of SCC‐DFTB‐based rather than B3LYP‐based geometries in the MD snapshots (red shift of ca. 0.1 eV), explicit inclusion of nuclear motion through the MD snapshots (red shift of ca. 0.1 eV), and intrinsic solvent effects (differences in the absorption maxima in the computed gas‐phase and solution spectra, typically ca. 0.1–0.3 eV). A detailed analysis of the results indicates that the intrinsic solvent effects arise both from solvent‐induced structural changes and from electrostatic solute–solvent interactions, the latter being dominant. © 2009 Wiley Periodicals, Inc. J Comput Chem 2010     

2,6-二巯基吡啶互变异构平衡体系溶剂效应的理论研究

在气相及甲苯、氯仿、乙腈和水等溶剂中对2,6-二巯基吡啶及其硫酮式互变异构体进行了HF/6-31G^**水平上的优化,其中溶液中的计算采用Onsager自洽反应场(SCRF)模型.探讨了溶剂对体系几何结构和能量的影响.结果表明:溶剂的存在与极性的增加有利于平衡体系中硫酮式异构体的存在.