Critical behavior of one-particle spectral weights in the transverse Ising model

We investigate the critical behavior of the spectral weight of a single quasiparticle, one of the key observables in experiment, for the particular case of the transverse Ising model. Series expansions are calculated for the linear chain and the square and simple cubic lattices. For the chain model, a conjectured exact result is discovered. For the square and simple cubic lattices, series analyses are used to estimate the critical exponents. The results agree with the general predictions of Sachdev [Quantum Phase Transitions (Cambridge University Press, Cambridge, England, 1999)].

Exploring synonymous codon usage preferences of disulfide-bonded and non-disulfide bonded cysteines in the E-coli genome

High-quality data about protein structures and their gene sequences are essential to the understanding of the relationship between protein folding and protein coding sequences. Firstly we constructed the EcoPDB database, which is a high-quality database of Escherichia coli genes and their corresponding PDB structures. Based on EcoPDB, we presented a novel approach based on information theory to investigate the correlation between cysteine synonymous codon usages and local amino acids flanking cysteines, the correlation between cysteine synonymous codon usages and synonymous codon usages of local amino acids flanking cysteines, as well as the correlation between cysteine synonymous codon usages and the disulfide bonding states of cysteines in the E. coli genome. The results indicate that the nearest neighboring residues and their synonymous codons of the C-terminus have the greatest influence on the usages of the synonymous codons of cysteines and the usage of the synonymous codons has a specific correlation with the disulfide bond formation of cysteines in proteins. The correlations may result from the regulation mechanism of protein structures at gene sequence level and reflect the biological function restriction that cysteines pair to form disulfide bonds. The results may also be helpful in identifying residues that are important for synonymous codon selection of cysteines to introduce disulfide bridges in protein engineering and molecular biology. The approach presented in this paper can also be utilized as a complementary computational method and be applicable to analyse the synonymous codon usages in other model organisms. (c) 2005 Elsevier Ltd. All rights reserved.

Influence of lattice interactions on the jahn-teller distortion of the [Cu(H2O)(6)](2+) ion: Dependence of the crystal structure of K-2[Cu(H2O)(6)](SO4)(2x)(SeO4)(2-2x) upon the sulfate/selenate ratio

The temperature dependence of the structure of the mixed-anion Tutton salt K-2[Cu(H2O)(6)](SO4)(2x)(SeO4)(2-2x) has been determined for crystals with 0, 17, 25, 68, 78, and 100% sulfate over the temperature range of 85-320 K. In every case, the [Cu(H2O)(6)](2+) ion adopts a tetragonally elongated coordination geometry with an orthorhombic distortion. However, for the compounds with 0, 17, and 25% sulfate, the long and intermediate bonds occur on a different pair of water molecules from those with 68, 78, and 100% sulfate. A thermal equilibrium between the two forms is observed for each crystal, with this developing more readily as the proportions of the two counterions become more similar. Attempts to prepare a crystal with approximately equal amounts of sulfate and selenate were unsuccessful. The temperature dependence of the bond lengths has been analyzed using a model in which the Jahn-Teller potential surface of the [Cu(H2O)(6)](2+) ion is perturbed by a lattice-strain interaction. The magnitude and sign of the orthorhombic component of this strain interaction depends on the proportion of sulfate to selenate. Significant deviations from Boltzmann statistics are observed for those crystals exhibiting a large temperature dependence of the average bond lengths, and this may be explained by cooperative interactions between neighboring complexes.

Understanding the crystallography of the eutectoid microstructure in a Zn-Al alloy using the edge-to-edge matching model

The orientation relationship (OR) between the beta(Zn) phase and the alpha(Al) phase and the corresponding habit planes in a Zn-Al eutectoid alloy were accurately determined using convergent beam Kikuchi line diffraction patterns. In addition to the previously reported OR. [11 (2) over bar0](beta)parallel to[110](alpha), (0002)(beta)parallel to ((1) over bar 11)alpha, two new ORs were observed. They are: [11 (2) over bar0](beta)parallel to [110], ((1) over bar 101)(beta) 0.82 degrees from (002)(alpha) and [(1) over bar 100](beta)parallel to[112](alpha), (0002)(beta) 4.5 degrees from (111)(alpha). These ORs can be explained and understood using the recently developed edge-to-edge matching model. (c) 2006 Acta Materialia Inc. Published by Elsevier Ltd. All rights reserved.

Error-correcting codes on a Bethe-like lattice

We analyse Gallager codes by employing a simple mean-field approximation that distorts the model geometry and preserves important interactions between sites. The method naturally recovers the probability propagation decoding algorithm as a minimization of a proper free-energy. We find a thermodynamical phase transition that coincides with information theoretical upper-bounds and explain the practical code performance in terms of the free-energy landscape.

Lattice approach to the dynamics of phase-coded soliton trains

We address the collective dynamics of a soliton train propagating in a medium described by the nonlinear Schrödinger equation. Our approach uses the reduction of train dynamics to the discrete complex Toda chain (CTC) model for the evolution of parameters for each train constituent: such a simplification allows one to carry out an approximate analysis of the dynamics of positions and phases of individual interacting pulses. Here, we employ the CTC model to the problem which has relevance to the field of fibre optics communications where each binary digit of transmitted information is encoded via the phase difference between the two adjacent solitons. Our goal is to elucidate different scenarios of the train distortions and the subsequent information garbling caused solely by the intersoliton interactions. First, we examine how the structure of a given phase pattern affects the initial stage of the train dynamics and explain the general mechanisms for the appearance of unstable collective soliton modes. Then we further discuss the nonlinear regime concentrating on the dependence of the Lax scattering matrix on the input phase distribution; this allows one to classify typical features of the train evolution and determine the distance where the soliton escapes from its slot. In both cases, we demonstrate deep mathematical analogies with the classical theory of crystal lattice dynamics.

Studies of the effects of molecular encounters on N.M.R. spin-lattice relaxation times

This thesis is concerned with investigations of the effects of molecular encounters on nuclear magnetic resonance spin-lattice relaxation times, with particular reference to mesitylene in mixtures with cyclohexane and TMS. The purpose of the work was to establish the best theoretical description of T1 and assess whether a recently identified mechanism (buffeting), that influences n.m.r. chemical shifts, governs Tl also. A set of experimental conditions are presented that allow reliable measurements of Tl and the N. O. E. for 1H and 13C using both C. W. and F.T. n.m.r. spectroscopy. Literature data for benzene, cyclohexane and chlorobenzene diluted by CC14 and CS2 are used to show that the Hill theory affords the best estimation of their correlation times but appears to be mass dependent. Evaluation of the T1 of the mesitylene protons indicates that a combined Hill-Bloembergen-Purcell-Pound model gives an accurate estimation of T1; subsequently this was shown to be due to cancellation of errors in the calculated intra and intemolecular components. Three experimental methods for the separation of the intra and intermolecular relaxation times are described. The relaxation times of the 13C proton satellite of neat bezene, 1,4 dioxane and mesitylene were measured. Theoretical analyses of the data allow the calculation of Tl intra. Studies of intermolecular NOE's were found to afford a general method of separating observed T1's into their intra and intermolecular components. The aryl 1H and corresponding 13C T1 values and the NOE for the ring carbon of mesitylene in CC14 and C6H12-TMS have been used in combination to determine T1intra and T1inter. The Hill and B.P.P. models are shown to predict similarly inaccurate values for T1linter. A buffeting contribution to T1inter is proposed which when applied to the BPP model and to the Gutowsky-Woessner expression for T1inter gives an inaccuracy of 12% and 6% respectively with respect to theexperimentally based T1inter.

Lattice approach to the dynamics of phase-coded soliton trains

We address the collective dynamics of a soliton train propagating in a medium described by the nonlinear Schrödinger equation. Our approach uses the reduction of train dynamics to the discrete complex Toda chain (CTC) model for the evolution of parameters for each train constituent: such a simplification allows one to carry out an approximate analysis of the dynamics of positions and phases of individual interacting pulses. Here, we employ the CTC model to the problem which has relevance to the field of fibre optics communications where each binary digit of transmitted information is encoded via the phase difference between the two adjacent solitons. Our goal is to elucidate different scenarios of the train distortions and the subsequent information garbling caused solely by the intersoliton interactions. First, we examine how the structure of a given phase pattern affects the initial stage of the train dynamics and explain the general mechanisms for the appearance of unstable collective soliton modes. Then we further discuss the nonlinear regime concentrating on the dependence of the Lax scattering matrix on the input phase distribution; this allows one to classify typical features of the train evolution and determine the distance where the soliton escapes from its slot. In both cases, we demonstrate deep mathematical analogies with the classical theory of crystal lattice dynamics.

Anisotropic model of kinetic roughening:he strong-coupling regime

We study the strong coupling (SC) limit of the anisotropic Kardar-Parisi-Zhang (KPZ) model. A systematic mapping of the continuum model to its lattice equivalent shows that in the SC limit, anisotropic perturbations destroy all spatial correlations but retain a temporal scaling which shows a remarkable crossover along one of the two spatial directions, the choice of direction depending on the relative strength of anisotropicity. The results agree with exact numerics and are expected to settle the long-standing SC problem of a KPZ model in the infinite range limit. © 2007 The American Physical Society.

Exotic meson decay widths using lattice quantum chromodynamics

Most experiments in particle physics are scattering experiments, the analysis of which leads to masses, scattering phases, decay widths and other properties of one or multi-particle systems. Until the advent of Lattice Quantum Chromodynamics (LQCD) it was difficult to compare experimental results on low energy hadron-hadron scattering processes to the predictions of QCD, the current theory of strong interactions. The reason being, at low energies the QCD coupling constant becomes large and the perturbation expansion for scattering; amplitudes does not converge. To overcome this, one puts the theory onto a lattice, imposes a momentum cutoff, and computes the integral numerically. For particle masses, predictions of LQCD agree with experiment, but the area of decay widths is largely unexplored. ^ LQCD provides ab initio access to unusual hadrons like exotic mesons that are predicted to contain real gluonic structure. To study decays of these type resonances the energy spectra of a two-particle decay state in a finite volume of dimension L can be related to the associated scattering phase shift δ(k) at momentum k through exact formulae derived by Lüscher. Because the spectra can be computed using numerical Monte Carlo techniques, the scattering phases can thus be determined using Lüscher's formulae, and the corresponding decay widths can be found by fitting Breit-Wigner functions. ^ Results of such a decay width calculation for an exotic hybrid( h) meson (JPC = 1-+) are presented for the decay channel h → πa 1. This calculation employed Lüscher's formulae and an approximation of LQCD called the quenched approximation. Energy spectra for the h and πa1 systems were extracted using eigenvalues of a correlation matrix, and the corresponding scattering phase shifts were determined for a discrete set of πa1 momenta. Although the number of phase shift data points was sparse, fits to a Breit-Wigner model were made, resulting in a decay width of about 60 MeV. ^

Lattice Boltzmann modeling of fluid flow and solute transport in karst aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Lattice Boltzmann Modeling of Fluid Flow and Solute Transport in Karst Aquifers

A novel modeling approach is applied to karst hydrology. Long-standing problems in karst hydrology and solute transport are addressed using Lattice Boltzmann methods (LBMs). These methods contrast with other modeling approaches that have been applied to karst hydrology. The motivation of this dissertation is to develop new computational models for solving ground water hydraulics and transport problems in karst aquifers, which are widespread around the globe. This research tests the viability of the LBM as a robust alternative numerical technique for solving large-scale hydrological problems. The LB models applied in this research are briefly reviewed and there is a discussion of implementation issues. The dissertation focuses on testing the LB models. The LBM is tested for two different types of inlet boundary conditions for solute transport in finite and effectively semi-infinite domains. The LBM solutions are verified against analytical solutions. Zero-diffusion transport and Taylor dispersion in slits are also simulated and compared against analytical solutions. These results demonstrate the LBM’s flexibility as a solute transport solver. The LBM is applied to simulate solute transport and fluid flow in porous media traversed by larger conduits. A LBM-based macroscopic flow solver (Darcy’s law-based) is linked with an anisotropic dispersion solver. Spatial breakthrough curves in one and two dimensions are fitted against the available analytical solutions. This provides a steady flow model with capabilities routinely found in ground water flow and transport models (e.g., the combination of MODFLOW and MT3D). However the new LBM-based model retains the ability to solve inertial flows that are characteristic of karst aquifer conduits. Transient flows in a confined aquifer are solved using two different LBM approaches. The analogy between Fick’s second law (diffusion equation) and the transient ground water flow equation is used to solve the transient head distribution. An altered-velocity flow solver with source/sink term is applied to simulate a drawdown curve. Hydraulic parameters like transmissivity and storage coefficient are linked with LB parameters. These capabilities complete the LBM’s effective treatment of the types of processes that are simulated by standard ground water models. The LB model is verified against field data for drawdown in a confined aquifer.

Numerical study of the enlarged O(5) symmetry of the 3D antiferromagnetic RP^(2) spin model

We investigate by means of Monte Carlo simulation and finite-size scaling analysis the critical properties of the three dimensional O (5) non-linear σ model and of the antiferromagnetic RP^(2) model, both of them regularized on a lattice. High accuracy estimates are obtained for the critical exponents, universal dimensionless quantities and critical couplings. It is concluded that both models belong to the same universality class, provided that rather non-standard identifications are made for the momentum-space propagator of the RP^(2) model. We have also investigated the phase diagram of the RP^(2) model extended by a second-neighbor interaction. A rich phase diagram is found, where most of the phase transitions are of the first order.

Phase diagram of the bosonic double-exchange model

The phase diagram of the simplest approximation to double-exchange systems, the bosonic double-exchange model with antiferromagnetic (AFM) superexchange coupling, is fully worked out by means of Monte Carlo simulations, large-N expansions, and variational mean-field calculations. We find a rich phase diagram, with no first-order phase transitions. The most surprising finding is the existence of a segmentlike ordered phase at low temperature for intermediate AFM coupling which cannot be detected in neutron-scattering experiments. This is signaled by a maximum (a cusp) in the specific heat. Below the phase transition, only short-range ordering would be found in neutron scattering. Researchers looking for a quantum critical point in manganites should be wary of this possibility. Finite-size scaling estimates of critical exponents are presented, although large scaling corrections are present in the reachable lattice sizes.

Monte Carlo determination of the phase diagram of the double-exchange model

We study the phase diagram of the double exchange model, with antiferromagnetic interactions, in a cubic lattice both at zero and finite temperature. There is a rich variety of magnetic phases, combined with regions where phase separation takes place. We identify phases, intrinsic to the cubic lattice, which are stable for realistic values of the interactions and dopings. Some of these phases break chiral symmetry, leading to unusual features.