Free energy approximations in simple lattice proteins

This work addresses the question of whether it is possible to define simple pairwise interaction terms to approximate free energies of proteins or polymers. Rather than ask how reliable a potential of mean force is, one can ask how reliable it could possibly be. In a two-dimensional, infinite lattice model system one can calculate exact free energies by exhaustive enumeration. A series of approximations were fitted to exact results to assess the feasibility and utility of pairwise free energy terms. Approximating the true free energy with pairwise interactions gives a poor fit with little transferability between systems of different size. Adding extra artificial terms to the approximation yields better fits, but does not improve the ability to generalize from one system size to another. Furthermore, one cannot distinguish folding from nonfolding sequences via the approximated free energies. Most usefully, the methodology shows how one can assess the utility of various terms in lattice protein/polymer models. (C) 2001 American Institute of Physics.

Determination of the stability constants of mercury(II) complexes of mixed donor macrobicyclic encapsulating ligands

The reactions of mercury(II) with the mixed donor encapsulating ligands 3,6,16-trithia-6,11,19-triazabicyclo[6.6.6]icosane (AMN(3)S(3)sar) and 1-amino-8-methyl-6,19-dithia-3,10,13,16-tetraazabicyclo[6.6.6]icosane (AMN(4)S(2)sar) have been studied. NMR ligand-ligand competition experiments with the ligands 1,4,8,11-tetraazaeyclotetradecane ([14]aneN(4)), 1-thia-4,7,10-triazacyclododecane ([12]aneN(3)S) and ethylenediaminetetraacetic acid (EDTA) with AMN(3)S(3)sar and Hg(II) indicated that [14]aneN(4) would be an appropriate competing ligand for the, determination of the Hg(II) stability constant. Calculations indicated the ratio of concentrations of AMN3S3sar, [14]aneN(4) and Hg(II) required for the determination of the stability constant ranged from 1:1:1 to 1:5:1. Refinement of the titration curves yielded log(10)K[Hg(AMN(3)S(3)sar)](2+) = 17.7. A similar competition titration resulted in the determination of the stability constant for the AMN(4)S(2)sar system as log(10)K[Hg(AMN(4)S(2)sar)](2+) = 19.5. The observed binding constants for the mixed N/S donor systems and the hexaaza analogues sar (3,6,10,13,16,19-hexaazabicyclo [6.6.6]icosane) and diamsar (1,8-diamino-3,6,10,13,16,19 -hexazabicyclo [6.6.6] icosane (log(10)K-[Hg(diamsar)](2+) = 26.4; log(10)K[Hg(sar)](2+) = 28.1) differ by approximately ten orders of magnitude. The difference is ascribed not to a cryptate effect but to a mismatch in the Hg-N and Hg-S bond lengths in the N/S systems.

Ladder operators for integrable one-dimensional lattice models

A generalised ladder operator is used to construct the conserved operators for any one-dimensional lattice model derived from the Yang-Baxter equation. As an example, the low order conserved operators for the XYh model are calculated explicitly.

Signature of mott-insulator transition with ultracold fermions in a one-dimensional optical lattice

Using the exact Bethe ansatz solution of the Hubbard model and Luttinger liquid theory, we investigate the density profiles and collective modes of one-dimensional ultracold fermions confined in an optical lattice with a harmonic trapping potential. We determine a generic phase diagram in terms of a characteristic filling factor and a dimensionless coupling constant. The collective oscillations of the atomic mass density, a technique that is commonly used in experiments, provide a signature of the quantum phase transition from the metallic phase to the Mott-insulator phase. A detailed experimental implementation is proposed.

Ferromagnetism, paramagnetism, and a Curie-Weiss metal in an electron-doped Hubbard model on a triangular lattice

Motivated by the unconventional properties and rich phase diagram of NaxCoO2 we consider the electronic and magnetic properties of a two-dimensional Hubbard model on an isotropic triangular lattice doped with electrons away from half-filling. Dynamical mean-field theory (DMFT) calculations predict that for negative intersite hopping amplitudes (t < 0) and an on-site Coulomb repulsion, U, comparable to the bandwidth, the system displays properties typical of a weakly correlated metal. In contrast, for t > 0 a large enhancement of the effective mass, itinerant ferromagnetism, and a metallic phase with a Curie-Weiss magnetic susceptibility are found in a broad electron doping range. The different behavior encountered is a consequence of the larger noninteracting density of states (DOS) at the Fermi level for t > 0 than for t < 0, which effectively enhances the mass and the scattering amplitude of the quasiparticles. The shape of the DOS is crucial for the occurrence of ferromagnetism as for t > 0 the energy cost of polarizing the system is much smaller than for t < 0. Our observation of Nagaoka ferromagnetism is consistent with the A-type antiferromagnetism (i.e., ferromagnetic layers stacked antiferromagnetically) observed in neutron scattering experiments on NaxCoO2. The transport and magnetic properties measured in NaxCoO2 are consistent with DMFT predictions of a metal close to the Mott insulator and we discuss the role of Na ordering in driving the system towards the Mott transition. We propose that the Curie-Weiss metal phase observed in NaxCoO2 is a consequence of the crossover from a bad metal with incoherent quasiparticles at temperatures T > T-* and Fermi liquid behavior with enhanced parameters below T-*, where T-* is a low energy coherence scale induced by strong local Coulomb electron correlations. Our analysis also shows that the one band Hubbard model on a triangular lattice is not enough to describe the unusual properties of NaxCoO2 and is used to identify the simplest relevant model that captures the essential physics in NaxCoO2. We propose a model which allows for the Na ordering phenomena observed in the system which, we propose, drives the system close to the Mott insulating phase even at large dopings.

Modeling the optical constants of solids using acceptance-probability-controlled simulated annealing with an adaptive move generation procedure

The acceptance-probability-controlled simulated annealing with an adaptive move generation procedure, an optimization technique derived from the simulated annealing algorithm, is presented. The adaptive move generation procedure was compared against the random move generation procedure on seven multiminima test functions, as well as on the synthetic data, resembling the optical constants of a metal. In all cases the algorithm proved to have faster convergence and superior escaping from local minima. This algorithm was then applied to fit the model dielectric function to data for platinum and aluminum.

Simulated-annealing-based genetic algorithm for modeling the optical constants of solids

We propose a simulated-annealing-based genetic algorithm for solving model parameter estimation problems. The algorithm incorporates advantages of both genetic algorithms and simulated annealing. Tests on computer-generated synthetic data that closely resemble optical constants of a metal were performed to compare the efficiency of plain genetic algorithms against the simulated-annealing-based genetic algorithms. These tests assess the ability of the algorithms to and the global minimum and the accuracy of values obtained for model parameters. Finally, the algorithm with the best performance is used to fit the model dielectric function to data for platinum and aluminum. (C) 1997 Optical Society of America.

Lattice-dome design using a knowledge-based system approach

In the design of lattice domes, design engineers need expertise in areas such as configuration processing, nonlinear analysis, and optimization. These are extensive numerical, iterative, and lime-consuming processes that are prone to error without an integrated design tool. This article presents the application of a knowledge-based system in solving lattice-dome design problems. An operational prototype knowledge-based system, LADOME, has been developed by employing the combined knowledge representation approach, which uses rules, procedural methods, and an object-oriented blackboard concept. The system's objective is to assist engineers in lattice-dome design by integrating all design tasks into a single computer-aided environment with implementation of the knowledge-based system approach. For system verification, results from design examples are presented.

Large-N solutions of the Heisenberg and Hubbard-Heisenberg models on the anisotropic triangular lattice: application to Cs2CuCl4 and to the layered organic superconductors kappa-(BEDT-TTF)(2)X (BEDT-TTF bis(ethylene-dithio)tetrathiafulvalene); X anion)

We solve the Sp(N) Heisenberg and SU(N) Hubbard-Heisenberg models on the anisotropic triangular lattice in the large-N limit. These two models may describe respectively the magnetic and electronic properties of the family of layered organic materials K-(BEDT-TTF)(2)X, The Heisenberg model is also relevant to the frustrated antiferromagnet, Cs2CuCl4. We find rich phase diagrams for each model. The Sp(N) :antiferromagnet is shown to have five different phases as a function of the size of the spin and the degree of anisotropy of the triangular lattice. The effects of fluctuations at finite N are also discussed. For parameters relevant to Cs2CuCl4 the ground state either exhibits incommensurate spin order, or is in a quantum disordered phase with deconfined spin-1/2 excitations and topological order. The SU(N) Hubbard-Heisenberg model exhibits an insulating dimer phase, an insulating box phase, a semi-metallic staggered flux phase (SFP), and a metallic uniform phase. The uniform and SFP phases exhibit a pseudogap, A metal-insulator transition occurs at intermediate values of the interaction strength.

Brief note on the variation of constants formula for fuzzy differential equations

This note gives a theory of state transition matrices for linear systems of fuzzy differential equations. This is used to give a fuzzy version of the classical variation of constants formula. A simple example of a time-independent control system is used to illustrate the methods. While similar problems to the crisp case arise for time-dependent systems, in time-independent cases the calculations are elementary solutions of eigenvalue-eigenvector problems. In particular, for nonnegative or nonpositive matrices, the problems at each level set, can easily be solved in MATLAB to give the level sets of the fuzzy solution. (C) 2002 Elsevier Science B.V. All rights reserved.

Gas permeation in silicon-oxide/polymer (SiOx/PET) barrier films: role of the oxide lattice, nano-defects and macro-defects

We propose a model for permeation in oxide coated gas barrier films. The model accounts for diffusion through the amorphous oxide lattice, nano-defects within the lattice, and macro-defects. The presence of nano-defects indicate the oxide layer is more similar to a nano-porous solid (such as zeolite) than silica glass with respect to permeation properties. This explains why the permeability of oxide coated polymers is much greater, and the activation energy of permeation much lower, than values expected for polymers coated with glass. We have used the model to interpret permeability and activation energies measured for the inert gases (He, Ne and Ar) in evaporated SiOx films of varying thickness (13-70 nm) coated on a polymer substrate. Atomic force and scanning electron microscopy were used to study the structure of the oxide layer. Although no defects could be detected by microscopy, the permeation data indicate that macro-defects (>1 nm), nano-defects (0.3-0.4 nm) and the lattice interstices (<0.3 nm) all contribute to the total permeation. (C) 2002 Elsevier Science B.V. All rights reserved.

Quasimode theory of quantum optical processes in photonic band gap materials

This paper deals with atomic systems coupled to a structured reservoir of quantum EM field modes, with particular relevance to atoms interacting with the field in photonic band gap materials. The case of high Q cavities has been treated elsewhere using Fano diagonalization based on a quasimode approach, showing that the cavity quasimodes are responsible for pseudomodes introduced to treat non-Markovian behaviour. The paper considers a simple model of a photonic band gap case, where the spatially dependent permittivity consists of a constant term plus a small spatially periodic term that leads to a narrow band gap in the spectrum of mode frequencies. Most treatments of photonic band gap materials are based on the true modes, obtained numerically by solving the Helmholtz equation for the actual spatially periodic permittivity. Here the field modes are first treated in terms of a simpler quasimode approach, in which the quasimodes are plane waves associated with the constant permittivity term. Couplings between the quasimodes occur owing to the small periodic term in the permittivity, with selection rules for the coupled modes being related to the reciprocal lattice vectors. This produces a field Hamiltonian in quasimode form. A matrix diagonalization method may be applied to relate true mode annihilation operators to those for quasimodes. The atomic transitions are coupled to all the quasimodes, and the true mode atom-EM field coupling constants (one-photon Rabi frequencies) are related to those for the quasimodes and also expressions are obtained for the true mode density. The results for the one-photon Rabi frequencies differ from those assumed in other work. Expressions for atomic decay rates are obtained using the Fermi Golden rule, although these are valid only well away from the band gaps.

Dynamic Behavior of the Jahn-Teller distorted Cu(H2O)(6)(2+) ion in Cu2+ doped Cs-2[Zn(H2O)(6)](ZrF6)(2) and the crystal structure of the host lattice

The temperature dependence of the X- and Q-band EPR spectra of Cs-2[Zn(H2O)(6)](ZrF6)(2) containing similar to1% Cu2+ is reported. All three molecular g-values vary with temperature, and their behavior is interpreted using a model in which the potential surface of the Jahn-Teller distorted Cu(H2O)(6)(2+) ion is perturbed by an orthorhombic strain induced by interactions with the surrounding lattice. The strain parameters are significantly smaller than those reported previously for the Cu(H2O)(6)(2+) ion in similar lattices. The temperature dependence of the two higher g-values suggests that in the present compound the lattice interactions change slightly with temperature. The crystal structure of the Cs-2[Zn(H2O)(6)](ZrF6)(2) host is reported, and the geometry of the Zn(H2O)(6)(2+) ion is correlated with lattice strain parameters derived from the EPR spectrum of the guest Cu2+ complex.