Liquid-Solid Phase Transition of the System with Two particles in a Rectangular Box

We study the statistical properties of two hard spheres in a two dimensional rectangular box. In this system, the relation like Van der Waals equation loop is obtained between the width of the box and the pressure working on side walls. The auto-correlation function of each particle's position is calculated numerically. By this calculation near the critical width, the time at which the correlation become zero gets longer according to the increase of the height of the box. Moreover, fast and slow relaxation processes like $\alpha$ and $\beta$ relaxations observed in supper cooled liquid are observed when the height of the box is sufficiently large. These relaxation processes are discussed with the probability distribution of relative position of two particles.Comment: 6 figure

Evaluation of a present-day climate simulation with a new coupled atmosphere-ocean model GENMOM

We present a new, non-flux corrected AOGCM, GENMOM, that combines the GENESIS version 3 atmospheric GCM (Global Environmental and Ecological Simulation of Interactive Systems) and MOM2 (Modular Ocean Model version 2) nominally at T31 resolution. We evaluate GENMOM by comparison with reanalysis products (e.g., NCEP2) and three models used in the IPCC AR4 assessment. GENMOM produces a global temperature bias of 0.6 °C. Atmospheric features such as the jet stream structure and major semi-permanent sea level pressure centers are well simulated as is the mean planetary-scale wind structure that is needed to produce the correct position of stormtracks. Most ocean surface currents are reproduced except where they are not resolvable at T31 resolution. Overall, GENMOM captures reasonably well the observed gradients and spatial distributions of annual surface temperature and precipitation and the simulations are on par with other AOGCMs. Deficiencies in the GENMOM simulations include a warm bias in the surface temperature over the southern oceans, a split in the ITCZ and weaker-than-observed overturning circulation

Magneto-Transport in the Two-Dimensional Lorentz Gas

We consider the two-dimensional Lorentz gas with Poisson distributed hard disk scatterers and a constant magnetic field perpendicular to the plane of motion. The velocity autocorrelation is computed numerically over the full range of densities and magnetic fields with particular attention to the percolation threshold between hopping transport and pure edge currents. The Ohmic and Hall conductance are compared with mode-coupling theory and a recent generalized kinetic equation valid for low densities and small fields. We argue that the long time tail as $t^{-2}$ persists for non-zero magnetic field.Comment: 7 pages, 14 figures. Uses RevTeX and epsfig.sty. Submitted to Physical Review

Global climate simulations at 3000-year intervals for the last 21,000 years with the GENMOM coupled atmosphere–ocean model

We apply GENMOM, a coupled atmosphere–ocean climate model, to simulate eight equilibrium time slices at 3000-year intervals for the past 21,000 years forced by changes in Earth–Sun geometry, atmospheric greenhouse gases (GHGs), continental ice sheets, and sea level. Simulated global cooling during the Last Glacial Maximum (LGM) is 3.8°C and the rate of post-glacial warming is in overall agreement with recently published temperature reconstructions. The greatest rate of warming occurs between 15 and 12 ka (2.4°C over land, 0.7°C over oceans, and 1.4°C globally) in response to changes in radiative forcing from the diminished extent of the Northern Hemisphere (NH) ice sheets and increases in GHGs and NH summer insolation. The modeled LGM and 6 ka temperature and precipitation climatologies are generally consistent with proxy reconstructions, the PMIP2 and PMIP3 simulations, and other paleoclimate data–model analyses. The model does not capture the mid-Holocene “thermal maximum” and gradual cooling to preindustrial (PI) global temperature found in the data. Simulated monsoonal precipitation in North Africa peaks between 12 and 9 ka at values ~50% greater than those of the PI, and Indian monsoonal precipitation peaks at 12 and 9 ka at values ~45% greater than the PI. GENMOM captures the reconstructed LGM extent of NH and Southern Hemisphere (SH) sea ice. The simulated present-day Antarctica Circumpolar Current (ACC) is ~48% weaker than the observed (62 versus 119 Sv). The simulated present-day Atlantic Meridional Overturning Circulation (AMOC) of 19:3 ± 1:4 Sv on the Bermuda Rise (33°N) is comparable with observed value of 18:7 ± 4:8 Sv. AMOC at 33°N is reduced by ~15% during the LGM, and the largest post-glacial increase (~11 %) occurs during the 15 ka time slice

Numerical Evidence for Divergent Burnett Coefficients

In previous papers [Phys. Rev. A {\bf 41}, 4501 (1990), Phys. Rev. E {\bf 18}, 3178 (1993)], simple equilibrium expressions were obtained for nonlinear Burnett coefficients. A preliminary calculation of a 32 particle Lennard-Jones fluid was presented in the previous paper. Now, sufficient resources have become available to address the question of whether nonlinear Burnett coefficients are finite for soft spheres. The hard sphere case is known to have infinite nonlinear Burnett coefficients (ie a nonanalytic constitutive relation) from mode coupling theory. This paper reports a molecular dynamics caclulation of the third order nonlinear Burnett coefficient of a Lennard-Jones fluid undergoing colour flow, which indicates that this term is diverges in the thermodynamic limit.Comment: 12 pages, 9 figure

Mode-coupling theory for multiple-time correlation functions of tagged particle densities and dynamical filters designed for glassy systems

The theoretical framework for higher-order correlation functions involving multiple times and multiple points in a classical, many-body system developed by Van Zon and Schofield [Phys. Rev. E 65, 011106 (2002)] is extended here to include tagged particle densities. Such densities have found an intriguing application as proposed measures of dynamical heterogeneities in structural glasses. The theoretical formalism is based upon projection operator techniques which are used to isolate the slow time evolution of dynamical variables by expanding the slowly-evolving component of arbitrary variables in an infinite basis composed of the products of slow variables of the system. The resulting formally exact mode-coupling expressions for multiple-point and multiple-time correlation functions are made tractable by applying the so-called N-ordering method. This theory is used to derive for moderate densities the leading mode coupling expressions for indicators of relaxation type and domain relaxation, which use dynamical filters that lead to multiple-time correlations of a tagged particle density. The mode coupling expressions for higher order correlation functions are also succesfully tested against simulations of a hard sphere fluid at relatively low density.Comment: 15 pages, 2 figure

Global Equation of State of two-dimensional hard sphere systems

Hard sphere systems in two dimensions are examined for arbitrary density. Simulation results are compared to the theoretical predictions for both the low and the high density limit, where the system is either disordered or ordered, respectively. The pressure in the system increases with the density, except for an intermediate range of volume fractions $0.65 \le \nu \le 0.75$, where a disorder-order phase transition occurs. The proposed {\em global equation of state} (which describes the pressure {\em for all densities}) is applied to the situation of an extremely dense hard sphere gas in a gravitational field and shows reasonable agreement with both experimental and numerical data.Comment: 4 pages, 2 figure

Internal Anisotropy of Collision Cascades

We investigate the internal anisotropy of collision cascades arising from the branching structure. We show that the global fractal dimension cannot give an adequate description of the geometrical structure of cascades because it is insensitive to the internal anisotropy. In order to give a more elaborate description we introduce an angular correlation function, which takes into account the direction of the local growth of the branches of the cascades. It is demonstrated that the angular correlation function gives a quantitative description of the directionality and the interrelation of branches. The power law decay of the angular correlation is evidenced and characterized by an exponent and an angular correlation length different from the radius of gyration. It is demonstrated that the overlapping of subcascades has a strong effect on the angular correlation.Comment: RevteX, 8 pages, 6 .eps figures include