Source author record

Chin-Kun Hu

Chin-Kun Hu appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

31works
22topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

31 published item(s)

preprint2019arXiv

Nematic and smectic liquid crystals modeling with stiff and free-joint Lennard-Jones chain molecules

Liquid crystals (LCs) composed of mesogens play important roles in various scientific and engineering problems. How a system with many mesogens can enter a LC state is an interesting and important problem. Using stiff and free-joint Lennard-Jones chain molecules as mesogens, we study the conditions under which the mesogens can enter various LC phases. The guideline is to eliminate the unwanted translational orders under a controlled fine-tuning procedure across a sequence of systems. Instead of monitoring the growth of order out of the disorder, we prepare a configuration of high orientation ordering and find out where it relaxes to. Such a procedure begins with a reference system, consisting of short chains of homogeneous soft spheres, in a liquid-vapor coexistence situation, at which the thermodynamic instability triggers a fast spontaneous growing process. By applying a short pulse of auxiliary field to align the dispersedly oriented clusters, followed by reducing the volume and, finally, changing the homogeneous molecules into heterogeneous chains, we are able to obtain a range of systems, including nematic and smectic LCs, at their stable ordered states. The model can be extended to study the influence of nanoparticles or external field on the LC structure.

preprint2019arXiv

Scaling relations and finite-size scaling in gravitationally correlated lattice percolation models

In some systems, the connecting probability (and thus the percolation process) between two sites depends on the geometric distance between them. To understand such process, we propose gravitationally correlated percolation models for link-adding networks on the two-dimensional lattice $G$ with two strategies $S_{\rm max}$ and $S_{\rm min}$, to add a link $l_{i,j}$ to connect site $i$ and site $j$ with mass $m_i$ and $m_j$, respectively; $m_i$ and $m_j$ are sizes of the clusters which contain site $i$ and site $j$, respectively. The probability to add the link $l_{i,j}$ is related to the generalized gravity $g_{ij} \equiv m_i m_j/r_{ij}^d$, where $r_{ij}$ is the geometric distance between $i$ and $j$, and $d$ is an adjustable decaying exponent. In the beginning of the simulation, all sites of $G$ are occupied and there is no link. In the simulation process, two inter-cluster links $l_{i,j}$ and $l_{k,n}$ are randomly chosen and the generalized gravities $g_{ij}$ and $g_{kn}$ are computed. In the strategy $S_{\rm max}$, the link with larger generalized gravity is added. In the strategy $S_{\rm min}$, the link with smaller generalized gravity is added, which include percolation on the Erd\H os-Rényi random graph and the Achlioptas process of explosive percolation as the limiting cases, $d \to \infty$ and $d \to 0$, respectively. Adjustable strategies facilitate or inhibit the network percolation in a generic view. We calculate percolation thresholds $T_c$ and critical exponents $β$ by numerical simulations. We also obtain various finite-size scaling functions for the node fractions in percolating clusters or arrival of saturation length with different intervening strategies.

preprint2018arXiv

Universality and scaling in human and social systems

The objective of statistical physics is to understand macroscopic behavior of a many-body system from the interactions of the constituents of that system. When many-body systems reach critical states, simple universal and scaling behaviors appear. In this talk, I first introduce the concepts of universality and scaling in critical physical systems, I then briefly review some examples of universal and scaling behaviors in human and social systems, e.g. universal crossover behavior of stock returns, universality and scaling in the statistical data of literary works, universal trend in the evolution of states or countries etc. Finally, I mention some interesting problems for further studies.

preprint2016arXiv

Finite-size corrections and scaling for the dimer model on the checkerboard lattice

Lattice models are useful for understanding behaviors of interacting complex many-body systems. The lattice dimer model has been proposed to study the adsorption of diatomic molecules on a substrate. Here we analyze the partition function of the dimer model on an $2 M \times 2 N$ checkerboard lattice wrapped on a torus and derive the exact asymptotic expansion of the logarithm of the partition function. We find that the internal energy at the critical point is equal to zero. We also derive the exact finite-size corrections for the free energy, the internal energy, and the specific heat. Using the exact partition function and finite-size corrections for the dimer model on finite checkerboard lattice we obtain finite-size scaling functions for the free energy, the internal energy, and the specific heat of the dimer model. We investigate the properties of the specific heat near the critical point and find that specific-heat pseudocritical point coincides with the critical point of the thermodynamic limit, which means that the specific-heat shift exponent $λ$ is equal to $\infty$. We have also considered the limit $N \to \infty$ for which we obtain the expansion of the free energy for the dimer model on the infinitely long cylinder. From a finite-size analysis we have found that two conformal field theories with the central charges $c = 1$ for the height function description and $c = -2$ for the construction using a mapping of spanning trees can be used to describe the dimer model on the checkerboard lattice.

preprint2015arXiv

Application of light pressure in the rheology of biopolymers

Precise measuring the deforming force and observation of the deformation dynamics along with difficulties of micro sample preparation are the main problems which make the micro-scale experiments difficult. Here we present a method for direct observation and measurements of mechanical parameters of elastic microfibrils in air using the pressure of optical light. For DNA microfibers, observation of the behavior dynamics of mechanical characteristics and the method sensitivity to changing the structure are demonstrated.

preprint2015arXiv

Mathematical model of influence of friction on the vortex motion

We study the influence of linear friction on the vortex motion in a non-viscous stratified compressible rotating media. Our method can be applied to describe the complex behavior of a tropical cyclone approaching land. In particular, we show that several features of the vortex in the atmosphere such as a significant track deflection, sudden decay and intensification, can be explained already by means of the simplest two dimensional barotropic model, which is a result of averaging over the height in the primitive equations of air motion in the atmosphere. Our theoretical considerations are in a good compliance with the experimental data. In contrast to other models, where first the additional physically reasonable simplifications are made, we deal with special solutions of the full system. Our method is able to explain the phenomenon of the cyclone attracting to the land and interaction of the cyclone with an island.

preprint2015arXiv

Nonlinear stability of two-dimensional axisymmetric vortices in compressible inviscid medium in a rotating reference frame

We study the stability of the vortex in a 2D model of continuous compressible media in a uniformly rotating reference frame. As it is known, the axisymmetric vortex in a fixed reference frame is stable with respect to asymmetric perturbations for the solution of the 2D incompressible Euler equations and basically instable for compressible Euler equations. We show that the situation is quite different for a compressible axisymmetric vortex in a rotating reference frame. First, we consider special solutions with linear profile of velocity (or with spatially-uniform velocity gradients), which are important because many real vortices have similar structure near their centers. We analyze both cyclonic and anticyclonic cases and show that the stability of the solution depends only on the ratio of the vorticity to the Coriolis parameter. Using a very delicate analysis along with computer aided proof, we show that the stability of solutions can take place only for a narrow range of this ratio. Our results imply that the rotation of the coordinate frame can stabilize the compressible vortex. Further, we perform both analytical and numerical analysis of stability for real-shaped vortices.

preprint2015arXiv

Thermal-induced proteinquake in oxyhemoglobin

Oxygen is released to living tissues via conformational changes of hemoglobin from R-state (oxyhemoglobin) to T-state (desoxyhemoglobin). The detailed mechanism of this process is not yet fully understood. We have carried out micromechanical experiments on oxyhemoglobin crystals to determine the behavior of the Young's modulus and the internal friction for temperatures between 20 C and 70 C. We have found that around 49 C oxyhemoglobin crystal samples undergo a sudden and strong increase of their Young's modulus, accompanied by a sudden decrease of the internal friction. This sudden mechanical change (proteinquake) takes place in a partially unfolded state and precedes the full denaturation transition at higher temperatures. The hemoglobin crystals after the proteinquake has the same mechanical properies as the initial state at room temperatures. We conjecture that it can be relevant for explaining the oxygen-releasing function of native oxyhemoglobin when the temperature is increased, e.g. due to active sport. The effect is specific for the quaternary structure of hemoglobin, and is absent for myoglobin with only one peptide sequence.

preprint2014arXiv

Non-equilibrium Lyapunov function and a fluctuation relation for stochastic systems: Poisson representation approach

We present a statistical physics framework for description of nonlinear non-equilibrium stochastic processes, modeled via chemical master equation, in the weak-noise limit. Using the Poisson representation approach and applying the large-deviation principle we first solve the master equation. Then we use the notion of the non-equilibrium free energy to derive an integral fluctuation relation for nonlinear non-equilibrium systems under feedback control. We point out that the free energy as well as some functionals can serve as non-equilibrium Lyapunov function which has an important property to decay to its minimal value monotonously at all times. The Poisson representation technique is illustrated via exact stochastic treatment of biophysical processes, such as bacterial chemosensing and molecular evolution.

preprint2013arXiv

A structure-based model fails to probe the mechanical unfolding pathways of the titin I27 domain

We discuss the use of a structure based C$α$-Go model and Langevin dynamics to study in detail the mechanical properties and unfolding pathway of the titin I27 domain. We show that a simple Go-model does detect correctly the origin of the mechanical stability of this domain. The unfolding free energy landscape parameters $x_u$ and $ΔG^{\ddagger}$, extracted from dependencies of unfolding forces on pulling speeds, are found to agree reasonably well with experiments. We predict that above $v=10^4$ nm/s the additional force-induced intermediate state is populated at an end-to-end extension of about $75 \mathring{A}$. The force-induced switch in the unfolding pathway occurs at the critical pulling speed $v_{crit} \approx 10^6-10^7$ nm/s. We argue that this critical pulling speed is an upper limit of the interval where Bell's theory works. However, our results suggest that the Go-model fails to reproduce the experimentally observed mechanical unfolding pathway properly, yielding an incomplete picture of the free energy landscape. Surprisingly, the experimentally observed intermediate state with the A strand detached is not populated in Go-model simulations over a wide range of pulling speeds. The discrepancy between simulation and experiment is clearly seen from the early stage of the unfolding process which shows the limitation of the Go model in reproducing unfolding pathways and deciphering the complete picture of the free energy landscape.

preprint2013arXiv

On diffusive stability of Eigen's quasispecies model

Eigen's quasispecies system with explicit space and global regulation is considered. Limit behavior and stability of the system in a functional space under perturbations of a diffusion matrix with nonnegative spectrum are investigated. It is proven that if the diffusion matrix has only positive eigenvalues then the solutions of the distributed system converge to the equilibrium solution of the corresponding local dynamical system. These results imply that the error threshold does not change if the spatial interactions under the principle of global regulation are taken into account.

preprint2012arXiv

Finite-size corrections for logarithmic representations in critical dense polymers

We study (analytic) finite-size corrections in the dense polymer model on the strip by perturbing the critical Hamiltonian with irrelevant operators belonging to the tower of the identity. We generalize the perturbation expansion to include Jordan cells, and examine whether the finite-size corrections are sensitive to the properties of indecomposable representations appearing in the conformal spectrum, in particular their indecomposability parameters. We find, at first order, that the corrections do not depend on these parameters nor even on the presence of Jordan cells. Though the corrections themselves are not universal, the ratios are universal and correctly reproduced by the conformal perturbative approach, to first order.

preprint2012arXiv

Phase diagram and universality of the Lennard-Jones gas-liquid system

The gas-liquid phase transition of the three-dimensional Lennard-Jones particles system is studied by molecular dynamics simulations. The gas and liquid densities in the coexisting state are determined with high accuracy. The critical point is determined by the block density analysis of the Binder parameter with the aid of the law of rectilinear diameter. From the critical behavior of the gas-liquid coexsisting density, the critical exponent of the order parameter is estimated to be $β= 0.3285(7)$. Surface tension is estimated from interface broadening behavior due to capillary waves. From the critical behavior of the surface tension, the critical exponent of the correlation length is estimated to be $ν= 0.63 (4)$. The obtained values of $β$ and $ν$ are consistent with those of the Ising universality class.

preprint2012arXiv

Universal geometrical factor of protein conformations as a consequence of energy minimization

The biological activity and functional specificity of proteins depend on their native three-dimensional structures determined by inter- and intra-molecular interactions. In this paper, we investigate the geometrical factor of protein conformation as a consequence of energy minimization in protein folding. Folding simulations of 10 polypeptides with chain length ranging from 183 to 548 residues manifest that the dimensionless ratio (V/(A<r>)) of the van der Waals volume V to the surface area A and average atomic radius <r> of the folded structures, calculated with atomic radii setting used in SMMP [Eisenmenger F., et. al., Comput. Phys. Commun., 138 (2001) 192], approach 0.49 quickly during the course of energy minimization. A large scale analysis of protein structures show that the ratio for real and well-designed proteins is universal and equal to 0.491\pm0.005. The fractional composition of hydrophobic and hydrophilic residues does not affect the ratio substantially. The ratio also holds for intrinsically disordered proteins, while it ceases to be universal for polypeptides with bad folding properties.

preprint2010arXiv

Boundary conditions and amplitude ratios for finite-size corrections of a one-dimensional quantum spin model

We study the influence of boundary conditions on the finite-size corrections of a one-dimensional (1D) quantum spin model by exact and perturbative theoretic calculations. We obtain two new infinite sets of universal amplitude ratios for the finite-size correction terms of the 1D quantum spin model of $N$ sites with free and antiperiodic boundary conditions. The results for the lowest two orders are in perfect agreement with a perturbative conformal field theory scenario proposed by Cardy [Nucl. Phys. B {\bf 270}, 186 (1986)].

preprint2010arXiv

DNA Bubble Life Time in Denaturation

We have investigated the denaturation bubble life time for a homogeneous as well as for a heterogeneous DNA within a Poland-Scheraga model. It is shown that at criticality the bubble life time for a homogeneous DNA is finite provided that the loop entropic exponent c>2 and has a scaling dependence on DNA length for c<2. Heterogeneity in the thermodynamical limit makes the bubble life time infinite for any entropic exponent.

preprint2010arXiv

Effect of time delay on the onset of synchronization of the stochastic Kuramoto model

We consider the Kuramoto model of globally coupled phase oscillators with time-delayed interactions, that is subject to the Ornstein-Uhlenbeck (Gaussian) colored or the non-Gaussian colored noise. We investigate numerically the interplay between the influences of the finite correlation time of noise $τ$ and the time delay $τ_{d}$ on the onset of the synchronization process. Both cases for identical and nonidentical oscillators had been considered. Among the obtained results for identical oscillators is a large increase of the synchronization threshold as a function of time delay for the colored non-Gaussian noise compared to the case of the colored Gaussian noise at low noise correlation time $τ$. However, the difference reduces remarkably for large noise correlation times. For the case of nonidentical oscillators, the incoherent state may become unstable around the maximum value of the threshold (as a function of time delay) even at lower coupling strength values in the presence of colored noise as compared to the noiseless case. We had studied the dependence of the critical value of the coupling strength (the threshold of synchronization) on given parameters of the stochastic Kuramoto model in great details and presented results for possible cases of colored Gaussian and non-Gaussian noises.

preprint2010arXiv

Glassy state of native collagen fibril

Our micromechanical experiments show that at physiological temperatures type I collagen fibril has several basic features of the glassy state. The transition out of this state [softening transition] essentially depends on the speed of heating v, e.g., for v=1 C/min it occurs around 70 C and is displayed by a peak of the internal friction and decreasing Young's modulus. The softening transition decreases by 45 C upon decreasing the heating speed to v=0.1 C/min. For temperatures 20-30 C the native collagen fibril demonstrates features of mechanical glassines at oscillation frequencies 0.1-3 kHz; in particular, the internal friction has a sharp maximum as a function of the frequency. This is the first example of biopolymer glassines at physiological temperatures, because well-known glassy features of DNA and globular proteins are seen only for much lower temperatures (around 200 K).

preprint2010arXiv

On trajectories of vortices in the compressible fluid on a two-dimensional manifold

For the model of a compressible barotropic fluid on a two dimensional rotating Riemmanian manifold we discuss a special class of smooth solutions having a form of a steady non-singular vortex moving with a bearing field. The model can be obtained from the system of primitive equations governing the motion of air over the Earth surface after averaging over the height and therefore the solution obtained can be interpreted as a tropical cyclone which is known as a long time existing stable vortex. We consider approximations of $l$- plane and $β$ - plane used in geophysics for modeling of middle scale processes and equations on the whole sphere as well. We show that the solutions of the mentioned form satisfy the equations of the model either exactly or with a discrepancy which is small in a neighborhood of the trajectory of the center of vortex. We perform a numeric study of the change of the shape of the vortex affected by the neglecting the discrepancy term.

preprint2010arXiv

Typhoon eye trajectory based on a mathematical model: comparing with observational data

We propose a model based on the primitive system of the Navier-Stokes equations in a bidimensional framework as the $l$ - plane approximation, which allows us to explain the variety of tracks of tropical cyclones (typhoons). Our idea is to construct special analytical solutions with a linear velocity profile for the Navier-Stokes systems. The evidence of the structure of linear velocity near the center of vortex can be proven by the observational data. We study solutions with the linear-velocity property for both barotropic and baroclinic cases and show that they follow the same equations in describing the trajectories of the typhoon eye at the equilibrium state (that relates to the conservative phase of the typhoon dynamics). Moreover, at the equilibrium state, the trajectories can be viewed as a superposition of two circular motions: one has period $2π/l,$ the other one has period $2π/b_0,$ where $l$ is the Coriolis parameter and $b_0$ is the height-averaged vorticity at the center of cyclone. Also, we compare our theoretical trajectories based on initial conditions from the flow with tracks obtained from the observational database. It is worth to mention that under certain conditions our results are still compatible with observational data although we did not truly consider the influence of steering effect. %Note that the %motion of the typhoon eye can not be totally determinated by initial %conditions due to the effect of ambient pressure field ("steering" %effect). Finally, we propose the parameter-adopting method so that one could correct the weather prediction in real time. Examples of our analysis and the use of parameter-adopting method for the historic trajectories are provided.

preprint2009arXiv

How adsorption influences DNA denaturation

The thermally induced denaturation of DNA in the presence of attractive solid surface is studied. The two strands of DNA are modeled via two coupled flexible chains without volume interactions. If the two strands are adsorbed on the surface, the denaturation phase-transition disappears. Instead, there is a smooth crossover to a weakly naturated state. Our second conclusion is that even when the inter-strand attraction alone is too weak for creating a naturated state at the given temperature, and also the surface-strand attraction alone is too weak for creating an adsorbed state, the combined effect of the two attractions can lead to a naturated and adsorbed state.

preprint2009arXiv

Protein mechanical unfolding: importance of non-native interactions

Mechanical unfolding of the fourth domain of Distyostelium discoideum filamin (DDFLN4) was studied by all-atom molecular dynamics simulations, using the GROMOS96 force field 43a1 and the simple point charge explicit water solvent. Our study reveals an important role of non-native interactions in the unfolding process. Namely, the existence of a peak centered at the end-to-end extension 22 nm in the force-extension curve, is associated with breaking of non-native hydrogen bonds. Such a peak has been observed in experiments but not in Go models, where non-native interactions are neglected. We predict that an additional peak occurs at 2 nm using not only GROMOS96 force field 43a1 but also Amber 94 and OPLS force fields. This result would stimulate further experimental studies on elastic properties of DDFLN4.

preprint2009arXiv

Replicators in Fine-grained Environment: Adaptation and Polymorphism

Selection in a time-periodic environment is modeled via the two-player replicator dynamics. For sufficiently fast environmental changes, this is reduced to a multi-player replicator dynamics in a constant environment. The two-player terms correspond to the time-averaged payoffs, while the three and four-player terms arise from the adaptation of the morphs to their varying environment. Such multi-player (adaptive) terms can induce a stable polymorphism. The establishment of the polymorphism in partnership games [genetic selection] is accompanied by decreasing mean fitness of the population.

preprint2009arXiv

The intermediate evolution phase in case of truncated selection

Using methods of statistical physics, we present rigorous theoretical calculations of Eigen's quasispecies theory with the truncated fitness landscape which dramatically limits the available sequence space of a reproducing quasispecies. Depending on the mutation rates, we observe three phases, a selective one, an intermediate one with some residual order and a completely randomized phase. Our results are applicable for the general case of fitness landscape.

preprint2006arXiv

Effect of finite size on cooperativity and rates of protein folding

We analyze the dependence of cooperativity of the thermal denaturation transition and folding rates of globular proteins on the number of amino acid residues, $N$, using lattice models with side chains,off-lattice Go models and the available experimental data. A dimensionless measure of cooperativity, $Ω_c$ ($0 < Ω_c < \infty$), scales as $Ω_c \sim N^ζ$. The results of simulations and the analysis of experimental data further confirm the earlier prediction that $ζ$ is universal with $ζ= 1 +γ$, where exponent $γ$ characterizes the susceptibility of a self-avoiding walk. This finding suggests that the structural characteristics in the denaturated state are manifested in the folding cooperativity at the transition temperature. The folding rates $k_F$ for the Go models and a dataset of 69 proteins can be fit using $k_F = k_F^0 \exp(-cN^β)$. Both $β= 1/2$ and 2/3 provide a good fit of the data. We find that $k_F = k_F^0 \exp(-cN^{1/2})$, with the average (over the dataset of proteins) $k_F^0 \approx (0.2μs)^{-1}$ and $c \approx 1.1$, can be used to estimate folding rates to within an order of magnitude in most cases. The minimal models give identical $N$ dependence with $c \approx 1$. The prefactor for off-lattice Go models is nearly four orders of magnitude larger than the experimental value.

preprint2006arXiv

Folding of the Protein Domain hbSBD

The folding of the alpha-helice domain hbSBD of the mammalian mitochondrial branched-chain alpha-ketoacid dehydrogenase (BCKD) complex is studied by the circular dichroism technique in absence of urea. Thermal denaturation is used to evaluate various thermodynamic parameters defining the equilibrium unfolding, which is well described by the two-state model with the folding temperature T_f = 317.8 K and the enthalpy change Delta H_g = 19.67 kcal/mol. The folding is also studied numerically using the off-lattice coarse-grained Go model and the Langevin dynamics. The obtained results, including the population of the native basin, the free energy landscape as a function of the number of native contacts and the folding kinetics, also suggest that the hbSBD domain is a two-state folder. These results are consistent with the biological function of hbSBD in BCKD.

preprint2002arXiv

Exact finite-size corrections for the square lattice Ising model with Brascamp-Kunz boundary conditions

Finite-size scaling, finite-size corrections, and boundary effects for critical systems have attracted much attention in recent years. Here we derive exact finite-size corrections for the free energy F and the specific heat C of the critical ferromagnetic Ising model on the M x 2N square lattice with Brascamp-Kunz (BK) boundary conditions [J. Math. Phys. 15, 66 (1974)] and compared such results with those under toroidal boundary conditions. In the limit N \to \infty we obtain the expansion of the free energy for infinitely long strip with BK boundary conditions. Our results are consistent with the conformal field theory prediction for the mixed boundary conditions by Cardy [Nucl. Phys. {\bf B 275}, 200 (1986)] although the definitions of boundary conditions in two cases are different in one side of the long strip.

preprint2002arXiv

Exact partition functions of the Ising model on MxN planar lattices with periodic-aperiodic boundary conditions

The Grassmann path integral approach is used to calculate exact partition functions of the Ising model on MxN square (sq), plane triangular (pt) and honeycomb (hc) lattices with periodic-periodic (pp), periodic-antiperiodic (pa), antiperiodic-periodic (ap) and antiperiodic-antiperiodic (aa) boundary conditions. The partition functions are used to calculate and plot the specific heat, $C/k_B$, as a function of the temperature, $θ=k_BT/J$. We find that for the NxN sq lattice, $C/k_B$ for pa and ap boundary conditions are different from those for aa boundary conditions, but for the NxN pt and hc lattices, $C/k_B$ for ap, pa, and aa boundary conditions have the same values. Our exact partition functions might also be useful for understanding the effects of lattice structures and boundary conditions on critical finite-size corrections of the Ising model.

preprint2001arXiv

Proteinlike behavior of a spin system near the transition between ferromagnet and spin glass

A simple spin system is studied as an analog for proteins. We investigate how the introduction of randomness and frustration into the system effects the designability and stability of ground state configurations. We observe that the spin system exhibits protein-like behavior in the vicinity of the transition between ferromagnet and spin glass. Our results illuminate some guiding principles in protein evolution.

preprint2001arXiv

Universal scaling functions for bond percolation on planar random and square lattices with multiple percolating clusters

Percolation models with multiple percolating clusters have attracted much attention in recent years. Here we use Monte Carlo simulations to study bond percolation on $L_{1}\times L_{2}$ planar random lattices, duals of random lattices, and square lattices with free and periodic boundary conditions, in vertical and horizontal directions, respectively, and with various aspect ratio $L_{1}/L_{2}$. We calculate the probability for the appearance of $n$ percolating clusters, $W_{n},$ the percolating probabilities, $P$, the average fraction of lattice bonds (sites) in the percolating clusters, $<c^{b}>_{n}$ ($<c^{s}>_{n}$), and the probability distribution function for the fraction $c$ of lattice bonds (sites), in percolating clusters of subgraphs with $n$ percolating clusters, $f_{n}(c^{b})$ ($f_{n}(c^{s})$). Using a small number of nonuniversal metric factors, we find that $W_{n}$, $P$, $<c^{b}>_{n}$ ($<c^{s}>_{n}$), and $f_{n}(c^{b})$ ($f_{n}(c^{s})$) for random lattices, duals of random lattices, and square lattices have the same universal finite-size scaling functions. We also find that nonuniversal metric factors are independent of boundary conditions and aspect ratios.

preprint1999arXiv

Exact phase diagrams for an Ising model on a two-layer Bethe lattice

Using an iteration technique, we obtain exact expressions for the free energy and the magnetization of an Ising model on a two - layer Bethe lattice with intralayer coupling constants J1 and J2 for the first and the second layer, respectively, and interlayer coupling constant J3 between the two layers; the Ising spins also couple with external magnetic fields, which are different in the two layers. We obtain exact phase diagrams for the system.