Statistical analysis of coherent structures in transitional pipe flow

Numerical and experimental studies of transitional pipe flow have shown the prevalence of coherent flow structures that are dominated by downstream vortices. They attract special attention because they contribute predominantly to the increase of the Reynolds stresses in turbulent flow. In the present study we introduce a convenient detector for these coherent states, calculate the fraction of time the structures appear in the flow, and present a Markov model for the transition between the structures. The fraction of states that show vortical structures exceeds 24% for a Reynolds number of about Re=2200, and it decreases to about 20% for Re=2500. The Markov model for the transition between these states is in good agreement with the observed fraction of states, and in reasonable agreement with the prediction for their persistence. It provides insight into dominant qualitative changes of the flow when increasing the Reynolds number.Comment: 11 pages, 26 (sub)figure

Collective Coordinate Control of Density Distributions

Real collective density variables $C(\boldsymbol{k})$ [c.f. Eq.\ref{Equation3})] in many-particle systems arise from non-linear transformations of particle positions, and determine the structure factor $S(\boldsymbol{k})$, where $\bf k$ denotes the wave vector. Our objective is to prescribe $C({\boldsymbol k})$ and then to find many-particle configurations that correspond to such a target $C({\bf k})$ using a numerical optimization technique. Numerical results reported here extend earlier one- and two-dimensional studies to include three dimensions. In addition, they demonstrate the capacity to control $S(\boldsymbol{k})$ in the neighborhood of $|\boldsymbol{k}| =$ 0. The optimization method employed generates multi-particle configurations for which $S(\boldsymbol{k}) \propto |\boldsymbol{k}|^{\alpha}$, $|\boldsymbol{k}| \leq K$, and $\alpha =$ 1, 2, 4, 6, 8, and 10. The case $\alpha =$ 1 is relevant for the Harrison-Zeldovich model of the early universe, for superfluid $^{4}{He}$, and for jammed amorphous sphere packings. The analysis also provides specific examples of interaction potentials whose classical ground state are configurationally degenerate and disordered.Comment: 26 pages, 8 figure

Competition of the connectivity with the local and the global order in polymer melts and crystals

The competition between the connectivity and the local or global order in model fully-flexible chain molecules is investigated by molecular-dynamics simulations. States with both missing (melts) and high (crystal) global order are considered. Local order is characterized within the first coordination shell (FCS) of a tagged monomer and found to be lower than in atomic systems in both melt and crystal. The role played by the bonds linking the tagged monomer to FCS monomers (radial bonds), and the bonds linking two FCS monomers (shell bonds) is investigated. The detailed analysis in terms of Steinhardt's orientation order parameters Q_l (l = 2 - 10) reveals that increasing the number of shell bonds decreases the FCS order in both melt and crystal. Differently, the FCS arrangements organize the radial bonds. Even if the molecular chains are fully flexible, the distribution of the angle formed by adjacent radial bonds exhibits sharp contributions at the characteristic angles {\theta} = 70{\deg}, 122{\deg}, 180{\deg}. The fractions of adjacent radial bonds with {\theta} = 122{\deg}, 180{\deg} are enhanced by the global order of the crystal, whereas the fraction with 70{\deg} < {\theta} < 110{\deg} is nearly unaffected by the crystallization. Kink defects, i.e. large lateral displacements of the chains, are evidenced in the crystalline state.Comment: J. Chem. Phys. in pres

Quantum computations with atoms in optical lattices: marker qubits and molecular interactions

We develop a scheme for quantum computation with neutral atoms, based on the concept of "marker" atoms, i.e., auxiliary atoms that can be efficiently transported in state-independent periodic external traps to operate quantum gates between physically distant qubits. This allows for relaxing a number of experimental constraints for quantum computation with neutral atoms in microscopic potential, including single-atom laser addressability. We discuss the advantages of this approach in a concrete physical scenario involving molecular interactions.Comment: 15 pages, 14 figure

Critical exponents of a three dimensional O(4) spin model

By Monte Carlo simulation we study the critical exponents governing the transition of the three-dimensional classical O(4) Heisenberg model, which is considered to be in the same universality class as the finite-temperature QCD with massless two flavors. We use the single cluster algorithm and the histogram reweighting technique to obtain observables at the critical temperature. After estimating an accurate value of the inverse critical temperature \Kc=0.9360(1), we make non-perturbative estimates for various critical exponents by finite-size scaling analysis. They are in excellent agreement with those obtained with the $4-\epsilon$ expansion method with errors reduced to about halves of them.Comment: 25 pages with 8 PS figures, LaTeX, UTHEP-28

Conditions for one-dimensional supersonic flow of quantum gases

One can use transsonic Bose-Einstein condensates of alkali atoms to establish the laboratory analog of the event horizon and to measure the acoustic version of Hawking radiation. We determine the conditions for supersonic flow and the Hawking temperature for realistic condensates on waveguides where an external potential plays the role of a supersonic nozzle. The transition to supersonic speed occurs at the potential maximum and the Hawking temperature is entirely determined by the curvature of the potential

Wind reversals in turbulent Rayleigh-Benard convection

The phenomenon of irregular cessation and subsequent reversal of the large-scale circulation in turbulent Rayleigh-B\'enard convection is theoretically analysed. The force and thermal balance on a single plume detached from the thermal boundary layer yields a set of coupled nonlinear equations, whose dynamics is related to the Lorenz equations. For Prandtl and Rayleigh numbers in the range $10^{-2} \leq \Pr \leq 10^{3}$ and 10^{7} \leq \Ra \leq 10^{12}, the model has the following features: (i) chaotic reversals may be exhibited at Ra $\geq 10^{7}$; (ii) the Reynolds number based on the root mean square velocity scales as \Re_{rms} \sim \Ra^{[0.41 ... 0.47]} (depending on Pr), and as $\Re_{rms} \sim \Pr^{-[0.66 ... 0.76]}$ (depending on Ra); and (iii) the mean reversal frequency follows an effective scaling law \omega / (\nu L^{-2}) \sim \Pr^{-(0.64 \pm 0.01)} \Ra^{0.44 \pm 0.01}. The phase diagram of the model is sketched, and the observed transitions are discussed.Comment: 4 pages, 5 figure

Infinite qubit rings with maximal nearest neighbor entanglement: the Bethe ansatz solution

We search for translationally invariant states of qubits on a ring that maximize the nearest neighbor entanglement. This problem was initially studied by O'Connor and Wootters [Phys. Rev. A {\bf 63}, 052302 (2001)]. We first map the problem to the search for the ground state of a spin 1/2 Heisenberg XXZ model. Using the exact Bethe ansatz solution in the limit of an infinite ring, we prove the correctness of the assumption of O'Connor and Wootters that the state of maximal entanglement does not have any pair of neighboring spins ``down'' (or, alternatively spins ``up''). For sufficiently small fixed magnetization, however, the assumption does not hold: we identify the region of magnetizations for which the states that maximize the nearest neighbor entanglement necessarily contain pairs of neighboring spins ``down''.Comment: 10 pages, 4 figures; Eq. (45) and Fig. 3 corrected, no qualitative change in conclusion

Numerical investigation of black hole interiors

Gravitational perturbations which are present in any realistic stellar collapse to a black hole, die off in the exterior of the hole, but experience an infinite blueshift in the interior. This is believed to lead to a slowly contracting lightlike scalar curvature singularity, characterized by a divergence of the hole's (quasi-local) mass function along the inner horizon. The region near the inner horizon is described to great accuracy by a plane wave spacetime. While Einstein's equations for this metric are still too complicated to be solved in closed form it is relatively simple to integrate them numerically. We find for generic regular initial data the predicted mass inflation type null singularity, rather than a spacelike singularity. It thus seems that mass inflation indeed represents a generic self-consistent picture of the black hole interior.Comment: 6 pages LaTeX, 3 eps figure

Kolmogorov spectrum of superfluid turbulence: numerical analysis of the Gross-Pitaevskii equation with the small scale dissipation

The energy spectrum of superfluid turbulence is studied numerically by solving the Gross-Pitaevskii equation. We introduce the dissipation term which works only in the scale smaller than the healing length, to remove short wavelength excitations which may hinder the cascade process of quantized vortices in the inertial range. The obtained energy spectrum is consistent with the Kolmogorov law.Comment: 4 pages, 4 figures and 1 table. Submitted to American Journal of Physic