Initial value representation for the SU(n) semiclassical propagator

The semiclassical propagator in the representation of SU(n) coherent states is characterized by isolated classical trajectories subjected to boundary conditions in a doubled phase space. In this paper we recast this expression in terms of an integral over a set of initial-valued trajectories. These trajectories are monitored by a filter that collects only the appropriate contributions to the semiclassical approximation. This framework is suitable for the study of bosonic dynamics in n modes with fixed total number of particles. We exemplify the method for a Bose-Einstein condensate trapped in a triple-well potential, providing a detailed discussion on the accuracy and efficiency of the procedure.Comment: 24 pages, 6 figure

Modelling Holling type II functional response in deterministic and stochastic food chain models with mass conservation

The Rosenzweig-MacArthur predator-prey model is the building block in modeling food chain, food webs and ecosystems. There are a number of hidden assumptions involved in the derivation. For instance the prey population growth is logistic without predation but also with predation. In order to reveal these we will start with modelling a resource-predator-prey system in a closed spatially homogeneous environment. This allows us to keep track of the nutrient flow. With an instantaneous remineralisation of the products excreted in the environment by the populations and dead body mass there is conservation of mass. This allows for a model dimension reduction and yields the mass balance predator-prey model. When furthermore the searching and handling processes are much faster that the population changing rates, the trophic interaction is described by a Holling type II functional response, also assumed in the Rosenzweig-MacArthur model. The derivation uses an extended deterministic model with number of searching and handling predators as model variables where the ratio of the predator/prey body masses is used as a mechanistic time-scale parameter. This extended model is also used as a starting point for the derivation of a stochastic model. We will investigate the stochastic effects of random switching between searching and handling of the predators and predator dying. Prey growth by consumption of ambient resources is still deterministic and therefore the stochastic model is hybrid. The transient dynamics is studied by numerical Monte Carlo simulations and also the quasi-equilibrium distribution for the population quantities is calculated. The body mass of the prey individual is the scaling parameter in the stochastic model formulation. This allows for a quantification of the mean-field approximation criterion for the justification of replacement of the stochastic by a deterministic model.Marie Skłodowska-Curie grant agreement No. 79249

Semiclassical coherent state propagator for systems with spin

We derive the semiclassical limit of the coherent state propagator for systems with two degrees of freedom of which one degree of freedom is canonical and the other a spin. Systems in this category include those involving spin-orbit interactions and the Jaynes-Cummings model in which a single electromagnetic mode interacts with many independent two-level atoms. We construct a path integral representation for the propagator of such systems and derive its semiclassical limit. As special cases we consider separable systems, the limit of very large spins and the case of spin 1/2.Comment: 19 pages, no figure

Semiclassical Propagation of Wavepackets with Real and Complex Trajectories

We consider a semiclassical approximation for the time evolution of an originally gaussian wave packet in terms of complex trajectories. We also derive additional approximations replacing the complex trajectories by real ones. These yield three different semiclassical formulae involving different real trajectories. One of these formulae is Heller's thawed gaussian approximation. The other approximations are non-gaussian and may involve several trajectories determined by mixed initial-final conditions. These different formulae are tested for the cases of scattering by a hard wall, scattering by an attractive gaussian potential, and bound motion in a quartic oscillator. The formula with complex trajectories gives good results in all cases. The non-gaussian approximations with real trajectories work well in some cases, whereas the thawed gaussian works only in very simple situations.Comment: revised text, 24 pages, 6 figure

Imaginary Phases in Two-Level Model with Spontaneous Decay

We study a two-level model coupled to the electromagnetic vacuum and to an external classic electric field with fixed frequency. The amplitude of the external electric field is supposed to vary very slow in time. Garrison and Wright [{\it Phys. Lett.} {\bf A128} (1988) 177] used the non-hermitian Hamiltonian approach to study the adiabatic limit of this model and obtained that the probability of this two-level system to be in its upper level has an imaginary geometric phase. Using the master equation for describing the time evolution of the two-level system we obtain that the imaginary phase due to dissipative effects is time dependent, in opposition to Garrison and Wright result. The present results show that the non-hermitian hamiltonian method should not be used to discuss the nature of the imaginary phases in open systems.Comment: 11 pages, new version, to appear in J. Phys.