Lattice Boltzmann scheme for relativistic fluids

A Lattice Boltzmann formulation for relativistic fluids is presented and numerically verified through quantitative comparison with recent hydrodynamic simulations of relativistic shock-wave propagation in viscous quark-gluon plasmas. This formulation opens up the possibility of exporting the main advantages of Lattice Boltzmann methods to the relativistic context, which seems particularly useful for the simulation of relativistic fluids in complicated geometries.Comment: Submitted to PR

Higher Order Methods for Simulations on Quantum Computers

To efficiently implement many-qubit gates for use in quantum simulations on quantum computers we develop and present methods reexpressing exp[-i (H_1 + H_2 + ...) \Delta t] as a product of factors exp[-i H_1 \Delta t], exp[-i H_2 \Delta t], ... which is accurate to 3rd or 4th order in \Delta t. The methods we derive are an extended form of symplectic method and can also be used for the integration of classical Hamiltonians on classical computers. We derive both integral and irrational methods, and find the most efficient methods in both cases.Comment: 21 pages, Latex, one figur

Derivation of the Lattice Boltzmann Model for Relativistic Hydrodynamics

A detailed derivation of the Lattice Boltzmann (LB) scheme for relativistic fluids recently proposed in Ref. [1], is presented. The method is numerically validated and applied to the case of two quite different relativistic fluid dynamic problems, namely shock-wave propagation in quark-gluon plasmas and the impact of a supernova blast-wave on massive interstellar clouds. Close to second order convergence with the grid resolution, as well as linear dependence of computational time on the number of grid points and time-steps, are reported

Jeans' gravitational instability and nonextensive kinetic theory

The concept of Jeans gravitational instability is rediscussed in the framework of nonextensive statistics and its associated kinetic theory. A simple analytical formula generalizing the Jeans criterion is derived by assuming that the unperturbed self- gravitating collisionless gas is kinetically described by the $q$-parameterized class of power law velocity distributions. It is found that the critical values of wavelength and mass depend explicitly on the nonextensive $q$-parameter. The standard Jeans wavelength derived for a Maxwellian distribution is recovered in the limiting case $q$=1. For power-law distributions with cutoff, the instability condition is weakened with the system becoming unstable even for wavelengths of the disturbance smaller than the standard Jeans length $\lambda_J$.Comment: 5 pages, including 3 figures. Accepted for publication in A&

A quantum algorithm providing exponential speed increase for finding eigenvalues and eigenvectors

We describe a new polynomial time quantum algorithm that uses the quantum fast fourier transform to find eigenvalues and eigenvectors of a Hamiltonian operator, and that can be applied in cases (commonly found in ab initio physics and chemistry problems) for which all known classical algorithms require exponential time. Applications of the algorithm to specific problems are considered, and we find that classically intractable and interesting problems from atomic physics may be solved with between 50 and 100 quantum bits.Comment: 10 page

Two-dimensional hydrodynamic lattice-gas simulations of binary immiscible and ternary amphiphilic fluid flow through porous media

The behaviour of two dimensional binary and ternary amphiphilic fluids under flow conditions is investigated using a hydrodynamic lattice gas model. After the validation of the model in simple cases (Poiseuille flow, Darcy's law for single component fluids), attention is focussed on the properties of binary immiscible fluids in porous media. An extension of Darcy's law which explicitly admits a viscous coupling between the fluids is verified, and evidence of capillary effects are described. The influence of a third component, namely surfactant, is studied in the same context. Invasion simulations have also been performed. The effect of the applied force on the invasion process is reported. As the forcing level increases, the invasion process becomes faster and the residual oil saturation decreases. The introduction of surfactant in the invading phase during imbibition produces new phenomena, including emulsification and micellisation. At very low fluid forcing levels, this leads to the production of a low-resistance gel, which then slows down the progress of the invading fluid. At long times (beyond the water percolation threshold), the concentration of remaining oil within the porous medium is lowered by the action of surfactant, thus enhancing oil recovery. On the other hand, the introduction of surfactant in the invading phase during drainage simulations slows down the invasion process -- the invading fluid takes a more tortuous path to invade the porous medium -- and reduces the oil recovery (the residual oil saturation increases).Comment: 48 pages, 26 figures. Phys. Rev. E (in press

Thermodynamics' 0-th-Law in a nonextensive scenario

Tsallis' thermostatistics is by now recognized as a new paradigm for statistical mechanical considerations. However, it is still affected by a serious hindrance: the generalization of thermodynamics' zero-th law to a nonextensive scenario is plagued by difficulties. Here we show how to overcome this problem.Comment: 4 pages, latex; added references for section

Correlations and Renormalization in Lattice Gases

A complete formulation is given of an exact kinetic theory for lattice gases. This kinetic theory makes possible the calculation of corrections to the usual Boltzmann / Chapman-Enskog analysis of lattice gases due to the buildup of correlations. It is shown that renormalized transport coefficients can be calculated perturbatively by summing terms in an infinite series. A diagrammatic notation for the terms in this series is given, in analogy with the diagrammatic expansions of continuum kinetic theory and quantum field theory. A closed-form expression for the coefficients associated with the vertices of these diagrams is given. This method is applied to several standard lattice gases, and the results are shown to correctly predict experimentally observed deviations from the Boltzmann analysis.Comment: 94 pages, pure LaTeX including all figure

Ideal gas in nonextensive optimal Lagrange multipliers formalism

Based on the prescription termed the optimal Lagrange multipliers formalism for extremizing the Tsallis entropy indexed by q, it is shown that key aspects of the treatment of the ideal gas problem are identical in both the nonextensive and extensive cases.Comment: 5 pages, no figure