Qualitative modelling of unknown interface behaviour

M. H. Lee, and S. M. Garrett, Qualitative modelling of unknown interface behaviour, International Journal of Human Computer Studies, Vol. 53, No. 4, pp. 493-515, 2000

First principles simulation of amorphous silicon bulk, interfaces, and nanowires for photovoltaics

Amorphous silicon has become the material of choice for many technologies, with major applications in large area electronics: displays, image sensing and thin film photovoltaic cells. This technology development has occurred because amorphous silicon is a thin film semiconductor that can be deposited on large, low cost substrates using low temperature. In this thesis, classical molecular dynamics and first principles DFT calculations have been performed to generate structural models of amorphous and hydrogenated amorphous silicon and interfaces of amorphous and crystalline silicon, with the ultimate aim of understanding the photovoltaic properties of core-shell crystalline amorphous Si nanowire structures. We have shown, unexpectedly, from the simulations, that our understanding of hydrogenated bulk a-Si needs to be revisited, with our robust finding that when fully saturated with hydrogen, bulk a-Si exhibits a constant optical energy gap, irrespective of the hydrogen concentration in the sample. Unsaturated a-Si:H, with a lower than optimum hydrogen content, shows a smaller optical gap, that increases with hydrogen content until saturation is reached. The mobility gaps obtained from an analysis of the electronic states show similar behavior. We also obtained that the optical and mobility gaps show a volcano curve as the H content is varied from 7% (undersaturation) to 18% (mild oversaturation). In the case of mild over saturation, the mid-gap states arise exclusively from an increase in the density of strained Si-Si bonds. Analysis of our structures shows the extra H atoms in this case form a bridge between neighboring silicon atoms which increases the corresponding Si-Si distance and promotes bond length disorder in the sample. That has the potential to enhance the Staebler-Wronski effect. Planar interface models of amorphous-crystalline silicon have been generated in Si (100), (110) and (111) surfaces. The interface models are characterized by structure, RDF, electronic density of states and optical absorption spectrum. We find that the least stable (100) surface will result in the formation of the thickest amorphous silicon layer, while the most stable (110) surface forms the smallest amorphous region. We calculated for the first time band offsets of a-Si:H/c-Si heterojunctions from first principles and examined the influence of different surface orientations and amorphous layer thickness on the offsets and implications for device performance. The band offsets depend on the amorphous layer thickness and increase with thickness. By controlling the amorphous layer thickness we can potentially optimise the solar cell parameters. Finally, we have successfully generated different amorphous layer thickness of the a-Si/c-Si and a-Si:H/c-Si 5 nm nanowires from heat and quench. We perform structural analysis of the a-Si-/c-Si nanowires. The RDF, Si-Si bond length distributions, and the coordination number distributions of amorphous regions of the nanowires reproduce similar behaviour compared to bulk amorphous silicon. In the final part of this thesis we examine different surface terminating chemical groups, -H, - OH and –NH2 in (001) GeNW. Our work shows that the diameter of Ge nanowires and the nature of surface terminating groups both play a significant role in both the magnitude and the nature of the nanowire band gaps, allowing tuning of the band gap by up to 1.1 eV. We also show for the first time how the nanowire diameter and surface termination shifts the absorption edge in the Ge nanowires to longer wavelengths. Thus, the combination of nanowire diameter and surface chemistry can be effectively utilised to tune the band gaps and thus light absorption properties of small diameter Ge nanowires.

Working at the interface of phylogenetics and population genetics: a biogeographical analysis of Triaenops spp. (Chiroptera: Hipposideridae).

New applications of genetic data to questions of historical biogeography have revolutionized our understanding of how organisms have come to occupy their present distributions. Phylogenetic methods in combination with divergence time estimation can reveal biogeographical centres of origin, differentiate between hypotheses of vicariance and dispersal, and reveal the directionality of dispersal events. Despite their power, however, phylogenetic methods can sometimes yield patterns that are compatible with multiple, equally well-supported biogeographical hypotheses. In such cases, additional approaches must be integrated to differentiate among conflicting dispersal hypotheses. Here, we use a synthetic approach that draws upon the analytical strengths of coalescent and population genetic methods to augment phylogenetic analyses in order to assess the biogeographical history of Madagascar's Triaenops bats (Chiroptera: Hipposideridae). Phylogenetic analyses of mitochondrial DNA sequence data for Malagasy and east African Triaenops reveal a pattern that equally supports two competing hypotheses. While the phylogeny cannot determine whether Africa or Madagascar was the centre of origin for the species investigated, it serves as the essential backbone for the application of coalescent and population genetic methods. From the application of these methods, we conclude that a hypothesis of two independent but unidirectional dispersal events from Africa to Madagascar is best supported by the data.

Suppression of conformational heterogeneity at a protein-protein interface.

Staphylococcal protein A (SpA) is an important virulence factor from Staphylococcus aureus responsible for the bacterium's evasion of the host immune system. SpA includes five small three-helix-bundle domains that can each bind with high affinity to many host proteins such as antibodies. The interaction between a SpA domain and the Fc fragment of IgG was partially elucidated previously in the crystal structure 1FC2. Although informative, the previous structure was not properly folded and left many substantial questions unanswered, such as a detailed description of the tertiary structure of SpA domains in complex with Fc and the structural changes that take place upon binding. Here we report the 2.3-Å structure of a fully folded SpA domain in complex with Fc. Our structure indicates that there are extensive structural rearrangements necessary for binding Fc, including a general reduction in SpA conformational heterogeneity, freezing out of polyrotameric interfacial residues, and displacement of a SpA side chain by an Fc side chain in a molecular-recognition pocket. Such a loss of conformational heterogeneity upon formation of the protein-protein interface may occur when SpA binds its multiple binding partners. Suppression of conformational heterogeneity may be an important structural paradigm in functionally plastic proteins.

VizMaps: A Bayesian Topic Modeling Based PubMed Search Interface

A common challenge that users of academic databases face is making sense of their query outputs for knowledge discovery. This is exacerbated by the size and growth of modern databases. PubMed, a central index of biomedical literature, contains over 25 million citations, and can output search results containing hundreds of thousands of citations. Under these conditions, efficient knowledge discovery requires a different data structure than a chronological list of articles. It requires a method of conveying what the important ideas are, where they are located, and how they are connected; a method of allowing users to see the underlying topical structure of their search. This paper presents VizMaps, a PubMed search interface that addresses some of these problems. Given search terms, our main backend pipeline extracts relevant words from the title and abstract, and clusters them into discovered topics using Bayesian topic models, in particular the Latent Dirichlet Allocation (LDA). It then outputs a visual, navigable map of the query results.

Models for magnetohydrodynamics of aluminium electrolysis cells

The electric current and the associated magnetic field in aluminium electrolysis cells create effects limiting the cell productivity and possibly cause instabilities: surface waving, ‘anode effects’, erosion of pot lining, feed material sedimentation, etc. The instructive analysis is presented via a step by step inclusion of different physical coupling factors affecting the magnetic field, electric current, velocity and wave development in the electrolysis cells. The full time dependent model couples the nonlinear turbulent fluid dynamics and the extended electromagnetic field in the cell, and the whole bus bar circuit with the ferromagnetic effects. Animated examples for the high amperage cells are presented. The theory and numerical model of the electrolysis cell is extended to the cases of variable cell bottom of aluminium layer and the variable thickness of the electrolyte due to the anode non-uniform burn-out process and the presence of the anode channels. The problem of the channel importance is well known Moreau-Evans model) for the stationary interface and the velocity field, and was validated against measurements in commercial cells, particularly with the recently published ‘benchmark’ test for the MHD models of aluminium cells [1]. The presence of electrolyte channels requires also to reconsider the previous magnetohydrodynamic instability theories and the dynamic wave development models. The results indicate the importance of a ‘sloshing’ parametrically excited MHD wave development in the aluminium production cells.

Transfer Across the Air-Sea Interface

The efficiency of transfer of gases and particles across the air-sea interface is controlled by several physical, biological and chemical processes in the atmosphere and water which are described here (including waves, large- and small-scale turbulence, bubbles, sea spray, rain and surface films). For a deeper understanding of relevant transport mechanisms, several models have been developed, ranging from conceptual models to numerical models. Most frequently the transfer is described by various functional dependencies of the wind speed, but more detailed descriptions need additional information. The study of gas transfer mechanisms uses a variety of experimental methods ranging from laboratory studies to carbon budgets, mass balance methods, micrometeorological techniques and thermographic techniques. Different methods resolve the transfer at different scales of time and space; this is important to take into account when comparing different results. Air-sea transfer is relevant in a wide range of applications, for example, local and regional fluxes, global models, remote sensing and computations of global inventories. The sensitivity of global models to the description of transfer velocity is limited; it is however likely that the formulations are more important when the resolution increases and other processes in models are improved. For global flux estimates using inventories or remote sensing products the accuracy of the transfer formulation as well as the accuracy of the wind field is crucial.

Absorption-based algorithm of primary production for total and size-fractionated phytoplankton in coastal waters

Most satellite models of production have been designed and calibrated for use in the open ocean. Coastal waters are optically more complex, and the use of chlorophyll a (chl a) as a first-order predictor of primary production may lead to substantial errors due to significant quantities of coloured dissolved organic matter (CDOM) and total suspended material (TSM) within the first optical depth. We demonstrate the use of phytoplankton absorption as a proxy to estimate primary production in the coastal waters of the North Sea and Western English Channel for both total, micro- and nano+pico-phytoplankton production. The method is implemented to extrapolate the absorption coefficient of phytoplankton and production at the sea surface to depth to give integrated fields of total and micro- and nano+pico-phytoplankton primary production using the peak in absorption coefficient at red wavelengths. The model is accurate to 8% in the Western English Channel and 22% in this region and the North Sea. By comparison, the accuracy of similar chl a based production models was >250%. The applicability of the method to autonomous optical sensors and remotely sensed aircraft data in both coastal and estuarine environments is discussed.

Effect of phytoplankton size classes on bio-optical properties of phytoplankton in the Western Iberian coast: Application of models

Chlorophyll-a satellite products are routinely used in oceanography, providing a synoptic and global view of phytoplankton abundance. However, these products lack information on the community structure of the phytoplankton, which is crucial for ecological modelling and ecosystem studies. To assess the usefulness of existing methods to differentiate phytoplankton functional types (PFT) or phytoplankton size classes from satellite data, in-situ phytoplankton samples collected in the Western Iberian coast, on the North-East Atlantic, were analysed for pigments and absorption spectra. Water samples were collected in five different locations, four of which were located near the shore and another in an open-ocean, seamount region. Three different modelling approaches for deriving phytoplankton size classes were applied to the in situ data. Approaches tested provide phytoplankton size class information based on the input of pigments data (Brewin et al., 2010), absorption spectra data (Ciotti et al., 2002) or both (Uitz et al., 2008). Following Uitz et al. (2008), results revealed high variability in microphytoplankton chlorophyll-specific absorption coefficients, ranging from 0.01 to 0.09 m2 (mg chl)− 1 between 400 and 500 nm. This spectral analysis suggested, in one of the regions, the existence of small cells (< 20 μm) in the fraction of phytoplankton presumed to be microphytoplankton (based on diagnostic pigments). Ciotti et al. (2002) approach yielded the highest differences between modelled and measured absorption spectra for the locations where samples had high variability in community structure and cell size. The Brewin et al. (2010) pigment-based model was adjusted and a set of model coefficients are presented and recommended for future studies in offshore water of the Western Iberian coast.

An assessment of phytoplankton primary productivity in the Arctic Ocean from satellite ocean color/in situ chlorophyll-abased models

We investigated 32 net primary productivity (NPP) models by assessing skills to reproduce integrated NPP in the Arctic Ocean. The models were provided with two sources each of surface chlorophyll-a concentration (chlorophyll), photosynthetically available radiation (PAR), sea surface temperature (SST), and mixed-layer depth (MLD). The models were most sensitive to uncertainties in surface chlorophyll, generally performing better with in situ chlorophyll than with satellite-derived values. They were much less sensitive to uncertainties in PAR, SST, and MLD, possibly due to relatively narrow ranges of input data and/or relatively little difference between input data sources. Regardless of type or complexity, most of the models were not able to fully reproduce the variability of in situ NPP, whereas some of them exhibited almost no bias (i.e., reproduced the mean of in situ NPP). The models performed relatively well in low-productivity seasons as well as in sea ice-covered/deep-water regions. Depth-resolved models correlated more with in situ NPP than other model types, but had a greater tendency to overestimate mean NPP whereas absorption-based models exhibited the lowest bias associated with weaker correlation. The models performed better when a subsurface chlorophyll-a maximum (SCM) was absent. As a group, the models overestimated mean NPP, however this was partly offset by some models underestimating NPP when a SCM was present. Our study suggests that NPP models need to be carefully tuned for the Arctic Ocean because most of the models performing relatively well were those that used Arctic-relevant parameters.

Adsorption isotherm models for basic dye adsorption by peat in single and binary component systems

Coloured effluents from textile industries are a problem in many rivers and waterways. Prediction of adsorption capacities of dyes by adsorbents is important in design considerations. The sorption of three basic dyes, namely Basic Blue 3, Basic Yellow 21 and Basic Red 22, onto peat is reported. Equilibrium sorption isotherms have been measured for the three single component systems. Equilibrium was achieved after twenty-one days. The experimental isotherm data were analysed using Langmuir, Freundlich, Redlich-Peterson, Temkin and Toth isotherm equations. A detailed error analysis has been undertaken to investigate the effect of using different error criteria for the determination of the single component isotherm parameters and hence obtain the best isotherm and isotherm parameters which describe the adsorption process. The linear transform model provided the highest R2 regression coefficient with the Redlich-Peterson model. The Redlich-Peterson model also yielded the best fit to experimental data for all three dyes using the non-linear error functions. An extended Langmuir model has been used to predict the isotherm data for the binary systems using the single component data. The correlation between theoretical and experimental data had only limited success due to competitive and interactive effects between the dyes and the dye-surface interactions.

Surface alloying and mixing at the Mn/Fe(001) interface: Real-time photoelectron spectroscopy and modified embedded atom simulations

Structural and magnetic properties of thin Mn films on the Fe(001) surface have been investigated by a combination of photoelectron spectroscopy and computer simulation in the temperature range 300 Kless than or equal toTless than or equal to750 K. Room-temperature as deposited Mn overlayers are found to be ferromagnetic up to 2.5-monolayer (ML) coverage, with a magnetic moment parallel to that of the iron substrate. The Mn atomic moment decreases with increasing coverage, and thicker samples (4-ML and 4.5-ML coverage) are antiferromagnetic. Photoemission measurements performed while the system temperature is rising at constant rate (dT/dtsimilar to0.5 K/s) detect the first signs of Mn-Fe interdiffusion at T=450 K, and reveal a broad temperature range (610 Kless than or equal toTless than or equal to680 K) in which the interface appears to be stable. Interdiffusion resumes at Tgreater than or equal to680 K. Molecular dynamics and Monte Carlo simulations allow us to attribute the stability plateau at 610 Kless than or equal toTless than or equal to680 K to the formation of a single-layer MnFe surface alloy with a 2x2 unit cell and a checkerboard distribution of Mn and Fe atoms. X-ray-absorption spectroscopy and analysis of the dichroic signal show that the alloy has a ferromagnetic spin structure, collinear with that of the substrate. The magnetic moments of Mn and Fe atoms in the alloy are estimated to be 0.8mu(B) and 1.1mu(B), respectively.

Mixed-dimensional coupling in finite element models

In many finite element analysis models it would be desirable to combine reduced- or lower-dimensional element types with higher-dimensional element types in a single model. In order to achieve compatibility of displacements and stress equilibrium at the junction or interface between the differing element types, it is important in such cases to integrate into the analysis some scheme for coupling the element types. A novel and effective scheme for establishing compatibility and equilibrium at the dimensional interface is described and its merits and capabilities are demonstrated. Copyright (C) 2000 John Wiley & Sons, Ltd.

Adsorption and photocatalysed destruction of cationic and anionic dyes on mesoporous titania films: Reactions at the air-solid interface

Cationic dyes, such as methylene blue (MB), Thionine (TH) and Basic Fuschin (BF), but not anionic dyes, such as Acid Orange 7 (AO7), Acid Blue 9 (AB9) and Acid Fuschin (AF), are readily adsorbed onto mesoporous titania films at high pH (pH 11), i.e. well above the pzc of titania (pH 6.5), due to electrostatic forces of attraction and repulsion, respectively. The same anionic dyes, but not the cationic dyes, are readily adsorbed on the same titania films at low pH (pH 3), i.e. well below titania's pzc. MB appears to adsorb on mesoporous titania films at pH 11 as the trimer (lambda(max) = 570 nm) but, upon drying, although the trimer still dominates, there is an absorption peak at 665 nm, especially notable at low [MB], which may be due to the monomer, but more likely MB J-aggregates. In contrast, the absorption spectrum of AO7 adsorbed onto the mesoporous titania film at low pH is very similar to the dye monomer. For both MB and AO7 the kinetics of adsorption are first order and yield high rate constants (3.71 and 1.481 g(-1) min(-1)), indicative of a strong adsorption process. Indeed, both MB and AO7 stained films retained much of their colour when left overnight in dye-free pH 11 and 3 solutions, respectively, indicating the strong nature of the adsorption. The kinetics of the photocatalytic bleaching of the MB-titania films at high pH are complex and not well-described by the Julson-Ollis kinetic model [A.J. Julson, D.F. Ollis, Appl. Catal. B. 65 (2006) 315]. Instead, there appears to be an initial fast but not simple demethylation step, followed by a zero-order bleaching and further demethylation steps. In contrast, the kinetics of photocatalytic bleaching of the AO7-titania film give a good fit to the Julson-Ollis kinetic model, yielding values for the various fitting parameters not too dissimilar to those reported for AO7 adsorbed on P25 titania powder. (C) 2008 Elsevier B.V. All rights reserved.

Validation of membrane protein topology models by oxidative labeling and mass spectrometry

Computer-assisted topology predictions are widely used to build low-resolution structural models of integral membrane proteins (IMPs). Experimental validation of these models by traditional methods is labor intensive and requires modifications that might alter the IMP native conformation. This work employs oxidative labeling coupled with mass spectrometry (MS) as a validation tool for computer-generated topology models. ·OH exposure introduces oxidative modifications in solvent-accessible regions, whereas buried segments (e.g., transmembrane helices) are non-oxidizable. The Escherichia coli protein WaaL (O-antigen ligase) is predicted to have 12 transmembrane helices and a large extramembrane domain (Pérez et al., Mol. Microbiol. 2008, 70, 1424). Tryptic digestion and LC-MS/MS were used to map the oxidative labeling behavior of WaaL. Met and Cys exhibit high intrinsic reactivities with ·OH, making them sensitive probes for solvent accessibility assays. Overall, the oxidation pattern of these residues is consistent with the originally proposed WaaL topology. One residue (M151), however, undergoes partial oxidation despite being predicted to reside within a transmembrane helix. Using an improved computer algorithm, a slightly modified topology model was generated that places M151 closer to the membrane interface. On the basis of the labeling data, it is concluded that the refined model more accurately reflects the actual topology of WaaL. We propose that the combination of oxidative labeling and MS represents a useful strategy for assessing the accuracy of IMP topology predictions, supplementing data obtained in traditional biochemical assays. In the future, it might be possible to incorporate oxidative labeling data directly as constraints in topology prediction algorithms.