Dynamic weighting in Monte Carlo and optimization

Dynamic importance weighting is proposed as a Monte Carlo method that has the capability to sample relevant parts of the configuration space even in the presence of many steep energy minima. The method relies on an additional dynamic variable (the importance weight) to help the system overcome steep barriers. A non-Metropolis theory is developed for the construction of such weighted samplers. Algorithms based on this method are designed for simulation and global optimization tasks arising from multimodal sampling, neural network training, and the traveling salesman problem. Numerical tests on these problems confirm the effectiveness of the method.

Cloning and characterization of KAP3: a novel kinesin superfamily-associated protein of KIF3A/3B.

We previously reported that KIF3A and KIF3B form a heterodimer that functions as a microtubule-based fast anterograde translocator of membranous organelles. We have also shown that this KIF3A/3B forms a complex with other associated polypeptides, named kinesin superfamily-associated protein 3 (KAP3). In the present study, we purified KAP3 protein by immunoprecipitation using anti-KIF3B antibody from mouse testis. Microsequencing was carried out, and we cloned the full-length KAP3 cDNA from a mouse brain cDNA library. Two isoforms of KAP3 exist [KAP3A (793 aa) and KAP3B (772 aa)], generated by alternative splicing in the carboxyl terminus region. Their amino acid sequences have no homology with those of any other known proteins, and prediction of their secondary structure indicated that almost the entire KAP3 molecule is alpha-helical. We produced recombinant KAP3 and KIF3A/3B using a baculovirus-Sf9 expression system. A reconstruction study in Sf9 cells revealed that KAP3 is a globular protein that binds to the tail domain of KIF3A/3B. The immunolocalization pattern of KAP3 was similar to that of KIF3A/3B in nerve cells. In addition, we found that KAP3 does not affect the motor activity of KIF3A/3B. KAP3 was associated with a membrane-bound form of KIF3A/3B in a fractional immunoprecipitation experiment, and since the KIF3 complex was found to bind to membranous organelles in an EM study, KAP3 may regulate membrane binding of the KIF3 complex.

GE models and algorithms for condensed phase equilibrium data regression in ternary systems: limitations and proposals

Phase equilibrium data regression is an unavoidable task necessary to obtain the appropriate values for any model to be used in separation equipment design for chemical process simulation and optimization. The accuracy of this process depends on different factors such as the experimental data quality, the selected model and the calculation algorithm. The present paper summarizes the results and conclusions achieved in our research on the capabilities and limitations of the existing GE models and about strategies that can be included in the correlation algorithms to improve the convergence and avoid inconsistencies. The NRTL model has been selected as a representative local composition model. New capabilities of this model, but also several relevant limitations, have been identified and some examples of the application of a modified NRTL equation have been discussed. Furthermore, a regression algorithm has been developed that allows for the advisable simultaneous regression of all the condensed phase equilibrium regions that are present in ternary systems at constant T and P. It includes specific strategies designed to avoid some of the pitfalls frequently found in commercial regression tools for phase equilibrium calculations. Most of the proposed strategies are based on the geometrical interpretation of the lowest common tangent plane equilibrium criterion, which allows an unambiguous comprehension of the behavior of the mixtures. The paper aims to show all the work as a whole in order to reveal the necessary efforts that must be devoted to overcome the difficulties that still exist in the phase equilibrium data regression problem.

Applying reactive models to column experiments to assess the hydrogeochemistry of seawater intrusion: optimising ACUAINTRUSION and selecting cation exchange coefficients with PHREEQC

Three sets of laboratory column experimental results concerning the hydrogeochemistry of seawater intrusion have been modelled using two codes: ACUAINTRUSION (Chemical Engineering Department, University of Alicante) and PHREEQC (U.S.G.S.). These reactive models utilise the hydrodynamic parameters determined using the ACUAINTRUSION TRANSPORT software and fit the chloride breakthrough curves perfectly. The ACUAINTRUSION code was improved, and the instabilities were studied relative to the discretisation. The relative square errors were obtained using different combinations of the spatial and temporal steps: the global error for the total experimental data and the partial error for each element. Good simulations for the three experiments were obtained using the ACUAINTRUSION software with slight variations in the selectivity coefficients for both sediments determined in batch experiments with fresh water. The cation exchange parameters included in ACUAINTRUSION are those reported by the Gapon convention with modified exponents for the Ca/Mg exchange. PHREEQC simulations performed using the Gains-Thomas convention were unsatisfactory, with the exchange coefficients from the database of PHREEQC (or range), but those determined with fresh water – natural sediment allowed only an approximation to be obtained. For the treated sediment, the adjusted exchange coefficients were determined to improve the simulation and are vastly different from those from the database of PHREEQC or batch experiment values; however, these values fall in an order similar to the others determined under dynamic conditions. Different cation concentrations were simulated using two different software packages; this disparity could be attributed to the defined selectivity coefficients that affect the gypsum equilibrium. Consequently, different calculated sulphate concentrations are obtained using each type of software; a smaller mismatch was predicted using ACUAINTRUSION. In general, the presented simulations by ACUAINTRUSION and PHREEQC produced similar results, making predictions consistent with the experimental data. However, the simulated results are not identical to the experimental data; sulphate (total S) is overpredicted by both models, most likely due to such factors as the kinetics of gypsum, the possible variations in the exchange coefficients due to salinity and the neglect of other processes.

Building bridges… and properties aplenty: cultural problems in Spanish real estate marketing for prospective British buyers

N.B. reproduced with permission of Peter Lang Verlag. For citation, please, use the original reference, that is Campos Pardillos, M.A. and Balteiro Fernández, I. 2009. “Building bridges… and properties aplenty: cultural problems in Spanish real estate marketing for prospective British buyers”. In: Guillén-Nieto, V., C. Marimón-Llorca and C. Vargas-Sierra. Eds. Intercultural Business Communication and Simulation and Gaming Methodology. Bern: Peter Lang. Pp. 155-174.

Oxygen isotope and paleomagnetic stratigraphy of sediment core V28-239

Core Vema 28-238 preserves an excellent oxygen isotope and magnetic stratigraphy and is shown to contain undisturbed sediments deposited continuously through the past 870,000 yr. Detailed correlation with sequences described by Emiliani in the Caribbean and Atlantic Ocean is demonstrated. The boundaries of 22 stages representing alternating times of high and low Northern Hemisphere ice volume are recognized and dated. The record is interpreted in terms of Northern Hemisphere ice accumulation, and is used to estimate the range of temperature variation in the Caribbean.

Development of an analytical approach to highway barrier design and evaluation.

"The Bureau of Physical Research has prepared this report primarily from two reports submitted by Cornell Aeronautical Laboratory... 'Summary report of highway barrier analysis and test program' July 1961, and 'Analysis and prediction of the performance of highway barrier' October 1962."

In honor of the ninetieth birthday of Charles Frederick Johnson, professor of English in Trinity college, 1883-1906; papers, essays, and stories by his former students,

CONTENTS. (cont.) --The Christianity of the fourth Gospel, by the Rev. C. B. Hedrick.--Correctness, by the Rev. M. B. Stewart.--An optimistic realist, Ernst Zahn, by B. Q. Morgan.--The hermit story, as used by Voltaire and Mark Twain, by F. A. G. Cowper.--Thumbnail classics. Jack Harker at Yalevard, by P. Curtiss.--The man who heard everything, by W. Trumbull.--"Whiskey, bullets, and..."--The novelist takes to real estate, by R. Wright.

Time-series Luminance Distribution Maps: implementation of annual daylight simulation methods for occupant visual comfort analysis

Thesis (Master's)--University of Washington, 2016-06

Influence apart from Adoption: How Interaction between Programming and Scientific Practices Shapes Modes of Inquiry in Four Oceanography Teams

Thesis (Ph.D.)--University of Washington, 2016-06

Refined method of potential enhancement in the equilibria characterization of activated carbon. Comparison with GCMC and DFT

In this paper, we present a technique for equilibria characterization of activated carbon having slit-shaped pores. This method was first developed by Do (Do, D. D. A new method for the characterisation of micro-mesoporous materials. Presented at the International Symposium on New Trends in Colloid and Interface Science, September 24-26, 1998 Chiba, Japan) and applied by his group and other groups for characterization of pore size distribution (PSD) as well as adsorption equilibria determination of a wide range of hydrocarbons. It is refined in this paper and compared with the grand canonical Monte Carlo (GCMG) simulation and density functional theory (DFT). The refined theory results in a good agreement between the pore filling pressure versus pore width and those obtained by GCMG and DFT. Furthermore, our local isotherms are qualitatively in good agreement with those obtained by the GCMC simulations. The main advantage of this method is that it is about 4 orders of magnitude faster than the GCMC simulations, making it suitable for optimization studies and design purposes. Finally, we apply our method and the GCMG in the derivation of the PSD of a commercial activated carbon. It was found that the PSD derived from our method is comparable to that derived from the GCMG simulations.

Gene duplication and hierarchical modularity in intracellular interaction networks

Networks of interactions evolve in many different domains. They tend to have topological characteristics in common, possibly due to common factors in the way the networks grow and develop. It has been recently suggested that one such common characteristic is the presence of a hierarchically modular organization. In this paper, we describe a new algorithm for the detection and quantification of hierarchical modularity, and demonstrate that the yeast protein-protein interaction network does have a hierarchically modular organization. We further show that such organization is evident in artificial networks produced by computational evolution using a gene duplication operator, but not in those developing via preferential attachment of new nodes to highly connected existing nodes. (C) 2004 Elsevier Ireland Ltd. All rights reserved.

Using the canonical modelling approach to simplify the simulation of function in functional-structural plant models

Functional-structural plant models that include detailed mechanistic representation of underlying physiological processes can be expensive to construct and the resulting models can also be extremely complicated. On the other hand, purely empirical models are not able to simulate plant adaptability and response to different conditions. In this paper, we present an intermediate approach to modelling plant function that can simulate plant response without requiring detailed knowledge of underlying physiology. Plant function is modelled using a 'canonical' modelling approach, which uses compartment models with flux functions of a standard mathematical form, while plant structure is modelled using L-systems. Two modelling examples are used to demonstrate that canonical modelling can be used in conjunction with L-systems to create functional-structural plant models where function is represented either in an accurate and descriptive way, or in a more mechanistic and explanatory way. We conclude that canonical modelling provides a useful, flexible and relatively simple approach to modelling plant function at an intermediate level of abstraction.

Predictive vaccinology: Optimisation of predictions using support vector machine classifiers

Promiscuous human leukocyte antigen (HLA) binding peptides are ideal targets for vaccine development. Existing computational models for prediction of promiscuous peptides used hidden Markov models and artificial neural networks as prediction algorithms. We report a system based on support vector machines that outperforms previously published methods. Preliminary testing showed that it can predict peptides binding to HLA-A2 and -A3 super-type molecules with excellent accuracy, even for molecules where no binding data are currently available.

GCMC-surface area of carbonaceous materials with N-2 and Ar adsorption as an alternative to the classical BET method

In this paper we apply a new method for the determination of surface area of carbonaceous materials, using the local surface excess isotherms obtained from the Grand Canonical Monte Carlo simulation and a concept of area distribution in terms of energy well-depth of solid–fluid interaction. The range of this well-depth considered in our GCMC simulation is from 10 to 100 K, which is wide enough to cover all carbon surfaces that we dealt with (for comparison, the well-depth for perfect graphite surface is about 58 K). Having the set of local surface excess isotherms and the differential area distribution, the overall adsorption isotherm can be obtained in an integral form. Thus, given the experimental data of nitrogen or argon adsorption on a carbon material, the differential area distribution can be obtained from the inversion process, using the regularization method. The total surface area is then obtained as the area of this distribution. We test this approach with a number of data in the literature, and compare our GCMC-surface area with that obtained from the classical BET method. In general, we find that the difference between these two surface areas is about 10%, indicating the need to reliably determine the surface area with a very consistent method. We, therefore, suggest the approach of this paper as an alternative to the BET method because of the long-recognized unrealistic assumptions used in the BET theory. Beside the surface area obtained by this method, it also provides information about the differential area distribution versus the well-depth. This information could be used as a microscopic finger-print of the carbon surface. It is expected that samples prepared from different precursors and different activation conditions will have distinct finger-prints. We illustrate this with Cabot BP120, 280 and 460 samples, and the differential area distributions obtained from the adsorption of argon at 77 K and nitrogen also at 77 K have exactly the same patterns, suggesting the characteristics of this carbon.