格点QCD中有限体积下形状因子的内插计算

在格点QCD的框架下,讨论了计算强子的形状因子的过程中,有限体积效应对所得结果的影响,并且给出了计算连续动量空间上形状因子数值的内插算法.本文以π介子的形状因子为例,介绍了从格点QCD中得到的三点关联函数出发计算π介子的形状因子的方法,并说明有限体积效应使得计算结果只能在分立动量上给出.本文探讨了一种保持旋转对称性的内...     

Numerical investigations of discrete scale invariance in fractals and multifractal measures

Fractals and multifractals and their associated scaling laws provide a quantification of the complexity of a variety of scale invariant complex systems. Here, we focus on lattice multifractals which exhibit complex exponents associated with observable log-periodicity. We perform detailed numerical analyses of lattice multifractals and explain the origin of three different scaling regions found in the moments. A novel numerical approach is proposed to extract the log-frequencies. In the non-lattice case, there is no visible log-periodicity, i.e., no preferred scaling ratio since the set of complex exponents spreads irregularly within the complex plane. A non-lattice multifractal can be approximated by a sequence of lattice multifractals so that the sets of complex exponents of the lattice sequence converge to the set of complex exponents of the non-lattice one. An algorithm for the construction of the lattice sequence is proposed explicitly.     

On the extendedness of eigenstates in a hierarchical lattice: A critical view

We take a critical view at the basic definition of extended single particle states in a non-translationally invariant system. For this, we present the case of a hierarchical lattice and incorporate long range interactions that are also distributed in a hierarchical fashion. We show that it is possible to explicitly construct eigenstates with constant amplitudes (normalized to unity) at every lattice point for special values of the electron-energy. However, the end-to-end transmission, corresponding to the above energy of the electron in such a hierarchical system depends strongly on a special correlation between the numerical values of the parameters of the Hamiltonian. Keeping the energy and the distribution of the amplitudes invariant, one can transform the lattice from conducting to insulating simply by tuning the numerical values of the long range interaction. The values of these interactions themselves display a fractal character.     

Asymptotics of the partition function of a general Markov random field on an infinite rectangular lattice

We investigate theoretically and numerically the asymptotics of the partition function of a general Markov random field (MRF) on an infinite rectangular lattice. We first propose the general local energy function (LEF)-parameterized MRF. Then we prove that the thermodynamic limit of the free energy of the MRF can be exactly characterized by the Perron root of the fundamental transfer matrix of a particular Markov additive process (MAP). This matrix possesses a special structure and many interesting properties that enable parallel computation of the Perron root and may be beneficial for deriving an analytical form of the free energy. We also develop another transfer matrix for numerical computation of the desired Perron root. Specifically, the former is a site-to-site transfer matrix on a twisted cylindrical lattice, while the latter is the one associated with a row-to-row transition on a vertical strip. Numerical results show that our methods exhibit consistent finite-size scaling behavior even for small values of the lattice width. This study reveals that the fundamental transfer matrix is an alternative direction of research on the analysis of the partition function of general MRFs within the scope of matrix algebra.     

From regular to growing small-world networks

We propose a growing model which interpolates between one-dimensional regular lattice and small-world networks. The model undergoes an interesting phase transition from large to small worlds. We investigate the structural properties by both theoretical predictions and numerical simulations. Our growing model is a complementarity for the important static Watts-Strogatz network model.     

Supercontinuum generation in a two-dimensional photonic kagome crystal

We demonstrate the generation of ultrabroad spectra in a photonic crystal fiber with a kagome-lattice transverse structure. This two-dimensional periodic photonic lattice allows for strong confinement of light without employing defect states nor using photonic bandgap guiding. Light guiding is mediated by total internal reflection in the intersections of the lattice structure, similar to tapered or micro-structured fibers. The kagome lattice structure is manufactured from a soft glass with a high nonlinearity. Using a Ti:sapphire oscillator as a pump source, we observe for the first time impressive supercontinuum generation in the guided modes of a 2D photonic lattice. Supercontinuum generation is caused by fission and radiation of higher-order solitons in the anomalous dispersion range. Our spectrum encompasses the spectral range from 200 to 1750 nm. The dependence of the continuum on coupling spot location, fiber length, and pump wavelength and power as well as on pulse duration and polarization state is investigated. Using a numerical simulation for the lattice structure, pulse propagation through this structure is theoretically studied. Our model reveals the mechanism of supercontinuuum generation in the 2D photonic structure and explains the essential experimental findings.Electronic supplementary material to this article is avaiable at and accessible for authorized users.     

Implementation of diffuse reflection boundary conditions in a thermal lattice Boltzmann model with flux limiters

We discuss three new implementation versions of diffuse reflection boundary conditions in a thermal lattice Boltzmann model. Their accuracy is investigated in the case of Couette flow by considering the slip regime. The best results are recovered with versions 2 and 3, which rely on outgoing fluxes to express the particle distribution functions in the ghost nodes outside the flow domain. Version 2 is found to be more economical since it involves no interpolation procedure. This version was thereafter used to investigate the temperature profile in Couette flow for various values of Prandtl number, as well as the capability of the thermal LB model to capture the Knudsen minimum in Poiseuille flow.     

Two-Dimensional Discrete Gap Breathers in a Two-Dimensional Diatomic β Fermi--Pasta--Ulam Lattice

We study the existence of two-dimensional discrete breathers in a two-dimensional face-centred square lattice consisting of alternating light and heavy atoms, with nearest-neighbour coupling containing quartic soft or hardnonlinearity. This study is focused on two-dimensional breathers with frequency in the gap that separates the acoustic and optical bands of the phonon spectrum. We demonstrate the possibility of existence of two-dimensional gap breathers by using the numerical method, the local anharmonicity approximation and the rotating wave approximation. We obtain six types of two-dimensional gap breathers, i.e., symmetric, mirror-symmetric and asymmetric, no matter whether the centre of the breather is on a light or a heavy atom.     

Lattice Boltzmann method for three-dimensional moving particles in a Newtonian fluid

A lattice Boltzmann method is developed to simulate three-dimensional solid particle motions in fluids. In the present model, a uniform grid is used and the exact spatial location of the physical boundary of the suspended particles is determined using an interpolation scheme. The numerical accuracy and efficiency of the proposed lattice Boltzmann method is demonstrated by simulating the sedimentation of a single sphere in a square cylinder. Highly accurate simulation results can be achieved with few meshes, compared with the previous lattice Boltzmann methods. The present method is expected to find applications on the flow systems with moving boundaries, such as the blood flow in distensible vessels, the particle－flow interaction and the solidification of alloys.     

Interaction versus dimerization in one-dimensional Fermi systems

In order to study the effect of interaction and lattice distortion on quantum coherence in one-dimensional Fermi systems, we calculate the ground state energy and the phase sensitivity of a ring of interacting spinless fermions on a dimerized lattice. Our numerical DMRG studies, in which we keep up to 1000 states for systems of about 100 sites, are supplemented by analytical considerations using bosonization techniques. We find a delocalized phase for an attractive interaction, which differs from that obtained for random lattice distortions. The extension of this delocalized phase depends strongly on the dimerization induced modification of the interaction. Taking into account the harmonic lattice energy, we find a dimerized ground state for a repulsive interaction only. The dimerization is suppressed at half filling, when the correlation gap becomes large. Received: 11 February 1998 / Revised: 1st April 1998 / Accepted: 30 April 1998     

Temperature and spatial diffusion of atoms cooled in a 3D lin ⊥lin bright optical lattice

We present a detailed experimental study of a three-dimensional lin⊥lin bright optical lattice. Measurements of the atomic temperature and spatial diffusion coefficients are reported for different angles between the lattice beams, i.e. for different lattice constants. The experimental findings are interpreted with the help of numerical simulations. In particular we show, both experimentally and theoretically, that the temperature is independent of the lattice constant. Received 5 July 2001 and Received in final form 13 August 2001     

Advances in multi-domain lattice Boltzmann grid refinement

Grid refinement has been addressed by different authors in the lattice Boltzmann method community. The information communication and reconstruction on grid transitions is of crucial importance from the accuracy and numerical stability point of view. While a decimation is performed when going from the fine to the coarse grid, a reconstruction must performed to pass form the coarse to the fine grid. In this context, we introduce a decimation technique for the copy from the fine to the coarse grid based on a filtering operation. We show this operation to be extremely important, because a simple copy of the information is not sufficient to guarantee the stability of the numerical scheme at high Reynolds numbers. Then we demonstrate that to reconstruct the information, a local cubic interpolation scheme is mandatory in order to get a precision compatible with the order of accuracy of the lattice Boltzmann method.These two fundamental extra-steps are validated on two classical 2D benchmarks, the 2D circular cylinder and the 2D dipole–wall collision. The latter is especially challenging from the numerical point of view since we allow strong gradients to cross the refinement interfaces at a relatively high Reynolds number of 5000. A very good agreement is found between the single grid and the refined grid cases.The proposed grid refinement strategy has been implemented in the parallel open-source library Palabos.     

Quasicontinuum approximation and iterative method for envelope solitons in anharmonic chains

For solitary waves on a monoatomic chain with nearest neighbor interactions the continuum approximation has a limited validity range and exhibits certein mathematical problems. For pulse solitons these problems are overcome by the Quasicontinuum Approach (QCA), and the validity range is considerably extended. We generalize the QCA to oscillatory excitations and derive analytic expressions for bright and dark envelope solitons, limiting ourselves to a polynomial interaction potential with harmonic, cubic and quartic terms. Moreover we describe and apply a numerical iteration procedure in Fourier space in order to take into account discreteness effects in a systematic way. This procedure yields envelope solitons with a width in the order of the lattice constant. In the case of zero velocity these solutions can be compared with intrinsic localized modes derived by other authors. The stability and accuracy of all our solutions are tested by numerical simulations.     

Non-perturbative renormalization-group approach to lattice models

The non-perturbative renormalization-group approach is extended to lattice models, considering as an example a φ⁴ theory defined on a d-dimensional hypercubic lattice. Within a simple approximation for the effective action, we solve the flow equations and obtain the renormalized dispersion epsilon(q) over the whole Brillouin zone of the reciprocal lattice. In the long-distance limit, where the lattice does not matter any more, we reproduce the usual flow equations of the continuum model. We show how the numerical solution of the flow equations can be simplified by expanding the dispersion in a finite number of circular harmonics.     

二维保单调保守恒插值算子

基于一个一维保单调保守恒插值算子,利用不完全双二次插值提出一个二维保单调保守恒插值算子.从插值逼近角度,通过几个数值实验验证该插值算子有效.用得到的二维插值算子作为结构网格自适应加密(structured adaptive mesh refinement,SAMR)算法中的细化插值算子,求解几个二维Euler方程数值例子,结果表明,提出的二维插值算子有效.     

The self-consistent diagram approximation for lattice systems

We apply the self-consistent diagram approximation to calculate equilibrium properties of lattice systems. The free energy of the system is represented by a diagram expansion in Mayer-like functions with averaging over states of a reference system. The latter is defined by one-particle mean potentials, which are calculated using the variational condition formulated. As an example, numerical computations for a two-dimensional lattice gas on a square lattice with attractive interaction between nearest neighbours were carried out. The critical temperature, the phase coexistence curve, the chemical potential and particle and vacancy distribution functions coincide within a few per cent with exact or with Monte Carlo data. Received 18 March 1999 and Received in final form 8 November 1999     

Cluster evaporation transition in finite system

Using Wang-Landau entropic sampling we study the Ising model in the framework of microcanonical ensemble (fixed magnetization). We are working for lattice size up to 1500×1500 in two dimensions and 100×100×100 in three dimensions. As we approach the coexistence curve from inside, varying temperature and keeping the magnetization constant, a first-order phase transition takes place for a temperature near the coexistence curve if the lattice size is large enough. We analyze various features of this transition as well as the scaling behavior of characteristic quantities and we compare our numerical results with existing theories.     

Convection-Diffusion Lattice Boltzmann Scheme for Irregular Lattices

In this paper, a lattice Boltzmann (LB) scheme for convection diffusion on irregular lattices is presented, which is free of any interpolation or coarse graining step. The scheme is derived using the axioma that the velocity moments of the equilibrium distribution equal those of the Maxwell–Boltzmann distribution. The axioma holds for both Bravais and irregular lattices, implying a single framework for LB schemes for all lattice types. By solving benchmark problems we have shown that the scheme is indeed consistent with convection diffusion. Furthermore, we have compared the performance of the LB schemes with that of finite difference and finite element schemes. The comparison shows that the LB scheme has a similar performance as the one-step second-order Lax–Wendroff scheme: it has little numerical diffusion, but has a slight dispersion error. By changing the relaxation parameter ω the dispersion error can be balanced by a small increase of the numerical diffusion.     

Characterisation of a three-dimensional Brownian motor in optical lattices

We present here a detailed study of the behaviour of a three dimensional Brownian motor based on cold atoms in a double optical lattice [P. Sj?lund et al., Phys. Rev. Lett. 96, 190602 (2006)]. This includes both experiments and numerical simulations of a Brownian particle. The potentials used are spatially and temporally symmetric, but combined spatiotemporal symmetry is broken by phase shifts and asymmetric transfer rates between potentials. The diffusion of atoms in the optical lattices is rectified and controlled both in direction and speed along three dimensions. We explore a large range of experimental parameters, where irradiances and detunings of the optical lattice lights are varied within the dissipative regime. Induced drift velocities in the order of one atomic recoil velocity have been achieved.     

Polymers on disordered hierarchical lattices: A nonlinear combination of random variables

The problem of directed polymers on disordered hierarchical and hypercubic lattices is considered. For the hierarchical lattices the problem can be reduced to the study of the stable laws for combining random variables in a nonlinear way. We present the results of numerical simulations of two hierarchical lattices, finding evidence of a phase transition in one case. For a limiting case we extend the perturbation theory developed by Derrida and Griffiths to nonzero temperature and to higher order and use this approach to calculate thermal and geometrical properties (overlaps) of the model. In this limit we obtain an interpolation formula, allowing one to obtain the noninteger moments of the partition function from the integer moments. We obtain bounds for the transition temperature for hierarchical and hypercubic lattices, and some similarities between the problem on the two different types of lattice are discussed.