Continuous and discrete Clebsch variational principles

The Clebsch method provides a unifying approach for deriving variational principles for continuous and discrete dynamical systems where elements of a vector space are used to control dynamics on the cotangent bundle of a Lie group \emph{via} a velocity map. This paper proves a reduction theorem which states that the canonical variables on the Lie group can be eliminated, if and only if the velocity map is a Lie algebra action, thereby producing the Euler-Poincar\'e (EP) equation for the vector space variables. In this case, the map from the canonical variables on the Lie group to the vector space is the standard momentum map defined using the diamond operator. We apply the Clebsch method in examples of the rotating rigid body and the incompressible Euler equations. Along the way, we explain how singular solutions of the EP equation for the diffeomorphism group (EPDiff) arise as momentum maps in the Clebsch approach. In the case of finite dimensional Lie groups, the Clebsch variational principle is discretised to produce a variational integrator for the dynamical system. We obtain a discrete map from which the variables on the cotangent bundle of a Lie group may be eliminated to produce a discrete EP equation for elements of the vector space. We give an integrator for the rotating rigid body as an example. We also briefly discuss how to discretise infinite-dimensional Clebsch systems, so as to produce conservative numerical methods for fluid dynamics

The Hamiltonian structure and Euler-Poincar\'{e} formulation of the Vlasov-Maxwell and gyrokinetic systems

We present a new variational principle for the gyrokinetic system, similar to the Maxwell-Vlasov action presented in Ref. 1. The variational principle is in the Eulerian frame and based on constrained variations of the phase space fluid velocity and particle distribution function. Using a Legendre transform, we explicitly derive the field theoretic Hamiltonian structure of the system. This is carried out with a modified Dirac theory of constraints, which is used to construct meaningful brackets from those obtained directly from Euler-Poincar\'{e} theory. Possible applications of these formulations include continuum geometric integration techniques, large-eddy simulation models and Casimir type stability methods. [1] H. Cendra et. al., Journal of Mathematical Physics 39, 3138 (1998)Comment: 36 pages, 1 figur

Ground state of two unlike charged colloids: An analogy with ionic bonding

In this letter, we study the ground state of two spherical macroions of identical radius, but asymmetric bare charge ((Q_{A}>Q_{B})). Electroneutrality of the system is insured by the presence of the surrounding divalent counterions. Using Molecular Dynamics simulations within the framework of the primitive model, we show that the ground state of such a system consists of an overcharged and an undercharged colloid. For a given macroion separation the stability of these ionized-like states is a function of the difference ((\sqrt{N_{A}}-\sqrt{N_{B}})) of neutralizing counterions (N_{A}) and (N_{B}). Furthermore the degree of ionization, or equivalently, the degree of overcharging, is also governed by the distance separation of the macroions. The natural analogy with ionic bonding is briefly discussed.Comment: published versio

The 1981 outburst of the old nova GK Persei

Old nova GK Per was observed in 1981 with the IUE, during its rise, maximum, and subsequent return to minimum. In outburst, GK Per is luminous but much redder than dwarf novae or standard model accretion disks. The observed spectrum can be explained qualitatively with the Ghosh and Lamb (1979) model for the interaction of an accretion disk with the magnetic field of the accreting white dwarf. The N V and He2 are enhanced relative to other emission lines during outburst. This can be understood with photoionization by very soft X-rays having a luminosity comparable to that of the hard X-rays

Random Hamiltonian in thermal equilibrium

A framework for the investigation of disordered quantum systems in thermal equilibrium is proposed. The approach is based on a dynamical model--which consists of a combination of a double-bracket gradient flow and a uniform Brownian fluctuation--that `equilibrates' the Hamiltonian into a canonical distribution. The resulting equilibrium state is used to calculate quenched and annealed averages of quantum observables.Comment: 8 pages, 4 figures. To appear in DICE 2008 conference proceeding

Competing superconducting and magnetic order parameters and field-induced magnetism in electron doped Ba(Fe1−x_{1-x}Cox_{x})2_{2}As2_{2}

We have studied the magnetic and superconducting properties of Ba(Fe$_{0.95}$Co$_{0.05}$)$_{2}$As$_{2}$ as a function of temperature and external magnetic field using neutron scattering and muon spin rotation. Below the superconducting transition temperature the magnetic and superconducting order parameters coexist and compete. A magnetic field can significantly enhance the magnetic scattering in the superconducting state, roughly doubling the Bragg intensity at 13.5 T. We perform a microscopic modelling of the data by use of a five-band Hamiltonian relevant to iron pnictides. In the superconducting state, vortices can slow down and freeze spin fluctuations locally. When such regions couple they result in a long-range ordered antiferromagnetic phase producing the enhanced magnetic elastic scattering in agreement with experiments.Comment: 9 pages, 6 figure

Onsager-Manning-Oosawa condensation phenomenon and the effect of salt

Making use of results pertaining to Painleve III type equations, we revisit the celebrated Onsager-Manning-Oosawa condensation phenomenon for charged stiff linear polymers, in the mean-field approximation with salt. We obtain analytically the associated critical line charge density, and show that it is severely affected by finite salt effects, whereas previous results focused on the no salt limit. In addition, we obtain explicit expressions for the condensate thickness and the electric potential. The case of asymmetric electrolytes is also briefly addressed.Comment: to appear in Phys. Rev. Let

Ion-ion correlations: an improved one-component plasma correction

Based on a Debye-Hueckel approach to the one-component plasma we propose a new free energy for incorporating ionic correlations into Poisson-Boltzmann like theories. Its derivation employs the exclusion of the charged background in the vicinity of the central ion, thereby yielding a thermodynamically stable free energy density, applicable within a local density approximation. This is an improvement over the existing Debye-Hueckel plus hole theory, which in this situation suffers from a "structuring catastrophe". For the simple example of a strongly charged stiff rod surrounded by its counterions we demonstrate that the Poisson-Boltzmann free energy functional augmented by our new correction accounts for the correlations present in this system when compared to molecular dynamics simulations.Comment: 5 pages, 2 figures, revtex styl

Leray and LANS-α\alpha modeling of turbulent mixing

Mathematical regularisation of the nonlinear terms in the Navier-Stokes equations provides a systematic approach to deriving subgrid closures for numerical simulations of turbulent flow. By construction, these subgrid closures imply existence and uniqueness of strong solutions to the corresponding modelled system of equations. We will consider the large eddy interpretation of two such mathematical regularisation principles, i.e., Leray and LANS$-\alpha$ regularisation. The Leray principle introduces a {\bfi smoothed transport velocity} as part of the regularised convective nonlinearity. The LANS$-\alpha$ principle extends the Leray formulation in a natural way in which a {\bfi filtered Kelvin circulation theorem}, incorporating the smoothed transport velocity, is explicitly satisfied. These regularisation principles give rise to implied subgrid closures which will be applied in large eddy simulation of turbulent mixing. Comparison with filtered direct numerical simulation data, and with predictions obtained from popular dynamic eddy-viscosity modelling, shows that these mathematical regularisation models are considerably more accurate, at a lower computational cost.Comment: 42 pages, 12 figure