Phase transition in inelastic disks

This letter investigates the molecular dynamics of inelastic disks without external forcing. By introducing a new observation frame with a rescaled time, we observe the virtual steady states converted from asymptotic energy dissipation processes. System behavior in the thermodynamic limit is carefully investigated. It is found that a phase transition with symmetry breaking occurs when the magnitude of dissipation is greater than a critical value.Comment: 9 pages, 6 figure

Long-Time Behavior of Velocity Autocorrelation Function for Interacting Particles in a Two-Dimensional Disordered System

The long-time behavior of the velocity autocorrelation function (VACF) is investigated by the molecular dynamics simulation of a two-dimensional system which has both a many-body interaction and a random potential. With strengthening the random potential by increasing the density of impurities, a crossover behavior of the VACF is observed from a positive tail, which is proportional to t^{-1}, to a negative tail, proportional to -t^{-2}. The latter tail exists even when the density of particles is the same order as the density of impurities. The behavior of the VACF in a nonequilibrium steady state is also studied. In the linear response regime the behavior is similar to that in the equilibrium state, whereas it changes drastically in the nonlinear response regime.Comment: 12 pages, 5 figure

The Semiclassical Coulomb Interaction

The semiclassical Coulomb excitation interaction is at times expressed in the Lorentz gauge in terms of the electromagnetic fields and a contribution from the scalar electric potential. We point out that the potential term can make spurious contributions to excitation cross sections, especially when the the decay of excited states is taken into account. We show that, through an appropriate gauge transformation, the excitation interaction can be expressed in terms of the electromagnetic fields alone.Comment: 12 pages. Phys. Rev. C, Rapid Communication, in pres

Velocity Correlations in Driven Two-Dimensional Granular Media

Simulations of volumetrically forced granular media in two dimensions produce s tates with nearly homogeneous density. In these states, long-range velocity correlations with a characteristic vortex structure develop; given sufficient time, the correlations fill the entire simulated area. These velocity correlations reduce the rate and violence of collisions, so that pressure is smaller for driven inelastic particles than for undriven elastic particles in the same thermodynamic state. As the simulation box size increases, the effects of veloc ity correlations on the pressure are enhanced rather than reduced.Comment: 12 pages, 6 figures, 21 reference

Intermediate energy Coulomb excitation as a probe of nuclear structure at radioactive beam facilities

The effects of retardation in the Coulomb excitation of radioactive nuclei in intermediate energy collisions (Elab ~100 MeV/nucleon) are investigated. We show that the excitation cross sections of low-lying states in 11Be, {38,40,42}S and {44,46}Ar projectiles incident on gold and lead targets are modified by as much as 20% due to these effects. The angular distributions of decaying gamma-rays are also appreciably modified.Comment: 21 pages, 3 figures, Phys. Rev. C, in pres

Stacking Entropy of Hard Sphere Crystals

Classical hard spheres crystallize at equilibrium at high enough density. Crystals made up of stackings of 2-dimensional hexagonal close-packed layers (e.g. fcc, hcp, etc.) differ in entropy by only about $10^{-3}k_B$ per sphere (all configurations are degenerate in energy). To readily resolve and study these small entropy differences, we have implemented two different multicanonical Monte Carlo algorithms that allow direct equilibration between crystals with different stacking sequences. Recent work had demonstrated that the fcc stacking has higher entropy than the hcp stacking. We have studied other stackings to demonstrate that the fcc stacking does indeed have the highest entropy of ALL possible stackings. The entropic interactions we could detect involve three, four and (although with less statistical certainty) five consecutive layers of spheres. These interlayer entropic interactions fall off in strength with increasing distance, as expected; this fall-off appears to be much slower near the melting density than at the maximum (close-packing) density. At maximum density the entropy difference between fcc and hcp stackings is $0.00115 +/- 0.00004 k_B$ per sphere, which is roughly 30% higher than the same quantity measured near the melting transition.Comment: 15 page

Event-Driven Molecular Dynamics in Parallel

Although event-driven algorithms have been shown to be far more efficient than time-driven methods such as conventional molecular dynamics, they have not become as popular. The main obstacle seems to be the difficulty of parallelizing event-driven molecular dynamics. Several basic ideas have been discussed in recent years, but to our knowledge no complete implementation has been published yet. In this paper we present a parallel event-driven algorithm including dynamic load-balancing, which can be easily implemented on any computer architecture. To simplify matters our explanations refer to a basic multi-particle system of hard spheres, but can be extended easily to a wide variety of possible models.Comment: 10 pages, 9 figure

Semiclassical treatment of fusion processes in collisions of weakly bound nuclei

We describe a semiclassical treatment of nuclear fusion reactions involving weakly bound nuclei. In this treatment, the complete fusion probabilities are approximated by products of two factors: a tunneling probability and the probability that the system is in its ground state at the strong absorption radius. We investigate the validity of the method in a schematic two-channel application, where the channels in the continuum are represented by a single resonant state. Comparisons with full coupled-channels calculations are performed. The agreement between semiclassical and quantal calculations isquite good, suggesting that the procedure may be extended to more sophisticated discretizations of the continuum.Comment: 11 pages, 5 figure

Effect of boundary conditions on diffusion in two-dimensional granular gases

We analyze the influence of boundary conditions on numerical simulations of the diffusive properties of a two dimensional granular gas. We show in particular that periodic boundary conditions introduce unphysical correlations in time which cause the coefficient of diffusion to be strongly dependent on the system size. On the other hand, in large enough systems with hard walls at the boundaries, diffusion is found to be independent of the system size. We compare the results obtained in this case with Langevin theory for an elastic gas. Good agreement is found. We then calculate the relaxation time and the influence of the mass for a particle of radius $R_s$ in a sea of particles of radius $R_b$. As granular gases are dissipative, we also study the influence of an external random force on the diffusion process in a forced dissipative system. In particular, we analyze differences in the mean square velocity and displacement between the elastic and inelastic cases.Comment: 15 figures eps figures, include

Possible Origin of RHIC R_{out}/R_{sid} HBT Results

The effects of opacity of the nuclei together with a blackbody type of emission along the system history are considered as a means to explain the ratio $R_{out}/R_{sid}$ observed by STAR and PHENIX collaborations at RHIC. Within our model, no flow is required to explain the data trend of this ratio for large surface emissivities.Comment: 4 pages, 2 Postscript figures, uses espcrc1.sty. Talk presented at Quark Matter 2002, Nantes, France, July 18-24, 2002; to appear in the proceedings (Nucl. Phys. A