Stochastic modeling of a billiard in a gravitational field: Power law behavior of Lyapunov exponents

We consider the motion of a point particle (billiard) in a uniform gravitational field constrained to move in a symmetric wedge-shaped region. The billiard is reflected at the wedge boundary. The phase space of the system naturally divides itself into two regions in which the tangent maps are respectively parabolic and hyperbolic. It is known that the system is integrable for two values of the wedge half-angle ₁ and ₂ and chaotic for ₁<< ₂. We study the system at three levels of approximation: first, where the deterministic dynamics is replaced by a random evolution; second, where, in addition, the tangent map in each region is, replaced by its average; and third, where the tangent map is replaced by a single global average. We show that at all three levels the Lyapunov exponent exhibits power law behavior near ₁ and ₂ with exponents 1/2 and 1, respectively. We indicate the origin of the exponent 1, which has not been observed in unaccelerated billiards.     

Restoration of universality for the rod-to-coil transition scaling in the infinite-dimensionality limit: Exact results for directed walks

We derive scaling forms for the thermodynamic and correlation quantities for the turn-weighted fully and partially directed self-avoiding walks on the hypercubic lattices ind2. In the grand canonical (fixed fugacity per step) ensemble, the conformational rod-to-coil transition sets up in the regimew¯N=O(1), wherew is the weight of each 90° turn and¯N is the (fugacity-dependent) average number of steps. Contrary to the conventional critical phenomena wisdom, the scaling functions for the two different walk models, directed and partially directed, become universal only in the limitd.     

Identifying potential binding modes and explaining partitioning behavior using flexible alignments and multidimensional scaling

A method is described for the rapid and automatic analysis of flexible molecular alignments using multidimensional scaling and a normalized scoring scheme. A projection scheme was devised to separate orientational and conformational effects. It is shown that the approach can be utilized for the identification of common binding orientations or to the study of differences in partioning behavior. It is suggested that the method can be employed as a novel approach exploring molecular similarity as a dynamic property, so that it includes aspects of motion (by way of mutual orientations), conformations and molecular properties.     

Dynamic percolation transition induced by phase separation: A Monte Carlo analysis

The percolation transition of geometric clusters in the three-dimensional, simple cubic, nearest neighbor Ising lattice gas model is investigated in the temperature and concentration region inside the coexistence curve. We consider quenching experiments, where the system starts from an initially completely random configuration (corresponding to equilibrium at infinite temperature), letting the system evolve at the considered temperature according to the Kawasaki spinexchange dynamics. Analyzing the distributionn _l(t) of clusters of sizel at timet, we find that after a time of the order of about 100 Monte Carlo steps per site a percolation transition occurs at a concentration distinctly lower than the percolation concentration of the initial random state. This dynamic percolation transition is analyzed with finite-size scaling methods. While at zero temperature, where the system settles down at a frozen-in cluster distribution and further phase separation stops, the critical exponents associated with this percolation transition are consistent with the universality class of random percolation, the critical behavior of the transient time-dependent percolation occurring at nonzero temperature possibly belongs to a different, new universality class.     

Modeling high-pressure densities at wide temperature range with volume scaling: Cyclohexane + n-hexadecane mixtures

High-pressure density data for cyclohexane + n-hexadecane mixtures at a wide temperature range was modeled with several classical equations of state (EOS) and correlative models. A modification for softening the co-volume and another for a volume scaling of the Peng–Robinson EOS (VS-PR) were proposed. The VS-PR model is able to correlate the pure component experimental data employing only five adjustable parameters, with root-mean-square deviation (RMSD) between calculated and experimental densities essentially within the experimental error. This result is superior to widely used approaches, i.e., a six parameter Tait model and six parameter volume translations (temperature and pressure dependent) for Peng–Robinson and Patel–Teja EOS. The VS-PR model also represents well the isobaric thermal expansion and the isothermal compressibility coefficients of the pure cyclohexane, a small naphthenic substance as well as a long chain n-alkane hydrocarbon, n-hexadecane. When modeling the mixture data, the use of VS-PR model of pure components along with the Redlich–Kister expansion, truncated at the first term, the density was correlated within a RMSD only 60% greater than the experimental error. The proposed model is able to accurately represent all the tested mixture data with a relatively small number of parameters.     

Linear scaling local correlation approach for solving the coupled cluster equations of large systems

A linear scaling local correlation approach is proposed for approximately solving the coupled cluster doubles (CCD) equations of large systems in a basis of orthogonal localized molecular orbitals (LMOs). By restricting double excitations from spatially close occupied LMOs into their associated virtual LMOs, the number of significant excitation amplitudes scales only linearly with molecular size in large molecules. Significant amplitudes are obtained to a very good approximation by solving the CCD equations of various subsystems, each of which is made up of a cluster associated with the orbital indices of a subset of significant amplitudes and the local environmental domain of the cluster. The combined effect of these two approximations leads to a linear scaling algorithm for large systems. By using typical thresholds, which are designed to target an energy accuracy, our numerical calculations for a wide range of molecules using the 6-31G or 6-31G* basis set demonstrate that the present local correlation approach recovers more than 98.5% of the conventional CCD correlation energy.     

Phase separation of symmetrical polymer mixtures in thin-film geometry

Monte Carlo simulations of the bond fluctuation model of symmetrical polymer blends confined between two neutral repulsive walls are presented for chain lengthN _A=N _B=32 and a wide range of film thicknessD (fromD=8 toD=48 in units of the lattice spacing). The critical temperaturesT _c (D) of unmixing are located by finite-size scaling methods, and it is shown that , wherev ₃0.63 is the correlation length exponent of the three-dimensional Ising model universality class. Contrary to this result, it is argued that the critical behavior of the films is ruled by two-dimensional exponents, e.g., the coexistence curve (difference in volume fraction of A-rich and A-poor phases) scales as , where ₂ is the critical exponent of the two-dimensional Ising universality class ( ₂=1/8). Since for largeD this asymptotic critical behavior is confined to an extremely narrow vicinity ofT _c (D), one observes in practice effective exponents which gradually cross over from ₂ to ₃ with increasing film thickness. This anomalous flattening of the coexistence curve should be observable experimentally.     

Multifractal analysis of Brownian zero set

We present a method for the derivation of the generating function and computation of critical exponents for several cluster models (staircase, bar-graph, and directed column-convex polygons, as well as partially directed self-avoiding walks), starting with nonlinear functional equations for the generating function. By linearizing these equations, we first give a derivation of the generating functions. The nonlinear equations are further used to compute the thermodynamic critical exponents via a formal perturbation ansatz. Alternatively, taking the continuum limit leads to nonlinear differential equations, from which one can extract the scaling function. We find that all the above models are in the same universality class with exponents _u =-1/2, _i =-1/3, and =2/3. All models have as their scaling function the logarithmic derivative of the Airy function.     

Schrödinger invariance and strongly anisotropic critical systems

The extension of strongly anisotropic or dynamical scaling to local scale invariance is investigated. For the special case of an anisotropy or dynamical exponent =z=2, the group of local scale transformation considered is the Schrödinger group, which can be obtained as the nonrelativistic limit of the conformal group. The requirement of Schrödinger invariance determines the two-point function in the bulk and reduces the three-point function to a scaling form of a single variable. Scaling forms are also derived for the two-point function close to a free surface which can be either spacelike or timelike. These results are reproduced in several exactly solvable statistical systems, namely the kinetic Ising model with Glauber dynamics, lattice diffusion, Lifshitz points in the spherical model, and critical dynamics of the spherical model with a nonconserved order parameter. For generic values of , evidence from higher-order Lifshitz points in the spherical model and from directed percolation suggests a simple scaling form of the two-point function.     

The covariance matrix of the Potts model: A random cluster analysis

We consider the covariance matrix,G ^mm =q ²<(_x,m);(_y,m)>, of thed-dimensionalq-states Potts model, rewriting it in the random cluster representation of Fortuin and Kasteleyn. In any of theq ordered phases, we identify the eigenvalues of this matrix both in terms of representations of the unbroken symmetry group of the model and in terms of random cluster connectivities and covariances, thereby attributing algebraic significance to these stochastic geometric quantities. We also show that the correlation length corresponding to the decay rate of one of the eigenvalues is the same as the inverse decay rate of the diameter of finite clusers. For dimensiond=2, we show that this correlation length and the correlation length of the two-point function with free boundary conditions at the corresponding dual temperature are equal up to a factor of two. For systems with first-order transitions, this relation helps to resolve certain inconsistencies between recent exact and numerical work on correlation lengths at the self-dual point _o. For systems with second order transitions, this relation implies the equality of the correlation length exponents from above and below threshold, as well as an amplitude ratio of two. In the course of proving the above results, we establish several properties of independent interest, including left continuity of the inverse correlation length with free boundary conditions and upper semicontinuity of the decay rate for finite clusters in all dimensions, and left continuity of the two-dimensional free boundary condition percolation probability at _o. We also introduce DLR equations for the random cluster model and use them to establish ergodicity of the free measure. In order to prove these results, we introduce a new class of events which we call decoupling events and two inequalities for these events. The first is similar to the FKG inequality, but holds for events which are neither increasing nor decreasing; the second is similar to the van den Berg-Kesten inequality in standard percolation. Both inequalities hold for an arbitrary FKG measure.