First order in the plaquette expansion

The plaquette expansion, a general non-perturbative method for calculating the properties of lattice Hamiltonian systems, is analytically investigated at the first non-trivial order for an arbitrary system. At this level the approximation describes systems with either a bounded or an unbounded spectrum, depending on simple inequalities formed from the first four cumulants. The analysis yields analytic forms for the ground state energy, mass gaps and density of states in the thermodynamic limit. An exact form for the resolvent oeprator is given for a finite number of sites, as well as the asymptotic form in the thermodynamic limit.     

Curvature corrections to the capillary wave Hamiltonian

The fluctuations of interfaces between fluid phases are governed by an effective nonlocal interface Hamiltonian whose gradient expansion leads to the Helfrich Hamiltonian. Based on density functional theory we discuss the validity of this expansion and study the expansion coefficients multiplying the corresponding curvature terms.     

Stability and instability in a class of car following model on a closed loop

A velocity-matching car following model is modified to represent the motion of n vehicles travelling on a closed loop. Each vehicle is given a preferred velocity profile, which it attempts to achieve while also attempting to maintain a zero relative velocity between itself and the vehicle in front. The crucial distinctive of the looped model, as opposed to ‘non-looped’ models, is that the last vehicle in the stream is itself being followed by the lead (first) vehicle. The model gives rise to a system of n coupled time delay differential equations which are solved approximately (using a Taylor series expansion in time delay) and numerically using a fourth-order Runge-Kutta routine.The stability of the model is considered and an analytic form of the stable region in parameter space is found in the limit as n approaches infinity.     

Persistent currents in superconducting quantum interference devices

Starting from the reduced dynamical model of a two-junction quantum interference device, it shown that a quantum analog of the system can be exhibited. This quantum model extends the well-known properties of the device when its characteristic dimensions are of the order of mesoscopic length scales. By finding eigenvalues of the corresponding Hamiltonian operator, the persistent currents flowing in the ring have been obtained. The resulting quantum analog of the overdamped two-junction quantum interference device can be seen as a supercurrent qubit operating in the limit of negligible capacitance and finite inductance.     

Spin glasses on Bethe lattices for large coordination number

We study spin glasses on random lattices with finite connectivity. In the infinite connectivity limit they reduce to the Sherrington Kirkpatrick model. In this paper we investigate the expansion around the high connectivity limit. Within the replica symmetry breaking scheme at two steps, we compute the free energy at the first order in the expansion in inverse powers of the average connectivity (z), both for the fixed connectivity and for the fluctuating connectivity random lattices. It is well known that the coefficient of the 1/z correction for the free energy is divergent at low temperatures if computed in the one step approximation. We find that this annoying divergence becomes much smaller if computed in the framework of the more accurate two steps breaking. Comparing the temperature dependance of the coefficients of this divergence in the replica symmetric, one step and two steps replica symmetry breaking, we conclude that this divergence is an artefact due to the use of a finite number of steps of replica symmetry breaking. The 1/z expansion is well defined also in the zero temperature limit. Received 15 July 2002 Published online 31 December 2002     

The ground state and the thermodynamics of an extended BCS model of superconductivity

The thermodynamics of an extended BCS model of superconductivity is investigated. A physical system is described by a Hamiltonian containing the BCS interaction and an attractive four-fermion interaction. The four-fermion potential is caused by attractions between Cooper pairs mediated by the phonon field. The weakness of this potential allows the use of perturbation theory. The perturbation expansion was restricted to the first order because in the ground state the second order terms are not larger than 0.5 percent of first order correction for parameters used for calculations. The BCS Hamiltonian is an unperturbed one. The ground state and the thermal properties are examined. As a result the jump in the specific heat is higher than that in the BCS case. Moreover, the squared critical field is larger than the corresponding one in the BCS theory. Additionally, we show connections with the Bogolyubov's mean field approach used earlier in order to investigate general physical consequences of the model.     

The eigenspectrum of the relativistic two-center Hamiltonian in the diatomic basis of bound hydrogenlike Dirac-spinors

During heavy ion collisions below the Coulomb barrier adiabatic molecular orbitals are formed for the inner electrons. For systems with total changes exceending Z=Z₁ + Z₂ 60 relativistic effects are important and we, therefore, use the two-center Dirac Hamiltonian to describe the dynamics. We investigate the method to approximate these molecular orbitals by an expansion in a di-atomic basis of Dirac wave functions and determine the binding energies by a variation of the expansion coefficients. Deviations from adiabaticity are taken into account to first order in the relative velocityv/c 0.1) of the nuclei which is considerably smaller than the electron velocity.     

ANALYTIC STUDY OF su(3) LATTICE GAUGE THEORY

The variational-cumulant expansion method has been extended to the case of lattice SU(3) Wilson model.The plaquette energy as an order parameter has calculated to the 2nd order expansion.No 1st order phase transition in the d=4 case is found which is in agreement with the monte Carlo results,and the 1st order phase transition in the d=5 case is clearly seen.The method can be used in the study of problems is LGT with SU(3) gauge group.     

Renormalization group dependence of the QCD coupling

The general relation between the standard expansion coefficients and the beta function for the QCD coupling is exactly derived in a mathematically strict way. It is accordingly found that an infinite number of logarithmic terms are lost in the standard expansion with a finite order, and these lost terms can be given in a closed form. Numerical calculations, by a new matching-invariant coupling with the corresponding beta function to four-loop level, show that the new expansion converges much faster.     

The multilevel pairing Hamiltonian versus the degenerate case

We study the pairing Hamiltonian in a set of non-degenerate levels. First, we review in the path integral framework the spontaneous breaking of the U(1) symmetry occurring in such a system for the degenerate situation. Then the behaviors with the coupling constant of the ground state energy in the multilevel and in the degenerate case are compared. Next we discuss, in the multilevel case, an exact strong coupling expansion for the ground state energy which introduces the moments of the single particle level distribution. The domain of validity of the expansion, which is known in the macroscopic limit, is explored for finite systems and its implications for the energy of the latter is discussed. Finally the seniority and Gaudin excitations of the pairing Hamiltonian are addressed and shown to display the same gap in leading order.     

Dual‐fermion approach to non‐equilibrium strongly correlated problems

We present a generalization of the recently developed dual‐fermion approach introduced for correlated lattices to non‐equilibrium problems. In its local limit, the approach has been used to devise an efficient impurity solver, the superperturbation solver for the Anderson impurity model (AIM). Here we show that the general dual perturbation theory can be formulated on the Keldysh contour. Starting from a reference Hamiltonian system, in which the time‐dependent solution is found by exact diagonalization, we make a dual perturbation expansion in order to account for the relaxation effects from the fermionic bath. Simple test results for closed as well as open quantum systems in a fermionic bath are presented.     

Variational pair theory of the attractive and extended Hubbard model ind-dimensions

We present a variational approach for treating the Hubbard Hamiltonian in one, two and three dimensions. It is based on 2M-fermion wavefunctions which are allowed to form correlated spin-singlet pairs. Expressions for the ground state energy and correlation functions are derived in terms of general pair coefficient functions. The presented approach offers a convenient starting point for improved variational treatments that allow to include different specific types of pair correlations. We present first applications to the attractive and to the extended Hubbard model using a very simple ansatz for the pair coefficient functions. The ground state energy, chemical potential, order parameter, momentum distribution as well as spin-spin and density-density correlation functions follow from a system of coupled nonlinear equations that has to be solved selfconsistently. All quantities are given for arbitrary band-filling in one, two and three dimensions. Our results are compared with those of other approximations and for the one-dimensional case with the exact results of Krivnov and Ovchinnikov.     

Double-exchange model for molecule-based magnets

We report a detailed study of a double-exchange model proposed for the molecule-based magnets. The model is applied to a two-dimensional periodic complex made of a transition metal and an organic molecule in which the electronic structure is described by effective d orbitals of the transition metal ion at infinite Hund's coupling limit and the lowest unoccupied molecular orbital of the organic molecule, π^∗. Depending on the average electron density of the organic molecules and various superexchange couplings between metal ions' core spins, magnetic states of the complex are investigated. Performing Monte Carlo calculations on a model Hamiltonian for various electron densities of the organic molecule, the average magnetization and critical magnetic ordering temperatures are determined.     

Modeling of Lamb wave propagation in plate with two-dimensional phononic crystal layer coated on uniform substrate using plane-wave-expansion method

We show that the conversional three-dimensional plane wave expansion method can be revised to investigate the lamb wave propagation in the plate with two-dimensional phononic crystal layer coated on uniform substrate. We find that an imaginary three-dimensional periodic system can be constructed by stacking the studied plates and vacuum layers alternately, and then the Fourier series expansion can be performed. The difference between our imaginary periodic system and the true three-dimensional one is that, in our system, the Bloch feature of the wave along the thickness direction is broken. Three different systems are investigated by the proposed method as examples. The principle and reliability of the method are also discussed.     

Local analysis of frustration based on Kagomé lattices

Magnetic frustration in the framework of the Edwards-Anderson model is studied for the ground level (T=0) of the Kagomé lattice (KL). A sample consists of a realization of a random distribution of ferromagnetic (F) and antiferromagnetic (AF) interactions of the same strength along lines connecting nearest neighbors. Our goal is to compare two methods to calculate the following parameters: ground state energy, frustration and average frustration segment as functions of x, the concentration of F bonds. In doing so we make use of topological concepts such as plaquettes and frustration segments. The probability of a plaquette being unfrustrated (or flat) is ℘_p(x), while the probability of a plaquette being frustrated (or curved) is ℘_c(x). The analysis is done locally on a representative portion of the lattice which is called cell; cells of two different sizes are used in the present work. One method (which is simpler) is based on the probability function of any plaquette configuration ?(℘_p,℘_c). The other method is more exact but also more complex and increasingly difficult to use for large cells; it is based on the probability of bond configurations ψ(x). These methods are compared between themselves noting that the simpler method can be enough for most of the range for x. In addition, numerical simulations for many random samples at different concentrations x for a size given by 75 spins with periodic boundary conditions were done. This provides reference lines to compare with the properties under study. The local frustration analysis to obtain both ?(℘_p,℘_c) and ψ(x) is done over two cells of different size. Robustness of the criteria used in the local frustration analysis is also investigated.     

Algebraic methods applied to the study of energy transfer in anharmonic systems

We make use of two different methodologies to study the transition probabilities in a molecular anharmonic system in the presence of an external perturbation. For the first method, we use a series expansion of the displacement coordinate keeping up to fourth order terms; for the second method we use a deformed algebra to approximate the anharmonic Hamiltonian via a harmonic oscillator's Hamiltonian written in terms of deformed operators. We evaluate vibrational transition probabilities as a function of the collision energy and compare the results obtained with the two approaches.     

Optical conductivity by simple non-adiabatic model Hamiltonian

Two simple nonadiabatic model Hamiltonian, a generalized Rabi Hamiltonian and a caricature of the Fröhlich Hamiltonian are presented. The eigenvalue problem defined by the first is solved in Bargmann Hilbert space of analytical functions to demonstrate the methods. Results for the second are reported. Eigenvalues are given in terms of continued fractions, the eigenvectors as power series or Neumann series. The coefficients of the series are determined by simple three term recurrence relations. Using these results, the optical conductivity is calculated with Kubo's formulae. We correlate our results with experimental observations and the conventional small polaron theory.     

The neutrino response of low-density neutron matter from the virial expansion

We generalize our virial approach to study spin-polarized neutron matter and the consistent neutrino response at low densities. In the long-wavelength limit, the virial expansion makes model-independent predictions for the density and spin response, based only on nucleon–nucleon scattering data. Our results for the neutrino response provide constraints for random-phase approximation or other model calculations, and we compare the virial vector and axial response to response functions used in supernova simulations. The virial expansion is suitable to describe matter near the supernova neutrinosphere, and this work extends the virial equation of state to predict neutrino interactions in neutron matter.