Some experiments with high order compact methods using a computer algebra software - Part II (non-uniform grid)

We generalize a procedure proposed by Mancera and Hunt [P.F.A. Mancera, R. Hunt, Some experiments with high order compact methods using a computer algebra software-Part 1, Appl. Math. Comput., in press, doi: 10.1016/j.amc.2005.05.015] for obtaining a compact fourth-order method to the steady 2D Navier-Stokes equations in the streamfunction formulation-vorticity using the computer algebra system Maple, which includes conformal mappings and non-uniform grids. To analyse the procedure we have solved a constricted stepped channel problem, where a fine grid is placed near the re-entrant corner by transformation of the independent variables. (c) 2006 Elsevier B.V. All rights reserved.

A study of a numerical solution of the steady two dimensions Navier-Stokes equations in a constricted channel problem by a compact fourth order method

We present a numerical solution for the steady 2D Navier-Stokes equations using a fourth order compact-type method. The geometry of the problem is a constricted symmetric channel, where the boundary can be varied, via a parameter, from a smooth constriction to one possessing a very sharp but smooth corner allowing us to analyse the behaviour of the errors when the solution is smooth or near singular. The set of non-linear equations is solved by the Newton method. Results have been obtained for Reynolds number up to 500. Estimates of the errors incurred have shown that the results are accurate and better than those of the corresponding second order method. (C) 2002 Elsevier B.V. All rights reserved.

s-wave approximation for asymmetry in nonmesonic decay of finite hypernuclei

We establish the bridge between the commonly used Nabetani-Ogaito-Sato-Kishimoto (NOSK) formula for the asymmetry parameter a(Lambda) in the Lambda p -> np emission of polarized hypernuclei, and the shell-model (SM) formalism for finite hypernuclei. We demonstrate that the s-wave approximation leads to a SM formula for a(Lambda) that is as simple as the NOSK one and that reproduces the exact results for (5)(Lambda)He and (12)(Lambda)C better than initially expected. The simplicity achieved here is indeed remarkable. The new formalism makes the theoretical evaluation of a(Lambda) more transparent and explains clearly why the one-meson exchange model is unable to account for the experimental data of (5)(Lambda)He.

Low-energy direct muon transfer from H to Ne10+, S16+ and Ar18+ using the two-state close-coupling approximation to the Faddeev-Hahn-type equation

We perform a three-body calculation of direct muon-transfer rates from thermalized muonic hydrogen isotopes to bare nuclei Ne10+, S16+ and Ar18+ employing integro-differential Faddeev-Hahn-type equations in configuration space with a two-state close-coupling approximation scheme. All Coulomb potentials including the strong final-state Coulomb repulsion are treated exactly. A long-range polarization potential is included in the elastic channel to take into account the high polarizability of the muonic hydrogen. The transfer rates so-calculated are in good agreement with recent experiments. We find that the muon is captured predominantly in the n = 6, 9 and 10 states of muonic Ne10+, S16+ and Ar18+, respectively.

Charges and fluxes in Maxwell theory on compact manifolds with boundary

We investigate the charges and fluxes that can occur in higher-order Abelian gauge theories defined on compact space-time manifolds with boundary. The boundary is necessary to supply a destination to the electric lines of force emanating from brane sources, thus allowing non-zero net electric charges, but it also introduces new types of electric and magnetic flux. The resulting structure of currents, charges, and fluxes is studied and expressed in the language of relative homology and de Rham cohomology and the corresponding abelian groups. These can be organised in terms of a pair of exact sequences related by the Poincare-Lefschetz isomorphism and by a weaker flip symmetry exchanging the ends of the sequences. It is shown how all this structure is brought into play by the imposition of the appropriately generalised Maxwell's equations. The requirement that these equations be integrable restricts the world-volume of a permitted brane (assumed closed) to be homologous to a cycle on the boundary of space-time. All electric charges and magnetic fluxes are quantised and satisfy the Dirac quantisation condition. But through some boundary cycles there may be unquantised electric fluxes associated with quantised magnetic fluxes and so dyonic in nature.

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

Combined D0 measurements constraining the CP-violating phase and width difference in the B-s(0) system

We combine the D0 measurement of the width difference between the light and heavy B-s(0) mass eigenstates and of the CP-violating mixing phase determined from the time-dependent angular distributions in the B-s(0)-> J/psi phi decays along with the charge asymmetry in semileptonic decays also measured with the D0 detector. With the additional constraint from the world average of the flavor-specific B-s(0) lifetime, we obtain Delta Gamma(s)equivalent to(Gamma(L)-Gamma(H))=0.13 +/- 0.09 ps(-1) and vertical bar phi(s)vertical bar=0.70(-0.47)(+0.39) or Delta Gamma(s)=-0.13 +/- 0.09 ps(-1) and vertical bar phi(s)vertical bar=2.44(-0.39)(+0.47). The data sample corresponds to an integrated luminosity of 1.1 fb(-1) accumulated with the D0 detector at the Fermilab Tevatron Collider.

Finite-volume form factors in semi-classical approximation

A semi-classical approach is used to obtain Lorentz covariant expressions for the form factors between the kink states of a quantum field theory with degenerate vacua. Implemented on a cylinder geometry it provides an estimate of the spectral representation of correlation functions in a finite volume. Illustrative examples of the applicability of the method are provided by the sine-Gordon and the broken phi(4) theories in 1 + 1 dimensions. (C) 2003 Elsevier B.V. All rights reserved.

A dynamical gluon mass solution in Mandelstam's approximation

We discuss the pure gauge Schwinger-Dyson equation for the gluon propagator in the Landau gauge within an approximation proposed by Mandelstam many years ago. We show that a dynamical gluon mass arises as a solution. This solution is obtained numerically in the full range of momenta that we have considered without the introduction of any ansatz or asymptotic expression in the infrared region. The vertex function that we use follows a prescription formulated by Cornwall to determine the existence of a dynamical gluon mass in the light cone gauge. The renormalization procedure differs from the one proposed by Mandelstam and allows for the possibility of a dynamical gluon mass. Some of the properties of this solution, such as its dependence on A(QCD) and its perturbative scaling behavior are also discussed.

Ps-He scattering using the three-state positronium close-coupling approximation

A three-state target elastic positronium close-coupling approximation (CCA) is employed to investigate Ps-He scattering in the energy range 0-200 eV with and without electron exchange. Low-lying phase shifts below the first excitation threshold and the total integrated cross sections using both the models are reported. Estimation of integrated excitation cross sections for Ps(1s --> 2s) and Ps(1s --> 2p) using CCA are presented for the first time. The present total cross sections are in good agreement with the measured data in the incident Ps energy range 20-30 eV.

Measurement of the lifetime difference in the B-s(0) system

We present a study of the decay B-s(0)-> J/psi phi. We obtain the CP-odd fraction in the final state at time zero R-perpendicular to = 0.16 +/- 0.10 (stat) +/- 0.02(syst), the average lifetime of the (B-s(0), (B) over bar (0)(s)) system, (tau) over bar (B-s(0)) = 1.39(-0.16)(+0.13)(stat)(-0.02)(+0.01)(syst) ps, and the relative width difference between the heavy and light mass eigen-states, Delta Gamma/(Gamma) over bar = (Gamma(L) - Gamma(H))/(Gamma) over bar = 0.24(-0.38)(+0.28)(stat)(-0.04)(+0.03)(syst). With the additional constraint from the world average of the B-s(0) lifetime measurements using semileptonic decays, we find (tau) over bar (B-s(0)) = 1.39 +/- 0.06 ps and Delta Gamma/(Gamma) over bar = 0.25(-0.15)(+0.14). For the ratio of the B-s(0) and B-0 lifetimes we obtain (tau) over bar (B-s(0))/tau(B-s(0)) = 0.91 +/- 0.09(stat) +/- 0.003(syst).

Miscibility in a degenerate fermionic mixture induced by linear coupling

We consider a one-dimensional mean-field-hydrodynamic model of a two-component degenerate Fermi gas in an external trap, each component representing a spin state of the same atom. We demonstrate that the interconversion between them (linear coupling), imposed by a resonant electromagnetic wave, transforms the immiscible binary gas into a miscible state, if the coupling constant, kappa, exceeds a critical value, kappa(cr). The effect is predicted in a variational approximation, and confirmed by numerical solutions. Unlike the recently studied model of a binary Bose-Einsten condensate with the linear coupling, the components in the immiscible phase of the binary fermion mixture never fill two separated domains with a wall between them, but rather form antilocked (pi-phase-shifted) density waves. Another difference from the bosonic mixture is spontaneous breaking of symmetry between the two components in terms of the numbers of atoms in them, N(1) and N(2). The latter effect is characterized by the parameter nu equivalent to(N(1)-N(2))/(N(1)+N(2)) (only N(1)+N(2) is a conserved quantity), the onset of miscibility at kappa >=kappa(cr) meaning a transition to nu equivalent to 0. At kappa difference between the two components. The relation between the imbalance in the spin population, induced by the linear coupling, and the developing spatial structure resembles the known Larkin-Ovchinnikov-Fulde-Ferrell states in the Fermi mixture. Dynamical states, when kappa is suddenly switched from zero to a value exceeding kappa(cr), are considered too. In the latter case, the system features oscillatory relaxation to the mixed state.

Lifetime difference and CP-violating phase in the B-s(0) system

From an analysis of the decay B-s(0)-> J/psi phi, we obtain the width difference between the light and heavy mass eigenstates Delta Gamma equivalent to(Gamma(L)-Gamma(H))=0.17 +/- 0.09(stat)+/- 0.02(syst) ps(-1) and the CP-violating phase phi(s)=-0.79 +/- 0.56(stat)(-0.01)(+0.14)(syst). Under the hypothesis of no CP violation (phi(s)equivalent to 0), we obtain 1/Gamma=tau/(B-s(0))=1.52 +/- 0.08(stat)(-0.03)(+0.01)(syst) ps and Delta Gamma=0.12(-0.10)(+0.08)(stat)+/- 0.02(syst) ps(-1). The data sample corresponds to an integrated luminosity of about 1.1 fb(-1) accumulated with the D0 detector at the Fermilab Tevatron collider. This is the first direct measurement of the CP-violating mixing phase in the B-s(0) system.

Gap solitons in a model of a superfluid fermion gas in optical lattices

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

Weak-field approximation of effective gravitational theory with local Galilean invariance

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