Computational mechanics of molecular systems:quantifying high-dimensional dynamics by distribution of Poincaré recurrence times

A framework that connects computational mechanics and molecular dynamics has been developed and described. As the key parts of the framework, the problem of symbolising molecular trajectory and the associated interrelation between microscopic phase space variables and macroscopic observables of the molecular system are considered. Following Shalizi and Moore, it is shown that causal states, the constituent parts of the main construct of computational mechanics, the e-machine, define areas of the phase space that are optimal in the sense of transferring information from the micro-variables to the macro-observables. We have demonstrated that, based on the decay of their Poincare´ return times, these areas can be divided into two classes that characterise the separation of the phase space into resonant and chaotic areas. The first class is characterised by predominantly short time returns, typical to quasi-periodic or periodic trajectories. This class includes a countable number of areas corresponding to resonances. The second class includes trajectories with chaotic behaviour characterised by the exponential decay of return times in accordance with the Poincare´ theorem.

Mechanisms of secondary dispersions separation in particulate beds

The mechanisms by which drops of secondary liquid dispersion ie. <100μ m, are collected, coalesced and transferred have been studied in particulate beds of different sizes and heights of glass ballotini. The apparatus facilitated different coalescer cell arrangements. The liquid-liquid system was toluene/de-ionised water. The inlet drop size distribution was measured by microscopy and using the Malvern Particle Size analyser; the outlet dispersion was sized by photography. The effect of packed height and packing size upon critical velocity, pressure drop and coalescence efficiency have been investigated. Single and two phase flow pressure drops across the packing were correlated by modified Blake-Kozeny equations. Two phase pressure drop was correlated by two equations, one for large ballotini sizes (267μm - 367μm), the other for small ballotini sizes (93μm- 147.5μm). The packings were efficient coalescers up to critical velocities of 3 x 10-2 m/s to 5 x 10-2 m/s. The saturation was measured across the bed using relative permeability and a mathematical model developed which related this profile to measured pressure drops. Filter coefficients for the range of packing studied were found to be accurately predicted from a modified queueing drop model. 

In situ imaging and height reconstruction of phase separation processes in polymer blends during spin coating

Spin coating polymer blend thin films provides a method to produce multiphase functional layers of high uniformity covering large surface areas. Applications for such layers include photovoltaics and light-emitting diodes where performance relies upon the nanoscale phase separation morphology of the spun film. Furthermore, at micrometer scales, phase separation provides a route to produce self-organized structures for templating applications. Understanding the factors that determine the final phase-separated morphology in these systems is consequently an important goal. However, it has to date proved problematic to fully test theoretical models for phase separation during spin coating, due to the high spin speeds, which has limited the spatial resolution of experimental data obtained during the coating process. Without this fundamental understanding, production of optimized micro- and nanoscale structures is hampered. Here, we have employed synchronized stroboscopic illumination together with the high light gathering sensitivity of an electron-multiplying charge-coupled device camera to optically observe structure evolution in such blends during spin coating. Furthermore the use of monochromatic illumination has allowed interference reconstruction of three-dimensional topographies of the spin-coated film as it dries and phase separates with nanometer precision. We have used this new method to directly observe the phase separation process during spinning for a polymer blend (PS-PI) for the first time, providing new insights into the spin-coating process and opening up a route to understand and control phase separation structures. © 2011 American Chemical Society.

A method of fundamental solutions for the one-dimensional inverse Stefan problem

We investigate an application of the method of fundamental solutions (MFS) to the one-dimensional inverse Stefan problem for the heat equation by extending the MFS proposed in [5] for the one-dimensional direct Stefan problem. The sources are placed outside the space domain of interest and in the time interval (-T, T). Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate and stable results can be obtained efficiently with small computational cost.

Numerical approximation of the one-dimensional inverse Cauchy–Stefan problem using a method of fundamental solutions

We investigate an application of the method of fundamental solutions (MFS) to the one-dimensional parabolic inverse Cauchy–Stefan problem, where boundary data and the initial condition are to be determined from the Cauchy data prescribed on a given moving interface. In [B.T. Johansson, D. Lesnic, and T. Reeve, A method of fundamental solutions for the one-dimensional inverse Stefan Problem, Appl. Math Model. 35 (2011), pp. 4367–4378], the inverse Stefan problem was considered, where only the boundary data is to be reconstructed on the fixed boundary. We extend the MFS proposed in Johansson et al. (2011) and show that the initial condition can also be simultaneously recovered, i.e. the MFS is appropriate for the inverse Cauchy-Stefan problem. Theoretical properties of the method, as well as numerical investigations, are included, showing that accurate results can be efficiently obtained with small computational cost.

Numerical solution of a Cauchy problem for Laplace equation in 3-dimensional domains by integral equations

A numerical method based on integral equations is proposed and investigated for the Cauchy problem for the Laplace equation in 3-dimensional smooth bounded doubly connected domains. To numerically reconstruct a harmonic function from knowledge of the function and its normal derivative on the outer of two closed boundary surfaces, the harmonic function is represented as a single-layer potential. Matching this representation against the given data, a system of boundary integral equations is obtained to be solved for two unknown densities. This system is rewritten over the unit sphere under the assumption that each of the two boundary surfaces can be mapped smoothly and one-to-one to the unit sphere. For the discretization of this system, Weinert’s method (PhD, Göttingen, 1990) is employed, which generates a Galerkin type procedure for the numerical solution, and the densities in the system of integral equations are expressed in terms of spherical harmonics. Tikhonov regularization is incorporated, and numerical results are included showing the efficiency of the proposed procedure.

Resonance properties of domain boundaries in quasi-two dimensional antiferromagnets

The so-called "internal modes" localized near the domain boundaries in quasi-two dimensional antiferromagnets are investigated. The possible localized states are classified and their frequency dependences on the system discreteness parameter λ=J/β, which describes the ratio of the magnitudes of the exchange interplane interaction and the magnetic anisotropy, are found. A sudden change in the spectrum of the local internal modes is observed at a critical value of this parameter, λ=λb=3/4, where the domain wall shifts from a collinear to a canted shape. When λ<λb there are one symmetric and two antisymmetric local modes, and when λ>λb the modes are two symmetric, one antisymmetric, and one shear. For discreteness parameters close to the critical value, the frequencies of some of the local modes lie deep inside the gap for the linear AFM magnon spectrum and can be observed experimentally. © 2010 American Institute of Physics.

One-dimensional and two-dimensional polyacrylamide gel electrophoresis: a tool for protein characterisation in aquatic samples

Dissolved organic nitrogen (DON) represents the least understood part of the nitrogen cycle. Due to recent methodological developments, proteins now represent a potentially characterisable fraction of DON at the macromolecular level. We have applied polyacrylamide gel electrophoresis to characterise proteins in samples from a range of aquatic environments in the Everglades National Park, Florida, USA. Sodium dodecyl sulphate polyacrylamide gel electrophoresis (SDS-PAGE) showed that each sample has a complex and characteristic protein distribution. Some proteins appeared to be common to more than one site, and these might derive from dominant higher plant vegetation. Two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) provided better resolution; however, strong background hindered interpretation. Our results suggest that the two techniques can be used in parallel as a tool for protein characterisation: SDS-PAGE to provide a sample-specific fingerprint and 2D-PAGE to focus on the characterisation of individual protein molecules.

Organic matter deposition/resuspension in a one-dimensional physical-biogeochemical model. A modelling study

The shallow water configuration of the gulf of Trieste allows the propagation of the stress due to wind and waves along the whole water column down to the bottom. When the stress overcomes a particular threshold it produces resuspension processes of the benthic detritus. The benthic sediments in the North Adriatic are rich of organic matter, transported here by many rivers. This biological active particulate, when remaining in the water, can be transported in all the Adriatic basin by the basin-wide circulation. In this work is presented a first implementation of a resuspension/deposition submodel in the oceanographic coupled physical-biogeochemical 1-dimensional numerical model POM-BFM. At first has been considered the only climatological wind stress forcing, next has been introduced, on the surface, an annual cycle of wave motion and finally have been imposed some exceptional wave event in different periods of the year. The results show a strong relationship between the efficiency of the resuspension process and the stratification of the water column. During summer the strong stratification can contained a great quantity of suspended matter near to the bottom, while during winter even a low concentration of particulate can reach the surface and remains into the water for several months without settling and influencing the biogeochemical system. Looking at the biologic effects, the organic particulate, injected in the water column, allow a sudden growth of the pelagic bacteria which competes with the phytoplankton for nutrients strongly inhibiting its growth. This happen especially during summer when the suspended benthic detritus concentration is greater.

Dynamics of spatial heterogeneity in a landfill with interacting phase densities:a stochastic analysis

A landfill represents a complex and dynamically evolving structure that can be stochastically perturbed by exogenous factors. Both thermodynamic (equilibrium) and time varying (non-steady state) properties of a landfill are affected by spatially heterogenous and nonlinear subprocesses that combine with constraining initial and boundary conditions arising from the associated surroundings. While multiple approaches have been made to model landfill statistics by incorporating spatially dependent parameters on the one hand (data based approach) and continuum dynamical mass-balance equations on the other (equation based modelling), practically no attempt has been made to amalgamate these two approaches while also incorporating inherent stochastically induced fluctuations affecting the process overall. In this article, we will implement a minimalist scheme of modelling the time evolution of a realistic three dimensional landfill through a reaction-diffusion based approach, focusing on the coupled interactions of four key variables - solid mass density, hydrolysed mass density, acetogenic mass density and methanogenic mass density, that themselves are stochastically affected by fluctuations, coupled with diffusive relaxation of the individual densities, in ambient surroundings. Our results indicate that close to the linearly stable limit, the large time steady state properties, arising out of a series of complex coupled interactions between the stochastically driven variables, are scarcely affected by the biochemical growth-decay statistics. Our results clearly show that an equilibrium landfill structure is primarily determined by the solid and hydrolysed mass densities only rendering the other variables as statistically "irrelevant" in this (large time) asymptotic limit. The other major implication of incorporation of stochasticity in the landfill evolution dynamics is in the hugely reduced production times of the plants that are now approximately 20-30 years instead of the previous deterministic model predictions of 50 years and above. The predictions from this stochastic model are in conformity with available experimental observations.

Degeneracy of velocity constraints in rigid body systems

The equations governing the dynamics of rigid body systems with velocity constraints are singular at degenerate configurations in the constraint distribution. In this report, we describe the causes of singularities in the constraint distribution of interconnected rigid body systems with smooth configuration manifolds. A convention of defining primary velocity constraints in terms of orthogonal complements of one-dimensional subspaces is introduced. Using this convention, linear maps are defined and used to describe the space of allowable velocities of a rigid body. Through the definition of these maps, we present a condition for non-degeneracy of velocity constraints in terms of the one dimensional subspaces defining the primary velocity constraints. A method for defining the constraint subspace and distribution in terms of linear maps is presented. Using these maps, the constraint distribution is shown to be singular at configuration where there is an increase in its dimension.

Asymptotic scattering and duality in the one-dimensional three-state quantum Potts model on a lattice

We determine numerically the single-particle and the two-particle spectrum of the three-state quantum Potts model on a lattice by using the density matrix renormalization group method, and extract information on the asymptotic (small momentum) S-matrix of the quasiparticles. The low energy part of the finite size spectrum can be understood in terms of a simple effective model introduced in a previous work, and is consistent with an asymptotic S-matrix of an exchange form below a momentum scale p*. This scale appears to vanish faster than the Compton scale, mc, as one approaches the critical point, suggesting that a dangerously irrelevant operator may be responsible for the behaviour observed on the lattice.

Separation and identification fresh water leeches of waters resources in Kordistan Province

In the primary studying, known leeches have included into free living and parasitic which the parasitic group, besides of direct impacts like: growth detraction, anemia, making wound in the connecting part with the skin, with making plat for secondary bacterial and molding infections in the toll place, is able to cause to transfer blood flagellate and virus to the fish too. Therefore, by accusing information related to fauna leeches of each area a risk factor from the viewpoint of the possibility of being or accusing transferred diseases by these leeches, one can predict and forehand about them. Freshwater leeches of Iran to present accurate recognition (morphological, molecular) have not been, and there are some limited reports from different parts of the country about them. One of the areas that its leeches have not been identified yet is Kurdistan, By having five latrines and big permanent rivers and 32000 springs and a lot of deep and semi-deep wells and this province is a convenient bed for growing aquatics in the country. Therefore, identifying risk factors for development of aquaculture on water resources is one important factor to access achieving development goals. For recognizing leeches of this province, some samples from 10 stands were token. Samples from under stones, sticking to the fish, turtles, plants and solid substances in the water were separated and after recording their physical characteristics, calming with 10% ethylic alcohol with 10% formalin become fixed and after painting with Carmen acetic acid by standard keys for 7 species of Helobdella stagnalis, Placobdella costata, Hemiclepsis marginata, Erpobdella octoculata, Hirudo medicinalis, Dina lineate lineata have been identified and described. Which Helobdella stagnalis has the highest distribution in the province and the minimum one is Hirudo medicinalis. However, that the data obtained in leeches in Kurdistan is a relatively complete collection in this research, recognizing fauna of these areas needs more studying. The Placobdella costata and Hemiclepsis marginata sticking to the fish were separated among identified species which showed that these are parasites for the fish. The sticking area of those leeches to the skin was accompanied with scales cast, damage to mucous membranes beneath the parasite and bleeding Was associated with Histopathology studying effect includes observing break and disconnection in the leech connecting place to the epithelial layer of epidermis in the skin, destroyed nucleus in skin Epithelial cells with observing necrosis in ulcerative place become of the leech and the sub acute inflammated penetration until acute necrosis with opening in Dermis layer is observable. Kidney of this fish have changes such as: proliferation, like proliferative kidney disease with increasing proliferative glomerular cells and increasing in membranous cells in Capillary corpuscle, observing necrotic cells in haematopoietic tissue of kidney along with increasing in infiltration of leukocyte's cells generally mono nucluars such as lymphocytes and less poly morpho nucluars such as neutrophiles that are symptoms of disorders causing anemia become of nourishing and sucking blood by the leech and creating a chronic kidney infection that originally root is in another place like the skin. Also Hemorrhagic anemia causes losing RBC's is because of using the host blood by the leech. (In this situation, one can see immature RBC red cells in Peripheral blood. To identify potential carriers of the leech to the viruses, after finding them in recorded stands and putting them in 75% ethanol for viruses cause IPN, VHS, IHN, they were tested by PCR that the conclusion of these experiments approved IPN virus in Hemiclepsis marginata and Hirudo medicinalis. This kind of leeches can act like a mechanical carrier and causing spreading the agent of this disease. It is worth mentioning that studying the pathogenicity of this virus for aquaculture sources, mentioned before needs more research. During the study of infected fish with leeches that was done after preparing bloody slides and staining them, no case blood parasites was observed. During a research about infecting fish experimentally to known leeches it become clear that 5 days after being in aquarium including leeches, samples of sticking Hirudo medicinalis leech to the golden carp with scales cast were observed. Including leeches to the fish started with molting the scales in the sticking area in the fish and fish become too uneasy and by rubbing themselves to the malls and things inside the aquarium, tried to separate them. Finally, after around 30 hours, leeches penetrate the skin, feeding from blood and tissue liquids and cause mortality the fish and then they become separated from them. If the corpse of these fish stayed in the aquarium, the Helobdella stagnalis and Erpobdella octoculata would start feeding them.