Finding Exponential Product Formulas of Higher Orders

In the present article, we review a continual effort on generalization of the Trotter formula to higher-order exponential product formulas. The exponential product formula is a good and useful approximant, particularly because it conserves important symmetries of the system dynamics. We focuse on two algorithms of constructing higher-order exponential product formulas. The first is the fractal decomposition, where we construct higher-order formulas recursively. The second is to make use of the quantum analysis, where we compute higher-order correction terms directly. As interludes, we also have described the decomposition of symplectic integrators, the approximation of time-ordered exponentials, and the perturbational composition.Comment: 22 pages, 9 figures. To be published in the conference proceedings ''Quantum Annealing and Other Optimization Methods," eds. B.K.Chakrabarti and A.Das (Springer, Heidelberg

Higher Order Force Gradient Symplectic Algorithms

We show that a recently discovered fourth order symplectic algorithm, which requires one evaluation of force gradient in addition to three evaluations of the force, when iterated to higher order, yielded algorithms that are far superior to similarly iterated higher order algorithms based on the standard Forest-Ruth algorithm. We gauge the accuracy of each algorithm by comparing the step-size independent error functions associated with energy conservation and the rotation of the Laplace-Runge-Lenz vector when solving a highly eccentric Kepler problem. For orders 6, 8, 10 and 12, the new algorithms are approximately a factor of $10^3$, $10^4$, $10^4$ and $10^5$ better.Comment: 23 pages, 10 figure

A Compact Gas Cerenkov Detector with Novel Optics

We discuss the design and performance of a threshold Cerenkov counter for identification of charged hadrons. The radiator is pressurized gas, which is contained in thin-walled cylindrical modules. A mirror system of novel design transports Cerenkov photons to photomultiplier tubes. This system is compact, contains relatively little material, and has a large fraction of active volume. A prototype of a module designed for the proposed CLEO III detector has been studied using cosmic rays. Results from these studies show good agreement with a detailed Monte Carlo simulation of the module and indicate that it should achieve separation of pions and kaons at the 2.5-3.0sigma level in the momentum range 0.8-2.8 GeV/c. We predict performance for specific physics analyses using a GEANT-based simulation package.Comment: Submitted to NIM. 23 pages, 11 postscript figures. Postscript file is also available at http://w4.lns.cornell.edu/public/CLNS/199

Causality and statistics on the Groenewold-Moyal plane

Quantum theories constructed on the noncommutative spacetime called the Groenewold-Moyal plane exhibit many interesting properties such as Lorentz and CPT noninvariance, causality violation and twisted statistics. We show that such violations lead to many striking features that may be tested experimentally. These theories predict Pauli forbidden transitions due to twisted statistics, anisotropies in the cosmic microwave background radiation due to correlations of observables in spacelike regions and Lorentz and CPT violations in scattering amplitudes.Comment: 12 pages, 1 figure. Based on the talk given by APB at the Workshop "Theoretical and Experimental Aspects of the Spin Statisics Connection and Related Symmetries", Stazione Marittima Conference Center, Trieste, Italy from the 21st to the 25th of October 200

QCD Sum Rules for the production of the X(3872) as a mixed molecule-charmonium state in B meson decay

We use QCD sum rules to calculate the branching ratio for the production of the meson X(3872) in the decay $B\to X(3872)K$, assumed to be a mixture between charmonium and exotic molecular $[c\bar{q}][q\bar{c}]$ states with $J^{PC}=1^{++}$. We find that in a small range for the values of the mixing angle, $5^\circ\leq\theta\leq13^\circ$, we get the branching ratio ${\mathcal{B}}(B\to XK)=(1.00\pm0.68)\times10^{-5}$, which is in agreement with the experimental upper limit. This result is compatible with the analysis of the mass and decay width of the mode $J/\psi(n\pi)$ and the radiative decay mode $J/\psi\gamma$ performed in the same approach.Comment: 6 pages, 3 figures; revised versions to appear on Phys. Lett.

Random Series and Discrete Path Integral methods: The Levy-Ciesielski implementation

We perform a thorough analysis of the relationship between discrete and series representation path integral methods, which are the main numerical techniques used in connection with the Feynman-Kac formula. First, a new interpretation of the so-called standard discrete path integral methods is derived by direct discretization of the Feynman-Kac formula. Second, we consider a particular random series technique based upon the Levy-Ciesielski representation of the Brownian bridge and analyze its main implementations, namely the primitive, the partial averaging, and the reweighted versions. It is shown that the n=2^k-1 subsequence of each of these methods can also be interpreted as a discrete path integral method with appropriate short-time approximations. We therefore establish a direct connection between the discrete and the random series approaches. In the end, we give sharp estimates on the rates of convergence of the partial averaging and the reweighted Levy-Ciesielski random series approach for sufficiently smooth potentials. The asymptotic rates of convergence are found to be O(1/n^2), in agreement with the rates of convergence of the best standard discrete path integral techniques.Comment: 20 pages, 4 figures; the two equations before Eq. 14 are corrected; other typos are remove

Promoter prediction using physico-chemical properties of DNA

The ability to locate promoters within a section of DNA is known to be a very difficult and very important task in DNA analysis. We document an approach that incorporates the concept of DNA as a complex molecule using several models of its physico-chemical properties. A support vector machine is trained to recognise promoters by their distinctive physical and chemical properties. We demonstrate that by combining models, we can improve upon the classification accuracy obtained with a single model. We also show that by examining how the predictive accuracy of these properties varies over the promoter, we can reduce the number of attributes needed. Finally, we apply this method to a real-world problem. The results demonstrate that such an approach has significant merit in its own right. Furthermore, they suggest better results from a planned combined approach to promoter prediction using both physicochemical and sequence based techniques

Experimental Studies of Magnetically Driven Plasma Jets

We present experimental results on the formation of supersonic, radiatively cooled jets driven by pressure due to the toroidal magnetic field generated by the 1.5 MA, 250 ns current from the MAGPIE generator. The morphology of the jet produced in the experiments is relevant to astrophysical jet scenarios in which a jet on the axis of a magnetic cavity is collimated by a toroidal magnetic field as it expands into the ambient medium. The jets in the experiments have similar Mach number, plasma beta and cooling parameter to those in protostellar jets. Additionally the Reynolds, magnetic Reynolds and Peclet numbers are much larger than unity, allowing the experiments to be scaled to astrophysical flows. The experimental configuration allows for the generation of episodic magnetic cavities, suggesting that periodic fluctuations near the source may be responsible for some of the variability observed in astrophysical jets. Preliminary measurements of kinetic, magnetic and Poynting energy of the jets in our experiments are presented and discussed, together with estimates of their temperature and trapped toroidal magnetic field.Comment: 7 pages, 6 figures, accepted for publication in Astrophysics & Space Scienc

From QFT to DCC

A quantum field theoretical model for the dynamics of the disoriented chiral condensate is presented. A unified approach to relate the quantum field theory directly to the formation, decay and signals of the DCC and its evolution is taken. We use a background field analysis of the O(4) sigma model keeping one-loop quantum corrections (quadratic order in the fluctuations). An evolution of the quantum fluctuations in an external, expanding metric which simulates the expansion of the plasma, is carried out. We examine, in detail, the amplification of the low momentum pion modes with two competing effects, the expansion rate of the plasma and the transition rate of the vacuum configuration from a metastable state into a stable state.We show the effect of DCC formation on the multiplicity distributions and the Bose-Einstein correlations.Comment: 34 pages, 10 figure

Determination of Omega_b From Big Bang Nucleosynthesis in the Presence of Regions of Antimatter

Production of regions of antimatter in the early universe is predicted in many baryogenesis models. Small scale antimatter regions would annihilate during or soon after nucleosynthesis, affecting the abundances of the light elements. In this paper we study how the acceptable range in Omega_b changes in the presence of antimatter regions, as compared to the standard big bang nucleosynthesis. It turns out that it is possible to produce at the same time both a low 4He value (Y_p < 0.240) and a low D/H value (D/H < 4e-5), but overproduction of 7Li is unavoidable at large Omega_b.Comment: 9 pages, PRD version, ref. 6 correcte