A coupling technique for non-matching finite element meshes

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

Position-dependent effective mass Dirac equations with PT-symmetric and non-PT-symmetric potentials

We present a new procedure to construct the one-dimensional non-Hermitian imaginary potential with a real energy spectrum in the context of the position-dependent effective mass Dirac equation with the vector-coupling scheme in 1 + 1 dimensions. In the first example, we consider a case for which the mass distribution combines linear and inversely linear forms, the Dirac problem with a PT-symmetric potential is mapped into the exactly solvable Schrodinger-like equation problem with the isotonic oscillator by using the local scaling of the wavefunction. In the second example, we take a mass distribution with smooth step shape, the Dirac problem with a non-PT-symmetric imaginary potential is mapped into the exactly solvable Schrodinger-like equation problem with the Rosen-Morse potential. The real relativistic energy levels and corresponding wavefunctions for the bound states are obtained in terms of the supersymmetric quantum mechanics approach and the function analysis method.

Extension of PT-symmetric quantum mechanics to the Dirac theory with position-dependent mass

We present a new method to construct the exactly solvable PT-symmetric potentials within the framework of the position-dependent effective mass Dirac equation with the vector potential coupling scheme in 1 + 1 dimensions. In order to illustrate the procedure, we produce three PT-symmetric potentials as examples, which are PT-symmetric harmonic oscillator-like potential, PT-symmetric potential with the form of a linear potential plus an inversely linear potential, and PT-symmetric kink-like potential, respectively. The real relativistic energy levels and corresponding spinor components for the bound states are obtained by using the basic concepts of the supersymmetric quantum mechanics formalism and function analysis method. (C) 2007 Elsevier B.V. All rights reserved.

Differential cross sections for target excitation and positronium formation in positron-helium scattering

Positronium formation and target excitation in positron-helium scattering have been investigated using the close-coupling approximation with realistic wave functions for the positronium and helium atoms. The following eight states have been used in the close-coupling scheme: He(1s1s), He(1s2(1)s), He(1s2(1)p), He(1s3(1)s), He(1s3(1)p), Ps(1s), Ps(2s), and Ps(2p), where Ps stands for the positronium atom. Calculations are reported of differential cross sections for elastic scatering,, inelastic target excitation to He(1s2(1)s) and He(1s2(1)p) slates, and rearrangement transition to Ps(1s), Ps(2s), and Ps(2p) states for incident positron energies between 40 and 200 eV. The coincidence parameters for the transition to the He(1s2(1)p) state of helium are also reported and briefly discussed. [S1050-2947(98)05101-4].

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.

QCD tests with an infrared finite gluon propagator and coupling constant

Recent progress in the solution of Schwinger-Dyson equations, as well as lattice simulation of pure glue QCD, indicate that the gluon propagator and coupling constant are infrared finite. Such non-perturbative information can be introduced in the QCD perturbative expansion in the scheme named Dynamical Perturbation Theory. We exemplify this procedure with the calculation of some two-body non-leptonic annihilation B meson decays, which show agreement with the experimental data in the case of a gluon propagator characterized by a dynamical gluon mass of 500MeV, compatible with the value found in several processes computed with this method. We give a. preliminary account of the application of this procedure at the loop level in the case of the Bjorken sum rule.

Extension of the Nambu-Jona-Lasinio model predictions at high densities and temperatures using an implicit regularization scheme

Traditional cutoff regularization schemes of the Nambu-Jona-Lasinio model limit the applicability of the model to energy-momentum scales much below the value of the regularizing cutoff. In particular, the model cannot be used to study quark matter with Fermi momenta larger than the cutoff. In the present work, an extension of the model to high temperatures and densities recently proposed by Casalbuoni, Gatto, Nardulli, and Ruggieri is used in connection with an implicit regularization scheme. This is done by making use of scaling relations of the divergent one-loop integrals that relate these integrals at different energy-momentum scales. Fixing the pion decay constant at the chiral symmetry breaking scale in the vacuum, the scaling relations predict a running coupling constant that decreases as the regularization scale increases, implementing in a schematic way the property of asymptotic freedom of quantum chromodynamics. If the regularization scale is allowed to increase with density and temperature, the coupling will decrease with density and temperature, extending in this way the applicability of the model to high densities and temperatures. These results are obtained without specifying an explicit regularization. As an illustration of the formalism, numerical results are obtained for the finite density and finite temperature quark condensate and applied to the problem of color superconductivity at high quark densities and finite temperature.

Left-right asymmetry and minimal coupling

In this paper we deal with an alternative approach to the description of massless particles of arbitrary spin. Within this scheme chiral components of a spinor field are regarded as fundamental quantities and treated as independent field variables. The free field Lagrangian is built up from the requirement of chiral invariance; This formulation is parallel to the neutrino theory and allows for a formulation that generalizes, to particles of arbitrary spin, the two-component neutrino theory. We achieve a spinor formulation of electrodynamics. In the case of the photon, the nonzero helicity components satisfy Weyl's equations and are associated to observables (electromagnetic fields) whereas the zero helicity components are related to nonobservables (electromagnetic potentials). Within the spinor formulation of electrodynamics the minimal coupling substitution follows as a consequence of the linearity of the interaction and the preference of nature for chiral components, that is, of the left-right asymmetry of nature. (C) 1996 American Institute of Physics.

The temperature dependence of the QCD running coupling

We study the running of the QCD coupling with the momentum squared (Q 2) and the temperature scales in the high temperature limit (T > Tc), using a mass dependent renormalization scheme to build the Renormalization Group Equations. The approach used guaranty gauge invariance, through the use of the Hard Thermal Loop approximation, and independence of the vertex chosen to renormalize the coupling. In general, the dependence of the coupling with the temperature is not logarithmical, although in the region Q2 ∼ T2 the logarithm approximation is reasonable. Finally, as known from Debye screening, color charge is screened in the coupling. The number of flavors, however, is anti-screened.

Development of a quality index method (QIM) sensory scheme and study of shelf-life of ice-stored blackspot seabream (Pagellus bogaraveo)

The aim of this work was to develop a quality index method (QIM) scheme for whole ice-boxed refrigerated blackspot seabream and to perform shelf-life evaluations, using sensory analysis, GR Torrymeter measurements and bacterial counts of specific spoilage organisms (SSO) during chilled storage. A QIM scheme based on a total of 30 demerit points was developed. Sensory, physical and microbiological data were integrated and used to determine the rejection point. Results indicated that the shelf-life of blackspot seabream is around 12-13 days. (C) 2011 Elsevier Ltd. All rights reserved.

Coupling Surface-Enhanced Resonance Raman Scattering and Electronic Tongue as Characterization Tools to Investigate Biological Membrane Mimetic Systems

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

A bounded upwinding scheme for computing convection-dominated transport problems

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

A continuously differentiable upwinding scheme for the simulation of fluid flow problems

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

Changing inter-molecular spin-orbital coupling for generating magnetic field effects in phosphorescent organic semiconductors

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

A fully adaptive multiresolution scheme for shock computations

The scheme is based on Ami Harten's ideas (Harten, 1994), the main tools coming from wavelet theory, in the framework of multiresolution analysis for cell averages. But instead of evolving cell averages on the finest uniform level, we propose to evolve just the cell averages on the grid determined by the significant wavelet coefficients. Typically, there are few cells in each time step, big cells on smooth regions, and smaller ones close to irregularities of the solution. For the numerical flux, we use a simple uniform central finite difference scheme, adapted to the size of each cell. If any of the required neighboring cell averages is not present, it is interpolated from coarser scales. But we switch to ENO scheme in the finest part of the grids. To show the feasibility and efficiency of the method, it is applied to a system arising in polymer-flooding of an oil reservoir. In terms of CPU time and memory requirements, it outperforms Harten's multiresolution algorithm.The proposed method applies to systems of conservation laws in 1Dpartial derivative(t)u(x, t) + partial derivative(x)f(u(x, t)) = 0, u(x, t) is an element of R-m. (1)In the spirit of finite volume methods, we shall consider the explicit schemeupsilon(mu)(n+1) = upsilon(mu)(n) - Deltat/hmu ((f) over bar (mu) - (f) over bar (mu)-) = [Dupsilon(n)](mu), (2)where mu is a point of an irregular grid Gamma, mu(-) is the left neighbor of A in Gamma, upsilon(mu)(n) approximate to 1/mu-mu(-) integral(mu-)(mu) u(x, t(n))dx are approximated cell averages of the solution, (f) over bar (mu) = (f) over bar (mu)(upsilon(n)) are the numerical fluxes, and D is the numerical evolution operator of the scheme.According to the definition of (f) over bar (mu), several schemes of this type have been proposed and successfully applied (LeVeque, 1990). Godunov, Lax-Wendroff, and ENO are some of the popular names. Godunov scheme resolves well the shocks, but accuracy (of first order) is poor in smooth regions. Lax-Wendroff is of second order, but produces dangerous oscillations close to shocks. ENO schemes are good alternatives, with high order and without serious oscillations. But the price is high computational cost.Ami Harten proposed in (Harten, 1994) a simple strategy to save expensive ENO flux calculations. The basic tools come from multiresolution analysis for cell averages on uniform grids, and the principle is that wavelet coefficients can be used for the characterization of local smoothness.. Typically, only few wavelet coefficients are significant. At the finest level, they indicate discontinuity points, where ENO numerical fluxes are computed exactly. Elsewhere, cheaper fluxes can be safely used, or just interpolated from coarser scales. Different applications of this principle have been explored by several authors, see for example (G-Muller and Muller, 1998).Our scheme also uses Ami Harten's ideas. But instead of evolving the cell averages on the finest uniform level, we propose to evolve the cell averages on sparse grids associated with the significant wavelet coefficients. This means that the total number of cells is small, with big cells in smooth regions and smaller ones close to irregularities. This task requires improved new tools, which are described next.