A Noisy Monte Carlo Algorithm

We propose a Monte Carlo algorithm to promote Kennedy and Kuti's linear accept/reject algorithm which accommodates unbiased stochastic estimates of the probability to an exact one. This is achieved by adopting the Metropolis accept/reject steps for both the dynamical and noise configurations. We test it on the five state model and obtain desirable results even for the case with large noise. We also discuss its application to lattice QCD with stochastically estimated fermion determinants.Comment: 10 pages, 1 tabl

The Kentucky Noisy Monte Carlo Algorithm for Wilson Dynamical Fermions

We develop an implementation for a recently proposed Noisy Monte Carlo approach to the simulation of lattice QCD with dynamical fermions by incorporating the full fermion determinant directly. Our algorithm uses a quenched gauge field update with a shifted gauge coupling to minimize fluctuations in the trace log of the Wilson Dirac matrix. The details of tuning the gauge coupling shift as well as results for the distribution of noisy estimators in our implementation are given. We present data for some basic observables from the noisy method, as well as acceptance rate information and discuss potential autocorrelation and sign violation effects. Both the results and the efficiency of the algorithm are compared against those of Hybrid Monte Carlo. PACS Numbers: 12.38.Gc, 11.15.Ha, 02.70.Uu Keywords: Noisy Monte Carlo, Lattice QCD, Determinant, Finite Density, QCDSPComment: 30 pages, 6 figure

Finite Density Algorithm in Lattice QCD -- a Canonical Ensemble Approach

I will review the finite density algorithm for lattice QCD based on finite chemical potential and summarize the associated difficulties. I will propose a canonical ensemble approach which projects out the finite baryon number sector from the fermion determinant. For this algorithm to work, it requires an efficient method for calculating the fermion determinant and a Monte Carlo algorithm which accommodates unbiased estimate of the probability. I shall report on the progress made along this direction with the Pad\'{e} - Z$_2$ estimator of the determinant and its implementation in the newly developed Noisy Monte Carlo algorithm.Comment: Invited talk at Nankai Symposium on Mathematical Physics, Tianjin, Oct. 2001, 18 pages, 3 figures; expanded and references adde

On the parabolic equation method in internal wave propagation

A parabolic equation for the propagation of periodic internal waves over varying bottom topography is derived using the multiple-scale perturbation method. Some computational aspects of the numerical implementation are discussed. The results of numerical experiments on propagation of an incident plane wave over a circular-type shoal are presented in comparison with the analytical result, based on Born approximation.Comment: Submitted to Coastal Engineering. 16 pages, 5 figures. One figure was excluded from article because of size problem

Observed photodetachment in parallel electric and magnetic fields

We investigate photodetachment from negative ions in a homogeneous 1.0-T magnetic field and a parallel electric field of approximately 10 V/cm. A theoretical model for detachment in combined fields is presented. Calculations show that a field of 10 V/cm or more should considerably diminish the Landau structure in the detachment cross section. The ions are produced and stored in a Penning ion trap and illuminated by a single-mode dye laser. We present preliminary results for detachment from S- showing qualitative agreement with the model. Future directions of the work are also discussed.Comment: Nine pages, five figures, minor revisions showing final publicatio

Rotational and Vibrational Dynamics of Interstitial Molecular Hydrogen

The calculation of the hindered roton-phonon energy levels of a hydrogen molecule in a confining potential with different symmetries is systematized for the case when the rotational angular momentum $J$ is a good quantum number. One goal of this program is to interpret the energy-resolved neutron time of flight spectrum previously obtained for H$_{2}$C$_{60}$. This spectrum gives direct information on the energy level spectrum of H$_2$ molecules confined to the octahedral interstitial sites of solid C$_{60}$. We treat this problem of coupled translational and orientational degrees of freedom a) by construction of an effective Hamiltonian to describe the splitting of the manifold of states characterized by a given value of $J$ and having a fixed total number of phonon excitations, b) by numerical solutions of the coupled translation-rotation problem on a discrete mesh of points in position space, and c) by a group theoretical symmetry analysis. Results obtained from these three different approaches are mutually consistent. The results of our calculations explain several hitherto uninterpreted aspects of the experimental observations, but show that a truly satisfactory orientational potential for the interaction of an H$_2$ molecule with a surrounding array of C atoms has not yet been developed.Comment: 53 pages, 9 figures, to appear in Phys. Rev B (in press). Phys. Rev. B (in press

Hydrodynamics of domain growth in nematic liquid crystals

We study the growth of aligned domains in nematic liquid crystals. Results are obtained solving the Beris-Edwards equations of motion using the lattice Boltzmann approach. Spatial anisotropy in the domain growth is shown to be a consequence of the flow induced by the changing order parameter field (backflow). The generalization of the results to the growth of a cylindrical domain, which involves the dynamics of a defect ring, is discussed.Comment: 12 revtex-style pages, including 12 figures; small changes before publicatio

Optical activity of neutrinos and antineutrinos

Using the one-loop helicity amplitudes for low-energy $\nu\gamma\to\nu\gamma$ and $\bar\nu\gamma\to\bar\nu\gamma$ scattering in the standard model with massless neutrinos, we study the optical activity of a sea of neutrinos and antineutrinos. In particular, we estimate the values of the index of refraction and rotary power of this medium in the absence of dispersion.Comment: Additional reference

Beauty is Attractive: Moduli Trapping at Enhanced Symmetry Points

We study quantum effects on moduli dynamics arising from the production of particles which are light at special points in moduli space. The resulting forces trap the moduli at these points, which often exhibit enhanced symmetry. Moduli trapping occurs in time-dependent quantum field theory, as well as in systems of moving D-branes, where it leads the branes to combine into stacks. Trapping also occurs in an expanding universe, though the range over which the moduli can roll is limited by Hubble friction. We observe that a scalar field trapped on a steep potential can induce a stage of acceleration of the universe, which we call trapped inflation. Moduli trapping ameliorates the cosmological moduli problem and may affect vacuum selection. In particular, rolling moduli are most powerfully attracted to the points with the largest number of light particles, which are often the points of greatest symmetry. Given suitable assumptions about the dynamics of the very early universe, this effect might help to explain why among the plethora of possible vacuum states of string theory, we appear to live in one with a large number of light particles and (spontaneously broken) symmetries. In other words, some of the surprising properties of our world might arise not through pure chance or miraculous cancellations, but through a natural selection mechanism during dynamical evolution.Comment: 50 pages, 4 figures; v2: added references and an appendix describing a related classical proces