Global dynamics of stationary solutions of the extended Fisher-Kolmogorov equation

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A Family of Stationary Solutions to the Euler Equations and Generalized Solutions

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

Theoretical study on the molecular mechanism for the reaction of VO2+ with C2H4

The complex reaction between VO2+ ((1)A(1)/(3)A) and C2H4 (Ag-1(g)/(3)A(1)) to yield VO+ ((1)Delta/(3)Sigma) and CH3CHO ('A'/(3)A) has been studied by means of B3LYP/6-31G* and B3LYP/6-311G(2d,p) calculations. The structures of all reactants, products, intermediates, and transition structures of this reaction have been optimized and characterized at the fundamental singlet and first excited triplet electronic states. Crossing points are localized, and possible spin inversion processes are discussed by means of the intrinsic reaction coordinate approach. Relevant stationary points along the most favorable reaction pathways have been studied at the CCSD/6-311G(2d,p)//B3LYP/6-311G(2d,p) calculation level. The theoretical results allow the development of thermodynamic and kinetic arguments about the reaction pathways of the title process. In the singlet state, the first step is the barrierless obtention of a reactant complex associated with the formation of a V-C bond, while in the triplet state a three-membered ring addition complex with the V bonded to the two C atoms is obtained. Similar behavior is found in the exit channels: the product complexes can be formed from isolated products without barriers. The reactant and product complexes are the most stable stationary points in the singlet and triplet electronic states. From the singlet state reactant complex, two reaction pathways are posssible to reach the triplet state product complex. (i) A mechanism in which a hydrogen transfer process is the first and rate limiting step and the second step is an oxygen transfer between vanadium and carbon atoms with a concomitant change in the spin state. The crossing point between singlet and triplet spin states is not kinetically relevant because it takes place at a later stage occurring in the exit channel. (ii) A mechanism in which the first stage renders a four-membered ring between vanadyl cation and the ethylene fragment and an oxygencarbon bond is formed; on going from this minimum to the second transition structure, associated with a carbon-vanadium bond breaking process, the crossing point between singlet and triplet spin states is reached. The final step is the hydrogen transfer between both carbon atoms to yield the product complex. In this case the spin change opens a lower barrier pathway. The transition structures with larger values of relative energies for both reactive channels of VO2+ ((1)A(1)) + C2H4 (Ag-1) --> VO+ ((3)Sigma) + CH3CHO ((1)A') present similar energies, and the two reaction pathways can be considered as competitive.

Energy-dependent point interaction: Self-adjointness

Recently, we constructed an energy-dependent point interaction (EDPI) in its most general form in one-dimensional quantum mechanics. In this paper, we show that stationary solutions of the Schrodinger equation with the EDPI form a complete set. Then any nonstationary solution of the time-dependent Schrodinger equation can be expressed as a linear combination of stationary solutions. This, however, does not necessarily mean that the EDPI is self-adjoint and the time-development of the nonstationary state is unitary. The EDPI is self-adjoint provided that the stationary solutions are all orthogonal to one another. We illustrate situations in which this orthogonality condition is not satisfied.

Fortran programs for the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

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

Self-trapping of a binary Bose-Einstein condensate induced by interspecies interaction

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

Dynamics in dumbbell domains I. Continuity of the set of equilibria

We analyze the dynamics of a reaction-diffusion equation with homogeneous Neumann boundary conditions in a dumbbell domain. We provide an appropriate functional setting to treat this problem and, as a first step, we show in this paper the continuity of the set of equilibria and of its linear unstable manifolds. (c) 2006 Elsevier B.V. All rights reserved.

Liquid-gas phase transition in Bose-Einstein condensates with time evolution

We study the effects of a repulsive three-body interaction on a system of trapped ultracold atoms in a Bose-Einstein condensed state. The stationary solutions of the corresponding s-wave nonlinear Schrödinger equation suggest a scenario of first-order liquid-gas phase transition in the condensed state up to a critical strength of the effective three-body force. The time evolution of the condensate with feeding process and three-body recombination losses has a different characteristic pattern. Also, the decay time of the dense (liquid) phase is longer than expected due to strong oscillations of the mean-squared radius.

PT-symmetric coupler with chi((2)) nonlinearity

We introduce the notion of a PT-symmetric dimer with a chi((2)) nonlinearity. Similarly to the Kerr case, we argue that such a nonlinearity should be accessible in a pair of optical waveguides with quadratic nonlinearity and gain and loss, respectively. An interesting feature of the problem is that because of the two harmonics, there exist in general two distinct gain and loss parameters, different values of which are considered herein. We find a number of traits that appear to be absent in the more standard cubic case. For instance, bifurcations of nonlinear modes from the linear solutions occur in two different ways depending on whether the first-or the second-harmonic amplitude is vanishing in the underlying linear eigenvector. Moreover, a host of interesting bifurcation phenomena appear to occur, including saddle-center and pitchfork bifurcations which our parametric variations elucidate. The existence and stability analysis of the stationary solutions is corroborated by numerical time-evolution simulations exploring the evolution of the different configurations, when unstable.

PT-symmetry management in oligomer systems

We study the effects of management of the PT-symmetric part of the potential within the setting of Schrodinger dimer and trimer oligomer systems. This is done by rapidly modulating in time the gain/loss profile. This gives rise to a number of interesting properties of the system, which are explored at the level of an averaged equation approach. Remarkably, this rapid modulation provides for a controllable expansion of the region of exact PT-symmetry, depending on the strength and frequency of the imposed modulation. The resulting averaged models are analysed theoretically and their exact stationary solutions are translated into time-periodic solutions through the averaging reduction. These are, in turn, compared with the exact periodic solutions of the full non-autonomous PT-symmetry managed problem and very good agreement is found between the two.

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

STABLE PIECEWISE POLYNOMIAL VECTOR FIELDS

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

Stable periodic waves in coupled Kuramoto-Sivashinsky-Korteweg-de Vries equations

Periodic waves are investigated in a system composed of a Kuramoto-Sivashinsky-Korteweg-de Vries (KS-KdV) equation linearly coupled to an extra linear dissipative one. The model describes, e.g., a two-layer liquid film flowing down an inclined plane. It has been recently shown that the system supports stable solitary pulses. We demonstrate that a perturbation analysis, based on the balance equation for the net field momentum, predicts the existence of stable cnoidal waves (CnWs) in the same system. It is found that the mean value u(0) of the wave field u in the main subsystem, but not the mean value of the extra field, affects the stability of the periodic waves. Three different areas can be distinguished inside the stability region in the parameter plane (L, u(0)), where L is the wave's period. In these areas, stable are, respectively, CnWs with positive velocity, constant solutions, and CnWs with negative velocity. Multistability, i.e., the coexistence of several attractors, including the waves with several maxima per period, appears at large value of L. The analytical predictions are completely confirmed by direct simulations. Stable waves are also found numerically in the limit of vanishing dispersion, when the KS-KdV equation goes over into the KS one.

Numerical and variational solutions of the dipolar Gross-Pitaevskii equation in reduced dimensions

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