New Analytical Solution of the Equilibrium Ampere’s Law Using the Walker’s Method: a Didactic Example

This work aims to demonstrate the analytical solution of the Grad-Shafranov (GS) equation or generalized Ampere’s law, which is important in the studies of self-consistent 2.5-D solution for current sheet structures. A detailed mathematical development is presented to obtain the generating function as shown by Walker (RSPSA 91, 410, 1915). Therefore, we study the general solution of the GS equation in terms of the Walker’s generating function in details without omitting any step. The Walker’s generating function g(ζ) is written in a new way as the tangent of an unspecified function K(ζ). In this trend, the general solution of the GS equation is expressed as exp(??2Ψ) =?4|K ^′(ζ)|²/cos²[K(ζ) ? K(ζ ^?)]. In order to investigate whether our proposal would simplify the mathematical effort to find new generating functions, we use Harris’s solution as a test, in this case K(ζ) = arctan(exp(i ζ)). In summary, one of the article purposes is to present a review of the Harris’s solution. In an attempt to find a simplified solution, we propose a new way to write the GS solution using g(ζ) = tan(K(ζ)). We also present a new analytical solution to the equilibrium Ampere’s law using g(ζ) = cosh(b ζ), which includes a generalization of the Harris model and presents isolated magnetic islands.     

The derivative nonlinear Schrödinger equation on the half-line

We analyze the derivative nonlinear Schrödinger equation on the half-line using the Fokas method. Assuming that the solution q(x,t) exists, we show that it can be represented in terms of the solution of a matrix Riemann-Hilbert problem formulated in the plane of the complex spectral parameter ζ. The jump matrix has explicit x,t dependence and is given in terms of the spectral functions a(ζ), b(ζ) (obtained from the initial data q₀(x)=q(x,0)) as well as A(ζ), B(ζ) (obtained from the boundary values g₀(t)=q(0,t) and g₁(t)=q_x(0,t)). The spectral functions are not independent, but related by a compatibility condition, the so-called global relation. Given initial and boundary values {q₀(x),g₀(t),g₁(t)} such that there exist spectral functions satisfying the global relation, we show that the function q(x,t) defined by the above Riemann-Hilbert problem exists globally and solves the derivative nonlinear Schrödinger equation with the prescribed initial and boundary values.     

Resolvent Expansion and Time Decay of the Wave Functions for Two-Dimensional Magnetic Schrödinger Operators

We consider two-dimensional Schrödinger operators H(B, V) given by Eq. (1.1) below. We prove that, under certain regularity and decay assumptions on B and V, the character of the expansion for the resolvent (H(B, V) ? λ)^?1 as λ → 0 is determined by the flux of the magnetic field B through \({\mathbb{R}^2}\) . Subsequently, we derive the leading term of the asymptotic expansion of the unitary group e ^{?i t H(B, V)} as t → ∞ and show how the magnetic field improves its decay in t with respect to the decay of the unitary group e ^{?i t H(0, V)}.     

Flux Fluctuations in the One Dimensional Nearest Neighbors Symmetric Simple Exclusion Process

Let J(t) be the the integrated flux of particles in the symmetric simple exclusion process starting with the product invariant measure ν _ρ with density ρ. We compute its rescaled asymptotic variance: $$\mathop {\lim }\limits_{t \to \infty } t^{ - 1/2} \mathbb{V}J(t) = \sqrt {2/\pi } (1 - \rho )\rho$$ Furthermore we show that t ^?1/4 J(t) converges weakly to a centered normal random variable with this variance. From these results we compute the asymptotic variance of a tagged particle in the nearest neighbor case and show the corresponding central limit theorem.     

Non-Markovian diffusion equations and processes: Analysis and simulations

In this paper we introduce and analyze a class of diffusion type equations related to certain non-Markovian stochastic processes. We start from the forward drift equation which is made non-local in time by the introduction of a suitable chosen memory kernel K(t). The resulting non-Markovian equation can be interpreted in a natural way as the evolution equation of the marginal density function of a random time process l(t). We then consider the subordinated process Y(t)=X(l(t)) where X(t) is a Markovian diffusion. The corresponding time evolution of the marginal density function of Y(t) is governed by a non-Markovian Fokker-Planck equation which involves the memory kernel K(t). We develop several applications and derive the exact solutions. We consider different stochastic models for the given equations providing path simulations.     

Nonlinearity associated with the motion of a Bloch wall-onset of a chaotic regime and analysis of the Barkhausen noise

The noise V(t) generated by the displacement of a magnetic domain wall is investigated. As the velocity of the wall increases. the trajectories in the (V, V) plane exhibit a transition from quasi-periodic to chaotic motion.     

Excitation of the Classical Electromagnetic Field in a Cavity Containing a Thin Slab with a Time-Dependent Conductivity

We derive an exact infinite set of coupled ordinary differential equations describing the evolution of the modes of the classical electromagnetic field inside an ideal cavity containing a thin slab with the time-dependent conductivity σ(t) and dielectric permittivity ε(t) for the dispersion-less media. We analyze this problem in connection with the attempts to simulate the so-called dynamical Casimir effect in three-dimensional electromagnetic cavities containing a thin semiconductor slab periodically illuminated by strong laser pulses. Therefore, we assume that functions σ(t) and δε(t) = ε(t) ? ε(0) are different from zero during short time intervals (pulses) only. Our main goal here is to find the conditions under which the initial nonzero classical field could be amplified after a single pulse (or a series of pulses). We obtain approximate solutions to the dynamical equations in the cases of “small” and “big” maximal values of the functions σ(t) and δε(t). We show that the single-mode approximation used in the previous studies can be justified in the case of “small” perturbations, but the initially excited field mode cannot be amplified in this case if the laser pulses generate free carriers inside the slab. The amplification could be possible, in principle, for extremely high maximum values of conductivity and the concentration of free carries (the model of an “almost ideal conductor”) created inside the slab under the crucial condition providing the negativity of the function δε(t). This result follows from a simple approximate analytical solution confirmed by exact numerical calculations. However, the evaluation shows that the necessary energy of laser pulses must be, probably, unrealistically high.     

CHAOTIC AND NON-LINEAR DYNAMIC ANALYSIS OF A TWO-AXIS RATE GYRO WITH FEEDBACK CONTROL MOUNTED ON A SPACE VEHICLE

An analysis is presented in this paper for a two-axis rate gyro subjected to linear feedback control mounted on a space vehicle, which is spinning with uncertain angular velocity ω_z(t) about its spin of the gyro. For the autonomous case in which ω_z(t) is steady, the stability analysis of the system is studied by Routh-Hurwitz theory. For the non-autonomous case in which ω_z(t) is sinusoidal function, this system is a strongly non-linear damped system subjected to parametric excitation. By varying the amplitude of sinusoidal motion, periodic and chaotic responses of this parametrically excited non-linear system are investigated using the numerical simulation. Some observations on symmetry-breaking bifurcations, period-doubling bifurcations, and chaotic behavior of the system are investigated by various numerical techniques such as phase portraits, Poincaré maps, average power spectra, and Lyapunov exponents. In addition, some discussions about chaotic motions of this system can be suppressed and changed into regular motions by a suitable constant motor torque are included.     

Detecting unknown paths on complex networks through random walks

In this article, we investigate the problem of detecting unknown paths on complex networks through random walks. To detect a given path on a network a random walker should pass through the path from its initial node to its terminal node in turn. We calculate probability ?(t) that a random walker detects a given path on a connected network in t steps when it starts out from source node s. We propose an iteration formula for calculating ?(t). Generating function of ?(t) is also derived. Major factors affecting ?(t), such as walking time t, path length l, starting point of the walker, structure of the path, and topological structure of the underlying network are further discussed. Among these factors, two most outstanding ones are walking time t and path length l. On the one hand, ?(t) increases as t increases, and ?(∞)=1, which shows that the longer the walking time, the higher the chance of detecting a given path, and the walker will discover the path sooner or later so long as it keeps wandering on the network. On the other hand, ?(t) drops substantially as path length l increases, which shows that the longer the path, the lower the chance for the walker to find it in a given time. Apart from path length, path structure also has obvious effect on ?(t). Starting point of the walker has only minor influence on ?(t), but topological structure of the underlying network has strong influence on ?(t). Simulations confirm our analytic results.     

Transient bi-fractional diffusion: space-time coupling inducing the coexistence of two fractional diffusions

Anomalous diffusion is researched within the framework of the coupled continuous time random walk model, in which the space-time coupling is considered through the correlated function g(t) ~ t ^γ , 0 ≤ γ< 2, and the probability density function ω(t) of a particle’s transition time t follows a power law for large t: ω(t) ~ t ^{? (1 + α)},1 <α< 2. The bi-fractional generalized master equation is derived analytically which can be applied to describe the transient bi-fractional diffusion phenomenon which is induced by the space-time coupling and the asymptotic behavior of ω(t). Numerical results show that for the transient bi-fractional diffusion, there is a transition from one fractional diffusion to another one in the diffusive process.     

Some exact solutions of the local induction equation for the motion of a vortex in a Bose–Einstein condensate with a Gaussian density profile

The dynamics of a vortex filament in a Bose–Einstein condensate whose equilibrium density in the reference frame rotating at the angular velocity Ω is Gaussian with the quadratic form r·D?r has been considered. It has been shown that the equation of motion of the filament in the local-induction approximation permits a class of exact solutions in the form R(β, t) = βM(t) + N(t) of a straight vortex, where β is the longitudinal parameter and is the time. The vortex slips over the surface of an ellipsoid, which follows from the conservation laws N · D?N=C ₁ and M · D?N=C ₀=0. The equation of the evolution of the tangential vector M(t) appears to be closed and has integrals of motion M ·D?M=C ₂ and (|M| ? M· G?Ω) = C, with the matrix G? = 2(I?TrD? ? D?)^?1. Crossing of the respective isosurfaces specifies trajectories in the phase space.     

Internal friction in alloys due to plasticity of the diffusion phase transformation

An analytic expression is obtained for the time dependence Q ^?1(t) of internal friction associated with plasticity of a phase transformation. Time dependences Q ^?1(t) of internal friction of the Pb-62Sn and Pb-1.9Sn alloys (wt.%) alloys were studied in the regime of continuous excitation of resonant flexural vibrations. The measurements of the Q ^?1(t) dependences for 1 h at room temperature and a fixed strain amplitude ε₀ ≈ 7 and 19 min) for the Pb-62Sn alloy. For the Pb-1.9Sn alloy under the same conditions, an exponential decrease followed by an internal friction peak (at t _m ≈ 7 min) is observed. It is shown numerically that the above singularities of internal friction are formed by processes of intermittent phase decomposition of Pb-Sn alloys in the cyclic stress field produced by an external load. Experimental data on Q ^?1(t) are used for reconstructing the kinetic curves describing the decomposition (conversion) ratio as a function of time and for calculating the corresponding values of parameters K and n of the Avrami kinetic equation for the Pb-62Sn alloy.     

The inverse problem for a general N × N spectral equation

The spectral transform ?u/?x = {A(ζ)+B(x,ζ)} · u, where u is an n-element column vector and A(ζ) and B(x,ζ) are nxn matrices, is considered. A set of spectral data is given and the problem of reconcstructing B(x, ζ) from this spectral data is solved for a large class of cases. A special case of this spectral transform is used to solve the Boussinesq equation.     

Limit behavior of saturated approximations of nonlinear Schrödinger equation

We consider the solutionu _?(t) of the saturated nonlinear Schrödinger equation (1) $$i\partial u/\partial t = - \Delta u - \left| u \right|^{4/N} u + \varepsilon \left| u \right|^{q - 1} uandu(0,.) = \varphi (.)$$ where \(N \geqslant 2,\varepsilon > 0,1 + 4/N< q< (N + 2)/(N - 2),u:\mathbb{R} \times \mathbb{R}^N \to \mathbb{C},\varphi \) , ? is a radially symmetric function inH ¹(R ^N ). We assume that the solution of the limit equation is not globally defined in time. There is aT>0 such that \(\mathop {\lim }\limits_{t \to T} \left\| {u(t)} \right\|_{H^1 } = + \infty \) , whereu(t) is solution of (1) $$i\partial u/\partial t = - \Delta u - \left| u \right|^{4/N} uandu(0,.) = \varphi (.)$$ For ?>0 fixed,u _?(t) is defined for all time. We are interested in the limit behavior as ?→0 ofu _?(t) fort≥T. In the case where there is no loss of mass inu _? at infinity in a sense to be made precise, we describe the behavior ofu _? as ? goes to zero and we derive an existence result for a solution of (1) after the blow-up timeT in a certain sense. Nonlinear Schrödinger equation with supercritical exponents are also considered.     

Dynamic model for alpha particle emission during fission

We investigate the motion of anα-particle in the average time depencent potentialV(R _α ,t) of a fissioning nucleus. The emission process is treated quantum mechanically via a numerical solution of the one-body Schroedinger equation withV(R _α ,t). This solution yields the distribution of initial conditions for classical trajectories describing theα-particles outside the Coulomb barrier. The time and shape dependence ofV(R _α,t) is shown to have significant influence on the observable angle and energy distribution of theα-particles emitted during fission.     

Current fluctuations in chain of nanoparticles

We have fabricated and analyzed devices where an electrode gap was bridged by a chain of citrate stabilized gold particles. The current-voltage characteristic I∼ζ(V−Vt) with voltage threshold Vt≈0 and scaling exponent ζ=0.91 is attributed to charge transport in a one-dimensional path. We observe current fluctuations in the frequency range 3-100 mHz as function of time as a result of conformational changes in the citrate molecules induced by charge transfer across the molecules.     

Studying the effect of biaxial anisotropy on nonlinear magnetization oscillations accompanying the 90° pulsed magnetization process in ferrite–garnet films with easy-plane anisotropy

Time dependences of the azimuthal component of the torque T _φ(t) acting on magnetization are calculated to understand the nature of the delayed magnetization acceleration effect observed during the 90° pulsed magnetization of real ferrite–garnet films, in which biaxial anisotropy exists alongside with in-plane anisotropy. A calculation technique based on analyzing an operating point trajectory is used. Calculations show that if the effective anisotropy field H _K2 is comparable to the magnetizing pulse amplitude H _ma, abruptly ascending regions at characteristic times t* in curves T _φ(t) arise, in the limit of which nonlinear magnetization oscillations formed. The shape of these regions depends weakly on the magnetizing pulse front duration τ_f. This explains the reason of the weak dependence of the nonlinear magnetization oscillations on duration of the magnetizing pulse front. Calculations also show that the main features of the delayed acceleration effect are less clear upon an increase of the pulse amplitude: the behavior of curves T _φ(t) becomes smoother near times t*, and an increase in the pulse front duration is accompanied by a stronger drop in the intensity of magnetization oscillations.     

Domain Number Distribution in the Nonequilibrium Ising Model

We study domain distributions in the one-dimensional Ising model subject to zero-temperature Glauber and Kawasaki dynamics. The survival probability of a domain, S(t)～t ^?ψ, and an unreacted domain, Q ₁(t)～t ^?δ, are characterized by two independent nontrivial exponents. We develop an independent interval approximation that provides close estimates for many characteristics of the domain length and number distributions including the scaling exponents.     

Kinetic Limit of N-Body Description of Wave–Particle Self-Consistent Interaction

A system of N particles $\xi ^N = x_1 ,\upsilon_1,...,x_N ,\upsilon _N )$ interacting self-consistently with one wave Z = A exp(iφ) is considered. Given initial data (Z ^(N)(0), ξ ^N (0)), it evolves according to Hamiltonian dynamics to (Z ^(N)(t), ξ ^N (t)). In the limit N → ∞, this generates a Vlasov-like kinetic equation for the distribution function f(x, v, t), abbreviated as f(t), coupled to the envelope equation for Z: initial data (Z ^(∞)(0), f(0)) evolve to (Z ^(∞)(t), f(t)). The solution (Z, f) exists and is unique for any initial data with finite energy. Moreover, for any time T>0, given a sequence of initial data with N particles distributed so that the particle distribution f ^N(0) → f(0) weakly and with Z ^(N)(0) → Z(0) as N → ∞, the states generated by the Hamiltonian dynamics at all times 0 ≤ t ≤ T are such that (Z ^(N)(t), f ^N(t)) converges weakly to (Z ^(∞)(t), f(t)).     

Chaotic time series prediction: From one to another

In this Letter, a new local linear prediction model is proposed to predict a chaotic time series of a component x(t) by using the chaotic time series of another component y(t) in the same system with x(t). Our approach is based on the phase space reconstruction coming from the Takens embedding theorem. To illustrate our results, we present an example of Lorenz system and compare with the performance of the original local linear prediction model.