Evaluation of bridge function for hard sphere fluid confined in a narrow slit-pore via TMMC Mayer-sampling

The bridge function required to yield a singlet integral equation (IE) up to the second order in density expansion for the hard sphere fluid confined in a slit-pore is evaluated. The slit-fluid bridge function can be divided into wall-particle bridge diagrams with h _b-bond, which were evaluated by recently proposed Transition Matrix Monte Carlo (TMMC) Mayer-sampling method. The bulk-fluid total correlation function h _b(r) used in cluster integrals is determined by solution of the bulk-fluid Ornstein–Zernike (OZ) equation with a hypernetted chain closure (HNC). The calculation is performed for the reduced density of bulk fluid in equilibrium with the fluid in slit-pores from 0.3 to 0.7 with narrow slit width of 3.0σ and 4.0σ. The quantity of the slit-fluid bridge function is assessed by comparison of the density profile obtained from the singlet IE theory and the grand canonical Monte Carlo (GCMC) simulation. Good agreement between the proposed approach and the GCMC data is observed. The reduced normal pressure is also calculated, and agrees well with the simulation data at low to medium densities but becomes a little larger at high density. It is expected that the result can be improved by adding higher order bridge coefficients. The direct evaluation of the slit-fluid bridge function seems to be practical since a great improvement of the quality of the singlet IE theory has been achieved for predicting the structural and thermodynamic properties of fluids confined in narrow slit pores.     

Molecular dynamics calculation of the dielectric constant without periodic boundary conditions. I

We outline the difficulties in obtaining a reliable value of the dielectric constant of a fluid using molecular dynamics calculations with periodic boundary conditions, and give some explanation of the observed asymptotic behaviour of h_D (r) and hΔ(r) in Monte Carlo simulations of dipolar hard spheres. An alternative method consisting in simulating a dielectric in vacuum is described. This is applied to two dimensional systems. The pertinent theoretical relations for a dielectric disc in vacuum are therefore derived. It is concluded that relations involving MC or MD computation of <m ²> must be carefully handled.     

Monte Carlo calculation of y(r) for the hard-sphere fluid

The function y(r) = exp {βu(r)}g(r) is calculated for hard spheres in the region r < σ using umbrella-sampling Monte Carlo techniques. The resulting values are found to be well represented over the entire range 0 < r < σ by a simple function proposed by Grundke and Henderson.     

E.S.R. evidence for the planarity of the t-butyl radical

Monte Carlo calculations are reported for the radial distribution function g ₂(r; λ) of a fluid in which the intermolecular pair potential is [u ^ref(r) + λu ^p(r)], u ^ref(r) being the Weeks-Chandler-Andersen (WCA) reference fluid, and [u ^ref(r) + u ^p(r)] being the Lennard-Jones (6, 12) fluid. The calculations are performed for λ values in the range 0 to 1, at the state condition ρσ³ = 0·80, kT/ε = 0·719. It is shown that at high densities the perturbation expansion of g ₂(r; λ = 1) about g ₂(r; λ = 0) is rapidly convergent, but that the corresponding expansion for y ₂(r; λ) = exp [βu(r; λ)] × g ₂(r; λ) is not. In addition Monte Carlo estimates of the individual terms that contribute to the first-order perturbation term, (?g ₂/?λ)_λ=0, are presented. It is shown that these terms are individually large, but that (?g ₂/?λ)_λ=0 is small because there is strong cancellation between the various terms. Consequently, the calculation of (?g ₂/?λ)_λ=0 is highly sensitive to the approximation used to evaluate the individual terms.     

Analytic solution of the RISM equation for symmetric diatomics with Yukawa closure

We present the site-site direct correlation function c(r) for a fluid of hard diatomic symmetric molecules obtained from Monte Carlo simulation data via the RISM integral equation. This c(r) ensures that the site-site correlation function given by the RISM equation is exact, and thus provides a basis for critically examining the usual closure for the RISM equation. As an example of an improved closure we present the analytic solution of the RISM integral equation with a Yukawa closure for c(r).     

The bridge function of hard spheres by direct inversion of computer simulation data

The bridge function of the hard sphere fluid has been calculated from our new highly accurate Monte Carlo and molecular dynamics simulation data on the radial distribution function using the (inverted) Ornstein-Zernike equation. Both the systematic errors (finite size, grid size, tail) and statistical errors are analysed in detail and ways to suppress them are proposed. Uncertainties in the resulting values of B(r) are about 0.001. In contrast with many previous findings the bridge function is both positive and negative.     

Bilinear quantum Monte Carlo: Expectations and energy differences

We propose a bilinear sampling algorithm in the Green's function Monte Carlo for expectation values of operators that do not commute with the Hamiltonian and for differences between eigenvalues of different Hamiltonians. The integral representations of the Schrödinger equations are transformed into two equations whose solution has the form _a(x) t(x, y)_b(y), where _a and _b are the wavefunctions for the two related systems andt(x, y) is a kernel chosen to couplex andy. The Monte Carlo process, with random walkers on the enlarged configuration spacex y, solves these equations by generating densities whose asymptotic form is the above bilinear distribution. With such a distribution, exact Monte Carlo estimators can be obtained for the expectation values of quantum operators and for energy differences. We present results of these methods applied to several test problems, including a model integral equation, and the hydrogen atom.     

具有D4h对称性构型的C2+4分子的Jahn-Teller效应与能级分裂

依据量子理论与配位场理论,利用群论和对称性分析的方法探讨了C²⁺₄分子在具有D_4h对称性构型时,E×(b_1g+b_2g)系统的Jahn-Teller效应中的相关问题.研究了C²⁺₄分子的电子态与声子态的对称性及其活跃声子态,讨论了系统声子间的耦合与CG系数,构建了E×(b_1g+b_{2g 关键词： 2+4分子')" href="#">C²⁺4分子 对称性 能级分裂 Jahn-Teller畸变}     

Density expansion of the radial distribution and bridge functions of the hard sphere fluid

The bridge function and the background correlation function (and consequently the radial distribution function) of the pure hard sphere fluid are expanded up to the sixth power in density. The calculations are based on the Ree–Hoover representation of the diagrams and Monte Carlo integration. The coefficients as functions of the particle–particle separation are fitted to splines taking into account discontinuities in higher derivatives up to the term of the order of (r- const)⁵.     

核物质状态方程中动量相关作用可能的探针

研究了动量相关作用对于中子-质子比动能谱rb(Ek)的效应,发现rb(Ek)灵敏的依赖于动量相关作用而弱的依赖于介质中核子-核子碰撞截面和对称势.因此rb(Ek)是提取重离子碰撞中动量相关作用信息的可能探针.同时,对于丰中子弹核和相同质量稳定弹核在相同入射道条件下,丰中子碰撞系统明显加强了动量相关作用对rb(Ek)的效应.故两个碰撞系统rb(Ek)结果的比较为在重离子碰撞中提取动量相关作用的知识提供了另一个重要的判据。     

Determining bottom price-levels after a speculative peak

During a stock market peak the price of a given stock (i) jumps from an initial level p ₁(i) to a peak level p ₂(i) before falling back to a bottom level p ₃(i). The ratios A(i) = p ₂(i)/p ₁(i) and B(i)= p ₃(i)/p ₁(i) are referred to as the peak- and bottom-amplitude respectively. The paper shows that for a sample of stocks there is a linear relationship between A(i) and B(i) of the form: B=0.4A+b. In words, this means that the higher the price of a stock climbs during a bull market the better it resists during the subsequent bear market. That rule, which we call the resilience pattern, also applies to other speculative markets. It provides a useful guiding line for Monte Carlo simulations. Received 9 June 2000     

Quantum virial coefficients of molecular nitrogen

We report virial coefficients up to third order in density for molecular nitrogen, investigating 103 temperatures in the range (15 K, 3000 K). All calculations are based on an ab initio-based potential taken from the literature. Path-integral Monte Carlo (PIMC) is applied to account for nuclear quantum effects, and these results are compared to a more approximate but faster semiclassical treatment. Additionally, we examine a PIMC approach that employs semiclassical beads for the path-integral images, but find that it offers marginal advantage. A recently developed orientation sampling algorithm is used in conjunction with Mayer sampling to compute precise virial coefficients. We find that, within the precision of our calculations of the second-order coefficient (B₂), semiclassical methods are adequate for temperatures greater than 250 K, and are needed to correct classical behaviour for temperatures as high as 800 K. For the third-order coefficient (B₃), the semiclassical methods are adequate above 150 K, and are required up to the highest temperature examined (3000 K) in order to correct the classical treatment within the precision of the calculations. However, three-body contributions to the potential are much more significant than nuclear quantum effects for the evaluation of B₃.     

Nonstandard GLh(n) quantum groups and contraction of covariant q-bosonic algebras

GL_h(n) × GL_h(m)-covariant h-bosonic algebras are built by contracting the GL_q(n) × GL_q(m)-covariant q-bosonic algebras considered by the present author some years ago. Their defining relations are written in terms of the corresponding R _h-matrices. Whenever n = 2, and m = 1 or 2, it is proved by using U_h(sl(2)) Clebsch-Gordan coefficients that they can also be expressed in terms of coupled commutators in a way entirely similar to the classical case. Some U_h(sl(2)) rank-(1/2) irreducible tensor operators, recently constructed by Aizawa in terms of standard bosonic operators, are shown to provide a realization of the h-bosonic algebra corresponding to n = 2 and m = 1.     

Dependence of structure factor and correlation energy on the width of electron wires

The structure factor and correlation energy of a quantum wire of thickness b ? a _B are studied in random phase approximation (RPA) and for the less investigated region r _s < 1. Using the single-loop approximation, analytical expressions of the structure factor are obtained. The exact expressions for the exchange energy are also derived for a cylindrical and harmonic wire. The correlation energy in RPA is found to be represented by ? _c (b, r _s ) = α(r _s )/b + β(r _s ) ln(b) + η(r _s ), for small b and high densities. For a pragmatic width of the wire, the correlation energy is in agreement with the quantum Monte Carlo simulation data.     

Asymptotic properties of the solutions of a differential equation appearing in QED and QCD

Stimulated by the study described in the preceding paper, we establish the asymptotic behaviour of the ratio h′(0)/h(0) for g → ∞, where h(r) is a solution, vanishing at infinity, of the differential equation h″(r) = igω(r)h(r) on the domain 0 ≤ r ≤ ∞ and ω(r) = (1 − √rK₁(√r))/r. Some results are valid for more general ω's.     

Wei Hua’s four-parameter potential: Comments and computation of molecular constants αe and εeXe

The value of adjustable parameterC in the four-parameter potentialU(r) =D _e [(1 - exp[-b(r -r _e)])/(1 -C exp[-b(r -r _e)])]² has been expressed in terms of molecular parameters and its significance has been brought out. The potential so constructed, withC derived from the molecular parameters, has been applied to ten electronic states in addition to the states studied by Wei Hua. Average mean deviation for these 25 states has been found to be 3.47 as compared to 6.93, 6.95 and 9.72 obtained from Levine, Varshni and Morse potentials, respectively. Also Dunham’s method has been used to express rotation-vibration interaction constant (α_e) and anharmonicity constant (ω_ex_e) in terms ofC and other molecular constants. These relations have been employed to determine these quantities for 37 electronic states. For α_e, the average mean deviation is 7.2% compared to 19.7% for Lippincott’s potential which is known to be the best to predict these values. Average mean deviation for (ω_ex_e) turns out to be 17.4% which is almost the same as found from Lippincott’s potential function.     

Smoothening of a random surface as it results from the edwards-wilkinson kinetic equation

The time evolution of a random surfacez=h(r, t) (r=x, y) formed by a deposition process of the Edwards-Wilkinson type is discussed. The discussion is based on the author’s former derivation of the autocorrelation functionA _h(|r − r′|,t, t′)=〈h(r,t)h(r′,t′)〉 of the height functionh(r,t) under the assumption of a stochastic initial condition [V. Bezák: Acta Physica Univ. Comenianae39 (1998) 135]. Under the assumption of a steady (non-zero) deposition rate, the varianceσ _h ² (t)=〈[h(r,t)]²〉 increases logarithmically in time whilst the correlation lengthl _h(t) of the height functionh(r,t) increases as ∼t ^1/2. Therefore, the ratioσ _h(t)/l _h (t) tends to zero and the surfacez=h(r,t) does always tend towards a smoothened appearance. This work has been supported by the Slovak Grant Agency VEGA under contract No. 1/4319/97.     

Bound state nature of the exotic Z<Subscript>b</Subscript> states

The assumption that the newly observed charged bottomonia states Z _b(10610) and Z _b(10650) are of molecular nature is confronted with the measured invariant-mass distributions for the transitions of the ϒ(5S) to the final states h _b π ⁺ π ⁻ and h _b(2P)π ⁺ π ⁻. It is shown that the assumption that the Z _b(10610) and Z _b(10650) are B [`(B)]*\bar B* + c.c. and B <sup>*</sup> [`(B)]*\bar B* bound states, respectively, with very small binding energies is consistent with the data. The calculation is based on a power counting for bottom meson loops, which is explicitly given up to two-loop in the framework of a nonrelativistic effective-field theory. We also show that if the Z_b states are of molecular nature, then the data should not be analyzed by using a Breit-Wigner parametrization.     

Determining Liquid Structure from the Tail of the Direct Correlation Function

In important early work, Stell showed that one can determine the pair correlation function h(r) of the hard-sphere fluid for all distances r by specifying only the tail of the direct correlation function c(r) at separations greater than the hard-core diameter. We extend this idea in a very natural way to potentials with a soft repulsive core of finite extent and a weaker and longer ranged tail. We introduce a new continuous function T(r) which reduces exactly to the tail of c(r) outside the (soft) core region and show that both h(r) and c(r) depend only on the out projection of T(r): i.e., the product of the Boltzmann factor of the repulsive core potential times T(r). Standard integral equation closures can thus be reinterpreted and assessed in terms of their predictions for the tail of c(r) and simple approximations for its form suggest new closures. A new and very efficient variational method is proposed for solving the Ornstein–Zernike equation given an approximation for the tail of c. Initial applications of these ideas to the Lennard-Jones and the hard-core Yukawa fluid are discussed.