Landau-Zener problem with waiting at the minimum gap and related quench dynamics of a many-body system

We discuss a technique for solving the Landau-Zener (LZ) problem of finding the probability of excitation in a two-level system. The idea of time reversal for the Schrodinger equation is employed to obtain the state reached at the final time and hence the excitation probability. Using this method, which can reproduce the well-known expression for the LZ transition probability, we solve a variant of the LZ problem, which involves waiting at the minimum gap for a time t(w); we find an exact expression for the excitation probability as a function of t(w). We provide numerical results to support our analytical expressions. We then discuss the problem of waiting at the quantum critical point of a many-body system and calculate the residual energy generated by the time-dependent Hamiltonian. Finally, we discuss possible experimental realizations of this work.

Exact free-surface flows for shallow-water equations .2. the compressible case

Exact free surface flows with shear in a compressible barotropic medium are found, extending the authors' earlier work for the incompressible medium. The barotropic medium is of finite extent in the vertical direction, while it is infinite in the horizontal direction. The ''shallow water'' equations for a compressible barotropic medium, subject to boundary conditions at the free surface and at the bottom, are solved in terms of double psi-series, Simple wave and time-dependent solutions are found; for the former the free surface is of arbitrary shape while for the latter it is a damping traveling wave in the horizontal direction, For other types of solutions, the height of the free surface is constant either on lines of constant acceleration or on lines of constant speed. In the case of an isothermal medium, when gamma = 1, we again find simple wave and time-dependent solutions.

On the solution of bivariate population balance equations for aggregation: Pivotwise expansion of solution domain

The solution of a bivariate population balance equation (PBE) for aggregation of particles necessitates a large 2-d domain to be covered. A correspondingly large number of discretized equations for particle populations on pivots (representative sizes for bins) are solved, although at the end only a relatively small number of pivots are found to participate in the evolution process. In the present work, we initiate solution of the governing PBE on a small set of pivots that can represent the initial size distribution. New pivots are added to expand the computational domain in directions in which the evolving size distribution advances. A self-sufficient set of rules is developed to automate the addition of pivots, taken from an underlying X-grid formed by intersection of the lines of constant composition and constant particle mass. In order to test the robustness of the rule-set, simulations carried out with pivotwise expansion of X-grid are compared with those obtained using sufficiently large fixed X-grids for a number of composition independent and composition dependent aggregation kernels and initial conditions. The two techniques lead to identical predictions, with the former requiring only a fraction of the computational effort. The rule-set automatically reduces aggregation of particles of same composition to a 1-d problem. A midway change in the direction of expansion of domain, effected by the addition of particles of different mean composition, is captured correctly by the rule-set. The evolving shape of a computational domain carries with it the signature of the aggregation process, which can be insightful in complex and time dependent aggregation conditions. (c) 2012 Elsevier Ltd. All rights reserved.

Existence, uniqueness, and stability of solutions of a class of nonlinear partial differential equations

In this study we investigate the existence, uniqueness and asymptotic stability of solutions of a class of nonlinear integral equations which are representations for some time dependent non- linear partial differential equations. Sufficient conditions are established which allow one to infer the stability of the nonlinear equations from the stability of the linearized equations. Improved estimates of the domain of stability are obtained using a Liapunov Functional approach. These results are applied to some nonlinear partial differential equations governing the behavior of nonlinear continuous dynamical systems.

Teoria de Landau-Ginzburg para o estado supercondutor nemático

O objetivo geral deste projeto é propor um modelo bidimensional que descreva o novo estado supercondutor, que apresenta simetria de cristal líquido, chamado de supercondutor nemático. O estudo começa com uma introdução sobre a teoria de Landau-Ginzburg das transições de fase, onde são discutidos conceitos como parâmetro de ordem e as ordens das transições de fase, que são essenciais para o desenvolvimento deste projeto. Em seguida, é feita uma discussão sobre as principais características dos supercondutores como a resistência zero, o efeito Meissner-Ochsenfeld, os tipos de supercondutores, o surgimento de vórtices e uma análise sobre a teoria de Landau-Ginzburg para transição de fase metal-supercondutor. Após isto, é feita uma abordagem sobre os principais tipos de cristais líquidos, com destaque ao cristal líquido nemático, onde é desenvolvida a teoria de Landau-Ginzburg para transição de fase isotrópica-nemática e um estudo sobre o surgimento de disclinações no cristal líquido nemático em duas dimensões. Por fim, é apresentado o modelo proposto para descrever o estado supercondutor nemático, com a construção da teoria de Landau-Ginzburg, o estudo do acoplamento entre as fases e os defeitos topológicos presentes nesse estado.

Development of realistic liquid equations of state for multi-phase CFD modelling

Several equations of state (EOS) have been incorporated into a novel algorithm to solve a system of multi-phase equations in which all phases are assumed to be compressible to varying degrees. The EOSs are used to both supply functional relationships to couple the conservative variables to the primitive variables and to calculate accurately thermodynamic quantities of interest, such as the speed of sound. Each EOS has a defined balance of accuracy, robustness and computational speed; selection of an appropriate EOS is generally problem-dependent. This work employs an AUSM+-up method for accurate discretisation of the convective flux terms with modified low-Mach number dissipation for added robustness of the solver. In this paper we show a newly-developed time-marching formulation for temporal discretisation of the governing equations with incorporated time-dependent source terms, as well as considering the system of eigenvalues that render the governing equations hyperbolic.

Uniformly convergent finite element methods for singularly perturbed parabolic partial differential equations

This thesis is concerned with uniformly convergent finite element methods for numerically solving singularly perturbed parabolic partial differential equations in one space variable. First, we use Petrov-Galerkin finite element methods to generate three schemes for such problems, each of these schemes uses exponentially fitted elements in space. Two of them are lumped and the other is non-lumped. On meshes which are either arbitrary or slightly restricted, we derive global energy norm and L2 norm error bounds, uniformly in the diffusion parameter. Under some reasonable global assumptions together with realistic local assumptions on the solution and its derivatives, we prove that these exponentially fitted schemes are locally uniformly convergent, with order one, in a discrete L∞norm both outside and inside the boundary layer. We next analyse a streamline diffusion scheme on a Shishkin mesh for a model singularly perturbed parabolic partial differential equation. The method with piecewise linear space-time elements is shown, under reasonable assumptions on the solution, to be convergent, independently of the diffusion parameter, with a pointwise accuracy of almost order 5/4 outside layers and almost order 3/4 inside the boundary layer. Numerical results for the above schemes are presented. Finally, we examine a cell vertex finite volume method which is applied to a model time-dependent convection-diffusion problem. Local errors away from all layers are obtained in the l2 seminorm by using techniques from finite element analysis.

Collective excitations of trapped Bose-Einstein condensates in the energy and time domains

A time-dependent method for calculating the collective excitation frequencies and densities of a trapped, inhomogeneous Bose-Einstein condensate with circulation is presented. The results are compared with time-independent solutions of the Bogoliubov-de Gennes equations. The method is based on time-dependent linear-response theory combined with spectral analysis of moments of the excitation modes of interest. The technique is straightforward to apply, extremely efficient in our implementation with parallel fast Fourier transform methods, and produces highly accurate results. For high dimensionality or low symmetry the time-dependent approach is a more practical computational scheme and produces accurate and reliable data. The method is suitable for general trap geometries, condensate flows and condensates permeated with defects and vortex structures.

Localized structures in complex plasmas in the presence of a magnetic field

In this work, the general framework in which fits our investigation is that of modeling the dynamics of dust grains therein dusty plasma (complex plasma) in the presence of electromagnetic fields. The generalized discrete complex Ginzburg-Landau equation (DCGLE) is thus obtained to model discrete dynamical structure in dusty plasma with Epstein friction. In the collisionless limit, the equation reduces to the modified discrete nonlinear Schrödinger equation (MDNLSE). The modulational instability phenomenon is studied and we present the criterion of instability in both cases and it is shown that high values of damping extend the instability region. Equations thus obtained highlight the presence of soliton-like excitation in dusty plasma. We studied the generation of soliton in a dusty plasma taking in account the effects of interaction between dust grains and theirs neighbours. Numerical simulations are carried out to show the validity of analytical approach.

A time-domain probe method for three-dimensional rough surface reconstructions

The task of this paper is to develop a Time-Domain Probe Method for the reconstruction of impenetrable scatterers. The basic idea of the method is to use pulses in the time domain and the time-dependent response of the scatterer to reconstruct its location and shape. The method is based on the basic causality principle of timedependent scattering. The method is independent of the boundary condition and is applicable for limited aperture scattering data. In particular, we discuss the reconstruction of the shape of a rough surface in three dimensions from time-domain measurements of the scattered field. In practise, measurement data is collected where the incident field is given by a pulse. We formulate the time-domain fieeld reconstruction problem equivalently via frequency-domain integral equations or via a retarded boundary integral equation based on results of Bamberger, Ha-Duong, Lubich. In contrast to pure frequency domain methods here we use a time-domain characterization of the unknown shape for its reconstruction. Our paper will describe the Time-Domain Probe Method and relate it to previous frequency-domain approaches on sampling and probe methods by Colton, Kirsch, Ikehata, Potthast, Luke, Sylvester et al. The approach significantly extends recent work of Chandler-Wilde and Lines (2005) and Luke and Potthast (2006) on the timedomain point source method. We provide a complete convergence analysis for the method for the rough surface scattering case and provide numerical simulations and examples.

A Trotter-Kato product formula for a class of non-autonomous evolution equations

In this article dedicated to Professor V. Lakshmikantham on the occasion of the celebration of his 84th birthday, we announce new results concerning the existence and various properties of an evolution system UA+B(t, s)(0 <= s <= t <= T) generated by the sum -(A(t)+B(t)) of two linear, time-dependent and generally unbounded operators defined on time-dependent domains in a complex and separable Banach space B. In particular, writing G(B) for the algebra of all linear bounded operators on B, we can express UA+B(t, s)(0 <= s <= t <= T) as the strong limit in L(B) of a product of the holomorphic contraction semigroups generated by -A(t) and -B(t), thereby getting a product formula of the Trotter-Kato type under very general conditions which allow the domain D(A(t)+B(t)) to evolve with time provided there exists a fixed set D subset of boolean AND D-t epsilon[0,D-T](A(t)+B(t)) everywhere dense in B. We then mention several possible applications of our product formula to various classes of non-autonomous parabolic initial-boundary value problems, as well as to evolution problems of Schrodinger type related to the theory of time-dependent singular perturbations of self-adjoint operators in quantum mechanics. We defer all the proofs and all the details of the applications to a separate publication. (C) 2008 Elsevier Ltd. All rights reserved.

Mean-field equations for cigar- and disc-shaped Bose and Fermi superfluids

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Noise and ultraviolet divergences in the dynamics of the chiral condensate in QCD

The time evolution of the matter produced in high energy heavy-ion collisions seems to be well described by relativistic viscous hydrodynamics. In addition to the hydrodynamic degrees of freedom related to energy-momentum conservation, degrees of freedom associated with order parameters of broken continuous symmetries must be considered because they are all coupled to each other. of particular interest is the coupling of degrees of freedom associated with the chiral symmetry of QCD. Quantum and thermal fluctuations of the chiral fields act as noise sources in the classical equations of motion, turning them into stochastic differential equations in the form of Ginzburg-Landau-Langevin (GLL) equations. Analytic solutions of GLL equations are attainable only in very special circumstances and extensive numerical simulations are necessary, usually by discretizing the equations on a spatial lattice. However, a not much appreciated issue in the numerical simulations of GLL equations is that ultraviolet divergences in the form of lattice-spacing dependence plague the solutions. The divergences are related to the well-known Rayleigh-Jeans catastrophe in classical field theory. In the present communication we present a systematic lattice renormalization method to control the catastrophe. We discuss the implementation of the method for a GLL equation derived in the context of a model for the QCD chiral phase transition and consider the nonequilibrium evolution of the chiral condensate during the hydrodynamic flow of the quark-gluon plasma.

Stability of trapped Bose-Einstein condensates

In three-dimensional trapped Bose-Einstein condensate (BEC), described by the time-dependent Gross-Pitaevskii-Ginzburg equation, we study the effect of initial conditions on stability using a Gaussian variational approach and exact numerical simulations. We also discuss the validity of the criterion for stability suggested by Vakhitov and Kolokolov. The maximum initial chirp (initial focusing defocusing of cloud) that can lead a stable condensate to collapse even before the number of atoms reaches its critical limit is obtained for several specific cases. When we consider two- and three-body nonlinear terms, with negative cubic and positive quintic terms, we have the conditions for the existence of two phases in the condensate. In this case, the magnitude of the oscillations between the two phases are studied considering sufficient large initial chirps. The occurrence of collapse in a BEC with repulsive two-body interaction is also shown to be possible.

Dynamics of the deconfinement transition

We estimate the dissipation coefficient Γ that appears in Ginzburg-Landau-Langevin equations that describe phenomenologically the deconfinement transition in QCD. This is done through the implementation of Glauber dynamics of pure SU(3) lattice gauge theory. The coefficient Γ is extracted from the short-time exponential growth of the equal time correlation function of the order parameter. Although the absolute determination of Γ is ambiguous due to the difficulties in relating real time and Monte Carlo time, its relative temperature dependence can be obtained with much less arbitrariness. © 2007 American Institute of Physics.