Characteristics of atomic vibrational motion in a one-component defective crystal

Vibrations of atoms in a defective argon crystal are considered. Frequencies are calculated in the harmonic approximation and Mie and Einstein approximations. The vibrations are calculated for a set of local clusters differing in the position of a vacancy at different distances from a selected atom. Probabilities for these clusters are calculated within a quasichemical approximation of the lattice gas model. It is shown that when combined contributions from lateral interactions and vibrational motions are allowed for in the free crystal energy, there is an increase in the lattice constant upon a rise in temperature in all approximations. It is found that the frequencies calculated using the Mie model become closer to the frequency distribution in the harmonic approximation as the degree of crystal defectiveness increases.     

Effects of hydration on scale factors for ab initio force constants

Experimentally measured vibrational frequencies from the polar groups of peptides in aqueous solutions do not agree with frequencies calculated from scaled quantum mechanical force fields (SQMFF) using differential scale factors developed for molecules in the vapor phase. Measured stretching frequencies for carbonyl groups are more than 50 wavenumbers lower than the calculated values. On the other hand, frequencies for non-polar groups calculated using these scale factors are relatively accurate. Our goal is to develop a SQMFF that yields accurate calculated frequencies for peptides in aqueous solutions. To this end, we have calculated scale factors for ab initio force constants for formic acid, acetic acid, and acetone using a least squares fit of calculated and experimental frequencies. We compare these scale factors with changes observed in the ab initio force constants calculated for these molecules at various states of hydration. These force constants are calculated using fully optimized geometries for these hydrated molecules using the 4-31G basis. We present a comparison of the experimental and calculated frequencies, along with their potential energy distributions, for both vapor and aqueous phases. The results indicate that scale factors can simulate the effects of solvation on molecular force constants to yield accurate scaled ab initio force fields.     

Benchmarking the performance of density functional theory based Green's function formalism utilizing different self-energy models in calculating electronic transmission through molecular systems

Electronic transmission through a metal-molecule-metal system is calculated by employing a Green's function formalism in the scattering based scheme. Self-energy models representing the bulk and the potential bias are used to describe electron transport through the molecular system. Different self-energies can be defined by varying the partition between device and bulk regions of the metal-molecule-metal model system. In addition, the self-energies are calculated with different representations of the bulk through its Green's function. In this work, the dependence of the calculated transmission on varying the self-energy subspaces is benchmarked. The calculated transmission is monitored with respect to the different choices defining the self-energy model. In this report, we focus on one-dimensional model systems with electronic structures calculated at the density functional level of theory.     

Vibrational contributions to indirect spin-spin coupling constants calculated via variational anharmonic approaches

Zero-point vibrational contributions to indirect spin-spin coupling constants for N2, CO, HF, H2O, C2H2, and CH4 are calculated via explicitly anharmonic approaches. Thermal averages of indirect spin-spin coupling constants are calculated for the same set of molecules and for C2X4, X = H, F, Cl. Potential energy surfaces have been calculated on a grid of points and analytic representations have been obtained by a linear least-squares fit in a direct product polynomial basis. Property surfaces have been represented by a fourth-order Taylor expansion around the equilibrium geometry. The electronic structure calculations employ density functional theory, and vibrational contributions to indirect spin-spin coupling constants are calculated employing vibrational self-consistent-field and vibrational configuration-interaction methods. The performance of vibrational perturbation theory and various approximate variational calculations are discussed. Thermal averages are computed by state-specific and virtual vibrational self-consistent-field methods.     

Theoretical prediction of emulsion color

The perceived quality of many commercial products that are based on emulsions is determined by their color. In this article, a theory is presented to relate the color of emulsions to their composition and microstructure. First, the scattering characteristics (Qs and g) of individual droplets are calculated using Mie Theory. Second, the scattering (S) and absorption (K) coefficients of a concentrated emulsion are calculated using radiative transfer theory. Third, the reflectance spectrum (R) of the emulsion is calculated using Kubelka-Munk Theory. Finally, the tristimulus coordinates (XYZ, or L*a*b*) of the emulsion are calculated using color theory. There is excellent agreement between theoretical predictions and experimental measurements of the influence of droplet and chromophore characteristics on the tristimulus coordinates of concentrated oil-in-water emulsions.     

Experimental and theoretical study on cation radicals of cyclopropane,cyclobutane and cyclopentane

Cation radicals of cycloalkanes have been produced for the first time in γ-irradiated rigid solutions using fluoromethane as a matrix. Observed ESR spectra are analyzed by ab initio MO calculations; molecular geometries of the cations are optimized for several low-lying doublet states by the energy gradient method. Based on the optimized geometry thus determined the hyperfine splitting constants are calculated by the pseudo-orbital theory. Large Jahn-Teller distortions are calculated for the carbon rings and the H-C-H frames of the cation radicals. The distortions are consistent with the nodal picture of singly occupied orbitals. The calculated hfs constants are sensitive to the change in geometry due to a large contribution of spin-delocalization. The average of the calculated proton hfs constants is compatible with the observed ESR spectra which indicates that all the protons of the individual cations have become equivalent owing to the dynamic Jahn-Teller distortion.     

NMR shielding constants for hydrogen guest molecules in structure II clathrates

Proton NMR shielding constants and chemical shifts for hydrogen guests in small and large cages of structure II clathrates are calculated using density-functional theory and the gauge-invariant atomic-orbital method. Shielding constants are calculated at the B3LYP level with the 6-311++G(d,p) basis set. The calculated chemical shifts are corrected with a linear regression to reproduce the experimental chemical shifts of a set of standard molecules. The calculated chemical shifts of single hydrogen molecules in the small and large structure II cages are 4.94 and 4.84 ppm, respectively, which show that within the error range of the method the H2 guest molecules in the small and large cages cannot be distinguished. Chemical shifts are also calculated for double occupancy of the hydrogen guests in small cages, and double, triple, and quadruple occupancy in large cages. Multiple occupancy changes the chemical shift of the hydrogen guests by approximately 0.2 ppm. The relative effects of other guest molecules and the cage on the chemical shift are studied for the cages with multiple occupancies.     

Ne-HCl势能面和振转光谱的理论研究

利用量子化学计算方法CCSD(T)和大基组aug-cc-pVTZ加键函数3_s3_p2_d对Ne-HCl体系的分子间势能面进行了理论研究.结果表明,势能面上有两个势阱,分别对应于线性Ne-ClH和Ne-HCl构型.通过精确求解核运动方程发现,该从头算势能面分别支持5个(对Ne-HCl)和7个(Ne-DCl)振动束缚态.计算得到的振转跃迁频率与实值值吻合.     

Configuration-dependent properties of polyoxyethylene chains in the third-order interaction approximation

Conformational energies of polyoxyethylene chains in the third-order interaction approximation are calculated using a semiempirical potential energy function in the usual manner. Configuration-dependent properties of POE in the third-order are also calculated, and more close to the experimental data.     

Molecular mechanics (MM3) calculations on conjugated hydrocarbons

The MM3 molecular mechanics program has been extended to conjugated systems. A VESCF method is applied to the pi-system to calculate bond orders, from which various stretching and torsional parameters are obtained. The procedure gives somewhat better results than the analogous MM2 calculations. It has been applied to a study of 81 compounds of aromatic and other conjugated hydrocarbons, as well as 45 alkenes and unconjugated polyenes. The structures calculated are generally in good agreement with experiment, and the heats of formation of these compounds can be calculated with a rms value of 0.62 kcal/mol, which may be compared with the average experimental error of 0.61 kcal/mol. In addition, vibrational frequencies for five representative conjugated model structures are calculated, with an rms value of 46 cm^?1, and from these, other properties such as entropy can be calculated.     

Calculation of relaxed orbital binding energies in atoms: Zinc and cadmium

Relaxed orbital binding energies of electrons in zinc and cadmium are calculated using a relativistic local density functional theory. The discrepancy between our calculated binding energies and those from gas phase experimental data are attributed to inhomogeniety and many body effects.     

The accuracy of diffusion quantum Monte Carlo simulations in the determination of molecular equilibrium structures

For a test set of 17 first-row small molecules, the equilibrium structures are calculated with Ornstein-Uhlenbeck diffusion quantum Monte Carlo simulations guiding by trial wave functions constructed from floating spherical Gaussian orbitals and spherical Gaussian geminals. To measure performance of the Monte Carlo calculations, the mean deviation, the mean absolute deviation, the maximum absolute deviation, and the standard deviation of Monte Carlo calculated equilibrium structures with respect to empirical equilibrium structures are given. This approach is found to yield results having a uniformly high quality, being consistent with empirical equilibrium structures and surpassing calculated values from the coupled cluster model with single, double, and noniterative triple excitations [CCSD(T)] with the basis sets of cc-pCVQZ and cc-pVQZ. The nonrelativistic equilibrium atomization energies are also presented to assess performance of the calculated methods. The mean absolute deviations regarding experimental atomization energy are 0.16 and 0.21 kcal/mol for the Monte Carlo and CCSD(T)/cc-pCV(56)Z calculations, respectively.     

Kinetics and thermodynamics of the exchange reactions of peroxy radicals with hydroperoxides

The enthalpies and equilibrium constants of the exchange reactions of peroxy radicals with hydroperoxides of various structures are calculated. The experimental data on the reactions of hydrogen atom abstraction by the peroxy radicals from the hydroperoxides are analyzed, and the kinetic parameters characterizing these reactions are calculated using the intersecting parabolas method. The activation energies and rate constants for nine reactions of H atom abstraction by a peroxy radical from the OOH group of a peroxide are calculated using the above parameters. The geometric parameters of the transition states for the reactions are calculated. The low triplet repulsion plays an important role in the fast occurrence of the reactions. The polar interaction in the transition state is manifested in the reactions of the peroxy radicals with hydroperoxides containing a polar group.     

A quantum wave packet dynamics study of the N(2D) + H2 reaction

We report a dynamics study of the reaction N((2)D) + H(2) (v=0, j=0-5) --> NH + H using the time-dependent quantum wave packet method and a recently reported single-sheeted double many-body expansion potential energy surface for NH(2)(1(2)A' ') which has been modeled from accurate ab initio multireference configuration-interaction calculations. The calculated probabilities for (v=0, j=0-5) are shown to display resonance structures, a feature also visible to some extent in the calculated total cross sections for (v=0, j=0). A comparison between the calculated centrifugal-sudden and coupled-channel reaction probabilities validate the former approximation for the title system. Rate constants calculated using a uniform J-shifting scheme and averaged over a Boltzmann distribution of rotational states are shown to be in good agreement with the available experimental values. Comparisons with other theoretical results are also made.     

Theoretical determination of the excited states and of g-factors of the Creutz-Taube ion, [(NH3)5-Ru-pyrazine-Ru-(NH3)5]5+

The Creutz-Taube complex [(NH3)5-Ru-pyrazine-Ru-(NH3)5]5+ is studied using wave function-based methods, namely the CASSCF/MS-CASPT2 method. Spin-orbit effects have been calculated with the SO-RASSI program. The nature of the ground state is analyzed, and all the excited states up to 50,000 cm(-1) are calculated. They form a quasi continuum from 25,000 cm(-1), theoretical bands are assigned to UV-visible spectra, and MCD bands are assigned by calculating transition moments from first principles for both spectroscopies. g-Factors are calculated from first principles and modeled by a model Hamiltonian: they compare well to experimental values and are shown to be the same as in the monomeric species.     

A polarizable ion model for the structure of molten CuI

The results are reported of the molecular dynamics simulations of the coherent static structure factor of molten CuI at 938 K using a polarizable ion model. This model is based on a rigid ion potential to which the many body interactions due to the anions induced polarization are added. The calculated structure factor reproduces the clear sharp prepeak observed in neutron diffraction data. The corresponding partial structure factors and the related radial distribution functions calculated by molecular dynamics are compared with those found in the literature derived from a combination of neutron and x-ray diffraction data with the aid of the reverse Monte Carlo simulation technique, as well as those calculated by ab initio MD simulations.     

Helmholtz energy,extended corresponding states and local composition model for fluid mixtures

We present a model for the calculation of mixture properties. This model is based on the mixture residual Helmholtz energy given by the sum of two terms: one is the residual Helmholtz energy calculated by an extended corresponding states (ECS) model while the other is a correction term. The ECS model is based on methane as the reference fluid and shape factors that carry out the scaling of properties between the fluid of interest and the reference fluid. These shape factors are given by correlations in terms of reduced temperature and density. The application to mixtures is carried out with the one-fluid van der Waals mixture model. The correction term is temperature- and density-dependent as is given by a local composition mixing rule. Local compositions were calculated from a coordination number model based on lattice gas theory. With respect to binary-mixture properties, densities were calculated with an average absolute deviation (AAD) of 0.12%; speeds of sound were calculated with an AAD of 0.16% and bubble pressures were calculated with an AAD of 1.77%. Also, natural gas densities were calculated with AAD of 0.03% and natural gas speeds of sound were calculated with AAD of 0.049%. All these results are very satisfactory when compared with those obtained by modern mixture models.     

Poisson-Boltzmann与Donnan模型计算压实膨润土孔隙水与外部溶液间离子平衡的差异性比较

压实膨润土孔隙水与外部溶液之间的离子平衡是影响离子在压实膨润土中扩散的影响因素之一, 表征这一平衡的离子平衡系数可用压实膨润土的宏观属性参数通过Donnan 模型计算得到. 通过对膨润土主体矿物蒙脱石的TOT层结构单元进行简化, 构建了一个压实膨润土的单类孔隙结构模型, 辅以一个尺度变量H, 用Poisson-Boltzmann (PB)理论模型计算上述离子平衡系数. 对比计算结果, 发现PB模型计算的离子平衡系数总是大于Donnan模型的结果, 而参数H是联系这两种模型之间的桥梁. 通过对参数H取极限H→0, 实现了从PB模型到Donnan 模型的数学变换, 并从机理上讨论了两种模型之间的差异及其在实际扩散问题中的应用.分析表明PB模型更符合离子在压实膨润土中扩散的实际情况, 更适于处理实际扩散问题.     

Resonant vibrational excitation and de-excitation of CO(v) by low energy electrons

Integral cross sections and rate coefficients for vibrational excitation of the excited carbon-monoxide molecule, via the (2)Pi shape resonance in the energy region from 0 to 5 eV have been calculated. Cross sections are calculated by using our recently measured cross sections for the ground level CO excitation and the most recent cross sections for elastic electron scattering, applying the principle of detailed balance. Rate coefficients are calculated for Maxwellian electron energy distribution, with mean electron energies below 5 eV. By using extended Monte Carlo simulations, electron energy distribution functions (EEDF) and rate coefficients are determined in nonequilibrium conditions, in the presence of homogeneous external electric field. Nonequilibrium rates are calculated for typical, moderate values of the electric field over gas number density ratios, E/N, from 1 to 220 Td. Maxwellian and nonequilibrium rate coefficients are compared and the difference is attributed to a specific shape of the electron energy distribution functions under considered conditions.