PubMed Health⌕ Search

Biomedical subjects

Seiji Tsuzuki

Publications and source records attributed to Seiji Tsuzuki.

At least 19 recordsLinked to original sources

Fluorescence spectroscopic properties and crystal structure of a series of donor-acceptor diphenylpolyenes.

A series of p-nitro-p'-alkoxy(OR)-substituted (E,E,E)-1,6-diphenyl-1,3,5-hexatrienes (1a, R = Me; 1b, R = Et; 1c, R = n-Pr; 1d, R = n-Bu) were prepared. The absorption and fluorescence spectra in solution were almost independent of the alkoxy chain length. The absorption maximum showed only a small dependence on the solvent polarity, whereas the fluorescence maximum red-shifted largely as the polarity increased. The solid-state absorption and fluorescence spectra were red-shifted relative to those in low polar solvents and were clearly dependent on the alkoxy chain length. The fluorescence maxima for the crystals of 1b and 1d were observed at 635-650 nm, which were red-shifted by 40-50 nm relative to those for 1a and 1c. The Stokes shifts were all relatively small (3000-3500 cm-1). For all four compounds, the fluorescence decay curves in the solid state were able to be analyzed by single-exponential fitting to give the lifetimes of 1.1-1.3 ns. This indicates that the emission of 1a-d is not originated from an excimer or molecular aggregates, but from only one emitting monomeric species. The fluorescence quantum yields of 1a-d were considerably high compared with the values for organic solids, which is consistent with their monomeric origin of emission. Single-crystal X-ray structure analyses of 1a, 1c, and 1d showed that the crystal packing was dependent on the alkoxy chain length. The crystals of 1a and 1c had herringbone structure, whereas that of 1d had pi-stacked structure. Strong pi-pi interaction in the crystal of 1d would be the cause of the spectral red shifts relative to those for 1a and 1c. No observation of excimer fluorescence from crystal 1d can be attributed to the limited overlap between the pi-planes of the molecules due to its "slipped-parallel" structure.

Journal Article↗

Alcohol and proton transport in perfluorinated ionomer membranes for fuel cells.

To clarify the transport mechanisms of alcohols and proton in perfluorosulfonated ionomer (PFSI) membranes for fuel cells, four membranes having different equivalent weight (EW) values were examined. Membranes were immersed in methanol, ethanol, and 2-propanol to prepare a total of 12 samples, and membrane swelling, mass (alcohol and proton) transports, and interactions between alcohols and proton were investigated systematically in the fully penetrated state. The membrane expansion fraction theta and alcohol content lambda increased with decreasing the EW value for all the samples. The self-diffusion coefficients (D's) of the alkyl group and of OH (including protons) were measured separately by the pulsed-gradient spin-echo (PGSE)-NMR method and the D's also increased with decreasing the EW value. These results implied that the alcohols penetrate into the hydrophilic regions of the PFSI membranes and diffuse through the space expanded by the alcohols. The ionic cluster regions formed by the alcohols resemble those induced by water in the water swollen membrane, where protons dissociated from sulfonic acid groups transport through the regions together with water molecules. The D values decreased with increasing the molecular weight of alcohols. This trend was supported by activation energies Ea estimated from the Arrhenius plots of D in the temperature range from 30 to -40 degrees C. The PGSE-NMR measurements also revealed that protons move faster than the alkyl groups in the membranes. The proton transport by the Grotthuss (hopping) mechanism was facilitated by the increase of the alcohol content and the decrease of the molecular weight. This result was also supported by the experimental results of proton conductivity kappa and mobility u(H(+)). Density functional theory (DFT) calculations of the interaction energy DeltaE(int) between proton and alcohol (including OH) showed that the /DeltaE(int)/ increases with increasing the molecular weight of alcohols, which is in a inverse relationship with the kappa and u(H(+)) values. The proton transport depends strongly on the DeltaE(int) in the membranes.

Journal Article↗

How ionic are room-temperature ionic liquids? An indicator of the physicochemical properties.

Room-temperature ionic liquids (RTILs) are liquids consisting entirely of ions, and their important properties, e.g., negligible vapor pressure, are considered to result from the ionic nature. However, we do not know how ionic the RTILs are. The ionic nature of the RTILs is defined in this study as the molar conductivity ratio (Lambda(imp)/Lambda(NMR)), calculated from the molar conductivity measured by the electrochemical impedance method (Lambda(imp)) and that estimated by use of pulse-field-gradient spin-echo NMR ionic self-diffusion coefficients and the Nernst-Einstein relation (Lambda(NMR)). This ratio is compared with solvatochromic polarity scales: anionic donor ability (Lewis basicity), E(T)(30), hydrogen bond donor acidity (alpha), and dipolarity/polarizability (pi), as well as NMR chemical shifts. The Lambda(imp)/Lambda(NMR) well illustrates the degree of cation-anion aggregation in the RTILs at equilibrium, which can be explained by the effects of anionic donor and cationic acceptor abilities for the RTILs having different anionic and cationic backbone structures with fixed counterparts, and by the inductive and dispersive forces for the various alkyl chain lengths in the cations. As a measure of the electrostatic interaction of the RTILs, the effective ionic concentration (C(eff)), which is a dominant parameter for the electrostatic forces of the RTILs, was introduced as the product of Lambda(imp)/Lambda(NMR) and the molar concentration and was compared with some physical properties, such as reported normal boiling points and distillation rates, glass transition temperature, and viscosity. A decrease in C(eff) of the RTILs is well correlated with the normal boiling point and distillation rate, whereas the liquid-state dynamics is controlled by a subtle balance between the electrostatic and other intermolecular forces.

Journal Article↗

Intermolecular interactions of nitrobenzene-benzene complex and nitrobenzene dimer: significant stabilization of slipped-parallel orientation by dispersion interaction.

The CCSD(T) level interaction energies of eight orientations of nitrobenzene-benzene complexes and nine orientations of nitrobenzene dimers at the basis set limit have been estimated. The calculated interaction energy of the most stable slipped-parallel (C(s)) nitrobenzene-benzene complex was -4.51 kcal/mol. That of the most stable slipped-parallel (antiparallel) (C(2h)) nitrobenzene dimer was -6.81 kcal/mol. The interaction energies of these complexes are significantly larger than that of the benzene dimer. The T-shaped complexes are substantially less stable. Although nitrobenzene has a polar nitro group, electrostatic interaction is always considerably weaker than the dispersion interaction. The dispersion interaction in these complexes is larger than that in the benzene dimer, which is the cause of the preference of the slipped-parallel orientation in these complexes.

Benzene↗

Origin of the attraction in aliphatic C-H/pi interactions: infrared spectroscopic and theoretical characterization of gas-phase clusters of aromatics with methane.

An attractive intermolecular interaction between an aliphatic C-H bond and a pi-electron system (C-H/pi interaction) was characterized on the basis of infrared spectroscopy and high level ab initio calculations. Infrared spectroscopy was applied to several isolated methane clusters with benzene, toluene, p-xylene, mesitylene, and naphthalene in the gas phase, and the spectral changes of the C-H stretch bands in the methane moiety upon the cluster formation were observed. In the theoretical approach, interaction energies of the clusters were evaluated by high-level ab initio calculations. The forbidden symmetric C-H stretch transition weakly appeared in the IR spectra of the clusters, and it confirmed the small deformation of the methane moiety from the T(d)() symmetry, which was predicted by the ab initio calculations. On the other hand, the degenerated asymmetric C-H stretch band showed complicated splitting, which is qualitatively interpreted by a hindered rotor model. Low-frequency shifts upon the cluster formation were seen in the symmetric C-H stretch frequency, though the magnitude of the shifts was extremely small and no clear correlation with the interaction energy was found. On the other hand, the size of the calculated interaction energy well correlates with the polarizability of aromatics. The S(1)-S(0) electronic transition of the aromatic moiety was also observed, and it showed low-frequency shifts upon cluster formation. These results support the dominance of the dispersion interaction over the electrostatic and charge-transfer terms in the aliphatic C-H/pi interaction.

Benzene↗

Magnitude and directionality of the interaction energy of the aliphatic CH/pi interaction: significant difference from hydrogen bond.

The CCSD(T) level interaction energies of CH/pi complexes at the basis set limit were estimated. The estimated interaction energies of the benzene complexes with CH(4), CH(3)CH(3), CH(2)CH(2), CHCH, CH(3)NH(2), CH(3)OH, CH(3)OCH(3), CH(3)F, CH(3)Cl, CH(3)ClNH(2), CH(3)ClOH, CH(2)Cl(2), CH(2)FCl, CH(2)F(2), CHCl(3), and CH(3)F(3) are -1.45, -1.82, -2.06, -2.83, -1.94, -1.98, -2.06, -2.31, -2.99, -3.57, -3.71, -4.54, -3.88, -3.22, -5.64, and -4.18 kcal/mol, respectively. Dispersion is the major source of attraction, even if substituents are attached to the carbon atom of the C-H bond. The dispersion interaction between benzene and chlorine atoms, which is not the CH/pi interaction, is the cause of the very large interaction energy of the CHCl(3) complex. Activated CH/pi interaction (acetylene and substituted methanes with two or three electron-withdrawing groups) is not very weak. The nature of the activated CH/pi interaction may be similar to the hydrogen bond. On the other hand, the nature of other typical (nonactivated) CH/pi interactions is completely different from that of the hydrogen bond. The typical CH/pi interaction is significantly weaker than the hydrogen bond. Dispersion interaction is mainly responsible for the attraction in the CH/pi interaction, whereas electrostatic interaction is the major source of attraction in the hydrogen bond. The orientation dependence of the interaction energy of the typical CH/pi interaction energy is very small, whereas the hydrogen bond has strong directionality. The weak directionality suggests that the hydrogen atom of the interacting C-H bond is not essential for the attraction and that the typical CH/pi interaction does not play critical roles in determining the molecular orientation in molecular assemblies.

Chlorine↗

Two-dimensional structure control by molecular width variation with metal coordination.

The self-assembled monolayer of bipyridine derivative 1, which has two alkyl chains on each end, at the HOPG/1-phenyloctane interface was studied by in situ scanning tunneling microscopy (STM). The detailed mechanism of a spontaneous change in the monolayer packing pattern by Pd coordination was studied. Uncomplexed 1 existed in a bent form in the monolayer, and the alkyl chains were interdigitated, whereas Pd-complexed 1 was in a straight form and the alkyl chains were not interdigitated. An intermediate state of 1 was successfully observed during metal coordination. The structure was the bent form with noninterdigitated alkyl chains. Equilibrium intermolecular distances reported from ab initio calculations indicate that the molecular width of the central aromatic part of uncomplexed 1 (7.5 A) is substantially smaller than that of the peripheral alkyl chain part (9.2 A). The bent form was suitable for covering up the surface to maximize the packing density. However, the molecular width of the aromatic unit of Pd-complexed 1 (9.1 A) was almost identical to that of the alkyl chain unit (9.2 A). Therefore, Pd-complexed 1 took the straight form in the monolayer. The observation of surface coverage by STM suggests that the bent form increases the packing density by as much as 16% compared with that of the straight form. These results indicate that the control of molecular width can be used to design molecular templates for nanostructure formation.

Journal Article↗

Magnitude of the CH/pi interaction in the gas phase: experimental and theoretical determination of the accurate interaction energy in benzene-methane.

The accurate CH/pi interaction energy of the benzene-methane model system was experimentally and theoretically determined. In the experiment, mass analyzed threshold ionization spectroscopy was applied to the benzene-methane cluster in the gas phase, prepared in a supersonic molecular beam. The binding energy in the neutral ground state of the cluster, which is regarded as the CH/pi interaction energy for this model system, was evaluated from the dissociation threshold measurements of the cluster cation. The experimentally determined binding energy (D(0)) was 1.03-1.13 kcal/mol. The interaction energy of the model system was calculated by ab initio molecular orbital methods. The estimated CCSD(T) interaction energy at the basis set limit (D(e)) was -1.43 kcal/mol. The calculated binding energy (D(0)) after the vibrational zero-point energy correction (1.13 kcal/mol) agrees well with the experimental value. The effects of basis set and electron correlation correction procedure on the calculated CH/pi interaction energy were evaluated. Accuracy of the calculated interaction energies by DFT methods using BLYP, B3LYP, PW91 and PBE functionals was also discussed.

Journal Article↗

Estimated MP2 and CCSD(T) interaction energies of n-alkane dimers at the basis set limit: comparison of the methods of Helgaker et al. and Feller.

The MP2 (the second-order Møller-Plesset calculation) and CCSD(T) (coupled cluster calculation with single and double substitutions with noniterative triple excitations) interaction energies of all-trans n-alkane dimers were calculated using Dunning's [J. Chem. Phys. 90, 1007 (1989)] correlation consistent basis sets. The estimated MP2 interaction energies of methane, ethane, and propane dimers at the basis set limit [EMP2(limit)] by the method of Helgaker et al. [J. Chem. Phys. 106, 9639 (1997)] from the MP2/aug-cc-pVXZ (X=D and T) level interaction energies are very close to those estimated from the MP2/aug-cc-pVXZ (X=T and Q) level interaction energies. The estimated EMP2(limit) values of n-butane to n-heptane dimers from the MP2/cc-pVXZ (X=D and T) level interaction energies are very close to those from the MP2/aug-cc-pVXZ (X=D and T) ones. The EMP2(limit) values estimated by Feller's [J. Chem. Phys. 96, 6104 (1992)] method from the MP2/cc-pVXZ (X=D, T, and Q) level interaction energies are close to those estimated by the method of Helgaker et al. from the MP2/cc-pVXZ (X=T and Q) ones. The estimated EMP2(limit) values by the method of Helgaker et al. using the aug-cc-pVXZ (X=D and T) are close to these values. The estimated EMP2(limit) of the methane, ethane, propane, n-butane, n-pentane, n-hexane, n-heptane, n-octane, n-nonane, and n-decane dimers by the method of Helgaker et al. are -0.48, -1.35, -2.08, -2.97, -3.92, -4.91, -5.96, -6.68, -7.75, and -8.75 kcal/mol, respectively. Effects of electron correlation beyond MP2 are not large. The estimated CCSD(T) interaction energies of the methane, ethane, propane, and n-butane dimers at the basis set limit by the method of Helgaker et al. (-0.41, -1.22, -1.87, and -2.74 kcal/mol, respectively) from the CCSD(T)/cc-pVXZ (X=D and T) level interaction energies are close to the EMP2(limit) obtained using the same basis sets. The estimated EMP2(limit) values of the ten dimers were fitted to the form m0+m1X (X is 1 for methane, 2 for ethane, etc.). The obtained m0 and m1 (0.595 and -0.926 kcal/mol) show that the interactions between long n-alkane chains are significant. Analysis of basis set effects shows that cc-pVXZ (X=T, Q, or 5), aug-cc-pVXZ (X=D, T, Q, or 5) basis set, or 6-311G** basis set augmented with diffuse polarization function is necessary for quantitative evaluation of the interaction energies between n-alkane chains.

Journal Article↗

Physicochemical properties and structures of room-temperature ionic liquids. 3. Variation of cationic structures.

A series of room-temperature ionic liquids (RTILs) were prepared with different cationic structures, 1-butyl-3-methylimidazolium ([bmim]), 1-butylpyridinium ([bpy]), N-butyl-N-methylpyrrolidinium, ([bmpro]), and N-butyl-N,N,N-trimethylammonium ([(n-C(4)H(9))(CH(3))(3)N]) combined with an anion, bis(trifluoromethane sulfonyl)imide ([(CF(3)SO(2))(2)N]), and the thermal property, density, self-diffusion coefficients of the cation and anion, viscosity, and ionic conductivity were measured over a wide temperature range. The self-diffusion coefficient, viscosity, ionic conductivity, and molar conductivity follow the Vogel-Fulcher-Tamman equation for temperature dependencies, and the best-fit parameters have been estimated, together with the linear fitting parameters for the density. The relative cationic and anionic self-diffusion coefficients for the RTILs, independently determined by the pulsed-field-gradient spin-echo NMR method, appear to be influenced by the shape of the cationic structure. A definite order of the summation of the cationic and anionic diffusion coefficients for the RTILs: [bmim][(CF(3)SO(2))(2)N] > [bpy][(CF(3)SO(2))(2)N] > [bmpro][(CF(3)SO(2))(2)N] > [(n-C(4)H(9))(CH(3))(3)N][(CF(3)SO(2))(2)N], has been observed, which coincides with the reverse order to the viscosity data. The ratio of molar conductivity obtained from the impedance measurements to that calculated by the ionic diffusivity using the Nernst-Einstein equation quantifies the active ions contributing to ionic conduction in the diffusion components and follows the order: [bmpro][(CF(3)SO(2))(2)N] > [(n-C(4)H(9))(CH(3))(3)N][(CF(3)SO(2))(2)N] > [bpy][(CF(3)SO(2))(2)N] > [bmim][(CF(3)SO(2))(2)N] at 30 degrees C.

Journal Article↗

Intermolecular interaction between hexafluorobenzene and benzene: Ab initio calculations including CCSD(T) level electron correlation correction.

The intermolecular interaction energy of hexafluorobenzene-benzene has been calculated with the ARS-E model (a model chemistry for the evaluation of the intermolecular interaction energy between aromatic systems using extrapolation), which was formerly called the AIMI model. The CCSD(T) interaction energy at the basis-set limit has been estimated from the MP2 interaction energy at the basis-set limit and the CCSD(T) correction term obtained using a medium-sized basis set. The slipped-parallel (Cs) complex has the largest (most negative) interaction energy (-5.38 kcal/mol). The sandwich (C6v) complex is slightly less stable (-5.07 kcal/mol). The interaction energies of two T-shaped (C2v) complexes are very small (-1.74 and -0.88 kcal/mol). The calculated interaction energy of the slipped-parallel complex is about twice as large as that of the benzene dimer. The dispersion interaction is found to be the major source of attraction in the complex, although electrostatic interaction also contributes to the attraction. The dispersion interaction increases the relative stability of the slipped-parallel benzene dimer and the hexafluorobenzene-benzene complex compared to T-shaped ones. The electrostatic interaction is repulsive in the slipped-parallel benzene dimer, whereas it stabilizes the slipped-parallel hexafluorobenzene-benzene complex. Both electrostatic and dispersion interactions stabilize the slipped-parallel hexafluorobenzene-benzene complex, which is the cause of the preference of the slipped-parallel orientation and the larger interaction energy of the complex compared to the benzene dimer.

Benzene↗

Crystalline-state Z,E-photoisomerization of a series of (Z,E,Z)-1,6-diphenylhexa-1,3,5-triene 4,4'-dicarboxylic acid dialkyl esters. Chain length effects on the crystal structure and photoreactivity.

[reaction: see text] Crystalline-state Z,E-photoisomerization of a series of (Z,E,Z)-1,6-diphenylhexa-1,3,5-triene 4,4'-dicarboxylic acid dialkyl (R) esters [(Z,E,Z)-1a, R = Me; (Z,E,Z)-1b, R = Et; (Z,E,Z)-1c, R = n-Pr; (Z,E,Z)-1d, R = n-Bu] was investigated. All Z,E,Z isomers underwent one-way isomerization to the corresponding E,E,E isomers. The reaction efficiency was strongly enhanced as the length of the alkyl chain increased. Single-crystal X-ray analyses of (Z,E,Z)-1a-d showed that the alkyl chain part of the crystals became larger as the chain length increased. The conformational flexibility of the alkyl chains made the large change in the triene geometry in the lattice possible, leading to the enhancement of the photoreactivity in the crystalline state.

Journal Article↗

Synthesis and structure of 16 pi octaalkyltetraphenylporphyrins.

The oxidized octaethyltetraphenylporphyrin (1, OETPP) and the corresponding newly prepared octaisobutyltetraphenylporphyrin (3, OisoBuTPP) could be isolated from the reaction of OETPPLi2 (or OisoBuTPPLi2) with SOCl2. The X-ray analysis and the characteristic UV-vis spectra of 1 and 3 revealed that these are the first examples of 16 pi nonaromatic porphyrins.

Crystallography, X-Ray↗

Regulation of growth-blocking peptide expression during embryogenesis of the cabbage armyworm.

Growth-blocking peptide (GBP) is an insect cytokine with diverse biological functions. Northern blot analysis revealed high heterogeneity in the size distribution of GBP mRNAs as well as in the tissues where they are detected. The spatio-temporal transcription pattern is dynamic, especially during embryogenesis. Gel shift assays demonstrated that the cabbage armyworm embryo nuclear extract specifically binds to a 178-bp element, at position +234 to +411 from the transcription start site of the 1.3 kb GBP transcript, in which two Drosophila Deformed (Dfd) binding sites are repeated in tandem. The specific binding between this element and Dfd was demonstrated using recombinant cabbage armyworm Dfd protein. Silencing the Dfd expression in embryos by treating with Dfd double-stranded RNA did not reduce the expression level of GBP, but ectopic GBP expression was observed in the lateral region of the embryo, suggesting that Dfd could serve as a transcriptional repressor for the GBP gene.

Amino Acid Sequence↗

Molecular dynamics study of the methanol effect on the membrane morphology of perfluorosulfonic ionomers.

The membranes of a perfluorosulfonic acid polymer swollen in 10-80 wt % methanol solution were investigated to elucidate the methanol effect on their morphologies, such as size of the solvent cluster, solvent location, and polymer structure, by using isothermal-isobaric molecular dynamics simulations. In higher methanol concentrations, we found less-spherical solvent aggregation and a more spread polymer structure because of the ampholytic nature of methanol. The partial radial distribution functions between solvent oxygen and fluorocarbons, which are composed of the main chain, clearly show that methanol is located closer to the polymer matrix than water. On the other hand, water is preferentially located in the vicinity of an acidic headgroup, SO(3)(-), compared with methanol, although both have similar attractive interaction energies to the acidic group. Furthermore, we discussed solvent dynamics and hydrogen bonding between sulfonic oxygen and solvent O-H groups.

Journal Article↗

Magnitude and directionality of interaction in ion pairs of ionic liquids: relationship with ionic conductivity.

The intermolecular interaction energies of nine ion pairs of room temperature ionic liquids were studied by MP2/6-311G level ab initio calculations. The magnitude of the interaction energies of 1-ethyl-3-methylimidazolium (emim) complexes follows the trend CF(3)CO(2)(-) > BF(4)(-) > CF(3)SO(3)(-) > (CF(3)SO(2))(2)N(-) approximately PF(6)(-) (-89.8, -85.2, -82.6, -78.8, and -78.4 kcal/mol, respectively). The interaction energies of BF(4)(-) complexes with emim, ethylpyridinium (epy), N-ethyl-N,N,N-trimethylammonium ((C(2)H(5))(CH(3))(3)N), and N-ethyl-N-methylpyrrolidinium (empro) are not very different (-85.2, -82.8, -84.6, and -84.4 kcal/mol, respectively), while the size of the orientation dependence of the interaction energies follows the trend emim > epy approximately (C(2)H(5))(CH(3))(3)N > empro. Comparison with the experimental ionic conductivities shows that the magnitude and directionality of the interaction energy of the ion pairs play a crucial role in determining the ionic dissociation/association dynamics in the ionic liquids. The electrostatic interaction is the major source of attraction between ions. The induction contribution is small but not negligible. The hydrogen bonding with the C(2)-H of imidazolium is not essential for the attraction in the ion pair. The interaction energy of the BF(4)(-) complex with 1-ethyl-2,3-dimethylimidazolium (em2im) (-81.8 kcal/mol) is only 4% smaller than that of the emim complex.

Journal Article↗

Lewis acidity/basicity of pi-electron systems: theoretical study of a molecular interaction between a pi system and a Lewis acid/base.

Molecular interactions between pi systems having different pi-electron character (benzene, hexafluorobenzene, and borazine), and a Lewis acid/base (borane and ammonia) were theoretically studied. An attractive interaction between benzene, the electron-rich pi system, and borane was observed. On the other hand, repulsive interactions between benzene and ammonia was observed when the lone pair of nitrogen points toward the benzene ring. In contrast, an attractive interaction between hexafluorobenzene, an electron-deficient pi system, and ammonia was observed. Unexpectedly, a weak attractive interaction between hexafluorobenzene and borane was also observed. Borazine shows an interaction both to borane and ammonia. The attraction between the nitrogen atom of borazine and borane was larger than that between the boron atom of borazine and ammonia.

Journal Article↗

Ab initio calculations of structures and interaction energies of toluene dimers including CCSD(T) level electron correlation correction.

The intermolecular interaction energy of the toluene dimer has been calculated with the ARS-F model (a model chemistry for the evaluation of intermolecular interaction energy between ARomatic Systems using Feller's method), which was formerly called as the AIMI model III. The CCSD(T) (coupled cluster calculations with single and double substitutions with noniterative triple excitations) interaction energy at the basis set limit has been estimated from the second-order Moller-Plesset perturbation interaction energy at the basis set limit obtained by Feller's method and the CCSD(T) correction term obtained using a medium-size basis set. The cross (C(2)) dimer has the largest (most negative) interaction energy (-4.08 kcal/mol). The antiparallel (C(2h)) and parallel (C(S)) dimers (-3.77 and -3.41 kcal/mol, respectively) are slightly less stable. The dispersion interaction is found to be the major source of attraction in the toluene dimer. The dispersion interaction mainly determines the relative stability of the stacked three dimers. The electrostatic interaction of the stacked three dimers is repulsive. Although the T-shaped and slipped-parallel benzene dimers are nearly isoenergetic, the stacked toluene dimers are substantially more stable than the T-shaped toluene dimer (-2.62 kcal/mol). The large dispersion interaction in the stacked toluene dimers is the cause of their enhanced stability.

Dimerization↗