Distributed Gaussian discrete variable representation

A discrete variable representation (DVR) made from distributed Gaussians g_n(x) = e, (n = ?∞, …, ∞) and its infinite grid limit is described. The infinite grid limit of the distributed Gaussian DVR (DGDVR) reduces to the sinc function DVR of Colbert and Miller in the limit c → 0. The numerical performance of both finite and infinite grid DGDVRs and the sinc function DVR is compared. If a small number of quadrature points are taken, the finite grid DGDVR performs much better than both infinite grid DGDVR and sinc function DVR. The infinite grid DVRs lose accuracy due to the truncation error. In contrast, the sinc function DVR is found to be superior to both finite and infinite grid DGDVRs if enough grid points are taken to eliminate the truncation error. In particular, the accuracy of DGDVRs does not get better than some limit when the distance between Gaussians d goes to zero with fixed c, whereas the accuracy of the sinc function DVR improves very quickly as d becomes smaller, and the results are exact in the limit d → 0. An analysis of the performance of distributed basis functions to represent a given function is presented in a recent publication. With this analysis, we explain why the sinc function DVR performs better than the infinite grid DGDVR. The analysis also traces the inability of Gaussians to yield exact results in the limit d → 0 to the incompleteness of this basis in this limit. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

稀有气体亚稳态原子与CCl4和CBr4的能量转移反应

CBr_4和CCl_4分子的解离反应前人已做了许多工作,他们分别采用射频放电、电子轰击、He~ 的传能反应等方法研究了CCl_4和CBr_4的解离反应,得到了CCI(A)、CCl~ 、CBr~ 等碎片的发射光谱。有关亚稳态原子与它们的传能反应,只有某些较简单的报道,对传能反应机理也未作深入探讨。本文研究了各种亚稳态原子He(2~3S)、Ne(~3P_(0.2))、Ar(~3P_(0.2))与CCl_4和CBr_4分子的传能反应,并对反应机理进行了初步的讨论和分析。     

Gaussian basis sets for calculation of spin densities in first-row atoms

Summary The suitability of Gaussian basis sets for ab initio calculation of Fermi contact spin densities is established by application to the prototype first-row atoms B-F having open shell p electrons. Small multiconfiguration self-consistent-field wave functions are used to describe relevant spin and orbital polarization effects. Basis sets are evaluated by comparing the results to highly precise numerical grid calculations previously carried out with the same wave function models. It is found that modest contracted Gaussian basis sets developed primarily for Hartree-Fock calculations can give semiquantitative results if augmented by diffuse functions and if further uncontracted in the outer core-inner valence region.     

Structure formation in heavy metal oxyhydrates at low rates of gel formation

Oxyhydrate gels have a hydrophilic surface, due to which they undergo destruction and secondary polymerization in aqueous media. Prolonged storage in aqueous solution gives rise to regions with selfsimilar helical ordering in gels. Structuring of this kind is also observed when synthesis is conducted under conditions that provide low gelation rates. Electromagnetic UV and visible radiation is another means to change the gel structure; it makes the oligomer species pass into the excited state, due to which one of the directions of structuring becomes dominant. This work summarizes the results of computer simulation of gel agglomerates. For oxyhydrate systems, helical ordering was found to be one of the local energy minima. The units of a macrohelix can lie at various angles relative to one another, and they can change, after absorption of energy, the helix pitch and the order of elements in the helix.     

Semi‐classical formulation of time‐dependent discrete variable representation method

In this article, we apply a novel time‐dependent discrete variable representation (TDDVR) method proposed by Barkakaty and Adhikari to investigate tunneling through an Eckart barrier. This semi‐classical method is theoretically rigorous and straightforward to implement. Among the TDDVR formulations, this report presents the first derivation of a rigorous form of quantum force (QF) for the present perspective. The validity of this semi‐classical approach is demanded based on the excellent agreement of the tunneling probability with the corresponding quantum results. © 2004 Wiley Periodicals, Inc. Int J Quantum Chem, 2004     

A six-dimensional quadrature procedure

An account is given of the use of Gaussian quadrature product formulae in the evaluation of certain six-dimensional, two-centre integrals involving one-electron Green's functions. These integrals occur in a new molecular variational principle recently proposed by Hall, Hyslop and Rees [1] from which an approximate energy may be derived which can be shown to be at least as good as that obtained from the Rayleigh-Ritz principle. Reductions in computing time are realized by removing certain singularities using a subtraction technique and also by using an empirically determined Richardson-type extrapolation formula.This paper was presented during the session on numerical integration methods for molecules of the 1970 Quantum Theory Conference in Nottingham. It has been revised in the light of the interesting discussion which followed.     

The Discrete Approximation of Continuous-state Nonlinear Branching Processes in Lévy Random Environments北大核心CSCD

In this paper, we establish the discrete approximation of continuous-state nonlinear branching processes in Lévy random environments by using tightness and convergence sequence in infinite dimensional product space via stochastic differential equations. Taking α-stable branching as an example, the conditions which are given to discretize continuous-state nonlinear branching processes in Lévy random environments are verified. © 2022 Chinese Academy of Sciences. All rights reserved.     

Damage spreading in a single-component irreversible reaction process: Dependence of the system's immunity on the Euclidean dimension

The spreading of a globally distributed damage, created in the stationary regime, is studied in a single-component irreversible reaction process, i.e., the BK model [Browne and Kleban,Phys. Rev. A 40, 1615 (1989)]. The BK model describes one variant of the A+AA₂ reaction process on a lattice in contact with a reservoir of A species. The BK model has a single parameter, namely the rate of arrival of A species to the lattice (Y). The model, exhibits an irreversible phase transition between a stationary reactive state with production of A₂ species and a poisoned state with the lattice fully covered by A species. The transition takes place at critical points (Y _C ) which solely depend on the Euclidean dimensiond. It is found that the system is immune ford=1 andd=2, in the sense that even 100% of initial damage is healed within a finite healing period (T _H ). Within the reactive regime,T _H diverges when approachingY _C according toT _H (Y _C –Y)^–, with 1.62 and 1.08 ford=1 andd=2, respectively. Ford=3 a frozen-chaotic transition is found close toY _s 0.4125, i.e., well inside the reactive regime 0YY _C 0.4985. Just atY _S the damageD(t) heals according toD(t) t ^–, with 0.71. For the frozen-chaotic transition atd=3 the order parameter critical exponent 0.997 is determined.     

Necessary conditions for characterization of laws via mixed sums

SupposeX andY are independent and identically distributed, and independent ofU which satisfies 0U1. Recent work has centered on finding the lawsL(X) for whichX U(X+Y) where denotes equality in law. We show that this equation corresponds to a certain projective invariance property under random rotations. Implicitly or explicitly, it has been assumed that the characteristic function ofX has an expansion property near the origin. We show that solutions may be admitted in the absence of this condition when –logU has a lattice law. A continuous version of the basic problem replaces sums with a Lévy process. Instead we consider self-similar processes, showing that a solution exists only whenU is constant, and then all processes of a given order are admitted.This research was in part supported by NSERC grant A-8466.     

Cramér theorem on symmetric spaces of noncompact type

We prove the Cramér theorem forK-invariant Gaussian measures on irreducible symmetric spacesX=G/K withG semisimple noncompact. To do this we use a kind of Abel transform ofK-invariant measures onX.This research is supported by KBN Grant.