Stability of pargasite during ultrahigh-temperature metamorphism:A consequence of titanium and REE partitioning

Orthopyroxene-clinopyroxene-plagioclase needles and symplectite along the cleavage planes and grain boundaries of fluorine-bearing titanian-ferroan pargasite from the Highland Complex, Sri Lanka, are interpreted as evidence for dehydration melting at ultrahigh-temperature conditions. High Ti (up to 0.4 pfu) and F (XF up to 0.56) content in pargasite extends its stability to higher temperatures, and the composition indicates the dehydration melting reaction may take place at ultrahigh-temperatures (~950 °C) at a pressure around 10 kbar, close to peak metamorphic conditions. The increase of Ti content close to the grain boundaries and cleavage planes in pargasite indicates titanium partitioning from the melt during dehydration melting enhanced the stability of the mineral toward ultrahigh-temperature conditions. The REE content in the pargasite shows a similar behavior to that of titanium. The cores with no breakdown assemblage consist of low and flat REE concentrations with respect to the high and Eu-depleted rim. Clinopyroxene in symplectite and needle-shaped lamellae within the pargasite porphyroblasts have similar REE patterns with slightly low-concentrations relative to that of pargasite. In the breakdown assemblage, LREEs are partitioned mainly into plagioclase while the HREEs are partitioned into orthopyroxene. The REE enrichment in the pargasite rims signals their relative partitioning between pargasite rims and melt. Modeling of the partitioning of Ti and REEs associated with pargasite breakdown demonstrates that its stability is greatly enhanced at UHT conditions. This investigation implies that the stability of hydrous minerals such as amphibole can be extended to UHT conditions, and expands our knowledge of metamorphism in the lower crust.

A Novel Operator-Splitting Technique for One-Dimensional Laminar Premixed Flames

This paper is concerned the calculation of flame structure of one-dimensional laminar premixed flames using the technique of operator-splitting. The technique utilizes an explicit method of solution with one step Euler for chemistry and a novel probabilistic scheme for diffusion. The relationship between diffusion phenomenon and Gauss-Markoff process is exploited to obtain an unconditionally stable explicit difference scheme for diffusion. The method has been applied to (a) a model problem, (b) hydrazine decomposition, (c) a hydrogen-oxygen system with 28 reactions with constant Dρ 2 approximation, and (d) a hydrogen-oxygen system (28 reactions) with trace diffusion approximation. Certain interesting aspects of behaviour of the solution with non-unity Lewis number are brought out in the case of hydrazine flame. The results of computation in the most complex case are shown to compare very favourably with those of Warnatz, both in terms of accuracy of results as well as computational time, thus showing that explicit methods can be effective in flame computations. Also computations using the Gear-Hindmarsh for chemistry and the present approach for diffusion have been carried out and comparison of the two methods is presented.

Structural Ordering of Disordered Ligand-Binding Loops of Biotin Protein Ligase into Active Conformations as a Consequence of Dehydration

Mycobacterium tuberculosis (Mtb), a dreaded pathogen, has a unique cell envelope composed of high fatty acid content that plays a crucial role in its pathogenesis. Acetyl Coenzyme A Carboxylase (ACC), an important enzyme that catalyzes the first reaction of fatty acid biosynthesis, is biotinylated by biotin acetyl-CoA carboxylase ligase (BirA). The ligand-binding loops in all known apo BirAs to date are disordered and attain an ordered structure only after undergoing a conformational change upon ligand-binding. Here, we report that dehydration of Mtb-BirA crystals traps both the apo and active conformations in its asymmetric unit, and for the first time provides structural evidence of such transformation. Recombinant Mtb-BirA was crystallized at room temperature, and diffraction data was collected at 295 K as well as at 120 K. Transfer of crystals to paraffin and paratone-N oil (cryoprotectants) prior to flash-freezing induced lattice shrinkage and enhancement in the resolution of the X-ray diffraction data. Intriguingly, the crystal lattice rearrangement due to shrinkage in the dehydrated Mtb-BirA crystals ensued structural order of otherwise flexible ligand-binding loops L4 and L8 in apo BirA. In addition, crystal dehydration resulted in a shift of similar to 3.5 angstrom in the flexible loop L6, a proline-rich loop unique to Mtb complex as well as around the L11 region. The shift in loop L11 in the C-terminal domain on dehydration emulates the action responsible for the complex formation with its protein ligand biotin carboxyl carrier protein (BCCP) domain of ACCA3. This is contrary to the involvement of loop L14 observed in Pyrococcus horikoshii BirA-BCCP complex. Another interesting feature that emerges from this dehydrated structure is that the two subunits A and B, though related by a noncrystallographic twofold symmetry, assemble into an asymmetric dimer representing the ligand-bound and ligand-free states of the protein, respectively. In-depth analyses of the sequence and the structure also provide answers to the reported lower affinities of Mtb-BirA toward ATP and biotin substrates. This dehydrated crystal structure not only provides key leads to the understanding of the structure/function relationships in the protein in the absence of any ligand-bound structure, but also demonstrates the merit of dehydration of crystals as an inimitable technique to have a glance at proteins in action.

A Wavelet based Teager Energy Operator for Spike Detection in Microelectrode Array Recordings

Spike detection in neural recordings is the initial step in the creation of brain machine interfaces. The Teager energy operator (TEO) treats a spike as an increase in the `local' energy and detects this increase. The performance of TEO in detecting action potential spikes suffers due to its sensitivity to the frequency of spikes in the presence of noise which is present in microelectrode array (MEA) recordings. The multiresolution TEO (mTEO) method overcomes this shortcoming of the TEO by tuning the parameter k to an optimal value m so as to match to frequency of the spike. In this paper, we present an algorithm for the mTEO using the multiresolution structure of wavelets along with inbuilt lowpass filtering of the subband signals. The algorithm is efficient and can be implemented for real-time processing of neural signals for spike detection. The performance of the algorithm is tested on a simulated neural signal with 10 spike templates obtained from [14]. The background noise is modeled as a colored Gaussian random process. Using the noise standard deviation and autocorrelation functions obtained from recorded data, background noise was simulated by an autoregressive (AR(5)) filter. The simulations show a spike detection accuracy of 90%and above with less than 5% false positives at an SNR of 2.35 dB as compared to 80% accuracy and 10% false positives reported [6] on simulated neural signals.

Time Dependent Operator-split and Unsplit Schemes for One Dimensional Premixed Flames

This paper is concerned with a study of an operator split scheme and unsplit scheme for the computation of adiabatic freely propagating one-dimensional premixed flames. The study uses unsteady method for both split and unsplit schemes employing implicit chemistry and explicit diffusion, a combination which is stable and convergent. Solution scheme is not sensitive to the initial starting estimate and provides steady state even with straight line profiles (far from steady state) in small number of time steps. Two systems H2-Air and H2-NO (involving complex nitrogen chemistry) are considered in presentinvestigation. Careful comparison shows that the operator split approach is slightly superior than the unsplit when chemistry becomes complex. Comparison of computational times with those of existing steady and unsteady methods seems to suggest that the method employing implicit-explicit algorithm is very efficient and robust.

Negative differential resistance in Gd0.5Sr0.5MnO3: A consequence of Joule heating

Negative differential resistance (NDR) in current-voltage (I-V) characteristics and apparent colossal electroresistance were observed in Gd0.5Sr0.5MnO3 single crystals at low temperatures. The continuous dc I-V measurements showed a marked thermal drift. In addition, temperature of the sample surface was found to be significantly higher than that of the base at high applied currents. Two different strategies namely estimation and diminution of the Joule heating (pulsed I-V measurements) were employed to investigate its role in the electric transport properties. Our experiments reveal that the NDR in Gd0.5Sr0.5MnO3 is a consequence of Joule heating rather than the melting of charge order. (C) 2010 American Institute of Physics. doi:10.1063/1.3486221]

Squeezed states, metaplectic group, and operator möbius transformations

We present an analysis, based on the metaplectic group Mp(2), of the recently introduced single-mode inverse creation and annihilation operators and of the associated eigenstates of different two-photon annihilation operators. We motivate and obtain a quantum operator form of the classical Mobius or fractional linear transformation. The subtle relation to the two unitary irreducible representations of Mp(2) is brought out. For problems involving inverse operators the usefulness of the Bargmann analytic function representation of quantum mechanics is demonstrated. Squeezing, bunching, and photon-number distributions of the four families of states that arise in this context are studied both analytically and numerically

Squeezed vacuum as an eigenstate of two-photon annihilation operator

We introduce the inverse annihilation and creation operators a-1 and a(dagger-1) by their actions on the number states. We show that the squeezed vacuum exp(1/2xia(dagger2)]\0] and squeezed first number state exp[1.2xia(dagger2)]\n = 1] are respectively the eigenstates of the operators (a(dagger-1)a) and (aa(dagger-1)) with the eigenvalue xi.

Polymorphism in Opto-Electronic Materials with a Benzothiazole-fluorene Core: A Consequence of High Conformational Flexibility of pi-Conjugated Backbone and Alkyl Side Chains

Fluorene and its derivatives are well-known organic semiconducting materials in the field of opto-electronic devices because of their charge transport properties. Three new organic semiconducting materials, namely, 2,2'-((9,9-butyl-9H-fluorene-2,7-diyl)bis(4,1 phenylene))bisbenzod]thiazole, C4; 2,2'-((octyl-9H-fluorene-2,7-diyl)bis(4,1 phenylene))bisbenzod]thiazole, C8; and 2,2'-((9,9-dodecayl-9H-fluorene-2,7-diyl)bis(4,1 phenylene))bisbenzod]thiazole, C12 with a benzothiazole-fluorene backbone, were synthesized and characterized for their photophysical properties. A phenomenon of concomitant polymorphism has been investigated in the first two derivatives (C4 and C8) and has been analyzed systematically in terms of the packing characteristics involving pi ... pi interactions. The conformational flexibility of the pi-conjugated 2,2'-(fluorene-2,7-diyl)bis(4,1 phenylene)bisbenzod]thiazole backbone coupled with orientational freedom of the terminal alkyl chains were found to be the key factors responsible for these polymorphic modifications. Attempts to grow suitable crystals for single crystal X-ray diffraction of compound C12 were unsuccessful.

The extended Enskog operator for simple fluids with continuous potentials: single particle and collective properties

We generalized the Enskog theory originally developed for the hard-sphere fluid to fluids with continuous potentials, such as the Lennard–Jones. We derived the expression for the k and ω dependent transport coefficient matrix which enables us to calculate the transport coefficients for arbitrary length and time scales. Our results reduce to the conventional Chapman–Enskog expression in the low density limit and to the conventional k dependent Enskog theory in the hard-sphere limit. As examples, the self-diffusion of a single atom, the vibrational energy relaxation, and the activated barrier crossing dynamics problem are discussed.

Wear of metals: a consequence of stable/unstable material response

Wear of metals in dry sliding is dictated by the material response to traction. This is demonstrated by considering the wear of aluminium and titanium alloys. In a regime of stable homogeneous deformation the material approaching the surface from the bulk passes through microprocessing zones of flow, fracture, comminution and compaction to generate a protective tribofilm that retains the interaction in the mild wear regime. If the response leads to microstructural instabilities such as adiabatic shear bands, the near-surface zone consists of stacks of 500 nm layers situated parallel to the sliding direction. Microcracks are generated below the surface to propagate normally away from the surface though microvoids situated in the layers, until it reaches a depth of 10-20 mum. A rectangular laminate debris consisting of a 20-40 layer stack is produced, The wear in this mode is severe.

An operator-splitting finite element method for the efficient parallel solution of multidimensional population balance systems

A finite element method for solving multidimensional population balance systems is proposed where the balance of fluid velocity, temperature and solute partial density is considered as a two-dimensional system and the balance of particle size distribution as a three-dimensional one. The method is based on a dimensional splitting into physical space and internal property variables. In addition, the operator splitting allows to decouple the equations for temperature, solute partial density and particle size distribution. Further, a nodal point based parallel finite element algorithm for multi-dimensional population balance systems is presented. The method is applied to study a crystallization process assuming, for simplicity, a size independent growth rate and neglecting agglomeration and breakage of particles. Simulations for different wall temperatures are performed to show the effect of cooling on the crystal growth. Although the method is described in detail only for the case of d=2 space and s=1 internal property variables it has the potential to be extendable to d+s variables, d=2, 3 and s >= 1. (C) 2011 Elsevier Ltd. All rights reserved.

Dilations of Gamma-contractions by solving operator equations

For a contraction P and a bounded commutant S of P. we seek a solution X of the operator equation S - S*P = (1 - P* P)(1/2) X (1 - P* P)(1/2) where X is a bounded operator on (Ran) over bar (1 - P* P)(1/2) with numerical radius of X being not greater than 1. A pair of bounded operators (S, P) which has the domain Gamma = {(z(1) + z(2), z(2)): vertical bar z(1)vertical bar < 1, vertical bar z(2)vertical bar <= 1} subset of C-2 as a spectral set, is called a P-contraction in the literature. We show the existence and uniqueness of solution to the operator equation above for a Gamma-contraction (S, P). This allows us to construct an explicit Gamma-isometric dilation of a Gamma-contraction (S, P). We prove the other way too, i.e., for a commuting pair (S, P) with parallel to P parallel to <= 1 and the spectral radius of S being not greater than 2, the existence of a solution to the above equation implies that (S, P) is a Gamma-contraction. We show that for a pure F-contraction (S, P), there is a bounded operator C with numerical radius not greater than 1, such that S = C + C* P. Any Gamma-isometry can be written in this form where P now is an isometry commuting with C and C. Any Gamma-unitary is of this form as well with P and C being commuting unitaries. Examples of Gamma-contractions on reproducing kernel Hilbert spaces and their Gamma-isometric dilations are discussed. (C) 2012 Elsevier Inc. All rights reserved.

An operator-splitting Galerkin/SUPG finite element method for population balance equations : stability and convergence

We present a heterogeneous finite element method for the solution of a high-dimensional population balance equation, which depends both the physical and the internal property coordinates. The proposed scheme tackles the two main difficulties in the finite element solution of population balance equation: (i) spatial discretization with the standard finite elements, when the dimension of the equation is more than three, (ii) spurious oscillations in the solution induced by standard Galerkin approximation due to pure advection in the internal property coordinates. The key idea is to split the high-dimensional population balance equation into two low-dimensional equations, and discretize the low-dimensional equations separately. In the proposed splitting scheme, the shape of the physical domain can be arbitrary, and different discretizations can be applied to the low-dimensional equations. In particular, we discretize the physical and internal spaces with the standard Galerkin and Streamline Upwind Petrov Galerkin (SUPG) finite elements, respectively. The stability and error estimates of the Galerkin/SUPG finite element discretization of the population balance equation are derived. It is shown that a slightly more regularity, i.e. the mixed partial derivatives of the solution has to be bounded, is necessary for the optimal order of convergence. Numerical results are presented to support the analysis.

Operator-splitting finite element algorithms for computations of high-dimensional parabolic problems

An operator-splitting finite element method for solving high-dimensional parabolic equations is presented. The stability and the error estimates are derived for the proposed numerical scheme. Furthermore, two variants of fully-practical operator-splitting finite element algorithms based on the quadrature points and the nodal points, respectively, are presented. Both the quadrature and the nodal point based operator-splitting algorithms are validated using a three-dimensional (3D) test problem. The numerical results obtained with the full 3D computations and the operator-split 2D + 1D computations are found to be in a good agreement with the analytical solution. Further, the optimal order of convergence is obtained in both variants of the operator-splitting algorithms. (C) 2012 Elsevier Inc. All rights reserved.