首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 260 毫秒
1.
Monte Carlo calculations are reported for the radial distribution function g 2(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 ρσ3 = 0·80, kT/ε = 0·719. It is shown that at high densities the perturbation expansion of g 2(r; λ = 1) about g 2(r; λ = 0) is rapidly convergent, but that the corresponding expansion for y 2(r; λ) = exp [βu(r; λ)] × g 2(r; λ) is not. In addition Monte Carlo estimates of the individual terms that contribute to the first-order perturbation term, (?g 2/?λ)λ=0, are presented. It is shown that these terms are individually large, but that (?g 2/?λ)λ=0 is small because there is strong cancellation between the various terms. Consequently, the calculation of (?g 2/?λ)λ=0 is highly sensitive to the approximation used to evaluate the individual terms.  相似文献   

2.
The Kirkwood-Buff statistical mechanical theory of surface tension γ for monatomic fluids is extended to molecular fluids. A rigorous expression for γ is derived in terms of the angular pair distribution function f(z 1 R 12θ1θ2) of an equilibrium fluid-fluid system (liquid-gas, liquid-liquid, or gas-gas). The Fowler approximation is applied for the liquid-gas case, and a simple expression for γ is derived in terms of the bulk liquid angular pair correlation function g(R 12θ1θ2). Thermodynamic perturbation theory for g(R 12θ1θ2) is also used to calculate γ theoretically. The theoretical results are compared with experimental values.  相似文献   

3.
Molecular dynamics (MD) calculations have been performed for the Lennard-Jones (12-6) potential function using 2048 particles. Using conventional parameters the results may be compared with those for liquid argon.

The dynamic structure factor S(k, ω) has been determined both by Fourier inversion of the intermediate scattering function F(k, t) and from the longitudinal current-current correlation function C (k, t). Particular attention was paid to the recurrence time of the system. The results for S(k, ω) by the two methods agree within 5 per cent for the whole region of small k-vectors considered. Double Fourier inversion of the van Hove function G(r, t) led to insufficiently accurate results for these small k-values. In view of the present data, the MD-results of Levesque et al. [1] for S(k, ω) have only a qualitative character. These latter data appear to contain truncation errors due to incomplete Fourier transformations.

Using a hydrodynamic assumption for F(k, t) we were able to extract the transport coefficients, the velocity of sound and the ratio of the specific heats in the limit of large wave lengths or small k. The velocity of sound was obtained by exploiting the MD generated anomalous dispersion curve of sound waves. Anomalous dispersion was found to set in for kσ ~ 0·25. A sound speed of 880 ms-1 has been determined which is in excellent agreement with experimental values for liquid argon. The total error for the MD value amounts to about 5 per cent. In contrast, the ratio of the specific heats γ and the transport coefficients D T and Γ (thermal diffusivity and sound attenuation) were determinable only with an accuracy of 15 per cent due to the need for larger extrapolations. Nevertheless, we found D T, Γ and γ in agreement with experimental values within 5-10 per cent.  相似文献   

4.
D. VIDUNA  W. R. SMITH 《Molecular physics》2013,111(17):2815-2821
An exact formula is derived relating the contact value of the solute-solvent radial distribution function for an additive binary hard-sphere (HS) mixture at infinite dilution, g 12(d 12), to the mixture equation of state (EOS) (1 denotes the solvent and 2 denotes the solute). This result can also be considered to be a consistency condition involving approximations for g 12(d 12) and for the mixture EOS. Employing three approximate HS mixture equations of state from the literature, we use our formula to derive corresponding analytical approximations for g 12(d 12) In addition, new computer simulations were performed to obtain accurate results for g 12(d 12) and for g 12(r 12) at the solute-solvent diameter ratios {1, 3, 5, 7, 10, 20} and the reduced solvent density ρ? = 0.8. We compare our results for g 12(d 12) with the simulation results and with the results of approximate analytical expressions for g 12(d 12) proposed by several workers. The results obtained from our formula in conjunction with two of the EOS expressions considered are more accurate than all previously proposed approximations, with the exception of the approximation of Matyushov and Ladanyi [1997, J. chem. Phys., 107, 5815], which is of comparable accuracy.  相似文献   

5.
A new type of light-scattering experiment, which should measure directly the triple static structure factor S (3) (k, q) of a fluid, is proposed. S (3)(k, q) is the full spatial Fourier transform of the equilibrium triplet distribution function g (3)(r 1, r 2, r 3). The experiment may also be used to study dynamic correlation functions of the form <ak (t)aq (t′)a_k_q(t″)> (where ak () is the kth spatial Fourier component of the density), thereby giving new information on mode-mode coupling. The method obtains its information from triple correlations in the arrival of scattered photons at three detectors. The detectors must be operated in the heterodyne mode (i.e. with a local oscillator); the scattering volume must be much larger than the volume over which molecular positions are correlated. Comparison is made with previous analyses of other multi-detector experiments.  相似文献   

6.
《Physics letters. A》1996,222(3):182-189
We propose an accurate approximate expression for the exact ℏ4 quantum correction to the pair distribution function g2q(r12) that we have derived recently in an OCP using Wigner-Kirkwood ℏ2 expansion. Our expression, depending only the classical pair distribution function g2c(r12), reproduces the behavior of Wigner-Kirkwood g2q(r12) at order ℏ4, at small, intermediate and large r12.  相似文献   

7.

Within the model of stable random matrices possessing translational invariance, a two-dimensional (on a square lattice) disordered oscillatory system with random strongly fluctuating bonds is considered. By a numerical analysis of the dynamic structure factor S(q, ω), it is shown that vibrations with frequencies below the Ioffe-Regel frequency ωIR are ordinary phonons with a linear dispersion law ω(q) ∝ q and a reciprocal lifetime б ~ q3. Vibrations with frequencies above ωIR, although being delocalized, cannot be described by plane waves with a definite dispersion law ω(q). They are characterized by a diffusion structure factor with a reciprocal lifetime б ~ q2, which is typical of a diffusion process. In the literature, they are often referred to as diffusons. It is shown that, as in the three-dimensional model, the boson peak at the frequency ωb in the reduced density of vibrational states g(ω)/ω is on the order of the frequency ωIR. It is located in the transition region between phonons and diffusons and is proportional to the Young’s modulus of the lattice, ω b E.

  相似文献   

8.
Details are given of the analytic features of this amplitudef g(l, k) in the g plane for a very large class of potentials that satisfy conditions (1). It is shown that there is a region (l, k2) around the point g = 0 that is free from singularities inf g(l, k), so the Mittag-Leffler method can be applied to findf g(l, k) and hence also the total amplitude Tg(k, t) for any g to any specified degree of accuracy with reference to the information contained in the coefficients of a finite number of terms of the series given by perturbation methods forf g(l, k).This paper was read in substance on 8th April 1963 at Uzhgood University [1] and on 29 June 1963 at the All-Union Symposium on the Theory of Elementary Particles at Novosibirsk.  相似文献   

9.
Perturbation theory is used to consider expansions for the radial distribution function, g 2(r), of a fluid with a soft core. We consider the Lennard-Jones (12, 6) potential and divide it into repulsive and attractive regions. In the repulsive region we expand the function exp (β u(r))g 2(r) about a hard sphere value. For the first-order contribution of the attractive region we consider a simple approximation to the exact analytical expression. The resulting g 2(r) is accurate at densities below about ρσ 3=0·5.  相似文献   

10.
The emission spectrum of BN has been investigated in the 1800–9000 cm−1region using a Fourier transform spectrometer. BN was formed in a microwave discharge of He with a trace of BCl3and N2. The bands observed in the 3000–7800 cm−1interval have been assigned as theb1Π–a1Σ+transition, with the 0–0 band at 3513.99040(43) cm−1. This transition is analogous to theA1Πu–X1Σ+g(Phillips) system of the isoelectronic C2molecule. The rotational analysis of the 0–0, 1–1, 1–0, 2–1, 3–2, 2–0, 3–1, 4–2, and 4–1 bands has been obtained and the molecular constants for theb1Π anda1Σ+states have been determined. A local perturbation has been observed in thev= 1 vibrational level of theb1Π state nearJ= 18 caused by the interaction with thev= 3 vibrational level of thea1Σ+state. The principal equilibrium constants for thea1Σ+state are: ωe= 1705.4032(11) cm−1, ωexe= 10.55338(52) cm−1,Be= 1.683771(10), αe= 0.013857(16) cm−1, andre= 1.2745081(37) Å. Although theb1Π–a1Σ+transition has recently been seen in emission from boron nitride trapped in solid neon matrices [J. Chem. Phys.104,3143–3146 (1996)], our work represents the first observation of this transition of BN in the gas phase.  相似文献   

11.
The spin‐polarized homogeneous electron gas with densities ρ and ρ for electrons with spin ‘up’ (↑) and spin ‘down’ (↓), respectively, is systematically analyzed with respect to its lowest‐order reduced densities and density matrices and their mutual relations. The three 2‐body reduced density matrices γ↑↑, γ↓↓, γa are 4‐point functions for electron pairs with spins ↑↑, ↓↓, and antiparallel, respectively. From them, three functions G↑↑(x,y), G↓↓(x,y), Ga(x,y), depending on only two variables, are derived. These functions contain not only the pair densities according to g↑↑(r) = G↑uarr;(0,r), g↓↓(r) = G↓↓(0,r), ga(r) = Ga(0,r) with r = | r 1 ‐ r 2|, but also the 1‐body reduced density matrices γ and γ being 2‐point functions according to γs = ρsfs and fs(r) = Gss(r, ∞) with s = ↑,↓ and r = | r 1 ‐ r 1|. The contraction properties of the 2‐body reduced density matrices lead to three sum rules to be obeyed by the three key functions Gss, Ga. These contraction sum rules contain corresponding normalization sum rules as special cases. The momentum distributions n(k) and n(k), following from f(r) and f(r) by Fourier transform, are correctly normalized through fs(0) = 1. In addition to the non‐negativity conditions ns(k),gss(r),ga(r) ≥ 0 [these quantities are probabilities], it holds ns(k) ≤ 1 and gss(0) = 0 due to the Pauli principle and ga(0) ≤ 1 due to the Coulomb repulsion. Recent parametrizations of the pair densities of the spin‐unpolarized homogeneous electron gas in terms of 2‐body wave functions (geminals) and corresponding occupancies are generalized (i) to the spin‐polarized case and (ii) to the 2‐body reduced density matrix giving thus its spectral resolutions.  相似文献   

12.
The g factors g // and g of the ground Γ6(2 T 2) and excited Γ4,5(2 E), Γ6(2 E) states for trigonal Cu2+ centres in ZnO crystals are calculated from three theoretical methods, the complete diagonalization (of the energy matrix) method, the second-order perturbation method (PTM-I) and the simplified second-order perturbation method (PTM-II, this method was described in an earlier paper). These methods are based on the cluster approach in which the spin-orbit coupling parameters ζ, ζ′ and the orbital reduction factors k, k′ are calculated from a semi-empirical molecular orbital method. The crystal-field parameters used in the calculations are obtained from the superposition model and so the defect structure of Cu2+ centres in ZnO can be acquired. The calculated g factors from the three methods are in reasonable agreement with the experimental values and the defect structure of Cu2+ centres in ZnO is acquired. It appears that in some cases the approximate PTM can be applied in the studies of g factors of various states. The conditions that the PTM are ineffective are discussed.  相似文献   

13.
In the frame work of TREDER 's gravitational theory we consider two classes of field equations which are derivable from two classes of LAGRANGE ian densities Ω(1)1, ω2), Ω(2)(s?1, s?2). ω1, ω2; s?1, s?2 are parameters. Ω(2)1, ω2) gives us field equations which are up to the post-NEWTON ian approximation in the sense of NORDTVEDT , THORNE and WILL equivalent to the field equations given by BRANS and DICKE . For ω2 = ?1 ?2ω1 field equations follow from Ω(1)1, ?1 ?2ω1) which are in the above mentioned sense of post-NEWTON ian approximation equivalent to EINSTEIN 's equations. The field equations following from Ω(1)1, ω2) have a cosmological model with the well known cosmological singularities for T → ± ∞ in case that ω1/(1 +3ω12) ? γ > 0. For ω1/(1 +3ω12) ≤ 0 cosmological models with no cosmological singularities exist. From Ω(2)(s?1, s?2) we obtain field equations which at the best give us perihelion rotation 7% above EINSTEIN 's value and light deflection 7% below the corresponding EINSTEIN 's value. But in that case we are able to show the existence of a cosmological model without any cosmological singularity.  相似文献   

14.
A mathematical method is presented for solving the Schr?dinger equation for a system of identical body forces. The N-body forces are more easily introduced and treated within the hyperspherical harmonics. The problem of the N-body potential has been used at the level of both classical and quantum mechanics. The hypercentral interacting potential is assumed to depend on the hyperradius x = (ξ12 + ξ22 + ⋯ + ξN−12)1/2 only, where ξ12,…,ξN−1 are Jacobi relative coordinates which are functions of N-particle relative positions r12,r23,…,rN1. The problem of the harmonic oscillator and the Coulomb-type potential has been widely studied in different contexts. Using the N-body potential V(x) = ax2 + bx − (c/x) as an example, and assuming an ansatz for the eigenfunction, an exact analytical solution of the Schr?dinger equation for an N-body system in three dimensions is obtained. This method is also applicable to some other types of potentials for N-identical interacting particles.  相似文献   

15.
Single crystals of nickel-doped lithium potassium sulphate were grown by slow evaporation method at room temperature. From the nature and position of the bands observed, a successful interpretation of all the bands could be made assumingO h symmetry for the Ni2+ ion in the crystal. The bands have been assigned transitions from the ground3A2g(F) state to the excited3T2g(F),1Eg(D),3T1g(F),1T2g(D) and3T1g(P) states. The crystal field parameters derived areDq=910cm–1,B=890cm–1 andC=3560cm–1.The authors wish to express their thanks to Prof. K. Sreerama Murthy for his constant encouragement throughout this investigation. The authors are also thankful to Prof. Mihir Chowdhury, Indian Association for the cultivation of Science, Calcutta for giving permission to take the spectra.  相似文献   

16.
Using the u-plane integral of Moore and Witten, we derive a simple expression for the Donaldson invariants of product ruled surfaces Σ g ×S 2, where Σ g is a Riemann surface of genus g. This expression generalizes a theorem of Morgan and Szabó for g=1 to any genus g. We give two applications of our results: (1) We derive Thaddeus' formulae for the intersection pairings on the moduli space of rank two stable bundles over a Riemann surface. (2) We derive the eigenvalue spectrum of the Fukaya–Floer cohomology of Σ g ×S 1. Received: 22 July 1999 / Accepted: 12 June 2000  相似文献   

17.
Kinetic processes taking place in the atomic-molecular system O-O2-O3 in the middle atmosphere with the participation of oxygen molecules in the excited electronic states O2(a 1Δg) and O2(b 1Σ g + ) are analyzed in detail. The possibility of increased ozone production under the influence of solar radiation during the laser excitation of O2 molecules in the a 1Δg state is demonstrated on the basis of numerical modeling. Upper and lower bounds are determined for the densities of O2(a 1Δg) molecules at which the ozone concentration increases in the irradiated zone. Zh. Tekh. Fiz. 68, 15–23 (August 1998)  相似文献   

18.
Raman spectra of liquid cyclohexane, C6H12, and deuterated cyclohexane, C6D12, were recorded with both parallel and perpendicular polarizations. The observed vibrational wavenumbers, depolarization ratios, and their intensities were measured and compared with the corresponding predicted values as well as the experimental values previously reported. The conformational energetics were obtained with the Møller–Plesset perturbation method to the second order [MP2(full)] as well as with density functional theory by the B3LYP method utilizing a variety of basis sets. The average ab initio predicted difference in energy between the more stable chair form (D3d) and the less stable twisted‐boat form (D2) is 2213 cm−1 (26.47 kJ/mol), with a similar value of 2223 cm−1 (26.59 kJ/mol) from the density function theory calculations. By using two dihedral angles as variables, we calculated the chair–boat interconversion pathway for cyclohexane at the MP2(full)/6‐31G(d) level. The harmonic force constants, Raman intensities, depolarization values, and the potential energy distribution were predicted from both MP2(full) and B3LYP calculations with the 6‐31G(d) basis set and compared with the experimental values for the chair form when available. The ‘adjusted’ r0 structural parameters were obtained from MP2/6‐311 + G(d,p) calculations and previously reported microwave rotational constants of five isotopomers of cyclohexane: i.e. 1,1‐d2, 13C‐1,1‐d2, 1,1,2,2,3,3‐d6, and d1 (equatorial and axial). The determined distances in Å are: r(CC) = 1.536(3), r(CH)ax = 1.098(1); r(CH)eq = 1.095(1); and the angles in degrees: ∠CCHax = 108.8(3); ∠CCHeq = 110.2(3); ∠CCC = 111.1(3); and ∠HCH = 107.6(3) with dihedral angle ∠CCCC = 55.7(3). These values are compared with those previously reported and it is found that the difference in the r0 distances (0.003 Å) between the two CH values is much smaller than the difference (0.008 Å) previously reported for the rs values. Copyright © 2008 John Wiley & Sons, Ltd.  相似文献   

19.
For a two-level system corresponding to a particle of spin ½ in a random field in the Z direction, the relaxation function

has been estimated, the magnitude Ω(t) being the sum of the isotropic interactions of the particle in question with particles j executing diffusional motion. Specifically Ω(t)=ω0 + Σ ω(Mj, rj), where ω0 = constant, Mj is a random time-independent parameter, ω(M, r) decreases with r faster than r -3 and r j = r j (t) is a diffusional, random function. From the expression for <σ+(t)>, we establish general features of the relaxation phenomenon for diffusional processes, and calculate the relaxation rate 1/T 2 and relaxation shift Δω to be 1/T 2-iΔω = 4π CAvM λ M , where C is the concentration of particles and λ M is the scattering length for an equation of the Schrödinger type with an imaginary potential -iω(M, r) instead of U/?, and diffusion coefficient instead of ?/2m. We also found that for the case of ‘external’ relaxation, the Redfield approach proved valid only under the simultaneous restrictions of low concentration and weak interaction.  相似文献   

20.
We present a direct method which allows an accurate and—at the same time—economical calculation of the pair distribution functions (PDFs) gij (r) of an additive ternary hard-sphere system within the Percus—Yevick approximation. The approach is based on the fact that for this approximation the Laplace transforms [?ij (s)] of the PDFs are known analytically, so that the inversion of the ?ij (s) into r space can be performed exactly. The expressions presented here allow the determination of the ?ij (r) for r values up to 8R 1, R 1 being the diameter of the smallest species; this range in r space should be sufficient for applications in standard algorithms of liquid state theory, such as thermodynamic perturbation theories or integral-equation approaches.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号