Hybrid Modelling of Titan's Interaction with the Magnetosphere of Saturn

In this dissertation we study the interaction between Saturn's moon Titan and the magnetospheric plasma and magnetic field. The method of research is a three-dimensional computer simulation model, that is used to simulate this interaction. The simulation model used is a hybrid model. Hybrid models enable individual tracking or tracing of ions and also take into account the particle motion in the propagation of the electromagnetic fields. The hybrid model has been developed at the Finnish Meteorological Institute. This thesis gives a general description of the effects that the solar wind has on Earth and other planets of our solar system. Planetary satellites can also have similar interactions with the solar wind but also with the plasma flows of planetary magnetospheres. Titan is clearly the largest among the satellites of Saturn and also the only known satellite with a dense atmosphere. It is the atmosphere that makes Titan's plasma interaction with the magnetosphere of Saturn so unique. Nevertheless, comparisons with the plasma interactions of other solar system bodies are valuable. Detecting charged plasma particles requires in situ measurements obtainable through scientific spacecraft. The Cassini mission has been one of the most remarkable international efforts in space science. Since 2004 the measurements and images obtained from instruments onboard the Cassini spacecraft have increased the scientific knowledge of Saturn as well as its satellites and magnetosphere in a way no one was probably able to predict. The current level of science on Titan is practically unthinkable without the Cassini mission. Many of the observations by Cassini instrument teams have influenced this research both the direct measurements of Titan as well as observations of its plasma environment. The theoretical principles of the hybrid modelling approach are presented in connection to the broader context of plasma simulations. The developed hybrid model is described in detail: e.g. the way the equations of the hybrid model are solved is shown explicitly. Several simulation techniques, such as the grid structure and various boundary conditions, are discussed in detail as well. The testing and monitoring of simulation runs is presented as an essential routine when running sophisticated and complex models. Several significant improvements of the model, that are in preparation, are also discussed. A main part of this dissertation are four scientific articles based on the results of the Titan model. The Titan model developed during the course of the Ph.D. research has been shown to be an important tool to understand Titan's plasma interaction. One reason for this is that the structures of the magnetic field around Titan are very much three-dimensional. The simulation results give a general picture of the magnetic fields in the vicinity of Titan. The magnetic fine structure of Titan's wake as seen in the simulations seems connected to Alfvén waves an important wave mode in space plasmas. The particle escape from Titan is also a major part of these studies. Our simulations show a bending or turning of Titan's ionotail that we have shown to be a direct result of the basic principles in plasma physics. Furthermore, the ion flux from the magnetosphere of Saturn into Titan's upper atmosphere has been studied. The modelled ion flux has asymmetries that would likely have a large impact in the heating in different parts of Titan's upper atmosphere.

Exploring the interaction energies for the binding of hydroxydiphenyl ethers to enoyl-acyl carrier protein reductases

It is now well established that the potent anti-microbial compound, triclosan, interrupts the type II fatty acid synthesis by inhibiting the enzyme enoyl-ACP reductase in a number of organisms. Existence of a high degree of similarity between the recently discovered enoyl-ACP reductase from R falciparum and B. napus enzyme permitted building of a satisfactory model for the former enzyme that explained some of the key aspects of the enzyme such as its specificity for binding to the cofactor and the inhibitor. We now report the interaction energies between triclosan and other hydroxydiphenyl ethers with the enzymes from B. napus, E. coli and R falciparum. Examination of the triclosan-enzyme interactions revealed that subtle differences exist in the ligand binding sites of the enzymes from different sources i.e., B. napus, E. coli and P falciparum. A comparison of their binding propensities thus determined should aid in the design of effective inhibitors for the respective enzymes.

Studies on the interaction of concanavalin A with glycoproteins

Lectins (phytohaemagglutinin) are known to have the unique property of binding with certain specific sugars, polysaccharides and glycoproteins. Although the kinetics of interaction between lectins and sugar have been extensively studied, the binding characteristics of the lectins with various glycoproteins are not well understood. In this laboratory a systematic study has been initiated in relation to the interaction of lectins with glycoproteins. Concanavalin A is known to bind alpha-glucosides, mannosides and biopolymers having these sugar configurations. A galactose binding protein from caster bean has been purified to homogeneity and was found to contain mannose. This lectin was used as the source of glycoprotein for studying its interaction with concanavalin A. This study showed that the interaction is temperature dependent and the dissociation is time and alpha-methyl glucoside concentration dependent. This has led to speculate a model for cell-lectin interaction. Using concanavalin A it has been shown that all the lysosomal enzymes from brain studied were glycoprotein in nature. Moreover, using Sepharose-bound concanavalin A it has been possible to devise a method by which these lysosomal enzymes could be purified considerably. With the knowledge that the interaction between lectin and glycoprotein is not only dependent on the specific sugar present in the glycoprotein, but also on the nature of the glycoprotein it was possible to develop a novel method for immobilizing various glycoprotein enzymes, such as arylsulphatase A, hyaluronidase and glucose oxidase.

Folding a Polymer via Two-Point Interaction with an External Folding Agent: Use of H-Bonding and Charge-Transfer Interactions

A polymer containing electron-rich aromatic donors (1,5-dialkoxynaphthalene (DAN)) was coerced into a folded state by an external folding agent that contained an electron-deficient aromatic acceptor (pyromellitic diimide (PM)) unit. The donor-containing polymer was designed to carry a tertiary amine moiety in the linking segment, which served as an H-bonding site for reinforcing the interaction with the acceptor containing folding agent that also bore a carboxylic acid group. The H-bonding interaction of the carboxylic acid and the tertiary amine brings the PDI unit between two adjacent DAN units along the polymer backbone to induce charge-transfer (C-T) interactions, and this in turn causes the polymer chain to form a pleated structure. Evidence for the formation of such a pleated structure was obtained from NMR titration studies and also by monitoring the C-T band in their UV-visible spectra. By varying the length of the segment that links the PDI acceptor to the carboxylic acid group, we showed that the most effective folding agent was the one that had a single carbon spacer, as evident from the highest value of the association constant. Control experiments with propionic acid clearly demonstrated the importance of the additional C-T interactions for venerating the folded structures. Further, solution viscosity measurements in the presence of varying amounts of the folding agent revealed a gradual stiffening of the chain in the case of the PDI carrying carboxylic acid, whereas no such affect was seen in the case of simple propionic acid. These observations were supported by D FT calculations of the interactions of a dimeric model of the polymer with the various folding agents; here too the stability of the complex was seen to be highest in the case of the single carbon spacer.

Studies on interaction of colloidal Ag nanoparticles with Bovine Serum Albumin (BSA)

Biofunctionalization of noble metal nanoparticles like Ag, Au is essential to obtain biocompatibility for specific biomedical applications. Silver nanciparticles are being increasingly used in bio-sensing applications owing to excellent optoelectronic properties. Among the serum albumins, the most abundant proteins in plasma, a wide range of physiological functions of Bovine Serum Albumin (BSA) has made it a model system for biofunctionalization. In absence of adequate prior reports, this study aims to investigate the interaction between silver nanoparticles and BSA. The interaction of BSA [0.05-0.85% concentrations] with Ag nanoparticles [50 ppm concentration] in aqueous dispersion was Studied through UV-vis spectral changes, morphological and surface structural changes. At pH 7, which is More than the isoelectric point of BSA, a decrease in absorbance at plasmon peak of uninteracted nanciparticles (425 mn) was noted till 0.45% BSA, beyond that a blue shift towards 410 urn was observed. The blue shift may be attributed to enhanced electron density on the particle surfaces. Increasing pH to 12 enhanced the blue shift further to 400 rim. The conformational changes in BSA at alkaline pH ranges and consequent hydrophobic interactions also played an important role. The equilibrium adsorption data fitted better to Freundlich isotherm compared to Langmuir Curve. The X-ray diffraction study revealed complete coverage of Ag nanoparticles by BSA. The scanning electron microscopic study of the interacted nanoparticles was also carried Out to decipher morphological changes. This study established that tailoring the concentration of BSA and pH of the interaction it was possible to reduce aggregation of nanoparticles. Biofunctionalized Ag nanoparticles with reduced aggregation will be more amenable towards bio-sensing applications. (C) 2009 Elsevier B.V. All rights reserved.

A neuro-adaptive augmented optimal dynamic inversion approach for effective and efficient treatment of chronic myelogenous leukemia

Combining the advanced techniques of optimal dynamic inversion and model-following neuro-adaptive control design, an efficient technique is presented for effective treatment of chronic myelogenous leukemia (CML). A recently developed nonlinear mathematical model for cell dynamics is used for the control (medication) synthesis. First, taking a set of nominal parameters, a nominal controller is designed based on the principle of optimal dynamic inversion. This controller can treat nominal patients (patients having same nominal parameters as used for the control design) effectively. However, since the parameters of an actual patient can be different from that of the ideal patient, to make the treatment strategy more effective and efficient, a model-following neuro-adaptive controller is augmented to the nominal controller. In this approach, a neural network trained online (based on Lyapunov stability theory) facilitates a new adaptive controller, computed online. From the simulation studies, this adaptive control design approach (treatment strategy) is found to be very effective to treat the CML disease for actual patients. Sufficient generality is retained in the theoretical developments in this paper, so that the techniques presented can be applied to other similar problem as well. Note that the technique presented is computationally non-intensive and all computations can be carried out online.

Multiple support excitations of open-plane frames by a filtered white noise and soil-structure interaction

Seismic structural design is essentially the estimation of structural response to a forced motion, which may be deterministic or stochastic, imposed on the ground. The assumption that the same ground motion acts at every point of the base of the structure (or at every support) is not always justifiable; particularly in case of very large structures when considerable spatial variability in ground motion can exist over significant distances example long span bridges. This variability is partly due to the delay in arrival of the excitation at different supports (which is called the wave passage effect) and due to heterogeneity in the ground medium which results in incoherency and local effects. The current study examines the influence of the wave passage effect (in terms of delay in arrival of horizontal ground excitation at different supports and neglecting transmission through the structure) on the response of a few open-plane frame building structures with soil-structure interaction. The ground acceleration has been modeled by a suitably filtered white noise. As a special case, the ground excitation at different supports has also been treated as statistically independent to model the extreme case of incoherence due to local effects and due to modifications to the ground motion resulting from wave reflections and refractions in heterogeneous soil media. The results indicate that, even for relatively short spanned building frames, wave passage effect can be significant. In the absence of soil-structure interaction, it can significantly increase the root mean square (rms) value of the shear in extreme end columns for the stiffer frames but has negligible effect on the flexible frames when total displacements are considered. It is seen that pseudo-static displacements increasingly contribute to the rms value of column shear as the time delay increases both for the stiffer and for the more flexible frames. When soil-structure interaction is considered, wave passage effect (in terms of total displacements) is significant only for low soil shear modulus, G. values (where soil-structure interaction significantly lowers the fundamental frequency) and for stiff frames. The contribution of pseudo-static displacement to these rms values is found to decrease with increase in G. In general, wave passage effect for most interactive frames is insignificant compared to the attenuating effect a decrease in G, has on the response of the interactive structure to uniform support excitation. When the excitations at different supports are statistically independent, it is seen that for both the stiff and flexible frames, the rms value of the column shear in extreme end columns is several times larger (more for the stiffer frames) than the value corresponding to uniform base excitation with the pseudo-static displacements contributing over 99% of the rms value of column shear. Soil-structure interaction has an attenuating effect on the rms value of the column shear, the effect decreasing with increase in G,. Here too, the pseudo-static displacements contribute very largely to the column shear. The influence of the wave passage effect on the response of three 2-bay frames with and without soil-structure interaction to a recorded horizontal accelerogram is also examined. (C) 2010 Elsevier Ltd. All rights reserved.

A soluble random-matrix model for relaxation in quantum systems

We study the relaxation of a degenerate two-level system interacting with a heat bath, assuming a random-matrix model for the system-bath interaction. For times larger than the duration of a collision and smaller than the Poincaré recurrence time, the survival probability of still finding the system at timet in the same state in which it was prepared att=0 is exactly calculated.

Aerodynamic noise from wave-turbulence interaction

Abstract. In order to estimate the acoustic energy scattered when a unit volume of free turbulence, such as in free jets, interacts with a plane steady sound wave, theoretical expressions are derived for two simple models of turbulence: eddy model and isotropic model. The effect of convection by mean motion of the energy-bearing eddies on the incident sound wave and on the sound generated from wave-turbulence interaction is taken into account. Finally, by means of a representative calculation,the directionality pattern and Mach number dependence of the noise so generated is discussed.

Interaction of synthetic analogs of distamycin and netropsin with nucleic acids. Does curvature of ligand play a role in distamycin-DNA interactions?

Distamycin and netropsin, a class of minor groove binding nonintercalating agents, are characterized by their B-DNA and A-T basespecific interactions. To understand the CQI I ~OIT~ ~ I ~ ~aOnMd ~c hemical basis of the above specificities, the DNA-binding characteristics of a novel synthetic analogue of distamycin have been studied. The analogue, mPD derivative, has the requisite charged end groups and a number of potential hydrogen-bonding loci equal to those of distamycin. The difference in the backbone curvatures of the ligands, distamycin, the mPD derivative, and NSC 101327 (another structurally analogous compound),is a major difference between these ligands. UV and CD spectrosoopic studies reported here show the following salient features: The mPD derivative recognizes only B-DNA, to which it binds via the minor groove. On the other hand, unlike distamycin, it binds with comparable affinities to A-T and G-C base pairs in a natural DNA. These DNA-binding properties are compared with those reported earlier for distamycin and NSC 101327 [Zimmer, Ch., & Wahnert, U. (1986) Prog. Biophys. Mol. Biol. 47, 31-1121. The backbone structures of these three ligands were compared to show the progressive decrease in curvatures in the order distamycin, mPD derivative, and NSC 101327. The plausible significance of the backbone curvature vis-&vis the characteristic B-DNA and AT-specific binding of distamycin is discussed. To our knowledge, this is the first attempt (with a model synthetic analogue) to probe the possible influence of backbone curvature upon the specificity of interactions of the distamycin class of groove-binding ligands with DNA.

Stabilization of the bio-membrane by small molecules: interaction of trehalose with the phospholipid bilayer

Anhydrobiotic organisms undergo periods of acute dehydration during their life cycle. It is of interest to understand how the biomembrane remains intact through such stress. A disaccharide, trehalose, which is metabolised during anhydrobiosis is found to prevent disruption of model membrane systems. Molecular modelling techniques are used to investigate the possible mode of interaction of trehalose with a model monolayer. The objective is to maximise hydrogen bonding between the two systems. A phospholipid matrix consisting of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) is chosen to represent the monolayer. The crystal structure of DMPC reveals that there are two distinct conformers designated as A and B. An expansion of the monolayer, coplanar with its surface, results in the trehalose molecule being accommodated in a pocket formed by four B conformers. One glucose ring of the sugar rests on the hydrophobic patch provided by the choline methyls of an A conformer. Five hydrogen bonds are formed involving the phosphate oxygens of three of the surrounding B conformers. The model will be discussed with reference to relevant experimental data on the interaction.

Theoretical Treatment of Interaction Parameters in Multicomponent Metallic Solutions

A new theoretical equation for interaction parameter in multicomponent metallic solutions is developed using the pseudopotential formalism coupled with the free energy of the hard sphere system. The approximate expression for the pseudopotential term is given in terms of the heat of solution at infinite dilution, to allow easy evaluation of the interaction parameter in various multicomponent systems. This theory has been applied to 23 non-ferrous alloys based on Pb, Sn, Bi and indium. Comparison with the results of previous theoretical calculations using only the hard sphere model suggests that the inclusion of the pseudopotential term yields a quantitatively more correct prediction of interaction parameters in multicomponent metallic solutions. Numerical calculations were also made for 320 Fe-base solutions relevant to steelmaking and the agreement between calculation and experimental data appears reasonable, with 90% reliability in predicting the correct sign.

Interaction of melittin with endotoxic lipid

Several amphipathic and cationic substances are known to bind lipid A, the toxic component of bacterial lipopolysaccharides. In this report, we have characterized, by fluorescence methods, the interaction of melittin, an amphipathic and basic 26-residue polypeptide isolated from bee venom, with lipid A. The stoichiometry of the complex appears to be two molecules of melittin to one of lipid A with a dissociation constant of 2.5 × 10 −6 M. The binding of melittin not only modifies the endotoxic properties of lipid A in a number of biological assays, but also results in abrogation of the hemolytic activity of melittin. A model of the complex is proposed based on the known structures of lipid A and melittin, and the observed stoichiometry of binding.

Metal-Insulator Transition in the Hubbard Model on a Triangular Lattice

We discuss the results of an extensive mean-field investigation of the half-filled Hubbard model on a triangular lattice at zero temperature. At intermediate U we find a first-order metal-insulator transition from an incommensurate spiral magnetic metal to a semiconducting state with a commensurate linear spin density wave ordering stabilized by the competition between the kinetic energy and the frustrated nature of the magnetic interaction. At large U the ground state is that of a classical triangular antiferromagnet within our approximation. In the incommensurate spiral metallic phase the Fermi surface has parts in which the wave function renormalization Z is extremely small. The evolution of the Fermi surface and the broadening of the quasi-particle band along with the variation of the plasma frequency and a charge stiffness constant with U/t are discussed.

A new model for the breakage frequency of drops in turbulent stirred dispersions

A model of drop breakage in turbulent stirred dispersions based on interaction of a drop with eddies of a length scale smaller than the drop diameter has been developed. It predicts that, unlike the equal breakage assumed by earlier models, a large drop reduces in size due to stripping of smaller segments off it through unequal breakage. It is only when the drop nears the value of the maximum stable drop diameter that it breaks into equal parts. This new model of drop breakage, coupled with the pattern of interaction of drops with eddies of different sizes existing in the vessel, has been used to evaluate not only the breakage frequency, but also the size distribution of the daughter droplets(which was hitherto assumed). The model has been incorporated in the population balance equation and the resulting cumulative size distributions compared with those availble in the literature.