Theoretical investigation of the binding energies of the iodide ion and xenon atom with decaborane

The interaction of decaborane (B(10)H(14)) with the I(-) ion and the (isoelectronic) Xe atom is investigated using a number of theoretical methods: MP2, CCSD(T), CCSD, spin-orbit CISD, and DFT using the B3LYP, B3PW91, PW91PW91, and PBE0 methods. All non-DFT and some DFT methods agree that B(10)H(14)I(-) is bound by charge-dipole electrostatic forces, charge- and dipole-induced-dipole forces, and dispersion forces, while B(10)H(14)Xe is bound by dipole-induced-dipole forces and dispersion forces. Counterpoise corrections are necessary to obtain reliable results. Relativistic effective core potentials were used for the I, Xe, and B atoms. Basis sets for use with these potentials are discussed as is the question of basis set balance in molecules. We find B(10)H(14)I(-) to be bound by 19.8 kcal/mol and B(10)H(14)Xe by 1.1 kcal/mol, indicating that the charge and polarizability of I(-) play the major role in the interaction energy.     

分子动力学模拟纳米尺寸限制体系下氩溶液中I2的振动能量弛豫

利用平衡态分子动力学方法(EMD)模拟了纳米尺寸限制球壳内I2在Ar溶液中的振动能量转移. 计算并讨论了I2振动能量弛豫时间T1随球壳半径、溶剂密度的变化规律. 通过分子间相互作用分析, 在原子、分子水平上, 揭示了随着球壳半径的减小, T1呈逐渐增大趋势的原因. 结果表明, 球壳的几何限制效应和表面作用对受限溶液密度分布的影响较大, 从而导致溶质振动弛豫的显著变化. 此外, 非限制体系模拟显示, 非平衡态分子动力学(NEMD)方法可以得到与平衡态分子动力学方法较一致的振动能量弛豫时间T1.     

Monte Carlo simulation of the influence of solvent on nucleic acid base associations

The results of Monte Carlo calculations of the association between nucleic acid bases in a nonpolar solvent (CCl4) are described. The influence of the solvent on planar and stacked associations of bases was examined by analyzing the total energy of the system, including solute-solute, solute-solvent, and solvent-solvent contributions. Good quantitative agreement with the available experimental data was obtained. Solute-solvent interactions are primarily determined by dispersion forces; consequently, solute-solvent interactions vertical to the solute plane that maximize dispersion interactions are most favored, and a rough proportionally between solute-solvent energy and the surface of the solute was observed. Analysis of solvent-solvent energy is not necessarily reduced when surface area decreases, contrary to the simple cavity concept. "Single molecule probe" calculations were performed to explain the differences in base associations in H2O and CCl4. In CCl4 dispersion forces dominate and planar complexes are stabilized by maximum exposure of molecular planes to the solvent. In H2O electrostatic forces dominate so that the most stable structures are stacked association that allow the maximum number of hydrophilic centers to be exposed to the solvent.     

Predictions of solvation and vaporization enthalpies. 1.n-alkane +n-alkane solutions

A simple method for the calculation of the enthalpy of solvation is presented and demonstrated for 35 n-alkane + n-alkane solutions at 25°C. There is a good agreement between the predicted and experimental values. The calculation was based on the separation of the solvation enthalpy into the cavity formation and solute-solvent interaction contributions. The former term was determined from the activation enthalpy of the solvent viscous flow and solute molar volume while the latter on the basis of the dispersion energy using van der Waals diameters for n-propyl group. The procedure was also successful in prediction of the vaporization enthalpy of C₅–C₁₇ n-alkanes.     

Vibrational dephasing and hydrodynamic effects on vibrational relaxation rates in acetophenone: Raman bandshape analysis

The isotropic component of Raman band for C=O stretching mode of acetophenone in solution was analyzed by estimating the correlation coefficient with reference to Lorentzian lineshape. In the intermediate region of solute/solvent concentration there is a sharp decrease in the correlation coefficient and there appears to be a transition from non-Lorentzian to Lorentzian lineshape. The vibrational relaxation rates have been estimated from the isotropic component of Raman band in different solvents. The rate is shown to be dependent on several macroscopic as well as microscopic properties of the solute-solvent system and intermolecular interactions. The hydrodynamic and dispersion forces appear to play a major role in determining the vibrational relaxation rate and the broadening of the bands.     

Molecular shape and orientational order. Effects in the energy of cavity formation in liquids

A recently proposed method for calculating the energy of cavity formation in liquids is presented in which the cavity formation process is described as work against the surface forces of the solvent, at the microscopic scale. The energy involved in the cavity formation process is, on the other hand, viewed as a strictly interaction potential energy and the reference cavity, which has the size and the shape of the space occupied by each molecule in the liquid, is considered as having short-range orientational order characteristic of the pure liquid. The method is successfully applied to binary alkane mixtures at infinite dilution whose components have different chemical structure (linear, cyclic and branched alkanes). The importance of the changes in the molecular order of the solute and the solvent occuring in the mixing process is emphasized.     

A response-function approach to the direct calculation of the transition-energy in a multiple-cluster expansion formalism

Approximate limiting collision induced dissociation (CID) probabilities are calculated as a function of the vibrational energy and temperature for HCl, HF and O₂, infinitely diluted in an inert gas, by assuming that only the molecules with an energy that is within kT below the threshold dissociate and do so with unit probability. Rotational energy is explicitly included in defining the threshold for dissociation. At the lowest temperatures of our calculations (2500 K for HCl and HF, 3500 K for O₂) the agreement between the calculated CID probabilities and those obtained by us earlier in fitting the experimental dissociation rates is quite remarkable for HCl and HF and satisfactory for O₂ at all vibrational energies. It is therefore argued that the ladder climbing model and the weak bias model for diatomic dissociation are not different in principle if rotational energy is properly accounted for.     

Low temperature luminescence spectra of the d10s2 complexes Cs2MX6 (M = Se,Te and X = Cl,Br). The Jahn—Teller effect in the Γ−4(3T1u) excited state

The emission spectra of the title compounds in microcrystalline form have been measured at 10 K. The extensive vibrational progression in the e_g mode is indicative of a tetragonal Jahn—Teller distortion in the Γ^?₄(³T_1u) excited state. The vibronic coupling of a threefold electronic state with a doubly degenerate e_g mode (T—e coupling), linear in the nuclear coordinates, has been reinvestigated considering spin—orbit coupling up to second order perturbation on energy levels which result from an a¹_1gt¹_1u electron configuration. For an estimation of Jahn—Teller coupling constants, the intensity distributions in the progressions were compared with the theoretical line shape functions which were obtained from a model which also permits the determination of potential energy minima and vibrational fundamentals of the excited state. The unusually large increase in the e_g vibrational frequency compared to the ground state is due to Jahn-Teller forces which distort the potential surface, yielding steeper excited state energy curves.     

On the correlation between bond-length change and vibrational frequency shift in hydrogen-bonded complexes: a computational study of Y...HCl dimers (Y = N2, CO, BF)

The H-Cl bond-length change and the harmonic vibrational frequency shift of the H-Cl stretch on formation of the linear isoelectronic Y...H-Cl complexes (Y = N(2), CO, BF) have been determined by ab initio computations at different levels of theory. These shifts are in agreement with predictions from a model based on perturbation theory and involving the first and second derivatives of the interaction energy with respect to displacement of the H-Cl bond length from its equilibrium value in the isolated monomer. At the highest level of theory, blue shifts were obtained for BF...HCl and CO...HCl, while red shifts were obtained for FB...HCl, OC...HCl, and N(2)...HCl. These vibrational characteristics are rationalized by considering the balance between the interaction energy derivatives obtained from the perturbative model. The widely believed correlation between the bond-length change and the sign of the frequency shift obtained on complexation is discussed and found to be unreliable.     

Theoretical study of the reaction H + ClCH3 → HCl + CH3

The reaction H + ClCH₃ has theoretically studied in a LEPS potential energy surface with a single-particle approximation for the methyl group. The LEPS adjustable parameters were selected to reach a good agreement with experimental values of activation energy and exothermicity. A wide set of quasi-classical trajectories for that system has been calculated within a energy range covering the significative values of relative velocities at temperatures between 300 and 1000 K. Calculated reactive cross sections increase with translational energy and with the initial vibrational level, but they are not influenced by rotational excitation of the reactants. Microscopic and total reaction rate constants have been obtained within the temperature range and agree quite well with available experimental results. Final energy distribution shows that most of the exoergicity is consumed in increasing the relative velocity of the products, while HCl molecules remain in their vibrational ground state.     

Controlling the Subtle Energy Balance in Protic Ionic Liquids: Dispersion Forces Compete with Hydrogen Bonds

The properties of ionic liquids are determined by the energy‐balance between Coulomb‐interaction, hydrogen‐bonding, and dispersion forces. Out of a set of protic ionic liquids (PILs), including trialkylammonium cations and methylsulfonate and triflate anions we could detect the transfer from hydrogen‐bonding to dispersion‐dominated interaction between cation and anion in the PIL [(C₆H₁₃)₃NH][CF₃SO₃]. The characteristic vibrational features for both ion‐pair species can be detected and assigned in the far‐infrared spectra. Our approach gives direct access to the relative strength of hydrogen‐bonding and dispersion forces in a Coulomb‐dominated system. Dispersion‐corrected density functional theory (DFT) calculations support the experimental findings. The dispersion forces could be quantified to contribute about 2.3 kJ mol^?1 per additional methylene group in the alkyl chains of the ammonium cation.     

Vibrational analysis of free and hydrogen bonded complexes of nicotinamide and picolinamide

Harmonic and anharmonic vibrations of free nicotinamide (NIA) and picolinamide (PIA) molecules together with their hydrogen bonded complexes H₂O–NIA and H₂O–PIA have been studied by means of density functional method. The calculation results of the vibrational spectra of free molecules have been investigated and are compared to the available experimental spectra. The vibrational wavenumbers of both molecules have also been calculated by polarizable continuum model (PCM) that represents the solvent as a polarizable continuum and places the solute in a cavity within the solvent (water is chosen as the solvent in this study). The results of PCM calculations and the H₂O–NIA, H₂O–PIA complexes, are used to investigate the H-bonding interactions of both molecules with the water molecule. The harmonic wavenumbers have been scaled by proper factors obtained by comparing the observed versus calculated wavenumbers and it is shown that anharmonic corrections on the vibrational spectra provided a better agreement between the observed and calculated wavenumbers compared to the results obtained by scaling factor method.     

Theoretical studies of solute vibrational energy relaxation at liquid interfaces

Recent advances in the theoretical understanding of solute vibrational energy relaxation at liquid interfaces and surfaces are described. Non-equilibrium molecular dynamics simulations of the relaxation of an initially excited solute molecule are combined with equilibrium force autocorrelation calculations to gain insight into the factors that influence the vibrational relaxation rate. Diatomic and triatomic nonpolar, polar, and ionic solute molecules adsorbed at the liquid/vapor interface of several liquids as well as at the water/CCl(4) liquid/liquid interface are considered. In general, the vibrational relaxation rate is significantly slower (a factor of 3 to 4) at the liquid/vapor and liquid/liquid interface than in the bulk due to the reduced density, which gives rise to a reduced contribution of the repulsive solvent-solute forces on the vibrational mode. The surface effects on the ionic solutes are much smaller (50% or less slower relaxation relative to the bulk). This is due to the fact that ionic solutes at the interface are able to keep part of their solvation shell to a degree that depends on their size. Thus, a significant portion of the repulsive forces is maintained. A high degree of correlation is found between the peak height of the solvent-solute radial distribution function and the vibrational relaxation rate. The relaxation rate at the liquid/liquid interface strongly depends on the location of the solute across the interface and correlates with the change in the density and polarity profile of the interface.     

Raman study of peak frequencies and linewidths of the vi mode of trichloroethane in solution

The vibrational displacements and isotropic Raman linewidths of the v₁ (CH₃ symmetric stretching) mode of 1,1,1-trichloroethane (TCE) have been measured in a number of solvents of diverse molecular properties. It was determined that the gas-solution frequency shifts are proportional to the solvent's polarizability, and thus could be interpreted on the basis of solution variations in the London dispersion forces. The isotropic linewidths, on the other hand, were found to be uncorrelated to solvent polarizability, indicating that dispersion interactions are not an important relaxation mechanism. On the other hand, a qualitative correlation was observed between bandwidth and the dipole moment of the solvent, suggesting that dipolar interactions may play a significant role in the line broadening of methyl group vibrations. An alternative explanation of the observed behavior was furnished by the Isolated Binary Collision model, which assumes that vibrational relaxation arises from collisionally modulated repulsive interactions in the liquid. A good correlation was observed between experimental linewidths and widths predicted by the IBC model.     

Hydrogen bond. Environmental effects on proton potential curves. An SCRF MO CNDO/2 calculation of a water dimer

A solvent and/or environmental effect has been introduced into the MO CNDO/2 calculation of a model hydrogen bonded system. Proton potential curves, potential energy, dipole moments, the polarizability component parallel to the reaction field, and the second order perturbation effects associated to the dispersion forces, have been studied as a function of a solute-solvent coupling parameter. This parameter may be related, through the self-consistent reaction field theory of solvent effects, to both the macroscopic dielectric properties of the solvent and to the local order (if any) around the solute. Numerical results corresponding to a water dimer are discussed.     

Competition between dissociation and exchange processes in a collinear A + BC collision. I. Exact quantum results

Exact quantum results for collision-induced dissociation on a reactive surface are presented. A modified LEPS potential-energy surface modeling the H + HD → H₂ + D system has been used. HD and H₂ bearing respectively 7 and 6 bound states. This system has been chosen because it displays significant reactive scattering for total energies above the dissociation threshold. Calculations have been performed using the time-dependent wavepacket method for two initial vibrational quantum numbers of the HD molecule (v = 0, 2). For each vibrational quantum number, two wavepackets with overlapping energy distributions have been run, covering a total energy range up to more than three times the dissociation energy. Comparison with previous collision-induced dissociation calculations shows that the dissociation is greatly enhanced by the presence of concomitant reactive scattering. A vibrational enhancement effect is also observed above the dissociation threshold; for higher energies the system exhibits a pronounced vibrational inhibition effect.     

A Monte Carlo simulation of energy transfer processes in thermal unimolecular reactions. II. An applications to polyatomic dissociation

The Monte Carlo method has been used to provide a numerical solution to the ro-vibrational master equation for the low pressure unimolecular decomposition of a polyatomic molecule. This type of solution is made possible through the use of a simple exponential transition probability function, that represents the efficiency with which energy transfer takes place between the reactant molecule and an unspecified heat bath gas. The Monte Carlo technique is used to generate random variables that are distributed in a manner prescribed by the transition probability function. In the case of the present simulation, these variables correspond to random energy jumps induced in the molecule through single collision events. In order to account for the energy dependence of the vibrational state densities, we have proposed that vibrational relaxation in the polyatomic takes place from a single vibrational mode. Under equilibrium conditions we are able to show that with this assumption, the Monte Carlo model is capable of reproducing molecular quantities, such as the average vibrational energy per molecule and the vibrational specific heat, that compare favourable with the corresponding values calculated from equilibrium statistical mechanics. The model has been applied to a study of the low pressure unimolecular decomposition of a series of polyatomics. For three of the molecules, CH₄, CD₄, and C₂H₆ the agreement between the calculated and the high temperature experimental rate constants is very good. The calculations indicate that a significant proportion of the molecules that dissociate are rotationally as well as vibrationally excited. Very few of the reactive molecules have a vibrational energy content equal to or greater than E₀, the dissociation energy. The extent of rotational excitation is found to be temperature dependent.     

Inclusion phenomena in a system designed for exclusion—A cavitary problem and a model

The Sephadex gels were introduced for their molecular-sieve-like properties. For some types of solute, however, they exhibit affinity for which the presence of water seems mandatory. Partitioning in these gels has been attributed to steric exclusion, hydrophobic interactions, dispersion interactions, etc. but may be due to a common structural feature, possibly a cavitary structure. Structural and partitioning similarities between the Sephadex gels and the cycloamyloses support this concept. A simple cavitary model vicinal water shell (CVS) is proposed in which the selectivity of the gel depends essentially on internal partitioning between a two phase system consisting of a vicinal-matrix-perturbed water shell lining the cavity and a core of normal (bulk) water. The CVS model cannot be used as yet to predict K_d values from the physico-chemical properties of the gels alone. The starting point is a known K_d value in one gel and the K_d value of the same solute in another gel is then calculated using the physical-chemical properties. The model was tested on several solutes and yielded promising results. Its implications are discussed.     

Product rotational distribution in the photodissociation of CH3CF2Cl

In the laser-induced unimolecular dissociation CH₃CF₂Cl an exceedingly narrow rotational distribution of the HCl product has been observed. The experimental evidence suggests that a few rotations around J = 16 are unusually open channels to either the vibrational relaxation of HCl or the dissociation.The dissociation is investigated in steady state as well as in a thermal beam by far-infrared emission.