Hybrid Monte Carlo algorithm for the double exchange model

The Hybrid Monte Carlo algorithm is adapted to the simulation of a system of classical degrees of freedom coupled to non self-interacting lattices fermions. The diagonalization of the Hamiltonian matrix is avoided by introducing a path-integral formulation of the problem, in d + 1 Euclidean space–time. A perfect action formulation allows to work on the continuum Euclidean time, without need for a Trotter–Suzuki extrapolation. To demonstrate the feasibility of the method we study the Double Exchange Model in three dimensions. The complexity of the algorithm grows only as the system volume, allowing to simulate in lattices as large as 163 on a personal computer. We conclude that the second order paramagnetic–ferromagnetic phase transition of Double Exchange Materials close to half-filling belongs to the Universality Class of the three-dimensional classical Heisenberg model.

Variational mean-field approach to the double-exchange model

It has been recently shown that the double exchange Hamiltonian, with weak antiferromagnetic interactions, has a richer variety of first- and second-order transitions than previously anticipated, and that such transitions are consistent with the magnetic properties of manganites. Here we present a thorough discussion of the variational mean-field approach that leads to these results. We also show that the effect of the Berry phase turns out to be crucial to produce first-order paramagnetic-ferromagnetic transitions near half filling with transition temperatures compatible with the experimental situation. The computation relies on two crucial facts: the use of a mean-field ansatz that retains the complexity of a system of electrons with off-diagonal disorder, not fully taken into account by the mean-field techniques, and the small but significant antiferromagnetic superexchange interaction between the localized spins.

Consistent nucleon-nucleon potentials and three-body forces

We construct microscopic three-nucleon forces consistent with the Bonn and Nijmegen two-nucleon potentials, and including , Roper, and nucleon-antinucleon excitations. Recent results for the choice of the meson parameters are discussed. The forces are used in Brueckner calculations and the saturation properties of nuclear matter are determined.

The reaction (p)over-barp -> (Lambda)over-bar(c)(-)Lambda(+)(c) close to threshold

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

Short-range repulsion and isospin dependence in the kaon-nucleon (KN) system

The short-range properties of the kaon-nucleon (KN) interaction are studied within the meson-exchange model of the Jülich group. Specifically, dynamical explanations for the phenomenological short-range repulsion, required in this model for achieving agreement with the empirical KN data, are explored. Evidence is found that contributions from the exchange of a heavy scalar-isovector meson [a0(980)] as well as from genuine quark-gluon exchange processes are needed. Taking both mechanisms into account, a satisfactory description of the KN phase shifts can be obtained without resorting to phenomenological pieces.

Influence of a Z+(1540) resonance on K+N scattering

The impact of a (I=0, JP=1/2+)Z+(1540) resonance with a width of 5 MeV or more on the K+N(I=0) elastic cross section and on the P01 phase shift is examined within the KN meson-exchange model of the Jülich group. It is shown that the rather strong enhancement of the cross section caused by the presence of a Z + with the above properties is not compatible with the existing empirical information on KN scattering. Only a much narrower Z+ state could be reconciled with the existing data - or, alternatively, the Z + state must lie at an energy much closer to the KN threshold.

Studies on Pseudoscalar Meson Bound States and Semileptonic Decays in a Relativistic Potential Model

In this thesis quark-antiquark bound states are considered using a relativistic two-body equation for Dirac particles. The mass spectrum of mesons includes bound states involving two heavy quarks or one heavy and one light quark. In order to analyse these states within a unified formalism, it is desirable to have a two-fermion equation that limits to one body Dirac equation with a static interaction for the light quark when the other particle's mass tends to infinity. A suitable two-body equation has been developed by Mandelzweig and Wallace. This equation is solved in momentum space and is used to describe the complete spectrum of mesons. The potential used in this work contains a short range one-gluon exchange interaction and a long range linear confining and constant potential terms. This model is used to investigate the decay processes of heavy mesons. Semileptonic decays are more tractable since there is no final state interactions between the leptons and hadrons that would otherwise complicate the situation. Studies on B and D meson decays are helpful to understand the nonperturbative strong interactions of heavy mesons, which in turn is useful to extract the details of weak interaction process. Calculation of form factors of these semileptonic decays of pseudo scalar mesons are also presented.

Quarkonia and heavy-light mesons in a covariant quark model

Preliminary calculations using the Covariant Spectator Theory (CST) employed a scalar linear confining interaction and an additional constant vector potential to compute the mesonic mass spectra. In this work we generalize the confining interaction to include more general structures, in particular a vector and also a pseudoscalar part, as suggested by a recent study. A one-gluon-exchange kernel is also implemented to describe the short-range part of the interaction. We solve the simplest CST approxima- tion to the complete Bethe-Salpeter equation, the one-channel spectator equation, using a numerical technique that eliminates all singularities from the kernel. The parameters of the model are determined through a fit to the experimental pseudoscalar meson spectra, with a good agreement for both quarkonia and heavy-light states.

Heavy-light mesons on a lattice

Several excited states of Ds and Bs mesons have been discovered in the last six years: BaBar, Cleo and Belle discovered the very narrow states D(s0)*(2317)+- and D(s1)(2460)+- in 2003, and CDF and DO Collaborations reported the observation of two narrow Bs resonances, B(s1)(5830)0 and B*(s2)(5840)0 in 2007. To keep up with experiment, meson excited states should be studied from the theoretical aspect as well. The theory that describes the interaction between quarks and gluons is quantum chromodynamics (QCD). In this thesis the properties of the meson states are studied using the discretized version of the theory - lattice QCD. This allows us to perform QCD calculations from first principles, and "measure" not just energies but also the radial distributions of the states on the lattice. This gives valuable theoretical information on the excited states, as we can extract the energy spectrum of a static-light meson up to D wave states (states with orbital angular momentum L=2). We are thus able to predict where some of the excited meson states should lie. We also pay special attention to the order of the states, to detect possible inverted spin multiplets in the meson spectrum, as predicted by H. Schnitzer in 1978. This inversion is connected to the confining potential of the strong interaction. The lattice simulations can also help us understand the strong interaction better, as the lattice data can be treated as "experimental" data and used in testing potential models. In this thesis an attempt is made to explain the energies and radial distributions in terms of a potential model based on a one-body Dirac equation. The aim is to get more information about the nature of the confining potential, as well as to test how well the one-gluon exchange potential explains the short range part of the interaction.

Carbon dioxide and methane exchange between a boreal pristine lake and the atmosphere

Lakes serve as sites for terrestrially fixed carbon to be remineralized and transferred back to the atmosphere. Their role in regional carbon cycling is especially important in the Boreal Zone, where lakes can cover up to 20% of the land area. Boreal lakes are often characterized by the presence of a brown water colour, which implies high levels of dissolved organic carbon from the surrounding terrestrial ecosystem, but the load of inorganic carbon from the catchment is largely unknown. Organic carbon is transformed to methane (CH4) and carbon dioxide (CO2) in biological processes that result in lake water gas concentrations that increase above atmospheric equilibrium, thus making boreal lakes as sources of these important greenhouse gases. However, flux estimates are often based on sporadic sampling and modelling and actual flux measurements are scarce. Thus, the detailed temporal flux dynamics of greenhouse gases are still largely unknown. ----- One aim here was to reveal the natural dynamics of CH4 and CO2 concentrations and fluxes in a small boreal lake. The other aim was to test the applicability of a measuring technique for CO2 flux, i.e. the eddy covariance (EC) technique, and a computational method for estimation of primary production and community respiration, both commonly used in terrestrial research, in this lake. Continuous surface water CO2 concentration measurements, also needed in free-water applications to estimate primary production and community respiration, were used over two open water periods in a study of CO2 concentration dynamics. Traditional methods were also used to measure gas concentration and fluxes. The study lake, Valkea-Kotinen, is a small, humic, headwater lake within an old-growth forest catchment with no local anthropogenic disturbance and thus possible changes in gas dynamics reflect the natural variability in lake ecosystems. CH4 accumulated under the ice and in the hypolimnion during summer stratification. The surface water CH4 concentration was always above atmospheric equilibrium and thus the lake was a continuous source of CH4 to the atmosphere. However, the annual CH4 fluxes were small, i.e. 0.11 mol m-2 yr-1, and the timing of fluxes differed from that of other published estimates. The highest fluxes are usually measured in spring after ice melt but in Lake Valkea-Kotinen CH4 was effectively oxidised in spring and highest effluxes occurred in autumn after summer stratification period. CO2 also accumulated under the ice and the hypolimnetic CO2 concentration increased steadily during stratification period. The surface water CO2 concentration was highest in spring and in autumn, whereas during the stable stratification it was sometimes under atmospheric equilibrium. It showed diel, daily and seasonal variation; the diel cycle was clearly driven by light and thus reflected the metabolism of the lacustrine ecosystem. However, the diel cycle was sometimes blurred by injection of hypolimnetic water rich in CO2 and the surface water CO2 concentration was thus controlled by stratification dynamics. The highest CO2 fluxes were measured in spring, autumn and during those hypolimnetic injections causing bursts of CO2 comparable with the spring and autumn fluxes. The annual fluxes averaged 77 (±11 SD) g C m-2 yr-1. In estimating the importance of the lake in recycling terrestrial carbon, the flux was normalized to the catchment area and this normalized flux was compared with net ecosystem production estimates of -50 to 200 g C m-2 yr-1 from unmanaged forests in corresponding temperature and precipitation regimes in the literature. Within this range the flux of Lake Valkea-Kotinen yielded from the increase in source of the surrounding forest by 20% to decrease in sink by 5%. The free water approach gave primary production and community respiration estimates of 5- and 16-fold, respectively, compared with traditional bottle incubations during a 5-day testing period in autumn. The results are in parallel with findings in the literature. Both methods adopted from the terrestrial community also proved useful in lake studies. A large percentage of the EC data was rejected, due to the unfulfilled prerequisites of the method. However, the amount of data accepted remained large compared with what would be feasible with traditional methods. Use of the EC method revealed underestimation of the widely used gas exchange model and suggests simultaneous measurements of actual turbulence at the water surface with comparison of the different gas flux methods to revise the parameterization of the gas transfer velocity used in the models.

Persistent currents in a one-dimensional ring for a disordered Hubbard model

We consider a one-dimensional Hubbard model in the presence of disorder. We compute the charge stiffness for a mesoscopic ring as a function of the size L, which is a measure of the persistent currents. We find that for finite disorder the persistent currents of the system with repulsive interactions are larger than those of the system with attractive interactions. This counterintuitive result is due to the fact that local-density fluctuations are reduced in the presence of repulsive interactions.

eta-production on the proton via electromagnetic and hadronic probes

The reactions pi(-)p -> eta n and gamma p -> eta p are investigated within a dynamical coupled-channels model of meson production reactions in the nucleon resonance region The meson-baryon channels included are pi N, pi Delta, sigma N, and rho N The direct eta-photoproduction process is studied within a formalism based on a chiral constituent quark model approach, complemented with a one-gluon-exchange mechanism, to take into account the breakdown of the SU(6)circle times O(3) symmetry In the models search, the following known nucleon resonances are embodied S-11(1535), S-11(1650), P-11(1440), P-11(1710), P-13(1720), D-13(1520), D-13(1700), D-15(1675), and F-15(1680). Data for the pi(-)p -> eta n reaction from threshold up to a total center-of-mass energy of W approximate to 2 GeV are satisfactorily reproduced For the photoproduction channel: two additional higher mass known resonances, P-13(1900) and F-15(2000), are also considered However, reproducing the data for gamma p -> eta p requires, within our approach, two new nucleon resonances, for which we extract, mass and width

Phenomenological analysis of the double-pion production in nucleon-nucleon collisions up to 2.2 GeV

With an effective Lagrangian approach, we analyze several NN -> NN pi pi channels by including various resonances with mass up to 1.72 GeV. For the channels with the pion pair of isospin zero, we confirm the dominance of N*(1440) -> N sigma in the near-threshold region. At higher energies and for channels with the final pion pair of isospin one, we find large contributions from N*(1440) -> Delta pi, double-Delta, Delta(1600) -> N*(1440)pi, Delta(1600) -> Delta pi and Delta(1620) -> Delta pi. There are also sizable contributions from Delta -> Delta pi, Delta -> N pi, N -> Delta pi, and nucleon pole at energies close to the threshold. We give a good reproduction to the total cross sections up to beam energies of 2.2 GeV except for the pp -> pp pi(0)pi(0) channel at energies around 1.1 GeV and our results agree with the existing data of differential cross sections of pp -> pp pi(+)p pi(-), pp -> nn pi(+)pi(+), and pp -> pp pi(0)pi(0) which are measured at CELSIUS and COSY.

Antikaon potential in hot and dense matter

The antikaon optical potential in hot and dense nuclear matter is studied within the framework of a coupled-channel self-consistent calculation taking, as bare meson-baryon interaction, the meson-exchange potential of the Jlich group. Typical conditions found in heavy-ion collisions at GSI are explored. As in the case of zero temperature, the angular momentum components larger than L=0 contribute significantly to the finite temperature antikaon optical potential at finite momentum. It is found that the particular treatment of the medium effects has a strong influence on the behavior of the antikaon potential with temperature. Our self-consistent model, in which antikaons and pions are dressed in the medium, gives a moderately temperature dependent antikaon potential which remains attractive at GSI temperatures, contrary to what one finds if only nuclear Pauli blocking effects are included.

A mass-transportation approach to a one dimensional fluid mechanics model with nonlocal velocity

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