Two-phase flow analogy as an effective boundary condition for modelling liquids at atomistic resolution

A hybrid Molecular Dynamics/Fluctuating Hydrodynamics framework based on the analogy with two-phase hydrodynamics has been extended to dynamically tracking the feature of interest at all-atom resolution. In the model, the hydrodynamics description is used as an effective boundary condition to close the molecular dynamics solution without resorting to standard periodic boundary conditions. The approach is implemented in a popular Molecular Dynamics package GROMACS and results for two biomolecular systems are reported. A small peptide dialanine and a complete capsid of a virus porcine circovirus 2 in water are considered and shown to reproduce the structural and dynamic properties compared to those obtained in theory, purely atomistic simulations, and experiment.

Immunological synapse: a mathematical model of the bond formation process:mathematical modelling of the immature synapse

The cell:cell bond between an immune cell and an antigen presenting cell is a necessary event in the activation of the adaptive immune response. At the juncture between the cells, cell surface molecules on the opposing cells form non-covalent bonds and a distinct patterning is observed that is termed the immunological synapse. An important binding molecule in the synapse is the T-cell receptor (TCR), that is responsible for antigen recognition through its binding with a major-histocompatibility complex with bound peptide (pMHC). This bond leads to intracellular signalling events that culminate in the activation of the T-cell, and ultimately leads to the expression of the immune eector function. The temporal analysis of the TCR bonds during the formation of the immunological synapse presents a problem to biologists, due to the spatio-temporal scales (nanometers and picoseconds) that compare with experimental uncertainty limits. In this study, a linear stochastic model, derived from a nonlinear model of the synapse, is used to analyse the temporal dynamics of the bond attachments for the TCR. Mathematical analysis and numerical methods are employed to analyse the qualitative dynamics of the nonequilibrium membrane dynamics, with the specic aim of calculating the average persistence time for the TCR:pMHC bond. A single-threshold method, that has been previously used to successfully calculate the TCR:pMHC contact path sizes in the synapse, is applied to produce results for the average contact times of the TCR:pMHC bonds. This method is extended through the development of a two-threshold method, that produces results suggesting the average time persistence for the TCR:pMHC bond is in the order of 2-4 seconds, values that agree with experimental evidence for TCR signalling. The study reveals two distinct scaling regimes in the time persistent survival probability density prole of these bonds, one dominated by thermal uctuations and the other associated with the TCR signalling. Analysis of the thermal fluctuation regime reveals a minimal contribution to the average time persistence calculation, that has an important biological implication when comparing the probabilistic models to experimental evidence. In cases where only a few statistics can be gathered from experimental conditions, the results are unlikely to match the probabilistic predictions. The results also identify a rescaling relationship between the thermal noise and the bond length, suggesting a recalibration of the experimental conditions, to adhere to this scaling relationship, will enable biologists to identify the start of the signalling regime for previously unobserved receptor:ligand bonds. Also, the regime associated with TCR signalling exhibits a universal decay rate for the persistence probability, that is independent of the bond length.

Structural flexibility and oxygen diffusion pathways in monomeric fluorescent proteins

Fluorescent proteins are valuable tools as biochemical markers for studying cellular processes. Red fluorescent proteins (RFPs) are highly desirable for in vivo applications because they absorb and emit light in the red region of the spectrum where cellular autofluorescence is low. The naturally occurring fluorescent proteins with emission peaks in this region of the spectrum occur in dimeric or tetrameric forms. The development of mutant monomeric variants of RFPs has resulted in several novel FPs known as mFruits. Though oxygen is required for maturation of the chromophore, it is known that photobleaching of FPs is oxygen sensitive, and oxygen-free conditions result in improved photostabilities. Therefore, understanding oxygen diffusion pathways in FPs is important for both photostabilites and maturation of the chromophores. We used molecular dynamics calculations to investigate the protein barrel fluctuations in mCherry, which is one of the most useful monomeric mFruit variants, and its GFP homolog citrine. We employed implicit ligand sampling and locally enhanced sampling to determine oxygen pathways from the bulk solvent into the mCherry chromophore in the interior of the protein. The pathway contains several oxygen hosting pockets, which were identified by the amino acid residues that form the pocket. We calculated the free-energy of an oxygen molecule at points along the path. We also investigated an RFP variant known to be significantly less photostable than mCherry and find much easier oxygen access in this variant. We showed that oxygen pathways can be blocked or altered, and barrel fluctuations can be reduced by strategic amino acid substitutions. The results provide a better understanding of the mechanism of molecular oxygen access into the fully folded mCherry protein barrel and provide insight into the photobleaching process in these proteins.

Tracing the behaviour of a bridge using time series modelling

Peer reviewed

Simulating Granite Emplacement and Deformation with Centrifuge Analogue Modelling

The Greater Himalayan leucogranites are a discontinuous suite of intrusions emplaced in a thickened crust during the Miocene southward ductile extrusion of the Himalayan metamorphic core. Melt-induced weakening is thought to have played a critical role in strain localization that facilitated the extrusion. Recent advancements in centrifuge analogue modelling techniques allow for the replication of a broader range of crustal deformation behaviors, enhancing our understanding of large hot orogens. Polydimethylsiloxane (PDMS) is commonly used in centrifuge experiments to model weak melt zones. Difficulties in handling PDMS had, until now, limited its emplacement in models prior to any deformation. A new modelling technique has been developed where PDMS is emplaced into models that have been subjected to some shortening. This technique aims to better understand the effects of melt on strain localization and potential decoupling between structural levels within an evolving orogenic system. Models are subjected to an early stage of shortening, followed by the introduction of PDMS, and then a final stage of shortening. Theoretical percentages of partial melt and their effect on rock strength are considered when adding a specific percentage of PDMS in each model. Due to the limited size of the models, only PDMS sheets of 3 mm thickness were used, which varied in length and width. Within undeformed packages, minimal surface and internal deformation occurred when PDMS is emplaced in the lower layer of the model, showing a vertical volume increase of ~20% within the package; whereas the emplacement of PDMS into the middle layer showed internal dragging of the middle laminations into the lower layer and a vertical volume increase ~30%. Emplacement of PDMS results in ~7% shortening for undeformed and deformed models. Deformed models undergo ~20% additional shortening after two rounds of deformation. Strain localization and decoupling between units occur in deformed models where the degree of deformation changes based on the amount of partial melt present. Surface deformation visible by the formation of a bulge, mode 1 extension cracks and varying surface strain ellipses varies depending if PDMS is present. Better control during emplacement is exhibited when PDMS is added into cooler models, resulting in reduced internal deformation within the middle layer.

Implementing a structural continuity constraint and a halting method for the topology optimization of energy absorbers

This study investigates topology optimization of energy absorbing structures in which material damage is accounted for in the optimization process. The optimization objective is to design the lightest structures that are able to absorb the required mechanical energy. A structural continuity constraint check is introduced that is able to detect when no feasible load path remains in the finite element model, usually as a result of large scale fracture. This assures that designs do not fail when loaded under the conditions prescribed in the design requirements. This continuity constraint check is automated and requires no intervention from the analyst once the optimization process is initiated. Consequently, the optimization algorithm proceeds towards evolving an energy absorbing structure with the minimum structural mass that is not susceptible to global structural failure. A method is also introduced to determine when the optimization process should halt. The method identifies when the optimization method has plateaued and is no longer likely to provide improved designs if continued for further iterations. This provides the designer with a rational method to determine the necessary time to run the optimization and avoid wasting computational resources on unnecessary iterations. A case study is presented to demonstrate the use of this method.

Application of GPCR Structures for Modelling of Free Fatty Acid Receptors

Five G protein-coupled receptors (GPCRs) have been identified to be activated by free fatty acids (FFA). Among them, FFA1 (GPR40) and FFA4 (GPR120) bind long-chain fatty acids, FFA2 (GPR43) and FFA3 (GPR41) bind short-chain fatty acids and GPR84 binds medium-chain fatty acids. Free fatty acid receptors have now emerged as potential targets for the treatment of diabetes, obesity and immune diseases. The recent progress in crystallography of GPCRs has now enabled the elucidation of the structure of FFA1 and provided reliable templates for homology modelling of other FFA receptors. Analysis of the crystal structure and improved homology models, along with mutagenesis data and structure activity, highlighted an unusual arginine charge pairing interaction in FFA1-3 for receptor modulation, distinct structural features for ligand binding to FFA1 and FFA4 and an arginine of the second extracellular loop as a possible anchoring point for FFA at GPR84. Structural data will be helpful for searching novel small molecule modulators at the FFA receptors.

Investigations on the structural, optical and magnetic properties of nanostructured cerium oxide in pure and doped forms and its polymer nanocomposites

In recent years, nanoscience and nanotechnology has emerged as one of the most important and exciting frontier areas of research interest in almost all fields of science and technology. This technology provides the path of many breakthrough changes in the near future in many areas of advanced technological applications. Nanotechnology is an interdisciplinary area of research and development. The advent of nanotechnology in the modern times and the beginning of its systematic study can be thought of to have begun with a lecture by the famous physicist Richard Feynman. In 1960 he presented a visionary and prophetic lecture at the meeting of the American Physical Society entitled “there is plenty of room at the bottom” where he speculated on the possibility and potential of nanosized materials. Synthesis of nanomaterials and nanostructures are the essential aspects of nanotechnology. Studies on new physical properties and applications of nanomaterials are possible only when materials are made available with desired size, morphology, crystal structure and chemical composition. Cerium oxide (ceria) is one of the important functional materials with high mechanical strength, thermal stability, excellent optical properties, appreciable oxygen ion conductivity and oxygen storage capacity. Ceria finds a variety of applications in mechanical polishing of microelectronic devices, as catalysts for three-way automatic exhaust systems and as additives in ceramics and phosphors. The doped ceria usually has enhanced catalytic and electrical properties, which depend on a series of factors such as the particle size, the structural characteristics, morphology etc. Ceria based solid solutions have been widely identified as promising electrolytes for intermediate temperature solid oxide fuel cells (SOFC). The success of many promising device technologies depends on the suitable powder synthesis techniques. The challenge for introducing new nanopowder synthesis techniques is to preserve high material quality while attaining the desired composition. The method adopted should give reproducible powder properties, high yield and must be time and energy effective. The use of a variety of new materials in many technological applications has been realized through the use of thin films of these materials. Thus the development of any new material will have good application potential if it can be deposited in thin film form with the same properties. The advantageous properties of thin films include the possibility of tailoring the properties according to film thickness, small mass of the materials involved and high surface to volume ratio. The synthesis of polymer nanocomposites is an integral aspect of polymer nanotechnology. By inserting the nanometric inorganic compounds, the properties of polymers can be improved and this has a lot of applications depending upon the inorganic filler material present in the polymer.

Modelling and estimating models of physician labour supply and productivity

Ma thèse s’intéresse aux politiques de santé conçues pour encourager l’offre de services de santé. L’accessibilité aux services de santé est un problème majeur qui mine le système de santé de la plupart des pays industrialisés. Au Québec, le temps médian d’attente entre une recommandation du médecin généraliste et un rendez-vous avec un médecin spécialiste était de 7,3 semaines en 2012, contre 2,9 semaines en 1993, et ceci malgré l’augmentation du nombre de médecins sur cette même période. Pour les décideurs politiques observant l’augmentation du temps d’attente pour des soins de santé, il est important de comprendre la structure de l’offre de travail des médecins et comment celle-ci affecte l’offre des services de santé. Dans ce contexte, je considère deux principales politiques. En premier lieu, j’estime comment les médecins réagissent aux incitatifs monétaires et j’utilise les paramètres estimés pour examiner comment les politiques de compensation peuvent être utilisées pour déterminer l’offre de services de santé de court terme. En second lieu, j’examine comment la productivité des médecins est affectée par leur expérience, à travers le mécanisme du "learning-by-doing", et j’utilise les paramètres estimés pour trouver le nombre de médecins inexpérimentés que l’on doit recruter pour remplacer un médecin expérimenté qui va à la retraite afin de garder l’offre des services de santé constant. Ma thèse développe et applique des méthodes économique et statistique afin de mesurer la réaction des médecins face aux incitatifs monétaires et estimer leur profil de productivité (en mesurant la variation de la productivité des médecins tout le long de leur carrière) en utilisant à la fois des données de panel sur les médecins québécois, provenant d’enquêtes et de l’administration. Les données contiennent des informations sur l’offre de travail de chaque médecin, les différents types de services offerts ainsi que leurs prix. Ces données couvrent une période pendant laquelle le gouvernement du Québec a changé les prix relatifs des services de santé. J’ai utilisé une approche basée sur la modélisation pour développer et estimer un modèle structurel d’offre de travail en permettant au médecin d’être multitâche. Dans mon modèle les médecins choisissent le nombre d’heures travaillées ainsi que l’allocation de ces heures à travers les différents services offerts, de plus les prix des services leurs sont imposés par le gouvernement. Le modèle génère une équation de revenu qui dépend des heures travaillées et d’un indice de prix représentant le rendement marginal des heures travaillées lorsque celles-ci sont allouées de façon optimale à travers les différents services. L’indice de prix dépend des prix des services offerts et des paramètres de la technologie de production des services qui déterminent comment les médecins réagissent aux changements des prix relatifs. J’ai appliqué le modèle aux données de panel sur la rémunération des médecins au Québec fusionnées à celles sur l’utilisation du temps de ces mêmes médecins. J’utilise le modèle pour examiner deux dimensions de l’offre des services de santé. En premierlieu, j’analyse l’utilisation des incitatifs monétaires pour amener les médecins à modifier leur production des différents services. Bien que les études antérieures ont souvent cherché à comparer le comportement des médecins à travers les différents systèmes de compensation,il y a relativement peu d’informations sur comment les médecins réagissent aux changementsdes prix des services de santé. Des débats actuels dans les milieux de politiques de santé au Canada se sont intéressés à l’importance des effets de revenu dans la détermination de la réponse des médecins face à l’augmentation des prix des services de santé. Mon travail contribue à alimenter ce débat en identifiant et en estimant les effets de substitution et de revenu résultant des changements des prix relatifs des services de santé. En second lieu, j’analyse comment l’expérience affecte la productivité des médecins. Cela a une importante implication sur le recrutement des médecins afin de satisfaire la demande croissante due à une population vieillissante, en particulier lorsque les médecins les plus expérimentés (les plus productifs) vont à la retraite. Dans le premier essai, j’ai estimé la fonction de revenu conditionnellement aux heures travaillées, en utilisant la méthode des variables instrumentales afin de contrôler pour une éventuelle endogeneité des heures travaillées. Comme instruments j’ai utilisé les variables indicatrices des âges des médecins, le taux marginal de taxation, le rendement sur le marché boursier, le carré et le cube de ce rendement. Je montre que cela donne la borne inférieure de l’élasticité-prix direct, permettant ainsi de tester si les médecins réagissent aux incitatifs monétaires. Les résultats montrent que les bornes inférieures des élasticités-prix de l’offre de services sont significativement positives, suggérant que les médecins répondent aux incitatifs. Un changement des prix relatifs conduit les médecins à allouer plus d’heures de travail au service dont le prix a augmenté. Dans le deuxième essai, j’estime le modèle en entier, de façon inconditionnelle aux heures travaillées, en analysant les variations des heures travaillées par les médecins, le volume des services offerts et le revenu des médecins. Pour ce faire, j’ai utilisé l’estimateur de la méthode des moments simulés. Les résultats montrent que les élasticités-prix direct de substitution sont élevées et significativement positives, représentant une tendance des médecins à accroitre le volume du service dont le prix a connu la plus forte augmentation. Les élasticitésprix croisées de substitution sont également élevées mais négatives. Par ailleurs, il existe un effet de revenu associé à l’augmentation des tarifs. J’ai utilisé les paramètres estimés du modèle structurel pour simuler une hausse générale de prix des services de 32%. Les résultats montrent que les médecins devraient réduire le nombre total d’heures travaillées (élasticité moyenne de -0,02) ainsi que les heures cliniques travaillées (élasticité moyenne de -0.07). Ils devraient aussi réduire le volume de services offerts (élasticité moyenne de -0.05). Troisièmement, j’ai exploité le lien naturel existant entre le revenu d’un médecin payé à l’acte et sa productivité afin d’établir le profil de productivité des médecins. Pour ce faire, j’ai modifié la spécification du modèle pour prendre en compte la relation entre la productivité d’un médecin et son expérience. J’estime l’équation de revenu en utilisant des données de panel asymétrique et en corrigeant le caractère non-aléatoire des observations manquantes à l’aide d’un modèle de sélection. Les résultats suggèrent que le profil de productivité est une fonction croissante et concave de l’expérience. Par ailleurs, ce profil est robuste à l’utilisation de l’expérience effective (la quantité de service produit) comme variable de contrôle et aussi à la suppression d’hypothèse paramétrique. De plus, si l’expérience du médecin augmente d’une année, il augmente la production de services de 1003 dollar CAN. J’ai utilisé les paramètres estimés du modèle pour calculer le ratio de remplacement : le nombre de médecins inexpérimentés qu’il faut pour remplacer un médecin expérimenté. Ce ratio de remplacement est de 1,2.

Random path sampling applied to vector models of magnetism

The recently reported Monte Carlo Random Path Sampling method (RPS) is here improved and its application is expanded to the study of the 2D and 3D Ising and discrete Heisenberg models. The methodology was implemented to allow use in both CPU-based high-performance computing infrastructures (C/MPI) and GPU-based (CUDA) parallel computation, with significant computational performance gains. Convergence is discussed, both in terms of free energy and magnetization dependence on field/temperature. From the calculated magnetization-energy joint density of states, fast calculations of field and temperature dependent thermodynamic properties are performed, including the effects of anisotropy on coercivity, and the magnetocaloric effect. The emergence of first-order magneto-volume transitions in the compressible Ising model is interpreted using the Landau theory of phase transitions. Using metallic Gadolinium as a real-world example, the possibility of using RPS as a tool for computational magnetic materials design is discussed. Experimental magnetic and structural properties of a Gadolinium single crystal are compared to RPS-based calculations using microscopic parameters obtained from Density Functional Theory.

Modelling the growth of Northern Cod

We develop a body size growth model of Northern cod (Gadus morhua) in Northwest Atlantic Fisheries Organization (NAFO) Divisions 2J3KL during 2009-2013. We use individual length-at-age data from the bottom trawl survey in these divisions during 2009–2013. We use the Von Bertalanffy (VonB) model extended to account for between-individual variations in growth, and variations that may be caused by the methods which fish are caught and sampled for length and age measurements. We assume between-individual variation in growth appears because individuals grow at a different rate (k), and they achieve different maximum sizes (l∞). We also included measurement error in length and age in our model since ignoring these errors can lead to biased estimates of the growth parameters. We use the structural errors-invariables (SEV) approach to estimate individual variation in growth, ageing error variation, and the true age distribution of the fish. Our results shows the existence of individual variation in growth and ME in age. According to the negative log likelihood ratio (NLLR) test, the best model indicated: 1) different growth patterns across divisions and years. 2) Between individual variation in growth is the same for the same division across years. 3) The ME in age and true age distribution are different for each year and division.

3D modelling of mango fruit skin blush in the tree canopy

The fruit of certain mango cultivars (e.g., 'Honey Gold') can develop blush on their skin. Skin blush due to red pigmentation is from the accumulation of anthocyanins. Anthocyanin biosynthesis is related to environmental determinants, including light received by the fruit. It has been observed that mango skin blush varies with position in the tree canopy. However, little investigation into this spatial relationship has been conducted. The objective of this preliminary study was to describe a 'Honey Gold' mango tree by capturing its three-dimensional (3D) architecture. A light path tracing model QuasiMC was then used to predict light received by fruit. The use of this 3D model was to better understand the relationship between mango fruit skin blush and fruit position in the canopy. The digitised mango tree mimicked the real tree at a high level of detail. Observations on mango skin blush distribution supported the proposition that sunlight exposure is an absolute requirement for anthocyanin development. No blush development occurred on shaded skin. It was affirmed that 3D mapping could allow for virtual experiments. For example, for virtual canopy thinning (e.g., 'window pruning') to admit more sunlight with a view to improve fruit blush. Improvements to 3D modelling of mango skin blush could focus on increasing accuracy, e.g., measurement of leaf light reflectance and transmission and the inclusion of the effect shading by branches.

Modelling the relationship between biomechanics and performance of young sprinting swimmers

The aim of this study was to compute a swimming performance confirmatory model based on biomechanical parameters. The sample included 100 young swimmers (overall: 12.3 ± 0.74 years; 49 boys: 12.5 ± 0.76 years; 51 girls: 12.2 ± 0.71 years; both genders in Tanner stages 1-2 by self-report) participating on a regular basis in regional and national-level events. The 100 m freestyle event was chosen as the performance indicator. Anthropometric (arm span), strength (throwing velocity), power output (power to overcome drag), kinematic (swimming velocity) and efficiency (propelling efficiency) parameters were measured and included in the model. The path-flow analysis procedure was used to design and compute the model. The anthropometric parameter (arm span) was excluded in the final model, increasing its goodness-of-fit. The final model included the throw velocity, power output, swimming velocity and propelling efficiency. All links were significant between the parameters included, but the throw velocity-power output. The final model was explained by 69% presenting a reasonable adjustment (model's goodness-of-fit; x(2)/df = 3.89). This model shows that strength and power output parameters do play a mediator and meaningful role in the young swimmers' performance.

A contribution towards the performance and structural understanding of copper and nickel salts as ammonia stores

The ongoing depletion of fossil fuels and the severe consequences of the greenhouse effect make the development of alternative energy systems crucially important. While hydrogen is, in principle, a promising alternative, releasing nothing but energy and pure water. Hydrogen storage is complicated and no completely viable technique has been proposed so far. This work is concerned with the study of one potential alternative to pure hydrogen: ammonia, and more specifically its storage in solids. Ammonia, NH3, can be regarded as a chemical hydrogen carrier with the advantages of strongly reduced flammability and explosiveness as compared to hydrogen. Furthermore, ammine metal salts presented here as promising ammonia stores easily store up to 50 wt.-% ammonia, giving them a volumetric energy density comparable to natural gas. The model system NiX2–NH3 ( X = Cl, Br, I) is studied thoroughly with respect to ammine salt formation, thermal decomposition, air stability and structural effects. The system CuX2–NH3 ( X = Cl, Br) has an adverse thermal decomposition behaviour, making it impractical for use as an ammonia store. This system is, however, most interesting from a structural point of view and some work concerning the study of the structural behaviour of this system is presented. Finally, close chemical relatives to the metal ammine halides, the metal ammine nitrates are studied. They exhibit interesting anion arrangements, which is an impressive showcase for the combination of diffraction and spectroscopic information. The characterisation techniques in this thesis range from powder diffraction over single crystal diffraction, spectroscopy, computational modelling, thermal analyses to gravimetric uptake experiments. Further highlights are the structure solutions and refinements from powder data of (NH4)2[NiCl4(H2O)(NH3)] and Ni(NH3)2(NO3)2, the combination of crystallographic and chemical information for the elucidation of the (NH4)2[NiCl4(H2O)(NH3)] formation reaction and the growth of single crystals under ammonia flow, a technique allowing the first documented successful growth and single crystal diffraction measurement for [Cu(NH3)6]Cl2.

Advanced parametrical modelling of 24 GHz radar sensor IC packaging components

This paper deals with the development of an advanced parametrical modelling concept for packaging components of a 24 GHz radar sensor IC used in automotive driver assistance systems. For fast and efficient design of packages for system-in-package modules (SiP), a simplified model for the description of parasitic electromagnetic effects within the package is desirable, as 3-D field computation becomes inefficient due to the high density of conductive elements of the various signal paths in the package. By using lumped element models for the characterization of the conductive components, a fast indication of the design's signal-quality can be gained, but so far does not offer enough flexibility to cover the whole range of geometric arrangements of signal paths in a contemporary package. This work pursues to meet the challenge of developing a flexible and fast package modelling concept by defining parametric lumped-element models for all basic signal path components, e.g. bond wires, vias, strip lines, bumps and balls. © Author(s) 2011. CC Attribution 3.0 License.