An introduction to numerical methods in low-dimensional quantum systems

This is an introductory course to the Lanczos Method and Density Matrix Renormalization Group Algorithms (DMRG), two among the leading numerical techniques applied in studies of low-dimensional quantum models. The idea of studying the models on clusters of a finite size in order to extract their physical properties is briefly discussed. The important role played by the model symmetries is also examined. Special emphasis is given to the DMRG.

Development and applications of three-dimensional gamma ray tomography system using ray casting volume rendering techniques

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

Study of stress distribution with finite element method in a root with prefabricated threaded split shaft post - Part I

A finite element analysis was carried out to study the role of prefabricated threaded split shaft post (Flexi-Post) on dentinal stress in pulpless tooth. Three dimensional plane strain model of mesio-distal section of a human maxillary central incisor without restoration was analysed with the MSC/NASTRAN (MacNeal/ Schwendler) general purpose finite analysis program was executed on a microcomputer. The model as discretized into 48.954 axisymmetric finite elements defined by 10.355 nodes. Each element was assigned unique elastic properties to represent the materials modeled. Homogeneity, isotropy and linear elasticity were assume for all material. A simulation of static load of 100N was applied to the incisal edge of the post; vertical. Maximal principal stresses and von Mises equivalent stress were calculated. Using the element analysis model employed in this study, the following can be concluded concerning threaded split shaft post (Flexi-Post): Maximum principal stresses in dentin were located at cervical place and at the post apex. The apical threads of the post not redirecting stresses away from the root.

Phase separation and gap bowing in zinc-blende InGaN, InAlN, BGaN, and BAlN alloy layers

We present first-principles calculations of the thermodynamic and electronic properties of the zinc-blende ternary InxGa1-xN. InxAl1-xN, BxGa1-xN, and BxAl1-xN alloys. They are based on a generalized quasi-chemical approximation and a pseudopotential-plane-wave method. T-x phase diagrams for the alloys are obtained, We show that due to the large difference in interatomic distances between the binary compounds a significant phase miscibility gap for the alloys is found. In particular for the InxGa1-xN alloy, we show also experimental results obtained from X-ray and resonant Raman scattering measurements, which indicate the presence of an In-rich phase with x approximate to 0.8. For the boron-containing alloy layers we found a very high value for the critical temperature for miscibility. similar to9000 K. providing an explanation for the difficulties encountered to grow these materials with higher boron content. The influence of a biaxial strain on phase diagrams, energy gaps and gap bowing of these alloys is also discussed. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Energy localization in the Phi(4) oscillator chain

We study energy localization in a finite one-dimensional Phi(4) oscillator chain with initial energy in a single oscillator of the chain. We numerically calculate the effective number of degrees of freedom sharing the energy on the lattice as a function of time. We find that for energies smaller than a critical value, energy equipartition among the oscillators is reached in a relatively short time. on the other hand, above the critical energy, a decreasing number of particles sharing the energy is observed. We give an estimate of the effective number of degrees of freedom as a function of the energy. Our results suggest that localization is due to the appearance, above threshold, of a breather-like structure. Analytic arguments are given, based on the averaging theory and the analysis of a discrete nonlinear Schrodinger equation approximating the dynamics, to support and explain the numerical results.

Spin, charge, and orbital correlations in the one-dimensional t(2g)-orbital Hubbard model

We present the zero-temperature phase diagram of the one-dimensional t(2g)-orbital Hubbard model, obtained using the density-matrix renormalization group and Lanczos techniques. Emphasis is given to the case of the electron density n=5 corresponding to five electrons per site, while several other cases for electron densities between n=3 and 6 are also studied. At n=5, our results indicate a first-order transition between a paramagnetic (PM) insulator phase, with power-law slowly decaying correlations, and a fully polarized ferromagnetic (FM) state by tuning the Hund's coupling. The results also suggest a transition from the n=5 PM insulator phase to a metallic regime by changing the electron density, either via hole or electron doping. The behavior of the spin, charge, and orbital correlation functions in the FM and PM states are also described in the text and discussed. The robustness of these two states against varying parameters suggests that they may be of relevance in quasi-one-dimensional Co-oxide materials, or even in higher dimensional cobaltite systems as well.

BERRY PHASE, AHARONOV-BOHM EFFECT AND TOPOLOGY

The Berry phase for an electron in a one-dimensional box rotated around a magnetic flux line has contributions from the geometry and the magnetic flux, which gives an Aharonov-Bohm effect. For a circular box enclosing the magnetic flux, the Berry phase depends on the boundary conditions.

NUMERICAL DETERMINATION OF THE ORDER OF PHASE-TRANSITIONS IN ISING SYSTEMS WITH MULTISPIN INTERACTIONS

The method of the fourth-order cumulant of Challa, Landau, and Binder is used together with the Monte Carlo histogram technique of Ferrenberg and Swendsen to study the order of the phase transitions of two-dimensional Ising systems with multispin interactions in the horizontal direction and two-body interactions in the vertical direction.

CLOSED SPACES IN COSMOLOGY

This paper deals with two aspects of relativistic cosmologies with closed spatial sections. These spacetimes are based on the theory of general relativity, and admit a foliation into space sections S(t), which are spacelike hypersurfaces satisfying the postulate of the closure of space: each S(t) is a three-dimensional closed Riemannian manifold. The topics discussed are: (i) a comparison, previously obtained, between Thurston geometries and Bianchi-Kantowski-Sachs metrics for such three-manifolds is here clarified and developed; and (ii) the implications of global inhomogeneity for locally homogeneous three-spaces of constant curvature are analyzed from an observational viewpoint.

CHIRAL SCHWINGER MODEL WITH A PODOLSKY TERM AT FINITE-TEMPERATURE

We study an exactly solvable two-dimensional model which mimics the basic features of the standard model. This model combines chiral coupling with an infrared behavior which resembles low energy QCD. This is done by adding a Podolsky higher-order derivative term in the gauge field to the Lagrangian of the usual chiral Schwinger model. We adopt a finite temperature regularization procedure in order to calculate the non-trivial fermionic Jacobian and obtain the photon and fermion propagators, first at zero temperature and then at finite temperature in the imaginary and real time formalisms. Both singular and non-singular cases, corresponding to the choice of the regularization parameter, are treated. In the nonsingular case there is a tachyonic mode as usual in a higher order derivative theory, however in the singular case there is no tachyonic excitation in the spectrum.

Designing spatial correlation of quantum dots: towards self-assembled three-dimensional structures

Buried two-dimensional arrays of InP dots were used as a template for the lateral ordering of self-assembled quantum dots. The template strain field can laterally organize compressive (InAs) as well as tensile (GaP) self-assembled nanostructures in a highly ordered square lattice. High-resolution transmission electron microscopy measurements show that the InAs dots are vertically correlated to the InP template, while the GaP dots are vertically anti-correlated, nucleating in the position between two buried InP dots. Finite InP dot size effects are observed to originate InAs clustering but do not affect GaP dot nucleation. The possibility of bilayer formation with different vertical correlations suggests a new path for obtaining three-dimensional pseudocrystals.

Numerical simulation of viscous three-dimensional flow over aerospace vehicles using Euler/boundary-layer zonal approach

This paper presents a viscous three-dimensional simulations coupling Euler and boundary layer codes for calculating flows over arbitrary surfaces. The governing equations are written in a general non orthogonal coordinate system. The Levy-Lees transformation generalized to three-dimensional flows is utilized. The inviscid properties are obtained from the Euler equations using the Beam and Warming implicit approximate factorization scheme. The resulting equations are discretized and approximated by a two-point fmitedifference numerical scheme. The code developed is validated and applied to the simulation of the flowfield over aerospace vehicle configurations. The results present good correlation with the available data.

Q-deformed fermionic Lipkin model at finite temperature

The interplay between temperature and q-deformation in the phase transition properties of many-body systems is studied in the particular framework of the collective q-deformed fermionic Lipkin model. It is shown that in phase transitions occuring in many-fermion systems described by su(2)q-like models are strongly influenced by the q-deformation.

PCM thermal insulation in buildings

This paper presents the results of a numerical and experimental study of phase change material (PCM) filled walls and roofs under real operational conditions to achieve passive thermal comfort. The numerical part of the study was based on a one-dimensional model for the phase change problem controlled by pure conduction. Real radiation data was used to determine the external face temperature. The numerical treatment was based upon using finite difference approximations and the ADI scheme. The results obtained were compared with field measurements. The experimental set-up consisted of a small room with movable roof and side wall. The roof was constructed in the traditional way but with the phase change material enclosed. Thermocouples were distributed across the cross section of the roof. Another roof, identical but without the PCM, was also used during comparative tests. The movable wall was also constructed as is done traditionally but with the PCM enclosed. Again, thermocouples were distributed across the wall thickness to enable measurement of the local temperatures. Another wall, identical but without the PCM, was also used during comparative tests. The PCM used in the numerical and experimental tests was composed of a mixture of two commercial grades of glycol in order to obtain the required fusion temperature range. Comparison between the simulation results and the experiments indicated good agreement. Field tests also indicated that the PCM used was adequate and that the concept was effective in maintaining the indoor temperature very close to the established comfort limits. Further economical analysis indicated that the concept could effectively help in reducing the electric energy consumption and improving the energy demand pattern. © 1997 by John Wiley & Sons, Ltd.

A two phase flow nonlinear parameter estimation for capillary tube-suction line heat exchanger

Here we present two-phase flow nonlinear parameter estimation for HFC's flow through capillary tube-suction line heat exchangers, commonly used as expansion devices in small refrigeration systems. The simplifying assumptions adopted are: steady state, pure refrigerant, one-dimensional flow, negligible axial heat conduction in the fluid, capillary tube and suction line walls. Additionally, it is considered that the refrigerant is free from oil and both phases are assumed to be at the same pressure, that is, surface tension effects are neglected. Metastable flow effects are also disregarded, and the vapor is assumed to be saturated at the local pressure. The so-called homogeneous model, involving three, first order, ordinary differential equations is applied to analyze the two-phase flow region. Comparison is done with experimental measurements of the mass flow rate and temperature distribution along capillary tubes working with refrigerant HFC-134a in different operating conditions.