Moment closure hierarchies for kinetic theories

This paper presents a systematicnonperturbative derivation of a hierarchy of closed systems of moment equations corresponding to any classical kinetic theory. The first member of the hierarchy is the Euler system, which is based on Maxwellian velocity distributions, while the second member is based on nonisotropic Gaussian velocity distributions. The closure proceeds in two steps. The first ensures that every member of the hierarchy is hyperbolic, has an entropy, and formally recovers the Euler limit. The second involves modifying the collisional terms so that members of the hierarchy beyound the second also recover the correct Navier-Stokes behavior. This is achieved through the introduction of a generalization of the BGK collision operator. The simplest such system in three spatial dimensions is a 14-moment closure, which also recovers the behavior of the Grad 13-moment system when the velocity distributions lie near local Maxwellians. The closure procedure can be applied to a general class of kinetic theories.     

Estimation of Multinomial Probabilities under a Model Constraint

In this paper estimation of the probabilities of a multinomial distribution has been studied. The five estimators considered are: unrestricted estimator (UE), restricted estimator (RE) (under model ), preliminary test estimator (PTE) based on a test of the model , shrinkage estimator (SE) and the positive-rule shrinkage estimator (PRSE). Asymptotic distributions of these estimators are given under Pitman alternatives and the asymptotic risk under a quadratic loss has been evaluated. The relative performance of the five estimators is then studied with respect to their asymptotic distributional risks (ADR). It is seen that neither of the preliminary test and shrinkage estimators dominates the other, though each fares well relative to the other estimators. However, the positive rule estimator is recommended for use for dimension 3 or more while the PTE is recommended for dimension less than 3.     

Measurement scheduling for recursive team estimation

We consider a decentralized LQG measurement scheduling problem in which every measurement is costly, no communication between observers is permitted, and the observers' estimation errors are coupled quadratically. This setup, motivated by considerations from organization theory, models measurement scheduling problems in which cost, bandwidth, or security constraints necessitate that estimates be decentralized, although their errors are coupled. We show that, unlike the centralized case, in the decentralized case the problem of optimizing the time integral of the measurement cost and the quadratic estimation error is fundamentally stochastic, and we characterize the -optimal open-loop schedules as chattering solutions of a deterministic Lagrange optimal control problem. Using a numerical example, we describe also how this deterministic optimal control problem can be solved by nonlinear programming.This research was supported in part by ARPA Grant N00174-91-C-0116 and NSF Grant NCR-92-04419.     

A note on smoothed estimating functions

The kernel estimate of regression function in likelihood based models has been studied in Staniswalis (1989,J. Amer. Statist. Assoc.,84, 276–283). The notion of optimal estimation for the nonparametric kernel estimation of semimartingale intensity (t) is proposed. The goal is to arrive at a nonparametric estimate of ₀=(t ₀) for a fixed pointt ₀ [0, 1]. We consider the estimator that is a solution of the smoothed optimal estimating equation is the optimal estimating function as in Thavaneswaran and Thompson (1986,J. Appl. Probab.,23, 409–417).     

AnL 1 smoothing spline algorithm with cross validation

We propose an algorithm for the computation ofL ₁ (LAD) smoothing splines in the spacesW _M (D), with . We assume one is given data of the formy _i =(f(t _i ) + _i , i=1,...,N with {itt_i} _i=1 ^{N D} , the _i are errors withE( _i )=0, andf is assumed to be inW _M . The LAD smoothing spline, for fixed smoothing parameter0, is defined as the solution,s , of the optimization problem (1/N) _i=1 ^N ¦y _i –g(t _i ¦+J _M (g), whereJ _M (g) is the seminorm consisting of the sum of the squaredL ₂ norms of theMth partial derivatives ofg. Such an LAD smoothing spline,s , would be expected to give robust smoothed estimates off in situations where the _i are from a distribution with heavy tails. The solution to such a problem is a thin plate spline of known form. An algorithm for computings is given which is based on considering a sequence of quadratic programming problems whose structure is guided by the optimality conditions for the above convex minimization problem, and which are solved readily, if a good initial point is available. The data driven selection of the smoothing parameter is achieved by minimizing aCV() score of the form .The combined LAD-CV smoothing spline algorithm is a continuation scheme in 0 taken on the above SQPs parametrized in, with the optimal smoothing parameter taken to be that value of at which theCV() score first begins to increase. The feasibility of constructing the LAD-CV smoothing spline is illustrated by an application to a problem in environment data interpretation.     

白藜芦醇分子的转动惯量和电偶极矩

通过分子轨道理论和杂化轨道理论推断出较稳定的白藜芦醇分子是平面型分子,然后根据白藜芦醇分子结构特点计算了该化合物的一种稳定异构体的转动惯量,用矢量合成法计算了其电偶极矩,为微波辅助白藜芦醇萃取理论研究提供转动惯量和电偶极矩的数据.     

New Monte Carlo method for the self-avoiding walk

We introduce a new Monte Carlo algorithm for the self-avoiding walk (SAW), and show that it is particularly efficient in the critical region (long chains). We also introduce new and more efficient statistical techniques. We employ these methods to extract numerical estimates for the critical parameters of the SAW on the square lattice. We find=2.63820 ± 0.00004 ± 0.00030=1.352 ± 0.006 ± 0.025v=0.7590 ± 0.0062 ± 0.0042 where the first error bar represents systematic error due to corrections to scaling (subjective 95% confidence limits) and the second bar represents statistical error (classical 95% confidence limits). These results are based on SAWs of average length 166, using 340 hours CPU time on a CDC Cyber 170–730. We compare our results to previous work and indicate some directions for future research.     

APPLICATION OF THE MOMENT METHODS TO ANALYSIS OF X-RAY DIFFRACTION LINE PROFILE FOR PA1010-BMI SYSTEM

Pure X-ray diffraction profiles have been analysed for polyamide 1010 and PA1010--BMI system by means of multipeak fitting resolution of X-ray diffraction.The methods ofvariance and fourth moment have been applied to determine the particle size and strainvalues for theparacrystalline materials.The results indicated that both variance and fourthmoment of X-ray diffraction line profile yielded approximately the same values of theparticle size and the strain.The particle sizes of(100)reflection have been found todecrease with increasing BMI content,whereas the strain values increased.     

三级相互作用对聚亚甲基链末端距四次矩的影响

本文用三级相互作用近似理论导出了高分子链末端距的四次矩的计算公式。计算聚亚甲基链的末端距四次矩,并同二级相互作用的结果进行了比较。对于聚亚甲基链,三级相互作用近似的四次矩小于二级相互作用近似的四次矩。     

Reduced-size polarized basis sets for calculations of molecular electric properties. IV. First-row transition metals

Recent studies of the perturbation-dependent basis sets have indicated the possibility of a significant reduction of the size of the usual CGTO sets without considerable loss of accuracy in calculations of molecular electric properties. The resulting (ZPolX) basis sets have been developed for several atoms of the first and second row of the Periodic Table. The same method of the ZPolX basis set generation is extended for the first-row transition metals and the corresponding contracted ZPolX basis sets of the size [6s5p3d1f] are determined for both nonrelativistic and scalar relativistic calculations. The performance of the ZPolX basis sets is verified in calculations on the first-row transition metal oxides at the level of the ROHF, ROHF/CASPT2, and ROHF/CCSD(T) approximations. Also the study of the dipole polarizability of TiCl₄ confirms the excellent features of these very compact basis sets. The ZPolX basis sets for nonrelativistic and relativistic calculations of molecular electric properties are available on the web page http://www.chem.uni.torun.pl/zchk/basis-sets.html.