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

Coulomb Excitation of Multi-Phonon Levels of the Giant Dipole Resonance

A closed expression is obtained for the cross-section for Coulomb excitation of levels of the giant dipole resonance of given angular momentum and phonon number. Applications are made to the Goldhaber-Teller and Steinwedel-Jensen descriptions of the resonance, at non-relativistic and relativistic bombarding energies.Comment: 16 pages, 5 figure

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

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

Quantal Brownian Motion - Dephasing and Dissipation

We analyze quantal Brownian motion in $d$ dimensions using the unified model for diffusion localization and dissipation, and Feynman-Vernon formalism. At high temperatures the propagator possess a Markovian property and we can write down an equivalent Master equation. Unlike the case of the Zwanzig-Caldeira-Leggett model, genuine quantum mechanical effects manifest themselves due to the disordered nature of the environment. Using Wigner picture of the dynamics we distinguish between two different mechanisms for destruction of coherence. The analysis of dephasing is extended to the low temperature regime by using a semiclassical strategy. Various results are derived for ballistic, chaotic, diffusive, both ergodic and non-ergodic motion. We also analyze loss of coherence at the limit of zero temperature and clarify the limitations of the semiclassical approach. The condition for having coherent effect due to scattering by low-frequency fluctuations is also pointed out. It is interesting that the dephasing rate can be either larger or smaller than the dissipation rate, depending on the physical circumstances.Comment: LaTex, 23 pages, 4 figures, published vesio

Free energies of crystalline solids: a lattice-switch Monte Carlo method

We present a method for the direct evaluation of the difference between the free energies of two crystalline structures, of different symmetry. The method rests on a Monte Carlo procedure which allows one to sample along a path, through atomic-displacement-space, leading from one structure to the other by way of an intervening transformation that switches one set of lattice vectors for another. The configurations of both structures can thus be sampled within a single Monte Carlo process, and the difference between their free energies evaluated directly from the ratio of the measured probabilities of each. The method is used to determine the difference between the free energies of the fcc and hcp crystalline phases of a system of hard spheres.Comment: 5 pages Revtex, 3 figure

Accuracy of B(E2; 0+ -> 2+) transition rates from intermediate-energy Coulomb excitation experiments

The method of intermediate-energy Coulomb excitation has been widely used to determine absolute B(E2; 0+ -> 2+) quadrupole excitation strengths in exotic nuclei with even numbers of protons and neutrons. Transition rates measured with intermediate-energy Coulomb excitation are compared to their respective adopted values and for the example of 26Mg to the B(E2; 0+ -> 2+) values obtained with a variety of standard methods. Intermediate-energy Coulomb excitation is found to have an accuracy comparable to those of long-established experimental techniques.Comment: to be published in Phys. Rev.

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