Logarithmic finite-size corrections in the three-dimensional mean spherical model

The finite-size scaling prediction about logarithmic corrections in the free energy arising from corners in the geometry of the system is tested on the three-dimensional mean spherical model. The general case of boundary conditions which are periodic ind 0 dimensions and free or fixed in the remaining 3 –d dimensions is considered. Logarithmic and double-logarithmic size corrections stemming from corners, edges, and surfaces are obtained.     

Finite-size scaling for the mean spherical model with inverse power law interaction

A new analytical technique based on integral transformations with Mittag-Leffler-type kernels is used to derive the finite-size scaling function for the free energy per particle of the mean spherical model with inverse power law asymptotics of the interaction potential. The asymptotic formation of the singularities in the specific heat and magnetic susceptibility at the bulk critical point is studied.     

On the finite-size scalling equation for the spherical model

The mean spherical model with an arbitrary interaction potential, the Fourier transform of which has a long-wavelength exponent , 0<2, is considered under periodic boundary conditions and fully finite geometry ind dimensions, when <d<2. A new form of the finite-size scaling equation for the spherical field in the critical region is derived, which relates the temperature shift to Madelung-type lattice constants. The method of derivation makes use of the Poisson summation formula and a Laplace transformation of the momentumspace correlation function.On leave of absence from Institute of Mechanics and Biomechanics, Bulgarian Academy of Sciences, 1113 Sofia, Bulgaria.     

Finite-size effects at first-order transitions

Finite-size rounding of a first-order phase transition is studied in “block”- and “cylinder”-shaped ferromagnetic scalar spin systems. Crossover in shape is investigated and the universal form of the rounded susceptibility peak is obtained. Scaling forms on the low-temperature side of the critical point are considered both above and below the borderline dimensionality,d _>=4. A method of phenomenological renormalization, applicable to both odd and even field derivatives, is suggested and used to estimate universal amplitudes for two-dimensional Ising models atT=T_c.     

Finite-size dependence of the helicity modulus within the mean spherical model

The validity of the finite-size scaling prediction about the existence of logarithmic corrections in the helicity modulus of three-dimensional O(n)-symmetric order parameter systems in confined geometries is studied for the three-dimensional mean spherical model of geometryL ^3/s-d×^d, 0d<3. For a fully finite geometry the general case ofd _p0 periodic,d _a0 antiperiodic,d ₀0 free, andd ₁0 fixed (d _p+d_a+d₀+d₁=d, d=3) boundary conditions is considered, whereas for film (d=2) and cylinder (d=1) geometries only the case of antiperiodic and/or periodic boundary conditions is investigated. The corresponding expressions for the finite-size scaling function of the helicity modulus and its asymptotics in the vicinity, below, and above the bulk critical temperatureT _c and the shifted critical temperatureT _c,L are derived. The obtained results are not in agreement with the hypothesis of the existence of a log(L) correction term to the finite-size behavior of the helicity modulus in the finite-size critical region if d=3. In the case of film and cylinder geometries there are no logarithmic corrections. In the case of a fully finite geometry a universal logarithmic correction term –[(d ₀ –d ₁)/4–2^da–1/²] lnL/L is obtained only for (T _c-T) LlnL.     

Interface formation and a structural phase transition for the spherical model of ferromagnetism

A detailed analysis is reported examining the local magnetic susceptibility (r), in relation to the correlation functionG(R) and correlation length , of a spherical model ferromagnet confined to geometry =L ^{d –d} × ^d ( ^d 2,d>2) under a continuous set oftwisted boundary conditions. The twist parameter in this problem may be interpreted as a measure of the geometry-dependent doping level of interfacial impurities (or antiferromagnetic seams) in theextended system at various temperatures. For _j 0, jd-d, no seams are present except at infinity, whereas if _j = 1/2, impurity saturation occurs. For 0 < _j < 1/2 the physical domain _phys =D ^{d –d} × ^d (D>L), defining the region between seams containing the origin, depends on temperature above a certain threshold (T>T ₀). Below that temperature (T>T ₀), seams are frozen at the same position (DL/2,d-d'=1), revealing a smoothly varying largescale structural phase transition.     

Sine-Gordon form factors in finite volume

We compare form factors in sine-Gordon theory, obtained via the bootstrap, to finite volume matrix elements computed using the truncated conformal space approach. For breather form factors, this is essentially a straightforward application of a previously developed formalism that describes the volume dependence of operator matrix elements up to corrections exponentially decaying with the volume. In the case of solitons, it is necessary to generalize the formalism to include effects of non-diagonal scattering. In some cases it is also necessary to take into account some of the exponential corrections (so-called μ-terms) to get agreement with the numerical data. For almost all matrix elements the comparison is a success, with the puzzling exception of some breather matrix elements that contain disconnected pieces. We also give a short discussion of the implications of the observed behavior of μ-terms on the determination of operator matrix elements from finite volume data, as occurs e.g. in the context of lattice field theory.     

The spectrum of a one-dimensional hierarchical model

The spectrum of a discrete Schrödinger operator with a hierarchically distributed potential is studied both by a renormalization group technique and by numerical analysis. A suitable choice of the potential makes it possible to reduce the original problem to a two-dimensional map. Scaling laws for the band-edge energyE _be and for the integrated density of states are predicted together with the global properties of the spectrum. Different scaling regimes are obtained depending on a hierarchy positive parameterR: for R<1/2 the usual scaling laws for the periodic case are obtained, while forR>1/2 the scaling behavior depends explicitly onR.     

The Yang-Lee edge singularity in spherical models

The density of Yang-Lee zeros in the thermodynamic limit is discussed for ferromagnetic spherical models of general dimensionalities and arbitrary range of interaction. In all cases the zeros lie on the imaginary axis in the complex magnetic field planeH=H+iH with a density (H) that exhibits a square root singularity(H) (H-H ₀), with=1/2, as the edge of the gap atH=H ₀(T) is approached forT>T _c. WhenTT _c one hasH ₀(T)(TT _c ) with critical exponent=+.Supported by the National Science Foundation in part through the Materials Science Center at Cornell University.     

An intermediate spherical model of a ferromagnet

A one-parameter family of partition functions is considered which for zero value of the parameter reduces to the spherical model of a ferromagnet. The model for > 0 is closer to the usual discrete lattice spin model of a ferromagnet than is the spherical model. The first four terms in of the limiting value of the partition function are calculated above and below the critical temperature for arbitrary interactions using the saddle point method to calculate certain correlation functions for the spherical model. These calculations indicate that the critical temperature is independent of for small and certain interactions.Part of this research appeared in the author's doctoral thesis.⁽³⁾     

A new fully nonisothermal coupled model for the simulation of ice sheet flow

A shallow ice thermocoupled model for the complex nonlinear polythermal ice sheet dynamics is proposed and solved by means of efficient numerical methods. A novelty is the obstacle problem formulation associated to a nonlinear integro-differential equation (with nonlocal temperature dependent coefficients) for the ice sheet profile. This formulation is motivated by the free boundary feature and the influence of the temperature on the profile (fully nonisothermal model). Concerning the temperature equation, a dynamically prescribed surface temperature, obtained from an Energy Balance model corrected by the altitude effect, is posed. As the profile and temperature equations are fully coupled, a nonlinear PDE system governing the upper ice sheet profile, the velocity field, the temperature and the basal stress is stated. In addition to the numerical difficulties associated to the new profile equation, several techniques have been considered for the numerical solution of the temperature, velocity and basal magnitudes. Discussions concerning the nonlinear dynamics of the different involved magnitudes and the improvement in their computed values with respect to previous works are also presented.     

Quantum spherical model with competing interactions

We analyse the phase diagram of a quantum mean spherical model in terms of the temperature T$T$, a quantum parameter g$g$, and the ratio p=−_J2/_J1$p=-J_2/J_1$, where _J1>0$J_1>0$ refers to ferromagnetic interactions between first-neighbour sites along the d$d$ directions of a hypercubic lattice, and _J2<0$J_2<0$ is associated with competing antiferromagnetic interactions between second neighbours along m≤d$m\le d$ directions. We regain a number of known results for the classical version of this model, including the topology of the critical line in the g=0$g=0$ space, with a Lifshitz point at p=1/4$p=1/4$, for d>2$d>2$, and closed-form expressions for the decay of the pair correlations in one dimension. In the T=0$T=0$ phase diagram, there is a critical border, _gc=_gc(p)$g_c=g_c\left(p\right)$ for d≥2$d\ge 2$, with a singularity at the Lifshitz point if d<(m+4)/2$d<(m+4)/2$. We also establish upper and lower critical dimensions, and analyse the quantum critical behavior in the neighborhood of p=1/4$p=1/4$.     

Trap-size scaling of finite Bose systems within an exact canonical ensemble

Based on an exact canonical partition function, we investigate the trap-size scaling for ideal Bose gases with a finite number of particles N confined in a cubic box or in a harmonic trap. We study the trap-size scaling behaviors of condensate fraction 〈n₀〉/N and specific heat C_N around the transition temperature T_c (i.e., t = T/T_c − 1 → 0) for the three different traps, where a trap exponent θ in dependence of the trapping potential and the universality class of transition are introduced. In the box trap with periodic and Dirichlet boundary conditions, where θ → 1, we find that the scaling functions governing the various critical behaviors are universal but respective of the boundary conditions. The calculated critical exponents are in nice agreement with analytical scaling predictions. The borders of universality validity are obtained numerically. In the case of the harmonic trap, the critical behavior of the system is also found to be universal, and the trap exponent is obtained as θ ? 0.068.     

球面旋涂光刻胶膜厚分布的数学模型

提出了球面旋涂微米级厚度光刻胶膜层薄化率公式及径向位置演变公式,并得到了膜厚分布的演变公式。与平面涂胶相比,球面涂胶离心力及重力分量是在不断的变化。根据平面旋涂运动方程及球面面形特征,给出了球面旋涂运动方程;结合流体层流的表面条件及不可压缩流体的质量连续方程,推导出了膜厚h及径向位置r对时间t的演变公式,并得到了在径向位置r处初始厚度为h0的膜厚演变的数学模型。通过对模型参数的分析可知,球面旋涂光刻胶应采用主从轴偏心旋涂,旋涂时工件的开口应朝向侧面(旋转轴水平)。     

Bethe ansatz calculations for the eight-vertex model on a finite strip

Bethe ansatz equations for the eigenvalues of the transfer matrix of the eight-vertex model are solved numerically to yield mass gap data on infinitely long strips of up to 512 sites in width. The finite-size corrections, at criticality, to the free energy per site and polarization gap are found to be in agreement with recent studies of theXXZ spin chain. The leading corrections to the finite-size scaling estimates of the critical line and thermal exponent are also found, providing an explanation of the poor convergence seen in earlier studies. Away from criticality, the linear scaling fields are derived exactly in the full parameter space of the spin system, allowing a thorough test of a recently proposed method of extracting linear scaling fields and related exponents from finite lattice data.     

球几何中子输运保正线性间断有限元格式

针对球几何中子输运方程线性间断有限元方法计算的负中子通量问题,构造了保正线性间断有限元格式,该格式保持中子角通量0阶矩和1阶矩。现有方法计算中子角通量非负时,采用传统的线性间断有限元方法,求解线性方程组;原方法计算出现负通量,则采用构造的保正格式,求解非线性方程组。编制了球几何中子输运问题保正格式程序模块,并集成到应用程序。数值算例表明构造的保正格式计算的中子通量非负,有效降低数值误差,提高数值计算的精度。     

Dimensional crossover and finite-size scaling belowT c

Using the formalism developed in earlier work, dimensional crossover on ad-dimensional layered Ising-type system satisfying periodic boundary conditions and of sizeL is considered belowT _c (L), T _c (L) being the critical temperature of the finite-size system. Effective critical exponents _eff and _eff are shown explicitly to crossover between theird- and (d–1)-dimensional values for _L in the limitsL/ _L andL/ _L 0, respectively, _L , being the correlation length in the layers. Using anL-dependent renormalization group, the effective exponents are shown to satisfy natural generalizations of the standard scaling laws. In addition,L-dependent global scaling fields which span the entire crossover are defined and a scaling form of the equation of state in terms of them derived. All the above assertions are verified explicitly to one loop in perturbation theory, in particular effective exponents and a universal crossover equation of state are obtained and shown in the above asymptotic limits to be in good agreement with known results.     

Finite-size scaling and the renormalization group

Renormalization group calculations ind = 4 andd = 4 – are performed for a system of finite size. A form of mean-field theory is used which yields a rounded transition for a finite system, and this allows a sensible expansion in fluctuations. A combination of Ewald and Poisson sum techniques is used to produce explicit numerical results for the specific heat ind = 4 which, with the setting of two nonuniversal metrical factors and the fourth-order coupling constant may be compared with simulations. The numerical visibility of logarithmic corrections is investigated. The universal scaling function for the specific heat to relativeO() is also evaluated numerically.     

The finite-size scaling study of the specific heat and the Binder parameter for the six-dimensional Ising model

The six-dimensional Ising model with nearest-neighbor pair interactions is simulated on the Creutz cellular automaton by using the finite-size lattices with the linear dimensions L=4,6,8,10. The temperature variations and the finite-size scaling plots of the specific heat and Binder parameter verify the theoretically predicted expression near the infinite lattice critical temperature. The approximate values for the critical temperature of the infinite lattice T_C=10.838(1), T_C=10.836(20) and T_C=10.835(1) are obtained from the intersection points of specific heat curves, Binder parameter curves and the straight line fit of specific heat maxima, respectively. These results are in agreement with the more precise value of T_C=10.835(5). The value obtained for the critical exponent of the specific heat, i.e., =0.012(2) is also in agreement with =0 predicted by the theory.