Higher order finite difference schemes for the magnetic induction equations

We describe high order accurate and stable finite difference schemes for the initial-boundary value problem associated with the magnetic induction equations. These equations model the evolution of a magnetic field due to a given velocity field. The finite difference schemes are based on Summation by Parts (SBP) operators for spatial derivatives and a Simultaneous Approximation Term (SAT) technique for imposing boundary conditions. We present various numerical experiments that demonstrate both the stability as well as high order of accuracy of the schemes.Comment: 20 page

AQFT from n-functorial QFT

There are essentially two different approaches to the axiomatization of quantum field theory (QFT): algebraic QFT, going back to Haag and Kastler, and functorial QFT, going back to Atiyah and Segal. More recently, based on ideas by Baez and Dolan, the latter is being refined to "extended" functorial QFT by Freed, Hopkins, Lurie and others. The first approach uses local nets of operator algebras which assign to each patch an algebra "of observables", the latter uses n-functors which assign to each patch a "propagator of states". In this note we present an observation about how these two axiom systems are naturally related: we demonstrate under mild assumptions that every 2-dimensional extended Minkowskian QFT 2-functor ("parallel surface transport") naturally yields a local net. This is obtained by postcomposing the propagation 2-functor with an operation that mimics the passage from the Schroedinger picture to the Heisenberg picture in quantum mechanics. The argument has a straightforward generalization to general pseudo-Riemannian structure and higher dimensions.Comment: 39 pages; further examples added: Hopf spin chains and asymptotic inclusion of subfactors; references adde

A lattice model for the kinetics of rupture of fluid bilayer membranes

We have constructed a model for the kinetics of rupture of membranes under tension, applying physical principles relevant to lipid bilayers held together by hydrophobic interactions. The membrane is characterized by the bulk compressibility (for expansion), the thickness of the hydrophobic part of the bilayer, the hydrophobicity and a parameter characterizing the tail rigidity of the lipids. The model is a lattice model which incorporates strain relaxation, and considers the nucleation of pores at constant area, constant temperature, and constant particle number. The particle number is conserved by allowing multiple occupancy of the sites. An equilibrium ``phase diagram'' is constructed as a function of temperature and strain with the total pore surface and distribution as the order parameters. A first order rupture line is found with increasing tension, and a continuous increase in proto-pore concentration with rising temperature till instability. The model explains current results on saturated and unsaturated PC lipid bilayers and thicker artificial bilayers made of diblock copolymers. Pore size distributions are presented for various values of area expansion and temperature, and the fractal dimension of the pore edge is evaluated.Comment: 15 pages, 8 figure

Equilibrium Properties of A Monomer-Monomer Catalytic Reaction on A One-Dimensional Chain

We study the equilibrium properties of a lattice-gas model of an $A + B \to 0$ catalytic reaction on a one-dimensional chain in contact with a reservoir for the particles. The particles of species $A$ and $B$ are in thermal contact with their vapor phases acting as reservoirs, i.e., they may adsorb onto empty lattice sites and may desorb from the lattice. If adsorbed $A$ and $B$ particles appear at neighboring lattice sites they instantaneously react and both desorb. For this model of a catalytic reaction in the adsorption-controlled limit, we derive analytically the expression of the pressure and present exact results for the mean densities of particles and for the compressibilities of the adsorbate as function of the chemical potentials of the two species.Comment: 19 pages, 5 figures, submitted to Phys. Rev.

Adsorption of Reactive Particles on a Random Catalytic Chain: An Exact Solution

We study equilibrium properties of a catalytically-activated annihilation $A + A \to 0$ reaction taking place on a one-dimensional chain of length $N$ ($N \to \infty$) in which some segments (placed at random, with mean concentration $p$) possess special, catalytic properties. Annihilation reaction takes place, as soon as any two $A$ particles land onto two vacant sites at the extremities of the catalytic segment, or when any $A$ particle lands onto a vacant site on a catalytic segment while the site at the other extremity of this segment is already occupied by another $A$ particle. Non-catalytic segments are inert with respect to reaction and here two adsorbed $A$ particles harmlessly coexist. For both "annealed" and "quenched" disorder in placement of the catalytic segments, we calculate exactly the disorder-average pressure per site. Explicit asymptotic formulae for the particle mean density and the compressibility are also presented.Comment: AMSTeX, 27 pages + 4 figure

Statistical distribution of quantum entanglement for a random bipartite state

We compute analytically the statistics of the Renyi and von Neumann entropies (standard measures of entanglement), for a random pure state in a large bipartite quantum system. The full probability distribution is computed by first mapping the problem to a random matrix model and then using a Coulomb gas method. We identify three different regimes in the entropy distribution, which correspond to two phase transitions in the associated Coulomb gas. The two critical points correspond to sudden changes in the shape of the Coulomb charge density: the appearance of an integrable singularity at the origin for the first critical point, and the detachement of the rightmost charge (largest eigenvalue) from the sea of the other charges at the second critical point. Analytical results are verified by Monte Carlo numerical simulations. A short account of some of these results appeared recently in Phys. Rev. Lett. {\bf 104}, 110501 (2010).Comment: 7 figure

Self-Similar Scalar Field Collapse: Naked Singularities and Critical Behaviour

Homothetic scalar field collapse is considered in this article. By making a suitable choice of variables the equations are reduced to an autonomous system. Then using a combination of numerical and analytic techniques it is shown that there are two classes of solutions. The first consists of solutions with a non-singular origin in which the scalar field collapses and disperses again. There is a singularity at one point of these solutions, however it is not visible to observers at finite radius. The second class of solutions includes both black holes and naked singularities with a critical evolution (which is neither) interpolating between these two extremes. The properties of these solutions are discussed in detail. The paper also contains some speculation about the significance of self-similarity in recent numerical studies.Comment: 27 pages including 5 encapsulated postcript figures in separate compressed file, report NCL94-TP1

Is weak temperature dependence of electron dephasing possible?

The first-principle theory of electron dephasing by disorder-induced two state fluctuators is developed. There exist two mechanisms of dephasing. First, dephasing occurs due to direct transitions between the defect levels caused by inelastic electron-defect scattering. The second mechanism is due to violation of the time reversal symmetry caused by time-dependent fluctuations of the scattering potential. These fluctuations originate from an interaction between the dynamic defects and conduction electrons forming a thermal bath. The first contribution to the dephasing rate saturates as temperature decreases. The second contribution does not saturate, although its temperature dependence is rather weak, $\propto T^{1/3}$. The quantitative estimates based on the experimental data show that these mechanisms considered can explain the weak temperature dependence of the dephasing rate in some temperature interval. However, below some temperature dependent on the model of dynamic defects the dephasing rate tends rapidly to zero. The relation to earlier studies of the dephasing caused by the dynamical defects is discussed.Comment: 14 pages, 6 figures, submitted to PR

Dark mammoth trunks in the merging galaxy NGC 1316 and a mechanism of cosmic double helices

NGC 1316 is a giant, elliptical galaxy containing a complex network of dark, dust features. The morphology of these features has been examined in some detail using a Hubble Space Telescope, Advanced Camera for Surveys image. It is found that most of the features are constituted of long filaments. There also exist a great number of dark structures protruding inwards from the filaments. Many of these structures are strikingly similar to elephant trunks in H II regions in the Milky Way Galaxy, although much larger. The structures, termed mammoth trunks, generally are filamentary and often have shapes resembling the letters V or Y. In some of the mammoth trunks the stem of the Y can be resolved into two or more filaments, many of which showing signs of being intertwined. A model of the mammoth trunks, related to a recent theory of elephant trunks, is proposed. Based on magnetized filaments, the model is capable of giving an account of the various shapes of the mammoth trunks observed, including the twined structures.Comment: Accepted for publication in Astrophysics & Space Scienc