Analytical Solitary Wave Solutions to a （3＋1）-Dimensional Gross-Pitaevskii Equation with Variable Coefficients

A （3＋1）-dimensional Gross-Pitaevskii （GP） equation with time variable coefficients is considered, and is transformed into a standard nonlinear Schrodinger （NLS） equation. Exact solutions of the （3＋1）D GP equation are constructed via those of the NLS equation. By applying specific time-modulated nonlinearities, dispersions, and potentials, the dynamics of the solutions can be controlled. Solitary and periodic wave solutions with snaking and breathing behavior are reported.     

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.     

110警车配置及巡逻方案

针对110警车调度问题,引入了图论中的最短路算法以及计算几何的相关理论,建立了车辆调配模型、巡逻路线模型以及基于模糊数学的评价指标模型.另外,用C++编写了一个可视化的软件,不仅实现了手动描点,自动求出覆盖线段集合的功能,同时利用计算机模拟警车的巡逻路线,最后通过计算机检验得到结果,其合理性和实用性都令人满意.针对问题一,通过人机结合,配置17辆警车就能实现D1的目标,很好地兼顾了警车巡逻的运行成本,减少公安部门车辆和人员等的投入.针对问题二,采用模糊数学相关理论使评价指标实现了从定性到定量的转变.针对问题三和六,在D1的基础上,兼顾了巡逻效果的显著性,采用最少被巡逻道路优先的贪心算法建立了动态巡逻模型,得到了合理的巡逻方案.在此方案中,我们动用了30辆警车完成了问题一的目标.另外还额外考虑了案发事件概率不均匀分布的情况,建立了改进模型.针对问题四,在完成问题三指标的基础上,为了尽可能提高巡逻车辆的隐蔽性和增强巡逻效果,采用轮盘赌算法来引入随机性.针对问题五,采用最远距离道路优先贪心策略,使模型尽可能满足条件D1,D2.针对问题七,提出了一些额外因素及其解决方案,进一步完善了模型,使模型更贴近现实.     

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.