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.     

Reflection of electrons and photons from solids bombarded by 0.1 - to 100-MeV electrons

Electron and photon reflection ratios (in number and energy) for absorbers bombarded by electrons have been computed with the ITS Monte Carlo system version 3. Electrons of energies from 0.1 to 100 MeV have been assumed normally incident on an effectively semi-infinite absorber. The absorbers considered are elemental solids of atomic numbers from 4 to 92. The data on the electron reflection ratios agree rather well with the experimental data collected from literature except some discrepancies when the number-reflection ratio is small. For photons, the number-reflection ratio increases with increasing energy, but the energy-reflection ratio shows a maximum around 10 MeV. Empirical equations for the electron reflection ratios and the photon energy-reflection ratio are given (for electrons, graphs only).     

Closer asymptotic approximations for the distributions of the power divergence goodness-of-fit statistics

Summary The members of the power divergence family of statistics all have an asymptotically equivalent χ² distribution (Cressie and Read [1]). An asymptotic expansion for the distribution function is derived which shows that the speed of convergence to this asymptotic limit is dependent on λ. Known results for Pearson'sX ² statistic and the log-likelihood ratio statistic then appear as special cases in a continuum rather than as separate (unrelated) expansions.     

A group additivity approach for the estimation of heat capacities of organic liquids and solids at 298 K

A group additivity method is described which provides heat capacity estimates of the condensed phase. The data base consists of 810 liquids and 446 solids. Group values for carbon in various common substitution and hybridization states and for 47 functional groups are provided. The standard error of estimation using this approach on this data base is 19.5 (liquids) and 26.9 J/ (mole K) (solids). This can be compared to typical experimental uncertainties of 8.12 and 23,4 J/ (mole K) associated with these measurements, respectively. Experimental uncertainties were estimated from the numerical differences obtained for a given substance from multiple independent literature reports.     

Evaluation of tocopherol recovery through simulation of molecular distillation process

DISMOL simulator was used to determine the best possible operating conditions to guide, in future studies, experimental works. This simulator needs several physical-chemical properties and often it is very difficult to determine them because of the complexity of the involved components. Their determinations must be made through correlations and/or predictions, in order to characterize the system and calculate it. The first try is to have simulation results of a system that later can be validated with experimental data. To implement, in the simulator, the necessary parameters of complex systems is a difficult task. In this work, we aimed to determe these properties in order to evaluate the tocopherol (vitamin E) recovery using a DISMOL simulator. The raw material used was the crude deodorizer distillate of soya oil. With this procedure, it is possible to determine the best operating conditions for experimental works and to evaluated the process in the separation of new systems, analyzing the profiles obtained from these simulations for the falling film molecular distillator.     

Accurate description of van der Waals complexes by density functional theory including empirical corrections

An empirical method to account for van der Waals interactions in practical calculations with the density functional theory (termed DFT-D) is tested for a wide variety of molecular complexes. As in previous schemes, the dispersive energy is described by damped interatomic potentials of the form C6R(-6). The use of pure, gradient-corrected density functionals (BLYP and PBE), together with the resolution-of-the-identity (RI) approximation for the Coulomb operator, allows very efficient computations for large systems. Opposed to previous work, extended AO basis sets of polarized TZV or QZV quality are employed, which reduces the basis set superposition error to a negligible extend. By using a global scaling factor for the atomic C6 coefficients, the functional dependence of the results could be strongly reduced. The "double counting" of correlation effects for strongly bound complexes is found to be insignificant if steep damping functions are employed. The method is applied to a total of 29 complexes of atoms and small molecules (Ne, CH4, NH3, H2O, CH3F, N2, F2, formic acid, ethene, and ethine) with each other and with benzene, to benzene, naphthalene, pyrene, and coronene dimers, the naphthalene trimer, coronene. H2O and four H-bonded and stacked DNA base pairs (AT and GC). In almost all cases, very good agreement with reliable theoretical or experimental results for binding energies and intermolecular distances is obtained. For stacked aromatic systems and the important base pairs, the DFT-D-BLYP model seems to be even superior to standard MP2 treatments that systematically overbind. The good results obtained suggest the approach as a practical tool to describe the properties of many important van der Waals systems in chemistry. Furthermore, the DFT-D data may either be used to calibrate much simpler (e.g., force-field) potentials or the optimized structures can be used as input for more accurate ab initio calculations of the interaction energies.     

Comparative methods in the determination of wax content and pour points of crude oils

In this research, differential scanning calorimetry (DSC) and gas chromatography is used to determine the wax content of fourteen crude oils of different sources. Different empirical equations were applied to compare the wax content of crude oils. For the fourteen crude oil samples with the wax content ranging from 7.5 to 43.8 mass%, it was observed that the results of empirical equations were in good agreement with those determined by DSC and GC. Accordingly, a correlation between ASTM pour point and the temperature at which 2 mass% of wax has precipitated out from crude oil is developed.     

Optimal convergence properties of kernel density estimators without differentiability conditions

Let X ₁, X ₂, ..., X _n be independent observations from an (unknown) absolutely continuous univariate distribution with density f and let % MathType!MTEF!2!1!+-% feaafeart1ev1aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLn% hiov2DGi1BTfMBaeXafv3ySLgzGmvETj2BSbqefm0B1jxALjhiov2D% aebbfv3ySLgzGueE0jxyaibaiiYdd9qrFfea0dXdf9vqai-hEir8Ve% ea0de9qq-hbrpepeea0db9q8as0-LqLs-Jirpepeea0-as0Fb9pgea% 0lrP0xe9Fve9Fve9qapdbaqaaeGacaGaaiaabeqaamaabaabcaGcba% GabmOzayaajaGaaiikaiaadIhacaGGPaGaeyypa0Jaaiikaiaad6ga% caWGObGaaiykamaaCaaaleqabaGaeyOeI0IaaGymaaaakmaaqadaba% Gaam4saiaacUfadaWcgaqaaiaacIcacaWG4bGaeyOeI0Iaamiwamaa% BaaaleaacaWGPbaabeaakiaacMcaaeaacaWGObGaaiyxaaaaaSqaai% aadMgacqGH9aqpcaaIXaaabaGaamOBaaqdcqGHris5aaaa!5356!\[\hat f(x) = (nh)^{ - 1} \sum\nolimits_{i = 1}^n {K[{{(x - X_i )} \mathord{\left/ {\vphantom {{(x - X_i )} {h]}}} \right. \kern-\nulldelimiterspace} {h]}}} \] be a kernel estimator of f(x) at the point x, \s-<x<, with h=h _n (h _n O and nh _n , as n) the bandwidth and K a kernel function of order r. Optimal rates of convergence to zero for the bias and mean square error of such estimators have been studied and established by several authors under varying conditions on K and f. These conditions, however, have invariably included the assumption of existence of the r-th order derivative for f at the point x. It is shown in this paper that these rates of convergence remain valid without any differentiability assumptions on f at x. Instead some simple regularity conditions are imposed on the density f at the point of interest. Our methods are based on certain results in the theory of semi-groups of linear operators and the notions and relations of calculus of finite differences.This research was supported in part by grants from the Natural Sciences and Engineering Research Council of Canada and the University of Alberta Central Research Fund.     

A review ofL=λW and extensions

A fundamental principle of queueing theory isL=W (Little's law), which states that the time-average or expected time-stationary number of customers in a system is equal to the product of the arrival rate and the customer-average or expected customer-stationary time each customer spends in the system. This principle is now well known and frequently applied. However, in recent years there have been extensions, such as H=G and the continuous, distributional, ordinal and central-limit-theorem versions, which show that theL=W relation, when viewed properly, has much more power than was previously realized. Moreover, connections have been established between H=G and other fundamental relations, such as the rate conservation law and PASTA (Poisson arrivals see time averages), which show that there is a much greater unity in the overall theory than was previously realized. This paper provides a review.This paper is dedicated to the memory of our colleague Professor Peter Franken (1937–1989), who contributed greatly to the subject of this paper and to queueing theory more generally.     

Bayes estimation with spherically symmetric,convex loss

Summary It is desired to estimate a parameter with the loss function of the formL(θ, a)=W(‖θ−a‖), where is convex, differentiable, and non-decreasing. With this structure a characterization of Bayes estimators is given. Also it is noted that if the sample space, , for the observation,X, is a complete separable metric space then a Bayes estimator exists.