Exact and Truncated Dynamics in Nonequilibrium Field Theory

Nonperturbative dynamics of quantum fields out of equilibrium is often described by the time evolution of a hierarchy of correlation functions, using approximation methods such as Hartree, large N, and nPI-effective action techniques. These truncation schemes can be implemented equally well in a classical statistical system, where results can be tested by comparison with the complete nonlinear evolution obtained by numerical methods. For a 1+1 dimensional scalar field we find that the early-time behaviour is reproduced qualitatively by the Hartree dynamics. The inclusion of direct scattering improves this to the quantitative level. We show that the emergence of nonthermal temperature profiles at intermediate times can be understood in terms of the fixed points of the evolution equations in the Hartree approximation. The form of the profile depends explicitly on the initial ensemble. While the truncated evolution equations do not seem to be able to get away from the fixed point, the full nonlinear evolution shows thermalization with a (surprisingly) slow relaxation.Comment: 30 pages with 12 eps figures, minor changes; to appear in Phys.Rev.

Covariant transport approach for strongly interacting partonic systems

The dynamics of partons, hadrons and strings in relativistic nucleus-nucleus collisions is analyzed within the novel Parton-Hadron-String Dynamics (PHSD) transport approach, which is based on a dynamical quasiparticle model for partons (DQPM) matched to reproduce recent lattice-QCD results - including the partonic equation of state - in thermodynamic equilibrium. Scalar- and vector-interaction densities are extracted from the DQPM as well as effective scalar- and vector-mean fields for the partons. The transition from partonic to hadronic degrees of freedom is described by covariant transition rates for the fusion of quark-antiquark pairs or three quarks (antiquarks), respectively, obeying flavor current-conservation, color neutrality as well as energy-momentum conservation. Since the dynamical quarks and antiquarks become very massive close to the phase transition, the formed resonant 'pre-hadronic' color-dipole states ($q\bar{q}$ or $qqq$) are of high invariant mass, too, and sequentially decay to the groundstate meson and baryon octets increasing the total entropy. When applying the PHSD approach to Pb+Pb colllisions at 158 A$\cdot$GeV we find a significant effect of the partonic phase on the production of multi-strange antibaryons due to a slightly enhanced $s{\bar s}$ pair production from massive time-like gluon decay and a larger formation of antibaryons in the hadronization process.Comment: 12 pages, 6 figures, to be published in the Proceedings of the 26th Winter Workshop on `Nuclear Dynamics', Ochto Rios, Jamaica, 2-9 January, 2010

Measurement of the spatial extent of inverse proximity in a Py/Nb/Py superconducting trilayer using low-energy muon-spin rotation

The authors acknowledge the financial support of the EPSRC (Grant No. EP/J01060X).Muon-spin rotation has been used to observe directly the spatial variation of the magnetic flux density near the ferromagnetic-superconducting interface in a permalloy-niobium trilayer. Above the superconducting transition temperature Tc the profile of the induced magnetic flux density within the niobium layer has been determined. Below Tc there is a significant reduction of the induced flux density, predominantly near the ferromagnetic-superconducting interfaces. We are uniquely able to determine the magnitude and spatial variation of this reduction in induced magnetization due to the presence of the Cooper pairs, yielding the magnitude and length scale associated with this phenomenon. Both are inconsistent with a simple Meissner screening and indicate the existence of another mechanism, the influence of which is localized within the vicinity of the ferromagnetic interface.Publisher PDFPeer reviewe

Continuous extremal optimization for Lennard-Jones Clusters

In this paper, we explore a general-purpose heuristic algorithm for finding high-quality solutions to continuous optimization problems. The method, called continuous extremal optimization(CEO), can be considered as an extension of extremal optimization(EO) and is consisted of two components, one is with responsibility for global searching and the other is with responsibility for local searching. With only one adjustable parameter, the CEO's performance proves competitive with more elaborate stochastic optimization procedures. We demonstrate it on a well known continuous optimization problem: the Lennerd-Jones clusters optimization problem.Comment: 5 pages and 3 figure

Convergence of simulated annealing by the generalized transition probability

We prove weak ergodicity of the inhomogeneous Markov process generated by the generalized transition probability of Tsallis and Stariolo under power-law decay of the temperature. We thus have a mathematical foundation to conjecture convergence of simulated annealing processes with the generalized transition probability to the minimum of the cost function. An explicitly solvable example in one dimension is analyzed in which the generalized transition probability leads to a fast convergence of the cost function to the optimal value. We also investigate how far our arguments depend upon the specific form of the generalized transition probability proposed by Tsallis and Stariolo. It is shown that a few requirements on analyticity of the transition probability are sufficient to assure fast convergence in the case of the solvable model in one dimension.Comment: 11 page

Convergence theorems for quantum annealing

We prove several theorems to give sufficient conditions for convergence of quantum annealing, which is a protocol to solve generic optimization problems by quantum dynamics. In particular the property of strong ergodicity is proved for the path-integral Monte Carlo implementation of quantum annealing for the transverse Ising model under a power decay of the transverse field. This result is to be compared with the much slower inverse-log decay of temperature in the conventional simulated annealing. Similar results are proved for the Green's function Monte Carlo approach. Optimization problems in continuous space of particle configurations are also discussed.Comment: 19 page

Optimal transport on supply-demand networks

Previously, transport networks are usually treated as homogeneous networks, that is, every node has the same function, simultaneously providing and requiring resources. However, some real networks, such as power grid and supply chain networks, show a far different scenario in which the nodes are classified into two categories: the supply nodes provide some kinds of services, while the demand nodes require them. In this paper, we propose a general transport model for those supply-demand networks, associated with a criterion to quantify their transport capacities. In a supply-demand network with heterogenous degree distribution, its transport capacity strongly depends on the locations of supply nodes. We therefore design a simulated annealing algorithm to find the optimal configuration of supply nodes, which remarkably enhances the transport capacity, and outperforms the degree target algorithm, the betweenness target algorithm, and the greedy method. This work provides a start point for systematically analyzing and optimizing transport dynamics on supply-demand networks.Comment: 5 pages, 1 table and 4 figure

Nonequilibrium evolution of Phi**4 theory in 1+1 dimensions in the 2PPI formalism

We consider the out-of-equilibrium evolution of a classical condensate field and its quantum fluctuations for a Phi**4 model in 1+1 dimensions with a symmetric and a double well potential. We use the 2PPI formalism and go beyond the Hartree approximation by including the sunset term. In addition to the mean field phi= the 2PPI formalism uses as variational parameter a time dependent mass M**2(t) which contains all local insertions into the Green function. We compare our results to those obtained in the Hartree approximation. In the symmetric Phi**4 theory we observe that the mean field shows a stronger dissipation than the one found in the Hartree approximation. The dissipation is roughly exponential in an intermediate time region. In the theory with spontaneous symmetry breaking, i.e., with a double well potential, the field amplitude tends to zero, i.e., to the symmetric configuration. This is expected on general grounds: in 1+1 dimensional quantum field theory there is no spontaneous symmetry breaking for T >0, and so there should be none at finite energy density (microcanonical ensemble), either. Within the time range of our simulations the momentum spectra do not thermalize and display parametric resonance bands.Comment: 14 pages, 18 encapsulated postscript figures; v2 minor changes, new appendix, accepted for publication in Phys.Rev.

An oil pipeline design problem

Copyright @ 2003 INFORMSWe consider a given set of offshore platforms and onshore wells producing known (or estimated) amounts of oil to be connected to a port. Connections may take place directly between platforms, well sites, and the port, or may go through connection points at given locations. The configuration of the network and sizes of pipes used must be chosen to minimize construction costs. This problem is expressed as a mixed-integer program, and solved both heuristically by Tabu Search and Variable Neighborhood Search methods and exactly by a branch-and-bound method. Two new types of valid inequalities are introduced. Tests are made with data from the South Gabon oil field and randomly generated problems.The work of the first author was supported by NSERC grant #OGP205041. The work of the second author was supported by FCAR (Fonds pour la Formation des Chercheurs et l’Aide à la Recherche) grant #95-ER-1048, and NSERC grant #GP0105574

Annealing schedule from population dynamics

We introduce a dynamical annealing schedule for population-based optimization algorithms with mutation. On the basis of a statistical mechanics formulation of the population dynamics, the mutation rate adapts to a value maximizing expected rewards at each time step. Thereby, the mutation rate is eliminated as a free parameter from the algorithm.Comment: 6 pages RevTeX, 4 figures PostScript; to be published in Phys. Rev.