Adiabatic hyperspherical study of triatomic helium systems

The 4He3 system is studied using the adiabatic hyperspherical representation. We adopt the current state-of-the-art helium interaction potential including retardation and the nonadditive three-body term, to calculate all low-energy properties of the triatomic 4He system. The bound state energies of the 4He trimer are computed as well as the 4He+4He2 elastic scattering cross sections, the three-body recombination and collision induced dissociation rates at finite temperatures. We also treat the system that consists of two 4He and one 3He atoms, and compute the spectrum of the isotopic trimer 4He2 3He, the 3He+4He2 elastic scattering cross sections, the rates for three-body recombination and the collision induced dissociation rate at finite temperatures. The effects of retardation and the nonadditive three-body term are investigated. Retardation is found to be significant in some cases, while the three-body term plays only a minor role for these systems.Comment: 24 pages 6 figures Submitted to Physical Review

An order parameter equation for the dynamic yield stress in dense colloidal suspensions

We study the dynamic yield stress in dense colloidal suspensions by analyzing the time evolution of the pair distribution function for colloidal particles interacting through a Lennard-Jones potential. We find that the equilibrium pair distribution function is unstable with respect to a certain anisotropic perturbation in the regime of low temperature and high density. By applying a bifurcation analysis to a system near the critical state at which the stability changes, we derive an amplitude equation for the critical mode. This equation is analogous to order parameter equations used to describe phase transitions. It is found that this amplitude equation describes the appearance of the dynamic yield stress, and it gives a value of 2/3 for the shear thinning exponent. This value is related to the mean field value of the critical exponent $\delta$ in the Ising model.Comment: 8 pages, 2 figure

Cold three-body collisions in hydrogen-hydrogen-alkali atomic system

We have studied hydrogen-hydrogen-alkali three-body systems in the adiabatic hyperspherical representation. For the spin-stretched case, there exists a single $X$H molecular state when $X$ is one of the bosonic alkali atoms: $^7$Li, $^{23}$Na, $^{39}$K, $^{87}$Rb and $^{133}$Cs. As a result, the {\em only} recombination process is the one that leads to formation of $X$H molecules, H+H+$X$$\rightarrow$$X$H+H, and such molecules will be stable against vibrational relaxation. We have calculated the collision rates for recombination and collision induced dissociation as well as the elastic cross-sections for H+$X$H collisions up to a temperature of 0.5 K, including the partial wave contributions from $J^\Pi$=$0^+$ to $5^-$. We have also found that there is just one three-body bound state for such systems for $J^\Pi$=$0^+$ and no bound states for higher angular momenta.Comment: 10 pages, 5 figures, 4 table

Diffusive transport in networks built of containers and tubes

We developed analytical and numerical methods to study a transport of non-interacting particles in large networks consisting of M d-dimensional containers C_1,...,C_M with radii R_i linked together by tubes of length l_{ij} and radii a_{ij} where i,j=1,2,...,M. Tubes may join directly with each other forming junctions. It is possible that some links are absent. Instead of solving the diffusion equation for the full problem we formulated an approach that is computationally more efficient. We derived a set of rate equations that govern the time dependence of the number of particles in each container N_1(t),N_2(t),...,N_M(t). In such a way the complicated transport problem is reduced to a set of M first order integro-differential equations in time, which can be solved efficiently by the algorithm presented here. The workings of the method have been demonstrated on a couple of examples: networks involving three, four and seven containers, and one network with a three-point junction. Already simple networks with relatively few containers exhibit interesting transport behavior. For example, we showed that it is possible to adjust the geometry of the networks so that the particle concentration varies in time in a wave-like manner. Such behavior deviates from simple exponential growth and decay occurring in the two container system.Comment: 21 pages, 18 figures, REVTEX4; new figure added, reduced emphasis on graph theory, additional discussion added (computational cost, one dimensional tubes

Highly accurate calculations of the rotationally excited bound states in three-body systems

An effective optimization strategy has been developed to construct highly accurate bound state wave functions in various three-body systems. Our procedure appears to be very effective for computations of weakly bound states and various excited states, including rotationally excited states, i.e. states with $L \ge 1$. The efficiency of our procedure is illustrated by computations of the excited $P^{*}(L = 1)-$states in the $dd\mu, dt\mu$ and $tt\mu$ muonic molecular ions, $P(L = 1)-$states in the non-symmetric $pd\mu, pt\mu$ and $dt\mu$ ions and $2^{1}P(L = 1)-$ and $2^{3}P(L = 1)-$states in He atom(s)

Noise-induced perturbations of dispersion-managed solitons

We study noise-induced perturbations of dispersion-managed solitons by developing soliton perturbation theory for the dispersion-managed nonlinear Schroedinger (DMNLS) equation, which governs the long-term behavior of optical fiber transmission systems and certain kinds of femtosecond lasers. We show that the eigenmodes and generalized eigenmodes of the linearized DMNLS equation around traveling-wave solutions can be generated from the invariances of the DMNLS equations, we quantify the perturbation-induced parameter changes of the solution in terms of the eigenmodes and the adjoint eigenmodes, and we obtain evolution equations for the solution parameters. We then apply these results to guide importance-sampled Monte-Carlo simulations and reconstruct the probability density functions of the solution parameters under the effect of noise.Comment: 12 pages, 6 figure

Approximate action-angle variables for the figure-eight and other periodic three-body orbits

We use the maximally permutation symmetric set of three-body coordinates, that consist of the "hyper-radius" $R = \sqrt{\rho^{2} + \lambda^{2}}$, the "rescaled area of the triangle" $\frac{\sqrt 3}{2 R^2} |{\bm \rho} \times {\bm \lambda}|$) and the (braiding) hyper-angle $\phi = \arctan(\frac{2{\bm \rho} \cdot {\bm \lambda}}{\lambda^2 - \rho^2})$, to analyze the "figure-eight" choreographic three-body motion discovered by Moore \cite{Moore1993} in the Newtonian three-body problem. Here ${\bm \rho}, {\bm \lambda}$ are the two Jacobi relative coordinate vectors. We show that the periodicity of this motion is closely related to the braiding hyper-angle $\phi$. We construct an approximate integral of motion ${\bar{G}}$ that together with the hyper-angle $\phi$ forms the action-angle pair of variables for this problem and show that it is the underlying cause of figure-eight motion's stability. We construct figure-eight orbits in two other attractive permutation-symmetric three-body potentials. We compare the figure-eight orbits in these three potentials and discuss their generic features, as well as their differences. We apply these variables to two new periodic, but non-choreographic orbits: One has a continuously rising $\phi$ in time $t$, just like the figure-eight motion, but with a different, more complex periodicity, whereas the other one has an oscillating $\phi(t)$ temporal behavior.Comment: 11 pages, 19 figure

Existence and Stability of Standing Pulses in Neural Networks : I Existence

We consider the existence of standing pulse solutions of a neural network integro-differential equation. These pulses are bistable with the zero state and may be an analogue for short term memory in the brain. The network consists of a single-layer of neurons synaptically connected by lateral inhibition. Our work extends the classic Amari result by considering a non-saturating gain function. We consider a specific connectivity function where the existence conditions for single-pulses can be reduced to the solution of an algebraic system. In addition to the two localized pulse solutions found by Amari, we find that three or more pulses can coexist. We also show the existence of nonconvex ``dimpled'' pulses and double pulses. We map out the pulse shapes and maximum firing rates for different connection weights and gain functions.Comment: 31 pages, 29 figures, submitted to SIAM Journal on Applied Dynamical System

Numerical Approximations Using Chebyshev Polynomial Expansions

We present numerical solutions for differential equations by expanding the unknown function in terms of Chebyshev polynomials and solving a system of linear equations directly for the values of the function at the extrema (or zeros) of the Chebyshev polynomial of order N (El-gendi's method). The solutions are exact at these points, apart from round-off computer errors and the convergence of other numerical methods used in connection to solving the linear system of equations. Applications to initial value problems in time-dependent quantum field theory, and second order boundary value problems in fluid dynamics are presented.Comment: minor wording changes, some typos have been eliminate

Geometry and symmetries of multi-particle systems

The quantum dynamical evolution of atomic and molecular aggregates, from their compact to their fragmented states, is parametrized by a single collective radial parameter. Treating all the remaining particle coordinates in d dimensions democratically, as a set of angles orthogonal to this collective radius or by equivalent variables, bypasses all independent-particle approximations. The invariance of the total kinetic energy under arbitrary d-dimensional transformations which preserve the radial parameter gives rise to novel quantum numbers and ladder operators interconnecting its eigenstates at each value of the radial parameter. We develop the systematics and technology of this approach, introducing the relevant mathematics tutorially, by analogy to the familiar theory of angular momentum in three dimensions. The angular basis functions so obtained are treated in a manifestly coordinate-free manner, thus serving as a flexible generalized basis for carrying out detailed studies of wavefunction evolution in multi-particle systems.Comment: 37 pages, 2 eps figure