Molecular potential-energy surfaces for chemical reaction dynamics

This paper reviews the construction of molecular potential-energy surfaces by an interpolation method which has been developed over the last several years. The method uses ab initio quantum chemistry calculations of the molecular electronic energy in an automated procedure to construct global potential- energy surfaces which can be used to simulate chemical reactions with either classical or quantum dynamics. The methodology is explained and several applications are presented to illustrate the approach. Received: 22 February 2002 / Accepted: 2 May 2002 / Published online: 6 November 2002 Correspondence to: M. A. Collins e-mail: collins@rsc.anu.edu.au Acknowledgements. The methods described in this overview are the result of collaborations with former members of my group, in particular with Josef Ischtwan, Meredith Jordon, Keiran Thompson and Ryan Bettens. I am also indebted for inspiration gained from many discussions with my colleagues Leo Radom and Donghui Zhang (National University of Singapore). This work has been supported by the Supercomputer Facility of the Australian National University and the Australian Partnership for Advanced Computing.     

Theoretical study of the effect of frit quality on chromatography using computational fluid dynamics

Summary Frits at both ends of a chromatographic column, especially for a preparative column, have significant influence on the flow distribution within the column and thus the column efficiency. However, frits have received little attention from chromatographers in the past. Here a theoretical study was conducted with the aid of CFD software FLUENT to investigate the effect of frits on the performance of homogeneous and heterogeneous chromatographic columns. A dimensionless number,FQ, was applied to characterize frit quality. This study visualized how frit quality affects the flow distribution and the concentration band, the shape of eluted pulse at the colum exit and column efficiency. Simulation results show that the development length of the flow distribution is related toFQ but has nothing to do with the packing heterogeneity. The curvature of the concentration band in a column depends onFQ and packing quality. This study shows column efficiency can be improved significantly by increasingFQ and/or frit permeability.     

Theoretical studies on the reaction path and dynamics of the reaction CH_3+H_2O→CH_4+OH

The reaction path, the dynamical properties along the reaction path and CVT rate constants are computed by the ab initio MO method, the reaction path Hamiltonian theory and the variational transition state theory. The results show that the effect of the electron correlation energy on activation barrier is large, the recrossing and tunneling effects exist in the reaction.     

On a class of least-element complementarity problems

The present paper studies the linear complementarity problem of finding vectorsx andy inR ₊ ⁿ such thatc + Dx + y 0,b – x 0 andx ^T (c + Dx + y) = y ^T (b – x) = 0 whereD is aZ-matrix andb > 0. Complementarity problems of this nature arise, for example, from the minimization of certain quadratic functions subject to upper and lower bounds on the variables. Two least-element characterizations of solutions to the above linear complementarity problem are established first. Next, a new and direct method to solve this class of problems, which depends on the idea of least-element solution is presented. Finally, applications and computational experience with its implementation are discussed.Research partially supported by the National Science Foundation Grant MCS 71-03341 A04 and the Air Force Office of Scientific Research Contract F 44620 14 C 0079.     

新型含氮螯合树脂的制备及其吸附性能

以三甘醇和苯磺酰氯为原料,二乙烯三胺为交联剂合成了新型含氮螯合树脂(5),其结构经IR表征.讨论了Ni2 浓度和pH对5吸附容量的影响.动力学研究表明,5对Ni2 的吸附速率符合鲛岛公式.     

Ion beam mixing of α-MoSi2

Abstract

Using molecular dynamics simulation, we study the mixing of an α-MoSi₂ crystal by 1 keV Ar ions. The observed order of magnitude of the mixing is compatible with a spike model of ion beam mixing. The influence of the target surface and of chemical effects are discussed.     

Short Communication: Application of WENO (Weighted Essentially Non-oscillatory) Approach to Navier–Stokes Computations

The direct implementation of the essentially non-oscillatory schemes for flow simulation over complex geometries sometimes results in insufficiently robust numerical algorithms. In order to overcome this difficulty, it is suggested to use the weighted essentially non-oscillatory approach for multidimensional Navier–Stokes computations. The results indicate a significant improvement in accuracy and robustness, especially for low Mach and high supersonic flows.     

Exclusion process and droplet shape

We use a mathematical isomorphism between the one-dimensional exclusion process and the two-dimensional stochastic Ising model in the low-temperature limit to describe the typical instantaneous shape of a supercritical droplet. We derive, specifically, the exact asymptotic distribution of the boundaries of a (+1) spin region, confined to Z ₊ ² and subjected to a positive magnetic field. In an appropriate scaling, the boundary distribution converges to a deterministic continuum limit.On leave at the EPFL, CH-1015 Lausanne, Switzerland.     

Some solvable models of nonuniform classical fluids

A model classical fluid is constructed by assuming that the direct correlation functionc(r – r) is independent of any applied external field. Thermodynamic consistency requires thatc(r – r) 0, and permits explicit representation of the model by a many-body interaction potential. In the canonical ensemble, the model shows a phase transition to an infinite density condensed phase, but in the grand canonical ensemble only an anomalous transition to zero density vapor is found to stably exist.     

Nonequilibrium dynamics of a diffusion-limited reaction driven by a cluster-memory mechanism

The diffusion-limited reaction A+AA+B is studied in general dimension. The asymptotic decay of the system is found to depend nontrivially upon the initial concentration of A particles for certain ranges of the diffusion constant, backward reaction rate, and total concentration of particles. This nonequilibrium behavior is due to the formation of clusters centered about the initial A particles. A perturbative analysis ind=1 shows that the transition to the nonequilibrium dynamics is sharp and is quite similar to another previously studied reaction A+AA. Ford>1, a scaling argument is presented which describes the dependence of the asymptotic decay on the initial concentration of A particles and the equilibrium concentration for large backward reaction rates. Monte Carlo data are shown which confirm the analytic work ind=1, 2, and 3.