Computational modeling of optical trapping of vaterite microspheres

The ability to grow microscopic spherical birefringent crystals of vaterite, a calcium carbonate mineral, has allowed the development of an optical microrheometer based on optical tweezers. However, since these crystals are birefringent, and worse, are expected to have non-uniform birefringence, computational modeling of the microrheometer is a highly challenging task. Modeling the microrheometer - and optical tweezers in general - typically requires large numbers of repeated calculations for the same trapped particle. This places strong demands on the efficiency of computational methods used. While our usual method of choice for computational modelling of optical tweezers - the T-matrix method - meets this requirement of efficiency, it is restricted to homogeneous isotropic particles. General methods that can model complex structures such as the vaterite particles, such as finite-difference time-domain (FDTD) or finite-difference frequency-domain (FDFD) methods, are inefficient. Therefore, we have developed a hybrid FDFD/T-matrix method that combines the generality of volume-discretisation methods such as FDFD with the efficiency of the T-matrix method. We have used this hybrid method to calculate optical forces and torques on model vaterite spheres in optical traps. We present and compare the results of computational modelling and experimental measurements.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.

Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Weibull model selection for reliability modelling

A large number of models have been derived from the two-parameter Weibull distribution and are referred to as Weibull models. They exhibit a wide range of shapes for the density and hazard functions, which makes them suitable for modelling complex failure data sets. The WPP and IWPP plot allows one to determine in a systematic manner if one or more of these models are suitable for modelling a given data set. This paper deals with this topic.

Synthesis of [Cu-2(Me-2[9]aneN(2)S)(2)(OH)(2)] (PF6)(2), a hydroxo-bridged copper(II) complex of N,N '-dimethyl-1-thia-4,7-diazacyclononane (Me-2[9]aneN(2)S). X-ray structural analysis, magnetic susceptibility, EPR and visible spectroscopy

The bis(mu-hydroxo) complex [Cu-2(Me-2[9]aneN(2)S)(2)(OH)(2)](PF6)(2) (Me-2[9]aneN(2)S = N,N'-dimethyl-1-thia-4,7-diazacyclononane) results after reaction of [Cu(Me-2[9]aneN(2)S)(MeCN)] (PF6) with dioxygen at -78 degrees C in acetonitrile. The complex has been characterized by X-ray crystallography: orthorhombic, space group Pnma, with a 18.710(3), b 16.758(2), c 9.593(2) Angstrom, and Z = 4. The structure refined to a final R value of 0.051. The complex contains two copper(II) ions bridged by two hydroxo groups with Cu ... Cu 2.866(1) Angstrom. The solid-state magnetic susceptibility study reveals ferromagnetic coupling, the fitting parameters being J = +46+/-5 cm(-1), g = 2.01+/-0.01 and theta = -0.58+/-0.03 K. The frozen-solution e.p.r. spectrum in dimethyl sulfoxide is characteristic of a monomeric copper(II) ion (g(parallel to) 2.300, g(perpendicular to) 2.063; A(parallel to) 156.2 x 10(-4) cm(-1), A(perpendicular to) 9.0 x 10(-4) cm(-1)) with an N2O2 donor set. Thioether coordination to the copper(II) in solution is supported by the presence of an intense absorption assigned to a sigma(S)-->Cu-II LMCT transition at c. 34000 cm(-1). The single-crystal spectrum of [Cu-2(Me-2[9]aneN(2)S)(2)(OH)(2)] (PF6)(2) (273 K) reveals d-->d transitions at 14500 and 18300 cm(-1) and a weak pi(S)-->Cu-II charge-transfer band at approximately 25000 cm(-1).

Reorganization of terminator DNA upon binding replication terminator protein: Implications for the functional replication fork arrest complex

Termination of DNA replication in Bacillus subtilis involves the polar arrest of replication forks by a specific complex formed between the replication terminator protein (RTP) and DNA terminator sites. While determination of the crystal structure of RTP has facilitated our understanding of how a single RTP dimer interacts with terminator DNA, additional information is required in order to understand the assembly of a functional fork arrest complex, which requires an interaction between two RTP dimers and the terminator site. In this study, we show that the conformation of the major B. subtilis DNA terminator, Terl, becomes considerably distorted upon binding RTP. Binding of the first dimer of RTP to the B site of Terl causes the DNA to become slightly unwound and bent by similar to 40 degrees. Binding of a second dimer of RTP to the A site causes the bend angle to increase to similar to 60 degrees. We have used this new data to construct two plausible models that might explain how the ternary terminator complex can block DNA replication in a polar manner, in the first model, polarity of action is a consequence of the two RTP-DNA half-sites having different conformations. These different conformations result from different RTP-DNA contacts at each half-site (due to the intrinsic asymmetry at the terminator DNA), as well as interactions (direct or indirect) between the RTP dimers on the DNA. In the second model, polar fork arrest activity is a consequence of the different affinities of RTP for the A and B sites of the terminator DNA, modulated significantly by direct or indirect interactions between the RTP dimers.

A peptide-binding motif for I-A(g7), the class II major histocompatibility complex (MHC) molecule of NOD and Biozzi AB/H mice

The class II major histocompatibility complex molecule I-A(g7) is strongly linked to the development of spontaneous insulin-dependent diabetes mellitus (IDDM) in non obese diabetic mice and to the induction of experimental allergic encephalomyelitis in Biozzi AB/H mice. Structurally, it resembles the HLA-DQ molecules associated with human IDDM, in having a non-Asp residue at position 57 in its beta chain. To identify the requirements for peptide binding to I-A(g7) and thereby potentially pathogenic T cell epitopes, we analyzed a known I-A(g7)-restricted T cell epitope, hen egg white lysozyme (HEL) amino acids 9-27. NH2- and COOH-terminal truncations demonstrated that the minimal epitope for activation of the T cell hybridoma 2D12.1 was M12-R21 and the minimum sequence for direct binding to purified I-A(g7) M12-Y20/K13-R21. Alanine (A) scanning revealed two primary anchors for binding at relative positions (p) 6 (L) and 9 (Y) in the HEL epitope. The critical role of both anchors was demonstrated by incorporating L and Y in poly(A) backbones at the same relative positions as in the HEL epitope. Well-tolerated, weakly tolerated, and nontolerated residues were identified by analyzing the binding of peptides containing multiple substitutions at individual positions. Optimally, p6 was a large, hydrophobic residue (L, I, V, M), whereas p9 was aromatic and hydrophobic (Y or F) or positively charged (K, R). Specific residues were not tolerated at these and some other positions. A motif for binding to I-A(g7) deduced from analysis of the model HEL epitope was present in 27/30 (90%) of peptides reported to be I-A(g7)-restricted T cell epitopes or eluted from I-A(g7). Scanning a set of overlapping peptides encompassing human proinsulin revealed the motif in 6/6 good binders (sensitivity = 100%) and 4/13 weak or non-binders (specificity = 70%). This motif should facilitate identification of autoantigenic epitopes relevant to the pathogenesis and immunotherapy of IDDM.

Analytical solution and numerical model for the interface in a stratified long wave system

For a two layered long wave propagation, linearized governing equations, which were derived earlier from the Euler equations of mass and momentum assuming negligible friction and interfacial mixing are solved analytically using Fourier transform. For the solution, variations of upper layer water level is assumed to be sinosoidal having known amplitude and variations of interface level is solved. As the governing equations are too complex to solve it analytically, density of upper layer fluid is assumed as very close to the density of lower layer fluid to simplify the lower layer equation. A numerical model is developed using the staggered leap-forg scheme for computation of water level and discharge in one dimensional propagation having known amplitude for the variations of upper layer water level and interface level to be solved. For the numerical model, water levels (upper layer and interface) at both the boundaries are assumed to be known from analytical solution. Results of numerical model are verified by comparing with the analytical solutions for different time period. Good agreements between analytical solution and numerical model are found for the stated boundary condition. The reliability of the developed numerical model is discussed, using it for different a (ratio of density of fluid in the upper layer to that in the lower layer) and p (ratio of water depth in the lower layer to that in the upper layer) values. It is found that as ‘CX’ increases amplification of interface also increases for same upper layer amplitude. Again for a constant lower layer depth, as ‘p’ increases amplification of interface. also increases for same upper layer amplitude.

Ligand-field analysis of an Er(III) complex with a heptadentate tripodal N4O3 ligand

Polarized absorption and emission spectra of trigonal single crystals of an Er(III) complex coordinated to a heptadentate tripodal ligand are reported at temperatures between 8 and 298 K. The assigned energy levels below the onset of ligand absorption (< 25 000 cm(-1)) are fitted to a parametrized electronic Hamiltonian. The C-3 site symmetry of the Er(HI) ion requires eight parameters for a full description of the ligand field within a one-electron operator description. This compound shows unusually large splittings of the multiplets, and the fitted parameters imply that this heptadentate ligand imparts the largest ligand field reported for an Er(III) complex. The ligand field was also interpreted within the angular overlap model (AOM). We derive the AOM matrix to include both sigma and anisotropic pi bonding and show that a useful description of the C-3 ligand field can be made using only five parameters. The success of the AOM description is encouraging for applications on isomorphous complexes within the lanthanide series and in describing the ligand field of low-symmetry complexes with less parameters than in the usual spherical harmonic expansion.

A formal model of real-time program compilation

Program compilation can be formally defined as a sequence of equivalence-preserving transformations, or refinements, from high-level language programs to assembler code, Recent models also incorporate timing properties, but the resulting formalisms are intimidatingly complex. Here we take advantage of a new, simple model of real-time refinement, based on predicate transformer semantics, to present a straightforward compilation formalism that incorporates real-time constraints. (C) 2002 Elsevier Science B.V. All rights reserved.