Testing predictions of macroscopic binary diffusion coefficients using lattice models with site heterogeneity

Quantitatively predicting mass transport rates for chemical mixtures in porous materials is important in applications of materials such as adsorbents, membranes, and catalysts. Because directly assessing mixture transport experimentally is challenging, theoretical models that can predict mixture diffusion coefficients using Only single-component information would have many uses. One such model was proposed by Skoulidas, Sholl, and Krishna (Langmuir, 2003, 19, 7977), and applications of this model to a variety of chemical mixtures in nanoporous materials have yielded promising results. In this paper, the accuracy of this model for predicting mixture diffusion coefficients in materials that exhibit a heterogeneous distribution of local binding energies is examined. To examine this issue, single-component and binary mixture diffusion coefficients are computed using kinetic Monte Carlo for a two-dimensional lattice model over a wide range of lattice occupancies and compositions. The approach suggested by Skoulidas, Sholl, and Krishna is found to be accurate in situations where the spatial distribution of binding site energies is relatively homogeneous, but is considerably less accurate for strongly heterogeneous energy distributions.

Hierarchical models for 2D presence/absence data having ambiguous zeroes: With a biogeographical case study on dingo behaviour

This dissertation is primarily an applied statistical modelling investigation, motivated by a case study comprising real data and real questions. Theoretical questions on modelling and computation of normalization constants arose from pursuit of these data analytic questions. The essence of the thesis can be described as follows. Consider binary data observed on a two-dimensional lattice. A common problem with such data is the ambiguity of zeroes recorded. These may represent zero response given some threshold (presence) or that the threshold has not been triggered (absence). Suppose that the researcher wishes to estimate the effects of covariates on the binary responses, whilst taking into account underlying spatial variation, which is itself of some interest. This situation arises in many contexts and the dingo, cypress and toad case studies described in the motivation chapter are examples of this. Two main approaches to modelling and inference are investigated in this thesis. The first is frequentist and based on generalized linear models, with spatial variation modelled by using a block structure or by smoothing the residuals spatially. The EM algorithm can be used to obtain point estimates, coupled with bootstrapping or asymptotic MLE estimates for standard errors. The second approach is Bayesian and based on a three- or four-tier hierarchical model, comprising a logistic regression with covariates for the data layer, a binary Markov Random field (MRF) for the underlying spatial process, and suitable priors for parameters in these main models. The three-parameter autologistic model is a particular MRF of interest. Markov chain Monte Carlo (MCMC) methods comprising hybrid Metropolis/Gibbs samplers is suitable for computation in this situation. Model performance can be gauged by MCMC diagnostics. Model choice can be assessed by incorporating another tier in the modelling hierarchy. This requires evaluation of a normalization constant, a notoriously difficult problem. Difficulty with estimating the normalization constant for the MRF can be overcome by using a path integral approach, although this is a highly computationally intensive method. Different methods of estimating ratios of normalization constants (N Cs) are investigated, including importance sampling Monte Carlo (ISMC), dependent Monte Carlo based on MCMC simulations (MCMC), and reverse logistic regression (RLR). I develop an idea present though not fully developed in the literature, and propose the Integrated mean canonical statistic (IMCS) method for estimating log NC ratios for binary MRFs. The IMCS method falls within the framework of the newly identified path sampling methods of Gelman & Meng (1998) and outperforms ISMC, MCMC and RLR. It also does not rely on simplifying assumptions, such as ignoring spatio-temporal dependence in the process. A thorough investigation is made of the application of IMCS to the three-parameter Autologistic model. This work introduces background computations required for the full implementation of the four-tier model in Chapter 7. Two different extensions of the three-tier model to a four-tier version are investigated. The first extension incorporates temporal dependence in the underlying spatio-temporal process. The second extensions allows the successes and failures in the data layer to depend on time. The MCMC computational method is extended to incorporate the extra layer. A major contribution of the thesis is the development of a fully Bayesian approach to inference for these hierarchical models for the first time. Note: The author of this thesis has agreed to make it open access but invites people downloading the thesis to send her an email via the 'Contact Author' function.

Design of experiments for bivariate binary responses modelled by Copula functions

Optimal design for generalized linear models has primarily focused on univariate data. Often experiments are performed that have multiple dependent responses described by regression type models, and it is of interest and of value to design the experiment for all these responses. This requires a multivariate distribution underlying a pre-chosen model for the data. Here, we consider the design of experiments for bivariate binary data which are dependent. We explore Copula functions which provide a rich and flexible class of structures to derive joint distributions for bivariate binary data. We present methods for deriving optimal experimental designs for dependent bivariate binary data using Copulas, and demonstrate that, by including the dependence between responses in the design process, more efficient parameter estimates are obtained than by the usual practice of simply designing for a single variable only. Further, we investigate the robustness of designs with respect to initial parameter estimates and Copula function, and also show the performance of compound criteria within this bivariate binary setting.

Lattice-based threshold changeability for standard shamir sectret-sharing schemes

We consider the problem of increasing the threshold parameter of a secret-sharing scheme after the setup (share distribution) phase, without further communication between the dealer and the shareholders. Previous solutions to this problem require one to start off with a nonstandard scheme designed specifically for this purpose, or to have communication between shareholders. In contrast, we show how to increase the threshold parameter of the standard Shamir secret-sharing scheme without communication between the shareholders. Our technique can thus be applied to existing Shamir schemes even if they were set up without consideration to future threshold increases. Our method is a new positive cryptographic application for lattice reduction algorithms, inspired by recent work on lattice-based list decoding of Reed-Solomon codes with noise bounded in the Lee norm. We use fundamental results from the theory of lattices (geometry of numbers) to prove quantitative statements about the information-theoretic security of our construction. These lattice-based security proof techniques may be of independent interest.

Lattice-based threshold-changeability for standard CRT secret-sharing schemes

We consider the problem of increasing the threshold parameter of a secret-sharing scheme after the setup (share distribution) phase, without further communication between the dealer and the shareholders. Previous solutions to this problem require one to start off with a non-standard scheme designed specifically for this purpose, or to have secure channels between shareholders. In contrast, we show how to increase the threshold parameter of the standard CRT secret-sharing scheme without secure channels between the shareholders. Our method can thus be applied to existing CRT schemes even if they were set up without consideration to future threshold increases. Our method is a positive cryptographic application for lattice reduction algorithms, and we also use techniques from lattice theory (geometry of numbers) to prove statements about the correctness and information-theoretic security of our constructions.

Lattice-based threshold-changeability for standard Shamir secret-sharing schemes

We consider the problem of increasing the threshold parameter of a secret-sharing scheme after the setup (share distribution) phase, without further communication between the dealer and the shareholders. Previous solutions to this problem require one to start off with a non-standard scheme designed specifically for this purpose, or to have communication between shareholders. In contrast, we show how to increase the threshold parameter of the standard Shamir secret-sharing scheme without communication between the shareholders. Our technique can thus be applied to existing Shamir schemes even if they were set up without consideration to future threshold increases. Our method is a new positive cryptographic application for lattice reduction algorithms, inspired by recent work on lattice-based list decoding of Reed-Solomon codes with noise bounded in the Lee norm. We use fundamental results from the theory of lattices (Geometry of Numbers) to prove quantitative statements about the information-theoretic security of our construction. These lattice-based security proof techniques may be of independent interest.

Efficient parameter estimation in longitudinal data analysis using a hybrid GEE method

The method of generalized estimating equations (GEEs) provides consistent estimates of the regression parameters in a marginal regression model for longitudinal data, even when the working correlation model is misspecified (Liang and Zeger, 1986). However, the efficiency of a GEE estimate can be seriously affected by the choice of the working correlation model. This study addresses this problem by proposing a hybrid method that combines multiple GEEs based on different working correlation models, using the empirical likelihood method (Qin and Lawless, 1994). Analyses show that this hybrid method is more efficient than a GEE using a misspecified working correlation model. Furthermore, if one of the working correlation structures correctly models the within-subject correlations, then this hybrid method provides the most efficient parameter estimates. In simulations, the hybrid method's finite-sample performance is superior to a GEE under any of the commonly used working correlation models and is almost fully efficient in all scenarios studied. The hybrid method is illustrated using data from a longitudinal study of the respiratory infection rates in 275 Indonesian children.

Mean-field results of a lattice-gas model of multilayer adsorption

A lattice-gas model of multilayer adsorption has been solved in the mean-field approximation by a different numerical method. Earlier workers obtained a single solution for all values of temperature and pressure. In the present work, multiple solutions have been obtained in certain regions of temperature and pressure which give rise to bysteresis in the adsorption isotherm. In addition, we have obtained a parameter which behaves like an order parameter for the transition. The potential-energy function shows a double minimum in the region of bysteresis and a single maximum elsewhere.

Critical point phenomena in the electrical resistivity of binary liquid systems: Cs2 + CH3CN and CS2 + CH3NO2

Electrical resistance measurements are reported on the binary liquid mixtures CS2 + CH3CN and CS2 + CH3NO2 with special reference to the critical region. Impurity conduction seems to be the dominant mechanism for charge transport. For the liquid mixture filled at the critical composition, the resistance of the system aboveT c follows the relationR=R c−A(T−T c) b withb=0·6±0·1. BelowT c the conductivities of the two phases obey a relation σ2−σ1=B(T c−T)β with β=0·34±0·02, the exponent of the transport coefficient being the same as the exponent of the order parameter, an equilibrium property.

Phase diagram and dielectric studies of binary organic materials

The phase equilibrium studies of organic system, involving resorcinol (R) and p-dimethylaminobenzaldehyde (DMAB), reveal the formation of a 1:1 molecular complex with two eutectics. The heat of mixing, entropy of fusion, roughness parameter, interfacial energy, and the excess thermodynamic functions were calculated based on enthalpy of fusion data determined via differential scanning calorimetric (DSC) method. X-ray powder diffraction studies confirm that the eutectics are not simple mechanical mixture of the components under investigation. The spectroscopic investigations (IR and NMR) suggest the occurrence of hydrogen bonding between the components forming the molecular complex. The dielectric measurements, carried out on hot-pressed addition compound (molecular complex), show higher dielectric constant at 320 K than that of individual components. The microstructural investigations of eutectic and addition compound indicate dendritic and faceted morphological features. (C) 2000 Elsevier Science B.V. All rights reserved.

A High Pressure Study of the Phase Diagram and the Critical Exponent β of Binary Liquid Systems

The pressure dependence of critical parameters xc, Tc, and β have been analysed in four systems namely cyclohexane + acetic anhydride, n-heptane + acetic anhydride, methanol + n-heptane, and carbon disulphide + acetonitrile. The separation temperature was found to increase linearly with pressure the value of dTc/dP being 28 mK, 11 mK, 22 mK, and 25 mK respectively. These are in fair agreement with earlier measurements available for two systems. For the methanol + n-heptane system dTc/dP is apparently not consistent with the value predicted from the specific heat and thermal expansion data.Die Druckabhängigkeit der kritischen Parameter xc, Tc und β ist in den vier Systemen Cyclohexan + Essigsäureanhydrid, n-Heptan + Essigsäureanhydrid, Methanol + n-Heptan und Schwefelkohlenstoff + Acetonitril analysiert worden. Es wurde gefunden, daß die kritische Temperatur linear mit dem Druck ansteigt. Die Werte für dTc/dP betragen 28 mK, 11 mK, 22 mK und 25 mK. Sie sind in guter überein-stimmung mit früheren Messungen an zweien dieser Systeme. Für Methanol + n-Heptan stimmt der Wert für dT/dP offensichtlich nicht mit Werten, die mit Hilfe von Daten für die spezifische Wärme und die thermische Ausdehnung vorhergesagt wurden, überein.

Kaon B parameter with Wilson fermions

We calculate the kaon B parameter in quenched lattice QCD at beta=6.0 using Wilson fermions at kappa=0.154 and 0.155. We use two kinds of nonlocal (''smeared'') sources for quark propagators to calculate the matrix elements between states of definite momentum. The use of smeared sources yields results with much smaller errors than obtained in previous calculations with Wilson fermions. By combining results for p=(0,0,0) and p=(0,0,1), we show that one can carry out the noperturbative subtraction necessary to remove the dominant lattice artifacts induced by the chiral-symmetry-breaking term in the Wilson action. Our final results are in good agreement with those obtained using staggered fermions. We also present results for B parameters of the DELTAI = 3/2 part of the electromagnetic penguin operators, and preliminary results for B(K) in the presence of two flavors of dynamical quarks.

Theoretical and Computational Analysis of Static and Dynamic Anomalies in Water-DMSO Binary Mixture at Low DMSO Concentrations

Experiments have repeatedly observed both thermodynamic and dynamic anomalies in aqueous binary mixtures, surprisingly at low solute concentration. Examples of such binary mixtures include water-DMSO, water-ethanol, water-tertiary butyl alcohol (TBA), and water-dioxane, to name a few. The anomalies have often been attributed to the onset of a structural transition, whose nature, however, has been left rather unclear. Here we study the origin of such anomalies using large scale computer simulations and theoretical analysis in water-DMSO binary mixture. At very low DMSO concentration (below 10%), small aggregates of DMSO are solvated by water through the formation of DMSO-(H2O)(2) moieties. As the concentration is increased beyond 10-12% of DMSO, spanning clusters comprising the same moieties appear in the system. Those clusters are formed and stabilized not only through H-bonding but also through the association of CH3 groups of DMSO. We attribute the experimentally observed anomalies to a continuum percolation-like transition at DMSO concentration X-DMSO approximate to 12-15%. The largest cluster size of CH3-CH3 aggregation clearly indicates the formation of such percolating clusters. As a result, a significant slowing down is observed in the decay of associated rotational auto time correlation functions (of the S = O bond vector of DMSO and O-H bond vector of water). Markedly unusual behavior in the mean square fluctuation of total dipole moment again suggests a structural transition around the same concentration range. Furthermore, we map our findings to an interacting lattice model which substantiates the continuum percolation model as the reason for low concentration anomalies in binary mixtures where the solutes involved have both hydrophilic and hydrophobic moieties.

Exact Solution of a One-Dimensional Multicomponent Lattice Gas with Hyperbolic Interaction

We present the exact solution to a one-dimensional multicomponent quantum lattice model interacting by an exchange operator which falls off as the inverse sinh square of the distance. This interaction contains a variable range as a parameter and can thus interpolate between the known solutions for the nearest-neighbor chain and the inverse-square chain. The energy, susceptibility, charge stiffness, and the dispersion relations for low-lying excitations are explicitly calculated for the absolute ground state, as a function of both the range of the interaction and the number of species of fermions.

Neutron diffraction study of the coupling between spin, lattice, and structural degrees of freedom in 0.8BiFeO(3)-0.2PbTiO(3)

A powder neutron diffraction study was carried out on 0.8BiFeO(3)-0.2PbTiO(3) in the temperature range 27-1000 degrees C. The system exhibits magnetic transition at similar to 300 degrees C and a rhombohedral (R3c)-cubic (Pm3m) ferroelectric phase transition at similar to 650 degrees C. Anomalous variation in the lattice parameters and the octahedral tilt angle is observed across the magnetic transition temperature. In the magnetic phase, the c parameter is contracted and the octahedral tilt angle is slightly increased. The results suggest coupling between the spin, lattice and structural degrees of freedom. (C) 2011 American Institute of Physics. doi:10.1063/1.3555093]