High-pressure studies on Ge-Te glasses. Evidence for a critical composition in IV-VI chalcogenide glassy systems

The electrical resistivity of bulk GexTe100-x glasses has been measured as a function of temperature and pressure. Under high pressure, all the glasses were found to undergo sharp discontinuous transitions from glassy semiconductors to crystalline metal. Several of the observed properties such as the transition pressure, conductivity activation energy and pre-exponential factor, exhibit anomalous trends at a composition x = 20. These results suggest that the x = 20 composition in the Ge-Te system should possess salient structural features. A model based on the unusual stability of structural units is proposed for explaining the anomaly at 20 at.% Ge concentration.

Elasto-Electrostatic Analysis of Circular Microplates Used in Capacitive Micromachined Ultrasonic Transducers

The active structural component of a capacitive micromachined ultrasonic transducer (CMUT) is the top plate which vibrates under the influence of a time-varying electrostatic force thereby producing ultrasound waves of the desired frequency in the surrounding medium. Analysis of MEMS devices which rely on electrostatic actuation is complicated due to the fact that the structural deformations alter the electrostatic forces, which redistribute and modify the applied loads. Hence, it becomes imperative to consider the electrostatics-structure coupling aspect in the design of these devices. This paper presents an approximate analytical solution for the static deflection of a thin, clamped circular plate caused by electrostatic forces which are inherently nonlinear. Traditionally, finite element simulations using some commercial software such as ANSYS are employed to determine the structural deflections caused by electrostatic forces. Since the structural deformation alters the electrostatic field, a coupled-field simulation is required wherein the electrostatic mesh is continuously updated to coincide with the deflection of the structure. Such simulations are extremely time consuming, in addition to being nontransparent and somewhat hard to implement. We employ the classical thin-plate theory which is adequate when the ratio of the diameter to thickness of the plate is very large, a situation commonly prevalent in many MEMS devices, especially the CMUTs. We solve the thin-plate electrostatic-elastic equation using the Galerkin-weighted residual technique, under the assumption that the deflections are small in comparison to the thickness of the plate. The evaluation of the electrostatic force between the two plates is simplified due to the fact that the electrostatic gap is much smaller than the lateral dimensions of the device. The results obtained are compared to those found from ANSYS simulations and an excellent agreement is observed between the two. The pull-in voltage predicted by our model is close to the value predicted by ANSYS simulations.

Temperature response of trehalase from a mesophilic (Neurospora crassa) and a thermophilic (Thermomyces lanuginosus) fungus

The purified trehalases of the mesophilic fungus, Neurospora crassa, and the thermophilic fungus, Thermomyces lanuginosus, had similar temperature and pH optima for activity, but differed in molecular weight, electrophoretic mobility and Michaelis constant. At lower concentration, trehalases from both fungi were inactivated to similar extent at 60°C. While purified trehalase of T. lanuginosus was afforded protection against heat-inactivation by proteinaceous protective factor(s) present in mycelial extracts, by bovine serum albumin and by casein, these did not afford protection to N. crassa trehalase against heat inactivation. Both trehalases exhibited discontinuous Arrhenius plots with temperature of discontinuity at 40°C. The activation energy calculated from the slope of the Arrhenius plot was higher for the T. lanuginosus enzyme. The plots of apparent K m versus 1/T for trehalases of N. crassa and T. lanuginosus were linear from 30° to 60°C. The results show that purified trehalases of the mesophilic and the thermophilic fungus are distinct. Although, these exhibit similar thermostability of their catalytic function at low concentration, distinctive thermal stability characteristics of thermophilic enzyme become apparent at high protein concentration. This could be brought about in the cell by the enzyme itself, or by other proteins.

Nucleation and growth of Al2O3/metal composites by oxidation of aluminum alloys

The nucleation and growth mechanisms during high temperature oxidation of liquid Al-3% Mg and Al-3% Mg-3% Si alloys were studied with the aim of enhancing our understanding of a new composite fabrication process. The typical oxidation sequence consists of an initial event of rapid but brief oxidation, followed by an incubation period of limited oxide growth after which bulk Al2O3/Al composite forms. A duplex oxide layer, MgO (upper) and MgAl2O4 (lower), forms on the alloy surface during initial oxidation and incubation. The spinel layer remains next to the liquid alloy during bulk oxide growth and is the eventual repository for most of the magnesium in the original alloy. Metal microchannels developed during incubation continuously supply alloy through the composite to the reaction interface. During the growth process, a layered structure exists at the upper extremity of the composite, consisting of MgO at the top surface, MgAl2O4 (probably discontinuous), Al alloy, and finally the bulk Al2O3 composite containing microchannels of the alloy. The bulk oxide growth mechanism appears to involve continuous formation and dissolution of the Mg-rich oxides at the surface, diffusion of oxygen through the underlying liquid metal, and epitaxial growth of Al2O3 on the existing composite body. The roles of Mg and Si in the composite growth process are discussed.

Nonlinear Transport in Hybrid Polypyrrole-Gold Nanostructures

We have studied charge transport in nanometer scale films of polypyrrole (PPy) that were grown electrochemically onto discontinuous ultrathin films of gold. The gold films consisted of 100 nm size islands, separated from each other by nanometer-size gaps. The thickness of PPy can be varied from 30 to 200 nm. The I-V characteristics of these hybrid PPy-Au nanostructures show strong non-linearity at low temperatures, and in particular for the more insulating samples. The hopping transport is further verified from the log / versus V-1/4 plots. Furthermore, the I-V data follow an empirical relation dlog//dV(1/4) similar to T-1/2.

A Smooth Finite Element Method Based on Reproducing Kernel DMS-Splines

The element-based piecewise smooth functional approximation in the conventional finite element method (FEM) results in discontinuous first and higher order derivatives across element boundaries Despite the significant advantages of the FEM in modelling complicated geometries, a motivation in developing mesh-free methods has been the ease with which higher order globally smooth shape functions can be derived via the reproduction of polynomials There is thus a case for combining these advantages in a so-called hybrid scheme or a `smooth FEM' that, whilst retaining the popular mesh-based discretization, obtains shape functions with uniform C-p (p >= 1) continuity One such recent attempt, a NURBS based parametric bridging method (Shaw et al 2008b), uses polynomial reproducing, tensor-product non-uniform rational B-splines (NURBS) over a typical FE mesh and relies upon a (possibly piecewise) bijective geometric map between the physical domain and a rectangular (cuboidal) parametric domain The present work aims at a significant extension and improvement of this concept by replacing NURBS with DMS-splines (say, of degree n > 0) that are defined over triangles and provide Cn-1 continuity across the triangle edges This relieves the need for a geometric map that could precipitate ill-conditioning of the discretized equations Delaunay triangulation is used to discretize the physical domain and shape functions are constructed via the polynomial reproduction condition, which quite remarkably relieves the solution of its sensitive dependence on the selected knotsets Derivatives of shape functions are also constructed based on the principle of reproduction of derivatives of polynomials (Shaw and Roy 2008a) Within the present scheme, the triangles also serve as background integration cells in weak formulations thereby overcoming non-conformability issues Numerical examples involving the evaluation of derivatives of targeted functions up to the fourth order and applications of the method to a few boundary value problems of general interest in solid mechanics over (non-simply connected) bounded domains in 2D are presented towards the end of the paper

Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

Analysis of reflection characteristics of a normal incidence plane wave on resonant sound absorbers: A finite element approach

Resonant sound absorbers are used widely as anechoic coatings in underwater applications. In this paper a finite element scheme based on the Galerkin technique is used to analyze the reflection characteristics of the resonant absorber when insonified by a normal incidence plane wave. A waveguide theory coupled with an impedance matching condition in the fluid is used to model the problem. It is shown in this paper that the fluid medium encompassing the absorber can be modeled as an elastic medium with equivalent Lamé constants. Quarter symmetry conditions within the periodic unit cell are exploited. The finite element results are compared with analytical results, and with results published elsewhere in the literature. It is shown in the process that meshing of the fluid domain can be obviated if the transmission coefficients or reflection coefficients only are desired as is often the case. Finally, some design curves for thin resonant absorbers with water closure are presented in this paper.

Antigenic determinants of human serum retinol binding protein as probed with monoclonal antibodies

Monoclonal antibodies raised against human serum retinol-binding protein (hRBP) were used as probes for the study of the antigenic determinants of hRBP and those shared with the same protein from other species. The antibodies could be classified into four distinct groups and react with the homologous proteins from the rat as well as the rabbit sera. Three of these antibodies recognize sequential or continuous epitopes while the remaining antibody is directed against a discontinuous or conformational epitope. By chemical cleavage with cyanogen bromide, the domains recognized by the monoclonal antibodies could be delineated. By solid-phase synthetic approach, the core sequences recognized by two of these monoclonal antibodies were identified to amino acid sequences 45–51 and 128–131 of the primary amino acid sequence of hRBP.

Electrical-Resistivity Studies on GE-SE Glasses at High-Pressures and Low-Temperatures

High pressure electrical resistivity measurements were carried out on GexSe100-x (0 less-than-or-equal-to x less-than-or-equal-to 40) glasses at ambient and low temperatures using the Bridgman anvil system. All the melt quenched glasses show a discontinuous glassy semiconductor to crystalline metal transition at high pressures. The high pressure phases of Ge-Se samples do not correspond to any of the equilibrium phases of the system. Additionally, the variation of transition pressure (P(T)), ambient resistivity (rho0) and the activation energy (DELTAE(t)) with composition, exhibit a change in behaviour at x = 20 and 33. The unusual variations observed in these glasses are discussed in the light of chemical and percolation thresholds occurring in the glassy system.

Pulsed NMR study of molecular motions and phase transitions in [N(CH3)4]PbX3 (X=Cl, Br, I)

Proton spin—lattice relaxation time (T1) is measured in [N(CH3)4]PbX3 (X=Cl, Br, I) from 300-77 K at 9.75 MHz. All the compounds show discontinuous changes in T1 values (at 256, 270 and 277 K, respectively), indicating phase transitions. Single T1 minimum is observed in all the cases and the T1 variation is explained in terms of [N(CH3)4] and CH3 group dynamics. The activation energy Eα decreases from chloride to iodide (from 4 to 2 kcal/mol). In bromide and iodide, T1 is found to decrease with increase in temperature at higher temperatures, indicating the presence of spin—rotation interaction.

Robust/Optimal Temperature Profile Control of a High-Speed Aerospace Vehicle Using Neural Networks

An approximate dynamic programming (ADP)-based suboptimal neurocontroller to obtain desired temperature for a high-speed aerospace vehicle is synthesized in this paper. A I-D distributed parameter model of a fin is developed from basic thermal physics principles. "Snapshot" solutions of the dynamics are generated with a simple dynamic inversion-based feedback controller. Empirical basis functions are designed using the "proper orthogonal decomposition" (POD) technique and the snapshot solutions. A low-order nonlinear lumped parameter system to characterize the infinite dimensional system is obtained by carrying out a Galerkin projection. An ADP-based neurocontroller with a dual heuristic programming (DHP) formulation is obtained with a single-network-adaptive-critic (SNAC) controller for this approximate nonlinear model. Actual control in the original domain is calculated with the same POD basis functions through a reverse mapping. Further contribution of this paper includes development of an online robust neurocontroller to account for unmodeled dynamics and parametric uncertainties inherent in such a complex dynamic system. A neural network (NN) weight update rule that guarantees boundedness of the weights and relaxes the need for persistence of excitation (PE) condition is presented. Simulation studies show that in a fairly extensive but compact domain, any desired temperature profile can be achieved starting from any initial temperature profile. Therefore, the ADP and NN-based controllers appear to have the potential to become controller synthesis tools for nonlinear distributed parameter systems.

A characteristic based numerical method with tracking for nonlinear wave equations

Many physical problems can be modeled by scalar, first-order, nonlinear, hyperbolic, partial differential equations (PDEs). The solutions to these PDEs often contain shock and rarefaction waves, where the solution becomes discontinuous or has a discontinuous derivative. One can encounter difficulties using traditional finite difference methods to solve these equations. In this paper, we introduce a numerical method for solving first-order scalar wave equations. The method involves solving ordinary differential equations (ODEs) to advance the solution along the characteristics and to propagate the characteristics in time. Shocks are created when characteristics cross, and the shocks are then propagated by applying analytical jump conditions. New characteristics are inserted in spreading rarefaction fans. New characteristics are also inserted when values on adjacent characteristics lie on opposite sides of an inflection point of a nonconvex flux function, Solutions along characteristics are propagated using a standard fourth-order Runge-Kutta ODE solver. Shocks waves are kept perfectly sharp. In addition, shock locations and velocities are determined without analyzing smeared profiles or taking numerical derivatives. In order to test the numerical method, we study analytically a particular class of nonlinear hyperbolic PDEs, deriving closed form solutions for certain special initial data. We also find bounded, smooth, self-similar solutions using group theoretic methods. The numerical method is validated against these analytical results. In addition, we compare the errors in our method with those using the Lax-Wendroff method for both convex and nonconvex flux functions. Finally, we apply the method to solve a PDE with a convex flux function describing the development of a thin liquid film on a horizontally rotating disk and a PDE with a nonconvex flux function, arising in a problem concerning flow in an underground reservoir.

Optimal management of beaver population using a reduced-order distributed parameter model and single network adaptive critics

Beavers are often found to be in conflict with human interests by creating nuisances like building dams on flowing water (leading to flooding), blocking irrigation canals, cutting down timbers, etc. At the same time they contribute to raising water tables, increased vegetation, etc. Consequently, maintaining an optimal beaver population is beneficial. Because of their diffusion externality (due to migratory nature), strategies based on lumped parameter models are often ineffective. Using a distributed parameter model for beaver population that accounts for their spatial and temporal behavior, an optimal control (trapping) strategy is presented in this paper that leads to a desired distribution of the animal density in a region in the long run. The optimal control solution presented, imbeds the solution for a large number of initial conditions (i.e., it has a feedback form), which is otherwise nontrivial to obtain. The solution obtained can be used in real-time by a nonexpert in control theory since it involves only using the neural networks trained offline. Proper orthogonal decomposition-based basis function design followed by their use in a Galerkin projection has been incorporated in the solution process as a model reduction technique. Optimal solutions are obtained through a "single network adaptive critic" (SNAC) neural-network architecture.

Reduced Order Suboptimal Control Design for a Class of Nonlinear Distributed Parameter Systems Using POD and /spl thetas/-D Techniques

A new computational tool is presented in this paper for suboptimal control design of a class of nonlinear distributed parameter systems. First proper orthogonal decomposition based problem-oriented basis functions are designed, which are then used in a Galerkin projection to come up with a low-order lumped parameter approximation. Next, a suboptimal controller is designed using the emerging /spl thetas/-D technique for lumped parameter systems. This time domain sub-optimal control solution is then mapped back to the distributed domain using the same basis functions, which essentially leads to a closed form solution for the controller in a state feedback form. Numerical results for a real-life nonlinear temperature control problem indicate that the proposed method holds promise as a good suboptimal control design technique for distributed parameter systems.