Error Estimates of Mixed Methods for Optimal Control Problems Governed by General Elliptic Equations

In this paper, we investigate the error estimates of mixed finite element methods for optimal control problems governed by general elliptic equations. The state and co-state are approximated by the lowest order Raviart-Thomas mixed finite element spaces and the control variable is approximated by piecewise constant functions. We derive $L^2$ and $H^{-1}$-error estimates both for the control variable and the state variables. Finally, a numerical example is given to demonstrate the theoretical results.     

Exact solutions of the Fokker-Planck equations with moving boundaries

By means of time-dependent similarity transformations, we derive exact solutions of the Fokker-Planck equations with moving boundaries in the presence of: (1) a time-dependent linear force and (2) a time-dependent nonlinear force. The method of similarity transformation is simple and can be easily applied to more general Fokker-Planck equations. Furthermore, the knowledge of the exact solutions in closed form can be useful as a benchmark to test approximate numerical or analytical procedures.     

Finite difference methods for second order in space,first order in time hyperbolic systems and the linear shifted wave equation as a model problem in numerical relativity

Motivated by the problem of solving the Einstein equations, we discuss high order finite difference discretizations of first order in time, second order in space hyperbolic systems. Particular attention is paid to the case when first order derivatives that can be identified with advection terms are approximated with non-centered finite difference operators. We first derive general properties of these discrete operators, then we extend a known result on numerical stability for such systems to general order of accuracy. As an application we analyze the shifted wave equation, including the behavior of the numerical phase and group speeds at different orders of approximations. Special attention is paid to when the use of off-centered schemes improves the accuracy over the centered schemes.     

Critical temperature and enhanced isotope effect in the presence of paramagnons in phonon-mediated superconductors

We reconsider the long-standing problem of the effect of spin fluctuations on the critical temperature and isotope effect in a phonon-mediated superconductor. Although the general physics of the interplay between phonons and paramagnons has been rather well understood, the existing approximate formulas fail to describe the correct behavior of Tc for general phonon and paramagnon spectra. Using a controllable approximation, we derive an analytical formula for Tc which agrees well with exact numerical solutions of the Eliashberg equations for a broad range of parameters. Based on both numerical and analytical results, we predict a strong enhancement of the isotope effect when the frequencies of spin fluctuation and phonons are of the same order. This effect may have important consequences for near-magnetic superconductors such as MgCNi3.     

Intriguing properties of kinks in a simple model with a two-component field

We present a summary of the topological and non-topological solitons of a two component field in 1+1 dimensions with application in field theory and in condensed matter physics. We note several intriguing analytical and numerical relationships between these solutions, which we believe to suggest that the relevant coupled nonlinear differential equations may be integrable exactly.     

射频场照射下多自旋体系弛豫的理论计算

根据Liouville-von Neumann方程，对射频场照射下多自旋体系的弛豫进行了理论描述，并用WBR理论推导出了体系的弛豫方程组，给出了各类弛豫速率的理论计算公式.在此基础上，编制了弛豫方程组数值解的计算程序，分别用此程序和Bloch方程计算了双自旋体系在不同情况下的稳态解，并对计算结果进行了简要的分析和讨论. 关键词： 核磁共振 弛豫 射频场 多自旋体系     

The Lorenz model for single-mode homogeneously broadened laser: analytical determination of the unpredictable zone

We have applied harmonic expansion to derive an analytical solution for the Lorenz-Haken equations. This method is used to describe the regular and periodic self-pulsing regime of the single mode homogeneously broadened laser. These periodic solutions emerge when the ratio of the population decay rate ? is smaller than 0:11. We have also demonstrated the tendency of the Lorenz-Haken dissipative system to behave periodic for a characteristic pumping rate “2C _P ”[7], close to the second laser threshold “2C _2th ”(threshold of instability). When the pumping parameter “2C” increases, the laser undergoes a period doubling sequence. This cascade of period doubling leads towards chaos. We study this type of solutions and indicate the zone of the control parameters for which the system undergoes irregular pulsing solutions. We had previously applied this analytical procedure to derive the amplitude of the first, third and fifth order harmonics for the laser-field expansion [7, 17]. In this work, we extend this method in the aim of obtaining the higher harmonics. We show that this iterative method is indeed limited to the fifth order, and that above, the obtained analytical solution diverges from the numerical direct resolution of the equations.     

Numerical Solutions of Variable Coefficient Higher-Order Partial Differential Equations Arising in Beam Models

In this work, an efficient and robust numerical scheme is proposed to solve the variable coefficients’ fourth-order partial differential equations (FOPDEs) that arise in Euler–Bernoulli beam models. When partial differential equations (PDEs) are of higher order and invoke variable coefficients, then the numerical solution is quite a tedious and challenging problem, which is our main concern in this paper. The current scheme is hybrid in nature in which the second-order finite difference is used for temporal discretization, while spatial derivatives and solutions are approximated via the Haar wavelet. Next, the integration and Haar matrices are used to convert partial differential equations (PDEs) to the system of linear equations, which can be handled easily. Besides this, we derive the theoretical result for stability via the Lax–Richtmyer criterion and verify it computationally. Moreover, we address the computational convergence rate, which is near order two. Several test problems are given to measure the accuracy of the suggested scheme. Computations validate that the present scheme works well for such problems. The calculated results are also compared with the earlier work and the exact solutions. The comparison shows that the outcomes are in good agreement with both the exact solutions and the available results in the literature.     

A mathematical solution for the parameters of three interfering resonances

The multiple-solution problem in determining the parameters of three interfering resonances from a fit to an experimentally measured distribution is considered from a mathematical viewpoint. It is shown that there are four numerical solutions for a fit with three coherent Breit-Wigner functions. Although explicit analytical formulae cannot be derived in this case, we provide some constraint equations between the four solutions. For the cases of nonrelativistic and relativistic Breit-Wigner forms of amplitude functions, a numerical method is provided to derive the other solutions from that already obtained, based on the obtained constraint equations. In real experimental measurements with more complicated amplitude forms similar to Breit-Wigner functions, the same method can be deduced and performed to get numerical solutions. The good agreement between the solutions found using this mathematical method and those directly from the fit verifies the correctness of the constraint equations and mathematical methodology used.     

Quasinormal modes of a Schwarzschild black hole immersed in an electromagnetic universe

We study the quasinormal modes(QNMs) of a Schwarzschild black hole immersed in an electromagnetic(EM) universe. The immersed Schwarzschild black hole(ISBH) originates from the metric of colliding EM waves with double polarization [Class. Quantum Grav. 12, 3013(1995)]. The perturbation equations of the scalar fields for the ISBH geometry are written in the form of separable equations. We show that these equations can be transformed to the confluent Heun's equations, for which we are able to use known techniques to perform analytical quasinormal(QNM) analysis of the solutions. Furthermore, we employ numerical methods(Mashhoon and 6~(th)-order Wentzel-Kramers-Brillouin(WKB)) to derive the QNMs. The results obtained are discussed and depicted with the appropriate plots.     

Covariant stringlike dynamics of scroll wave filaments in anisotropic cardiac tissue

It has been hypothesized that stationary scroll wave filaments in cardiac tissue describe a geodesic in a curved space whose metric is the inverse diffusion tensor. Several numerical studies support this hypothesis, but no analytical proof has been provided yet for general anisotropy. In this Letter, we derive dynamic equations for the filament in the case of general anisotropy. These equations are covariant under general spatial coordinate transformations and describe the motion of a stringlike object in a curved space whose metric tensor is the inverse diffusion tensor. Therefore the behavior of scroll wave filaments in excitable media with anisotropy is similar to the one of cosmic strings in a curved universe. Our dynamic equations are valid for thin filaments and for general anisotropy. We show that stationary filaments obey the geodesic equation.     

Current-induced conformational switching in single-molecule junctions

Current-induced conformational switching in single-molecule junctions constitutes a fundamental process in molecular electronics. Motivated by recent experiments on azobenzene derivatives, we study this process for molecules which exhibit two (meta)stable conformations in the neutral state but only a single stable conformation in the ionic state. We derive and analyze appropriate Fokker–Planck equations obtained from a density-matrix formalism starting from a generic model and present comprehensive analytical and numerical results for the switching dynamics in general and the quantum yield in particular.     

An Analysis of the Transition Zone Between the Various Scaling Regimes in the Small-World Model

We analyse the so-called small-world network model (originally devised by Strogatz and Watts), treating it, among other things, as a case study of non-linear coupled difference or differential equations. We derive a system of evolution equations containing more of the previously neglected (possibly relevant) non-linear terms. As an exact solution of this entangled system of equations is out of question we develop a (as we think, promising) method of enclosing the “exact” solutions for the expected quantities by upper and lower bounds, which represent solutions of a slightly simpler system of differential equation. Furthermore we discuss the relation between difference and differential equations and scrutinize the limits of the spreading idea for random graphs. We then show that there exists in fact a “broad” (with respect to scaling exponents) crossover zone, smoothly interpolating between linear and logarithmic scaling of the diameter or average distance. We are able to corroborate earlier findings in certain regions of phase or parameter space (as e.g. the finite size scaling ansatz) but find also deviations for other choices of the parameters. Our analysis is supplemented by a variety of numerical calculations, which, among other things, quantify the effect of various approximations being made. With the help of our analytical results we manage to calculate another important network characteristic, the (fractal) dimension, and provide numerical values for the case of the small-world network.     

Solitons in photovoltaic photorefractive media with an external electric field

We use the classical Lie-group method for studying the evolution equation in photovoltaic photorefractive media with an external electric field, reducing it to some similarity equations firstly, and then obtain some exact analytical solutions including the soliton solution, the period solution and the oscillatory solution. We also obtain the bright soliton, dark soliton, gray soliton from these similarity equations with the numerical method. Furthermore, we investigate what factors contribute to the beamwidth of these solitons with the numerical method and know the beamwidth of these solitons are associated with the external electric field, the photovoltaic field and the intensity ratio of the incident soliton.     

Localized polaritons and second-harmonic generation in a resonant medium with quadratic nonlinearity

We derive model equations for the propagation of ultrashort pulses in materials with resonant linear and quadratic nonlinear responses and find approximate soliton solutions describing all-bright and dark-bright polaritons. We report the specific phase matching condition for efficient 2nd harmonic generation, which involves detuning from the resonance. We also demonstrate that the 2nd harmonic emission by the polaritonic pulses can lead to reduction of their group velocity, having zero as a theoretical limit. Our analytical results are supported by numerical simulations.     

Transverse spectra of waves in random media

We consider the statistics of the transverse spectra of forward-propagating waves in a stationary random medium. A short-range perturbation solution is used to derive the difference equations that govern the long-range evolution of the ensemble-averaged transverse wave spectrum and coherence. The conditions under which these equations may be approximated by differential and integro-differential equations are given, and it is shown that the approximation is valid for the treatment of beam propagation provided that the transverse dimension of the beam is sufficiently large, and at ranges where the transverse coherence length of the beam remains larger than a wavelength. The equations that are derived are not limited by the parabolic approximation, and are amenable to numerical solution by marching techniques. We use the equation that governs the spectral density of the total energy flux, and also the propagation of waves which are statistically homogeneous in transverse planes, to show the conditions under which previously studied approximations derive from the present formulation, and to illustrate the numerical solution of the problem.     

Fractional diffusion-wave problem in cylindrical coordinates

In this Letter, we present analytical and numerical solutions for an axis-symmetric diffusion-wave equation. For problem formulation, the fractional time derivative is described in the sense of Riemann-Liouville. The analytical solution of the problem is determined by using the method of separation of variables. Eigenfunctions whose linear combination constitute the closed form of the solution are obtained. For numerical computation, the fractional derivative is approximated using the Grünwald-Letnikov scheme. Simulation results are given for different values of order of fractional derivative. We indicate the effectiveness of numerical scheme by comparing the numerical and the analytical results for α=1 which represents the order of derivative.     

Some theoretical and numerical results for delayed neural field equations

In this paper we study neural field models with delays which define a useful framework for modeling macroscopic parts of the cortex involving several populations of neurons. Nonlinear delayed integro-differential equations describe the spatio-temporal behavior of these fields. Using methods from the theory of delay differential equations, we show the existence and uniqueness of a solution of these equations. A Lyapunov analysis gives us sufficient conditions for the solutions to be asymptotically stable. We also present a fairly detailed study of the numerical computation of these solutions. This is, to our knowledge, the first time that a serious analysis of the problem of the existence and uniqueness of a solution of these equations has been performed. Another original contribution of ours is the definition of a Lyapunov functional and the result of stability it implies. We illustrate our numerical schemes on a variety of examples that are relevant to modeling in neuroscience.