Systematic Derivation of Amplitude Equations and Normal Forms for Dynamical Systems

We present a systematic approach to deriving normal forms and related amplitude equations for flows and discrete dynamics on the center manifold of a dynamical system at local bifurcations and unfoldings of these. We derive a general, explicit recurrence relation that completely determines the amplitude equation and the associated transformation from amplitudes to physical space. At any order, the relation provides explicit expressions for all the nonvanishing coefficients of the amplitude equation together with straightforward linear equations for the coefficients of the transformation. The recurrence relation therefore provides all the machinery needed to solve a given physical problem in physical terms through an amplitude equation. The new result applies to any local bifurcation of a flow or map for which all the critical eigenvalues are semisimple i.e. have Riesz index unity). The method is an efficient and rigorous alternative to more intuitive approaches in terms of multiple time scales. We illustrate the use of the method by deriving amplitude equations and associated transformations for the most common simple bifurcations in flows and iterated maps. The results are expressed in tables in a form that can be immediately applied to specific problems.Comment: 40 pages, 1 figure, 4 tables. Submitted to Chaos. Please address any correspondence by email to [email protected]

Radial sine-Gordon kinks as sources of fast breathers

We consider radial sine-Gordon kinks in two, three and higher dimensions. A full two dimensional simulation showing that azimuthal perturbations remain small allows to reduce the problem to the one dimensional radial sine-Gordon equation. We solve this equation on an interval $[r_0,r_1]$ and absorb all outgoing radiation. Before collision the kink is well described by a simple law derived from the conservation of energy. In two dimensions for $r_0 \le 2$, the collision disintegrates the kink into a fast breather while for $r_0 \ge 4$ we obtain a kink-breather meta-stable state where breathers are shed at each kink "return". In three and higher dimensions $d$ a kink-pulson state appears for small $r_0$. The three states then exist as shown by a study of the $(d,r_0)$ parameter space. On the application side, the kink disintegration opens the way for new types of terahertz microwave generators

Control of Multi-level Voltage States in a Hysteretic SQUID Ring-Resonator System

In this paper we study numerical solutions to the quasi-classical equations of motion for a SQUID ring-radio frequency (rf) resonator system in the regime where the ring is highly hysteretic. In line with experiment, we show that for a suitable choice of of ring circuit parameters the solutions to these equations of motion comprise sets of levels in the rf voltage-current dynamics of the coupled system. We further demonstrate that transitions, both up and down, between these levels can be controlled by voltage pulses applied to the system, thus opening up the possibility of high order (e.g. 10 state), multi-level logic and memory.Comment: 8 pages, 9 figure

Switching between dynamic states in intermediate-length Josephson junctions

The appearance of zero-field steps (ZFS’s) in the current-voltage characteristics of intermediate-length overlap-geometry Josephson tunnel junctions described by a perturbed sine-Gordon equation (PSGE) is associated with the growth of parametrically excited instabilities of the McCumber background curve (MCB). A linear stability analysis of a McCumber solution of the PSGE in the asymptotic linear region of the MCB and in the absence of magnetic field yields a Hill’s equation which predicts how the number, locations, and widths of the instability regions depend on the junction parameters. A numerical integration of the PSGE in terms of truncated series of time-dependent Fourier spatial modes verifies that the parametrically excited instabilities of the MCB evolve into the fluxon oscillations characteristic of the ZFS’s. An approximate analysis of the Fourier mode equations in the presence of a small magnetic field yields a field-dependent Hill’s equation which predicts that the major effect of such a field is to reduce the widths of the instability regions. Experimental measurements on Nb-NbxOy-Pb junctions of intermediate length, performed at different operating temperatures in order to vary the junction parameters and for various magnetic field values, verify the physical existence of switching from the MCB to the ZFS’s. Good qualitative, and in many cases quantitative, agreement between analytic, numerical, and experimental results is obtained

Close binary evolution. III. Impact of tides, wind magnetic braking, and internal angular momentum transport

Massive stars with solar metallicity lose important amounts of rotational angular momentum through their winds. When a magnetic field is present at the surface of a star, efficient angular momentum losses can still be achieved even when the mass-loss rate is very modest, at lower metallicities, or for lower-initial-mass stars. In a close binary system, the effect of wind magnetic braking also interacts with the influence of tides, resulting in a complex evolution of rotation. We study the interactions between the process of wind magnetic braking and tides in close binary systems. We discuss the evolution of a 10 M$_\odot$ star in a close binary system with a 7 M$_\odot$ companion using the Geneva stellar evolution code. The initial orbital period is 1.2 days. The 10 M$_\odot$ star has a surface magnetic field of 1 kG. Various initial rotations are considered. We use two different approaches for the internal angular momentum transport. In one of them, angular momentum is transported by shear and meridional currents. In the other, a strong internal magnetic field imposes nearly perfect solid-body rotation. The evolution of the primary is computed until the first mass-transfer episode occurs. The cases of different values for the magnetic fields and for various orbital periods and mass ratios are briefly discussed. We show that, independently of the initial rotation rate of the primary and the efficiency of the internal angular momentum transport, the surface rotation of the primary will converge, in a time that is short with respect to the main-sequence lifetime, towards a slowly evolving velocity that is different from the synchronization velocity. (abridged).Comment: 11 pages, 13 figures, accepted for publication in Astronomy and Astrophysic

Kink propagation in a two-dimensional curved Josephson junction

We consider the propagation of sine-Gordon kinks in a planar curved strip as a model of nonlinear wave propagation in curved wave guides. The homogeneous Neumann transverse boundary conditions, in the curvilinear coordinates, allow to assume a homogeneous kink solution. Using a simple collective variable approach based on the kink coordinate, we show that curved regions act as potential barriers for the wave and determine the threshold velocity for the kink to cross. The analysis is confirmed by numerical solution of the 2D sine-Gordon equation.Comment: 8 pages, 4 figures (2 in color

An Improved Global Model for Air-Sea Exchange of Mercury: High Concentrations over the North Atlantic

We develop an improved treatment of the surface ocean in the GEOS-Chem global 3-D biogeochemical model for mercury (Hg). We replace the globally uniform subsurface ocean Hg concentrations used in the original model with basin-specific values based on measurements. Updated chemical mechanisms for Hg0/HgII redox reactions in the surface ocean include both photochemical and biological processes, and we improved the parametrization of particle-associated Hg scavenging. Modeled aqueous Hg concentrations are consistent with limited surface water observations. Results more accurately reproduce high-observed MBL concentrations over the North Atlantic (NA) and the associated seasonal trends. High seasonal evasion in the NA is driven by inputs from Hg enriched subsurface waters through entrainment and Ekman pumping. Globally, subsurface waters account for 40% of Hg inputs to the ocean mixed layer, and 60% is from atmospheric deposition. Although globally the ocean is a net sink for 3.8 Mmol Hg y−1, the NA is a net source to the atmosphere, potentially due to enrichment of subsurface waters with legacy Hg from historical anthropogenic sources.Engineering and Applied Science

Multi-Decadal Decline of Mercury in the North Atlantic Atmosphere Explained by Changing Subsurface Seawater Concentrations

[1] We analyze 1977–2010 trends in atmospheric mercury (Hg) from 21 ship cruises over the North Atlantic (NA) and 15 over the South Atlantic (SA). We find a steep 1990–2009 decline of −0.046 ± 0.010 ng m−3 a−1 (−2.5% a−1) over the NA (steeper than at Northern Hemispheric land sites) but no significant decline over the SA. Surface water Hg0 measurements in the NA show a decline of −5.7% a−1since 1999, and limited subsurface ocean data show an ∼80% decline from 1980 to present. We use a coupled global atmosphere-ocean model to show that the decline in NA atmospheric concentrations can be explained by decreasing oceanic evasion from the NA driven by declining subsurface water Hg concentrations. We speculate that this large historical decline of Hg in the NA Ocean could have been caused by decreasing Hg inputs from rivers and wastewater and by changes in the oxidant chemistry of the atmospheric marine boundary layer.Engineering and Applied Science