Quasi-particle structure around a single vortex in high-Tc superconductors

For analyzing the checker-board like modulation of the local density of states (LDOS) around a vortex observed in the slightly overdoped Bi₂Sr₂CaCu₂O_x, we examined the effect of pseudogap state of high-T_c superconductors to the LDOS around the vortex. We first derived the Bogoliubov-de Gennes equation for d-wave superconductivity (d-SC) in the presence of d-spin density wave (d-SDW). Using the Fourier–Bessel expansion, we solved this equation for a single vortex state, numerically. We found that the peak of the bound states around E = 0 becomes small and modulation of the LDOS is observed for larger d-SDW order parameter.     

Classicalization times of parametrically amplified “Schrödinger cat” states coupled to phase-sensitive reservoirs

The exact Wigner function of a parametrically excited quantum oscillator in a phase-sensitive amplifying/attenuating reservoir is found for initial even/odd coherent states. Studying the evolution of negativity of the Wigner function we show the difference between the “initial positivization time” (IPT), which is inversely proportional to the square of the initial size of the superposition, and the “final positivization time” (FPT), which does not depend on this size. Both these times can be made arbitrarily long in maximally squeezed high-temperature reservoirs. Besides, we find the conditions when some (small) squeezing can exist even after the Wigner function becomes totally positive.     

Multi-moment advection scheme for Vlasov simulations

We present a new numerical scheme for solving the advection equation and its application to Vlasov simulations. The scheme treats not only point values of a profile but also its zeroth to second order piecewise moments as dependent variables, for better conservation of the information entropy. We have developed one-and two-dimensional schemes and show that they provide quite accurate solutions within reasonable usage of computational resources compared to other existing schemes. The two-dimensional scheme can accurately solve the solid body rotation problem of a gaussian profile for more than hundred rotation periods with little numerical diffusion. This is crucially important for Vlasov simulations of magnetized plasmas. Applications of the one- and two-dimensional schemes to electrostatic and electromagnetic Vlasov simulations are presented with some benchmark tests.     

ON THE EXPECTED DISCOUNTED PENALTY FUNCTION IN A MARKOV-DEPENDENT RISK MODEL WITH CONSTANT DIVIDEND BARRIER

This article considers a Markov-dependent risk model with a constant dividend barrier. A system of integro-differential equations with boundary conditions satisfied by the expected discounted penalty function, with given initial environment state, is derived and solved. Explicit formulas for the discounted penalty function are obtained when the initial surplus is zero or when all the claim amount distributions are from rational family. In two state model, numerical illustrations with exponential claim amounts are given.     

Optimal control for reinitialization in finite element level set methods

A new optimal control problem that incorporates the residual of the Eikonal equation into its objective is presented. The formulation of the state equation is based on the level set transport equation but extended by an additional source term, correcting the solution so as to minimize the objective functional. The method enforces the constraint so that the interface cannot be displaced at least in the continuous setting. The system of first‐order optimality conditions is derived, linearized, and solved numerically. The control also prevents numerical instabilities, so that no additional stabilization techniques are required. This approach offers the flexibility to include other desired design criteria into the objective functional. The methodology is evaluated numerically in three different examples and compared with other PDE‐based reinitialization techniques. Copyright © 2016 John Wiley & Sons, Ltd.     

Approximate wφ - Ωφ Relations in Quintessence Models

Quintessence field is a widely-studied candidate of dark energy. There is ＂tracker solution＂ in quintessence models, in which evolution of the field φ at present times is not sensitive to its initial conditions. When the energy density of dark energy is negleetable （Ωφ 〈〈 1）, evolution of the tracker solution can be well analysed from ＂tracker equation＂. In this paper, we try to study evolution of the quintessence field from ＂full tracker equation＂, which is valid for all spans of Ωφ. We get stable fixed points of we and wφ （noted as wφ and Ωφ） from the ＂full tracker equation＂, i.e., we and ωφ will always approach ωφ and Ωφ respectively. Since wφ and Ωφ are analytic functions of φ, analytic relation of φ can be obtained, which is a good approximation for the we φ relation and can be obtained for the most type of quintessence potentials. By using this approximation, we find that inequalities ωφ 〈 we and 〈ωφ are statisfied if the we （or ωφ） decreases with time. In this way, the potentiai U（φ） can be constrained directly from observations, by no need of solving the equations of motion numerically.     

Fast Gaussian wavepacket transforms and Gaussian beams for the Schrödinger equation

This paper introduces a wavepacket-transform-based Gaussian beam method for solving the Schrödinger equation. We focus on addressing two computational issues of the Gaussian beam method: how to generate a Gaussian beam representation for general initial conditions and how to perform long time propagation for any finite period of time. To address the first question, we introduce fast Gaussian wavepacket transforms and develop on top of them an efficient initialization algorithm for general initial conditions. Based on this new initialization algorithm, we address the second question by reinitializing the beam representation when the beams become too wide. Numerical examples in one, two, and three dimensions demonstrate the efficiency and accuracy of the proposed algorithms. The methodology can be readily generalized to deal with other semi-classical quantum mechanical problems.     

On the existence,uniqueness and stability of solutions for semi-linear generalized elasticity equation with general damping term

In this paper, we consider a semi-linear generalized hyperbolic boundary value problem associated to the linear elastic equations with general damping term and nonlinearities of variable exponent type. Under suitable conditions, local and global existence theorems are proved. The uniqueness of the solution have been gotten by eliminating some hypotheses that have been imposed by other authors for different particular problems. We show that any solution with nontrivial initial datum becomes stable.     

具交叉感染的流行病模型的数学分析

In this paper,an initial boundary value problem with homogeneous Neumann boundary condition is studied for a reaction-diffusion system which models the spread of infectious diseases within two population groups by means of self and criss-cross infection mechanism. Existence,uniqueness and boundedness of the nonnegative global solution (u1,u2,u3,u4),existence of the limit of the solution as t→+∞,namely,(u1,u2,u3,u4)→(u*1,0,u*3,0),are proved.