Statistical mechanics of energy transfer in gas–surface scattering: A dynamical Lie algebraic approach

The dynamical Lie algebraic (DLA) method is used to describe statistical mechanics of energy transfer in rotationally inelastic molecule–surface scattering. Statistical average values of an observable for the scattering system are calculated in terms of density operator formalism in statistical mechanics. Employing a cubic expansion procedure of molecule–surface interaction potential leads to generation of a dynamical Lie algebra. Thus these statistical average values as a function of the group parameters can be obtained analytically in this formulation. The group parameters can be found from solving a set of coupled nonlinear differential equations. The DLA method, which has no need for determination of transition probabilities in advance as made routinely in the calculation, offers an efficient alternative to the method for computing the statistical average values. This method is much less computationally intensive because most of calculations can be analytically carried out. The average final rotational energies and their dependence on the main dynamic variables and the average interaction potential are presented for the rotationally inelastic scattering of NO molecules from a flat, static Ag(111) surface. Direct comparison is made between the predictions of this model calculation and experiment. The model reproduces well the degree of rotational excitation and correlation between the average final translational and the average rotational energies. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

取代苯体系的二阶非线性光学性质:动力学李代数方法

提出了一种动力学李代数方法来研究取代苯体系的非线性光学性质. 对于给定的PPP模型(Pariser-Parr-Pople)哈密顿量, 生成了一个动力学李代数. 依据这些代数元构造出演化算子作为群参数的函数, 通过求解一组非线性微分方程能够得到这些群参数. 再按照统计力学中的密度算子公式给出取代苯分子体系偶极矩的统计平均值. 于是导出二阶极化率的表达式. 与其他量子力学计算结果比较, 表明这种动力学李代数方法在预言有机共轭分子的非线性光学性质上同样有用.     

Dynamic Lie algebraic formulation of second‐order optical nonlinearity for substituted benzenes

A dynamic Lie algebraic (DLA) formulation is applied to the study of nonlinear optical properties of the substituted benzenes. We have described the generation of the dynamic Lie algebra for the model Hamiltonian used in the present study. In terms of these elements of the dynamic algebra we express the evolution operator as a function of the group parameters, which can be determined by means of solving a system of coupled nonlinear differential equations. Thus, in terms of the density matrix operator formalism in statistical mechanics we obtain the statistical averages of the electric polarization and then derive readily an expression for the hyperpolarizability of the paradisubstituted benzenes. Comparisons with experimental observations and other quantum calculations of the hyperpolarizabilities for the paradisubstituted benzenes are made qualitatively and quantitatively. These results imply that the DLA method appears useful in describing the nonlinear optical phenomena in the substituted benzene molecules. © 2003 Wiley Periodicals, Inc. Int J Quantum Chem 93: 335–343, 2003     

Application of the Lie algebraic approach to diffractionally and rotationally inelastic molecule–surface scattering

The Lie algebraic approach of Alhassid and Levine [Phys. Rev. A 18 , 89 (1978)] is applied to the molecule–surface scattering. Specially, the diffractionally and rotationally inelastic scattering of a diatomic molecule from a solid surface is dealt with. Within the framework of the close-coupling method, we construct a Hamiltonian for the scattering system and use it to generate a dynamical algebra h₆. By solving equations of motion for the group parameters, the scattering wave functions near the surface are obtained. Computed transition probabilities of diffractively and rotationally inelastic scattering of H₂ from LiF(001) surface with the use of Lie algebraic method are seen to agree well with the coupled-channel calculations. The Lie algebraic method thus appears to have a wide range of validity for describing the dynamics of gas–surface scattering. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 981–989, 1997     

Rotationally inelastic molecule–surface scattering: dynamical lie algebraic method

The rotationally inelastic molecule–surface scattering is analyzed using dynamical Lie algebraic method. We treat, by example, the simple model of the scattering of NO from a rigid, flat Ag(111) surface. The explicit expressions of transition probability and the probability current density are obtained. It is proved that dynamical Lie algebraic method can be useful for describing the scattering problems. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 500–510, 2000     

The temperature dependence of the hard sphere diameter

The temperature dependence of the equivalent hard sphere diameter is analyzed using statistical mechanics. The analysis shows that algebraic equations for the hard sphere diameter which were developed by previous workers are not consistent with the statistical mechanics at high temperature. An algebraic equation which is consistent with the statistical mechanics is developed.     

Time-dependent wave packet propagation using quantum hydrodynamics

A new approach for propagating time-dependent quantum wave packets is presented based on the direct numerical solution of the quantum hydrodynamic equations of motion associated with the de Broglie–Bohm formulation of quantum mechanics. A generalized iterative finite difference method (IFDM) is used to solve the resulting set of non-linear coupled equations. The IFDM is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a “smart” Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is generalized to treat higher-dimensional problems and anharmonic potentials. The method is applied to a one-dimensional Gaussian wave packet scattering from an Eckart barrier, a one-dimensional Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. The 2D scattering results represent the first successful application of an accurate direct numerical solution of the quantum hydrodynamic equations to an anharmonic potential energy surface.     

分形表面上CO氧化反应的蒙特卡罗模拟

众所周知,人们在理想完好的单晶表面上的动力学研究已取得了很大的成功[1],然而对许多实际催化剂来说,传统催化理论和实验之间不符的情况仍然俯拾皆是．现在人们对造成这种情况的原因已经有了更加深入的认识：大量的实验事实表明,在表面催化反应体系中,不但催化剂表面具有复杂的分形结构,而且催化剂表面上的活性中心分布也具有复杂的分形特性[2]．“分形”是指那些具有分数维数的几何对象．这些对象往往是不规则的,不能用通常的欧氏几何来描述[3]．将分形引入催化科学中最早的是Pfeifer和Avnir等[4]人,至今已有近二十年的历史,现…     

On the determination of absolute vibrational excitation probabilities in molecule-surface scattering: Case study of NO on Au(111)

We describe a method to obtain absolute vibrational excitation probabilities of molecules scattering from a surface based on measurements of the rotational state, scattering angle, and temporal distributions of the scattered molecules and apply this method to the vibrational excitation of NO scattering from Au(111). We report the absolute excitation probabilities to the v = 1 and v = 2 vibrational states, rotational excitation distributions, and final scattering angle distributions for a wide range of incidence energies and surface temperatures. In addition to demonstrating the methodology for obtaining absolute scattering probabilities, these results provide an excellent benchmark for theoretical calculations of molecule-surface scattering.     

A+BC共线碰撞能量传递的动力学李代数结合中间绘景研究

在半经典近似下将动力学李代数方法和中间绘景结合应用于原子-双原子分子(非谐振子)共线碰撞中的平动-振动传能的研究。在群参量的一级近似下求解群参量的运动方程,进而确定时间演化算符。并以散射体系H₂+He为例,计算了H₂分子的振动跃迁几率,计算结果与精确量子力学计算结果符合得较好。     

Quantum hydrodynamics: capturing a reactive scattering resonance

The hydrodynamic equations of motion associated with the de Broglie-Bohm formulation of quantum mechanics are solved using a meshless method based upon a moving least-squares approach. An arbitrary Lagrangian-Eulerian frame of reference and a regridding algorithm which adds and deletes computational points are used to maintain a uniform and nearly constant interparticle spacing. The methodology also uses averaged fields to maintain unitary time evolution. The numerical instabilities associated with the formation of nodes in the reflected portion of the wave packet are avoided by adding artificial viscosity to the equations of motion. A new and more robust artificial viscosity algorithm is presented which gives accurate scattering results and is capable of capturing quantum resonances. The methodology is applied to a one-dimensional model chemical reaction that is known to exhibit a quantum resonance. The correlation function approach is used to compute the reactive scattering matrix, reaction probability, and time delay as a function of energy. Excellent agreement is obtained between the scattering results based upon the quantum hydrodynamic approach and those based upon standard quantum mechanics. This is the first clear demonstration of the ability of moving grid approaches to accurately and robustly reproduce resonance structures in a scattering system.     

State-to-state dynamics at the gas-liquid metal interface: rotationally and electronically inelastic scattering of NO[2Π(1/2)(0.5)] from molten gallium

Jet cooled NO molecules are scattered at 45° with respect to the surface normal from a liquid gallium surface at E(inc) from 1.0(3) to 20(6) kcal/mol to probe rotationally and electronically inelastic scattering from a gas-molten metal interface (numbers in parenthesis represent 1σ uncertainty in the corresponding final digits). Scattered populations are detected at 45° by confocal laser induced fluorescence (LIF) on the γ(0-0) and γ(1-1) A(2)Σ ← X(2)Π(Ω) bands, yielding rotational, spin-orbit, and λ-doublet population distributions. Scattering of low speed NO molecules results in Boltzmann distributions with effective temperatures considerably lower than that of the surface, in respectable agreement with the Bowman-Gossage rotational cooling model [J. M. Bowman and J. L. Gossage, Chem. Phys. Lett. 96, 481 (1983)] for desorption from a restricted surface rotor state. Increasing collision energy results in a stronger increase in scattered NO rotational energy than spin-orbit excitation, with an opposite trend noted for changes in surface temperature. The difference between electronic and rotational dynamics is discussed in terms of the possible influence of electron hole pair excitations in the conducting metal. While such electronically non-adiabatic processes can also influence vibrational dynamics, the γ(1-1) band indicates <2.6 × 10(-4) probability for collisional formation of NO(v = 1) at surface temperatures up to 580 K. Average translational to rotational energy transfer is compared from a hard cube model perspective with previous studies of NO scattering from single crystal solid surfaces. Despite a lighter atomic mass (70 amu), the liquid Ga surface is found to promote translational to rotational excitation more efficiently than Ag(111) (108 amu) and nearly as effectively as Au(111) (197 amu). The enhanced propensity for Ga(l) to transform incident translational energy into rotation is discussed in terms of temperature-dependent capillary wave excitation of the gas-liquid metal interface.     

Theoretical modeling of energy redistribution and stereodynamics in CF scattering from Si(100) under grazing incidence

We have simulated CF scattering from Si(100) using the molecular dynamics method. Translational energy loss spectra are presented. The shape of the energy loss distribution as a result of internal energy release is analyzed. At the classical turning point, the internal energy of the molecule is mainly in the form of rotational energy. The strong rotational excitation results in additional molecule-surfaces interactions during the latter half of the collision. These additional collisions permit some molecules that initially gain internal energy exceeding the bond strength to ultimately survive the collision process via rotational de-excitation. The rotational motion exhibited by surviving molecules is determined by the combination of the molecular axis orientation and the local surface structure during the collision process. The rotation planes of the surviving molecules are preferentially aligned with the surface normal (cartwheel-like and propeller-like motions). In this study, propeller-like motion of the surviving molecules is predicted. The majority of surviving molecules exhibit a cartwheel-like motion. However, molecules that gain a propeller-like rotation exhibit a much better alignment of their planes-of-rotation compared with molecules exhibiting cartwheel-like motion.     

Translational to rotational energy transfer in molecule-surface collisions

A theoretical approach that combines classical mechanics for treating translational and rotational degrees of freedom and quantum mechanics for describing the excitation of internal molecular modes is applied to the scattering of diatomic molecules from metal surfaces. Calculations are carried out for determining the extent of energy transfer to the rotational degrees of freedom of the projectile molecule. For the case of observed spectra of intensity versus final rotational energy, quantitative agreement with available experimental data for the scattering of NO and N(2) from close packed metal surfaces is obtained. It is shown that such measurements can be used to determine the average rotational energy of the incident molecular beam. Measurements of the exchange of energy between translational and rotational degrees of freedom upon collision are also described by calculations for these same systems.     

A stochastic, local mode study of neon-liquid surface collision dynamics

Equations of motion for a fast, light rare gas atom passing over a liquid surface are derived and used to infer the dynamics of neon collisions with squalane and perfluorinated polyether surfaces from experimental data. The equations incorporate the local mode model of a liquid surface via a stochastic process and explicitly account for impulsive collisional energy loss to the surface. The equations predict angular distributions for scattering of neon that are in good quantitative agreement with experimental data. Our key dynamical conclusions are that experimental angular distributions derive mainly from local mode surface topography rather than from structural features of individual surface molecules, and that the available data for these systems can be accounted for almost exclusively by single collisions between neon atoms and the liquid surface.     

Molecular theory of barycentric velocity: monatomic fluids

The notion of barycentric velocity appears in irreversible thermodynamics and fluid mechanics, in which it is a field variable obeying the hydrodynamic equations or, more specifically, the momentum balance equation, which is coupled to the rest of hydrodynamic equations. Therefore, its behavior is not known until the hydrodynamic equations are solved for the flow problem of interest. Unlike diffusion fluxes, heat fluxes, or stresses, it does not have its own constitutive relation similar to Fick's law, Fourier's law, and Newton's law of viscosity. In this work, the constitutive equation is derived for it. In parallel to the phenomenological notion of barycentric velocity, the notion of mean fluid velocity appears in statistical mechanics of irreversible dynamic processes according to the theory of Irving and Kirkwood [J. Chem. Phys. 18, 817 (1950)], and plays the same role of the phenomenological counterpart. In this work, we investigate the statistical mechanical meanings of the mean fluid velocity of a fluid in flow beyond its formal connection with the barycentric velocity. We show that it consists of two components; the center-of-gravity velocity of the packet of fluid molecules, which may be identified with the barycentric velocity in the phenomenological theory, and the diffusive contribution of its collective modes relative to the center of gravity. If the fluid is uniform in space or if the packet of fluid mass is rigid, the diffusive component vanishes. The statistical mechanical (molecular theory) formula for the mean fluid velocity provides the constitutive relation for it in terms of density and temperature gradients present in the fluid in flow. The constitutive relation obtained for the mean fluid velocity can be an important component in the theory of transport processes in liquids. Its significance to fluid mechanics is briefly discussed.     

Distribution functions and dynamical properties of stiff macromolecules

An analytically tractable model for chain molecules with bending stiffness is presented and the dynamical properties of such chains are investigated. The partition function is derived via the maximum entropy principle taking into account the chain connectivity as well as the bending restrictions in form of constraints. We demonstrate that second moments agree exactly with those known from the Kratky-Porod wormlike chain. Moreover, various distribution functions are calculated. In particular, the static structure factor is shown to be proportional to 1/q at large scattering vectors q. The equations of motion for a chain in a melt as well as in dilute solution are presented. In the latter case the hydrodynamic interaction is taken into account via the Rotne-Prager tensor. The dynamical equations are solved by a normal mode analysis. In the limit of a flexible chain the model reproduces the well-known Rouse and Zimm dynamics, respectively, on large length scales, whereas in the rod limit the eigenfunctions correspond to bending motion only. In addition, the coherent and incoherent dynamic structure factor is discussed. For melts we show that at large scattering vectors the incoherent dynamic structure factor is a universal function of only the combination q^8/3tp^1/3, where 1/(2p) is the persistence length of the macromolecules. The comparison of the theoretical results with quasielastic neutron and light scattering experiments of various polymers in solution and melt exhibits good agreement. Our investigations show that local stiffness strongly influences the dynamics of macromolecules on small length scales even for long and flexible chains.