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.

A MATLAB Code for Three Dimensional Linear Elastostatics using Constant Boundary Elements

Present work presents a code written in the very simple programming language MATLAB, for three dimensional linear elastostatics, using constant boundary elements. The code, in full or in part, is not a translation or a copy of any of the existing codes. Present paper explains how the code is written, and lists all the formulae used. Code is verified by using the code to solve a simple problem which has the well known approximate analytical solution. Of course, present work does not make any contribution to research on boundary elements, in terms of theory. But the work is justified by the fact that, to the best of author’s knowledge, as of now, one cannot find an open access MATLAB code for three dimensional linear elastostatics using constant boundary elements. Author hopes this paper to be of help to beginners who wish to understand how a simple but complete boundary element code works, so that they can build upon and modify the present open access code to solve complex engineering problems quickly and easily. The code is available online for open access (as supplementary file for the present paper), and may be downloaded from the website for the present journal.

Inter-domain interactions influence the stability and catalytic activity of the bi-domain protein tyrosine phosphatase PTP99A

The two protein tyrosine phosphatase (PTP) domains in bi-domain PTPs share high sequence and structural similarity. However, only one of the two PIP domains is catalytically active. Here we describe biochemical studies on the two tandem PTP domains of the bi-domain PTP, PTP99A. Phosphatase activity, monitored using small molecule as well as peptide substrates, revealed that the inactive (D2) domain activates the catalytic (D1) domain. Thermodynamic measurements suggest that the inactive D2 domain stabilizes the bi-domain (D1-D2) protein. The mechanism by which the D2 domain activates and stabilizes the bi-domain protein is governed by few interactions at the inter-domain interface. In particular, mutating Lys990 at the interface attenuates inter-domain communication. This residue is located at a structurally equivalent location to the so-called allosteric site of the canonical single domain PIP, PTP1B. These observations suggest functional optimization in bi-domain PTPs whereby the inactive PTP domain modulates the catalytic activity of the bi-domain enzyme. (C) 2012 Elsevier B.V. All rights reserved.

Enhanced J-protein interaction and compromised protein stability of mtHsp70 variants lead to mitochondrial dysfunction in Parkinsons disease

Parkinsons disease (PD) is the second most prevalent progressive neurological disorder commonly associated with impaired mitochondrial function in dopaminergic neurons. Although familial PD is multifactorial in nature, a recent genetic screen involving PD patients identified two mitochondrial Hsp70 variants (P509S and R126W) that are suggested in PD pathogenesis. However, molecular mechanisms underlying how mtHsp70 PD variants are centrally involved in PD progression is totally elusive. In this article, we provide mechanistic insights into the mitochondrial dysfunction associated with human mtHsp70 PD variants. Biochemically, the R126W variant showed severely compromised protein stability and was found highly susceptible to aggregation at physiological conditions. Strikingly, on the other hand, the P509S variant exhibits significantly enhanced interaction with J-protein cochaperones involved in folding and import machinery, thus altering the overall regulation of chaperone-mediated folding cycle and protein homeostasis. To assess the impact of mtHsp70 PD mutations at the cellular level, we developed yeast as a model system by making analogous mutations in Ssc1 ortholog. Interestingly, PD mutations in yeast (R103W and P486S) exhibit multiple in vivo phenotypes, which are associated with omitochondrial dysfunction', including compromised growth, impairment in protein translocation, reduced functional mitochondrial mass, mitochondrial DNA loss, respiratory incompetency and increased susceptibility to oxidative stress. In addition to that, R103W protein is prone to aggregate in vivo due to reduced stability, whereas P486S showed enhanced interaction with J-proteins, thus remarkably recapitulating the cellular defects that are observed in human PD variants. Taken together, our findings provide evidence in favor of direct involvement of mtHsp70 as a susceptibility factor in PD.

Evolution of quantum discord and its stability in two-qubit NMR systems

We investigate evolution of quantum correlations in ensembles of two-qubit nuclear spin systems via nuclear magnetic resonance techniques. We use discord as a measure of quantum correlations and the Werner state as an explicit example. We, first, introduce different ways of measuring discord and geometric discord in two-qubit systems and then describe the following experimental studies: (a) We quantitatively measure discord for Werner-like states prepared using an entangling pulse sequence. An initial thermal state with zero discord is gradually and periodically transformed into a mixed state with maximum discord. The experimental and simulated behavior of rise and fall of discord agree fairly well. (b) We examine the efficiency of dynamical decoupling sequences in preserving quantum correlations. In our experimental setup, the dynamical decoupling sequences preserved the traceless parts of the density matrices at high fidelity. But they could not maintain the purity of the quantum states and so were unable to keep the discord from decaying. (c) We observe the evolution of discord for a singlet-triplet mixed state during a radio-frequency spin-lock. A simple relaxation model describes the evolution of discord, and the accompanying evolution of fidelity of the long-lived singlet state, reasonably well.

Synchrotron Small-Angle X-ray Scattering Studies of Hemoglobin Nonaggregation Confined inside Polymer Capsules

The effect of confinement on the structure of hemoglobin (Hb) within polymer capsules was investigated here. Hemoglobin transformed from an aggregated state in solution to a nonaggregated state when confined inside the polymer capsules. This was directly confirmed using synchrotron small-angle X-ray scattering (SAXS) studies. The radius of gyration (R-g) and polydispersity (p) of the proteins in the confined state were smaller compared to those in solution. In fact, the R-g value is very similar to theoretical values obtained using protein structures generated from the Protein Databank. In the temperature range (25-85 degrees C, Tm 59 degrees C), the R-g values for the confined Hb remained constant. This observation is in contrary to the increasing R-g values obtained for the bare Hb in solution. This suggested higher thermal stability of Hb when confined inside the polymer capsule than when in solution. Changes in protein configuration were also reflected in the protein function. Confinement resulted in a beneficial enhancement of the electroactivity of Hb. While Hb in solution showed dominance of the cathodic process (Fe3+ -> Fe2+), efficient reversible Fe3+/Fe2+ redox response is observed in the case of the confined Hb. This has important protein functional implications. Confinement allows the electroactive heme to take up positions favorable for various biochemical activities such as sensing of analytes of various sizes from small to macromolecules and controlled delivery of drugs.

Monsoon sensitivity to aerosol direct radiative forcing in the community atmosphere model

Aerosol forcing remains a dominant uncertainty in climate studies. The impact of aerosol direct radiative forcing on Indian monsoon is extremely complex and is strongly dependent on the model, aerosol distribution and characteristics specified in the model, modelling strategy employed as well as on spatial and temporal scales. The present study investigates (i) the aerosol direct radiative forcing impact on mean Indian summer monsoon when a combination of quasi-realistic mean annual cycles of scattering and absorbing aerosols derived from an aerosol transport model constrained with satellite observed Aerosol Optical Depth (AOD) is prescribed, (ii) the dominant feedback mechanism behind the simulated impact of all-aerosol direct radiative forcing on monsoon and (iii) the relative impacts of absorbing and scattering aerosols on mean Indian summer monsoon. We have used CAM3, an atmospheric GCM (AGCM) that has a comprehensive treatment of the aerosol-radiation interaction. This AGCM has been used to perform climate simulations with three different representations of aerosol direct radiative forcing due to the total, scattering aerosols and black carbon aerosols. We have also conducted experiments without any aerosol forcing. Aerosol direct impact due to scattering aerosols causes significant reduction in summer monsoon precipitation over India with a tendency for southward shift of Tropical Convergence Zones (TCZs) over the Indian region. Aerosol forcing reduces surface solar absorption over the primary rainbelt region of India and reduces the surface and lower tropospheric temperatures. Concurrent warming of the lower atmosphere over the warm oceanic region in the south reduces the land-ocean temperature contrast and weakens the monsoon overturning circulation and the advection of moisture into the landmass. This increases atmospheric convective stability, and decreases convection, clouds, precipitation and associated latent heat release. Our analysis reveals a defining negative moisture-advection feedback that acts as an internal damping mechanism spinning down the regional hydrological cycle and leading to significant circulation changes in response to external radiative forcing perturbations. When total aerosol loading (both absorbing and scattering aerosols) is prescribed, dust and black carbon aerosols are found to cause significant atmospheric heating over the monsoon region but the aerosol-induced weakening of meridional lower tropospheric temperature gradient (leading to weaker summer monsoon rainfall) more than offsets the increase in summer-time rainfall resulting from the atmospheric heating effect of absorbing aerosols, leading to a net decrease of summer monsoon rainfall. Further, we have carried out climate simulations with globally constant AODs and also with the constant AODs over the extended Indian region replaced by realistic AODs. Regional aerosol radiative forcing perturbations over the Indian region is found to have impact not only over the region of loading but over remote tropical regions as well. This warrants the need to prescribe realistic aerosol properties in strategic regions such as India in order to accurately assess the aerosol impact.

Expanding universe: thermodynamical aspects from different models

The pivotal point of the paper is to discuss the behavior of temperature, pressure, energy density as a function of volume along with determination of caloric EoS from following two model: w(z)=w (0)+w (1)ln(1+z) & . The time scale of instability for this two models is discussed. In the paper we then generalize our result and arrive at general expression for energy density irrespective of the model. The thermodynamical stability for both of the model and the general case is discussed from this viewpoint. We also arrive at a condition on the limiting behavior of thermodynamic parameter to validate the third law of thermodynamics and interpret the general mathematical expression of integration constant U (0) (what we get while integrating energy conservation equation) physically relating it to number of micro states. The constraint on the allowed values of the parameters of the models is discussed which ascertains stability of universe. The validity of thermodynamical laws within apparent and event horizon is discussed.

Thermo-mechanical stability of a cellular assembly of carbon nanotubes in air

Carbon nanotubes (CNT) in bulk form offer outstanding structural and functional properties, and are shown to remain viscoelastic over a wide temperature range (77-1273 K) under inert conditions. We examine the quasi-static and dynamic compressive mechanical response of these cellular CNT materials in ambient air up to a temperature of 773 K. In uniaxial quasi-static compression, several displacement bursts are noted at large strains. These are results of the slippage and zipping of the CNT, and lead to significant mechanical energy absorption. Results of the dynamic mechanical analysis experiments show no degradation in storage modulus and loss coefficient for up to 20 h at 673 K. Hence, these stable cellular CNT structures can be utilized up to a maximum temperature of 673 K in air, which is much higher than the best polymers. (C) 2012 Elsevier Ltd. All rights reserved.

Seismic stability of a long unsupported circular tunnel

The stability of a long unsupported circular tunnel (opening) in a cohesive frictional soil has been assessed with the inclusion of pseudo-static horizontal earthquake body forces. The analysis has been performed under plane strain conditions by using upper bound finite element limit analysis in combination with a linear optimization procedure. The results have been presented in the form of a non-dimensional stability number (gamma H-max/c); where H = tunnel cover, c refers to soil cohesion and gamma(max) is the maximum unit weight of soil mass which the tunnel can support without collapse. The results have been obtained for various values of H/D (D = diameter of the tunnel), internal friction angle (phi) of soil, and the horizontal earthquake acceleration coefficient (alpha(h)). The computations reveal that the values of the stability numbers (i) decrease quite significantly with an increase in alpha(h), and (ii) become continuously higher for greater values of H/D and phi. As expected, the failure zones around the periphery of the tunnel becomes always asymmetrical with an inclusion of horizontal seismic body forces. (c) 2012 Elsevier Ltd. All rights reserved.

Role of Specific Cations and Water Entropy on the Stability of Branched DNA Motif Structures

DNA three-way junctions (TWJs) are important intermediates in various cellular processes and are the simplest of a family of branched nucleic acids being considered as scaffolds for biomolecular nanotechnology. Branched nucleic acids are stabilized by divalent cations such as Mg2+, presumably due to condensation and neutralization of the negatively charged DNA backbone. However, electrostatic screening effects point to more complex solvation dynamics and a large role of interfacial waters in thermodynamic stability. Here, we report extensive computer simulations in explicit water and salt on a model TWJ and use free energy calculations to quantify the role of ionic character and strength on stability. We find that enthalpic stabilization of the first and second hydration shells by Mg2+ accounts for 1/3 and all of the free energy gain in 50% and pure MgCl2 solutions, respectively. The more distorted DNA molecule is actually destabilized in pure MgCl2 compared to pure NaCl. Notably, the first shell, interfacial waters have very low translational and rotational entropy (i.e., mobility) compared to the bulk, an entropic loss that is overcompensated by increased enthalpy from additional electrostatic interactions with Mg2+. In contrast, the second hydration shell has anomalously high entropy as it is trapped between an immobile and bulklike layer. The nonmonotonic entropic signature and long-range perturbations of the hydration shells to Mg2+ may have implications in the molecular recognition of these motifs. For example, we find that low salt stabilizes the parallel configuration of the three-way junction, whereas at normal salt we find antiparallel configurations deduced from the NMR. We use the 2PT analysis to follow the thermodynamics of this transition and find that the free energy barrier is dominated by entropic effects that result from the decreased surface area of the antiparallel form which has a smaller number of low entropy waters in the first monolayer.

Seismic Stability and Liquefaction Analysis of Tailings Dam

Proper analysis for safe design of tailings earthen dam is necessary under static loading and more so under earthquake conditions to reduce damages of important geotechnical structure. This paper presents both static and seismic analyses of a typical section of tailings earthen dam constructed by downstream method and located at a site in eastern part India to store non-radioactive nuclear waste material. The entire analysis is performed using geotechnical softwares FLAC(3D) and TALREN 4. Results are obtained for various possible conditions of the reservoir to investigate the stability under both static and seismic loading condition using 1989 Loma Prieta earthquake acceleration-time history. FLAC(3D) analyses indicate the critical maximum displacement at crest of the proposed tailings dam section is 5.5 cm under the static loading but it increases to about 16.24 cm under seismic loading. The slope stability analyses provide the minimum value of factor of safety for seismic loading as 1.5 as compared to 2.31 for static loading. Amplification of base seismic acceleration is also observed. The liquefaction potential analysis in FLAC(3D) indicates considerable loss of shear strength in the tailings portion of the proposed earthen dam section with significantly high values of pore pressure ratio.

Precision and long-term stability of clumped-isotope analysis of CO2 using a small-sector isotope ratio mass spectrometer

RATIONALE The ratio of the measured abundance of 13C18O bonding CO2 to its stochastic abundance, prescribed by the delta 13C and delta 18O values from a carbonate mineral, is sensitive to its growth temperature. Recently, clumped-isotope thermometry, which uses this ratio, has been adopted as a new tool to elucidate paleotemperatures quantitatively. METHODS Clumped isotopes in CO2 were measured with a small-sector isotope ratio mass spectrometer. CO2 samples digested from several kinds of calcium carbonates by phosphoric acid at 25 degrees C were purified using both cryogenic and gas-chromatographic separations, and their isotopic composition (delta 13C, delta 18O, Delta 47, Delta 48 and Delta 49 values) were then determined using a dual-inlet Delta XP mass spectrometer. RESULTS The internal precisions of the single gas Delta 47 measurements were 0.005 and 0.02 parts per thousand (1 SE) for the optimum and the routine analytical conditions, respectively, which are comparable with those obtained using a MAT 253 mass spectrometer. The long-term variations in the Delta 47 values for the in-house working standard and the heated CO2 gases since 2007 were close to the routine, single gas uncertainty while showing seasonal-like periodicities with a decreasing trend. Unlike the MAT 253, the Delta XP did not show any significant relationship between the Delta 47 and delta 47 values. CONCLUSIONS The Delta XP gave results that were approximately as precise as those of the MAT 253 for clumped-isotope analysis. The temporal stability of the Delta XP seemed to be lower, although an advantage of the Delta XP was that no dependency of delta 47 on Delta 47 was found. Copyright (c) 2012 John Wiley & Sons, Ltd.

Stability of Adult Emergence and Activity/Rest Rhythms in Fruit Flies Drosophila melanogaster under Semi-Natural Condition

Here we report the results of a study aimed at examining stability of adult emergence and activity/rest rhythms under seminatural conditions (henceforth SN), in four large outbred fruit fly Drosophila melanogaster populations, selected for emergence in a narrow window of time under laboratory (henceforth LAB) light/dark (LD) cycles. When assessed under LAB, selected flies display enhanced stability in terms of higher amplitude, synchrony and accuracy in emergence and activity rhythms compared to controls. The present study was conducted to assess whether such differences in stability between selected and control populations, persist under SN where several gradually changing time-cues are present in their strongest form. The study revealed that under SN, emergence waveform of selected flies was modified, with even more enhanced peak and narrower gate-width compared to those observed in the LAB and compared to control populations in SN. Furthermore, flies from selected populations continued to exhibit enhanced synchrony and accuracy in their emergence and activity rhythms under SN compared to controls. Further analysis of zeitgeber effects revealed that enhanced stability in the rhythmicity of selected flies under SN was primarily due to increased sensitivity to light because emergence and activity rhythms of selected flies were as stable as controls under temperature cycles. These results thus suggest that stability of circadian rhythms in fruit flies D. melanogaster, which evolved as a consequence of selection for emergence in a narrow window of time under weak zeitgeber condition of LAB, persists robustly in the face of day-to-day variations in cycling environmental factors of nature.

Structure of SAICAR synthetase from Pyrococcus horikoshii OT3: Insights into thermal stability

The first native crystal structure of Phosphoribosylaminoimidazole-succinocarboxamide synthetase (SAICAR synthetase) from a hyperthermophilic organism Pyrococcus horikoshii OT3 was determined in two space groups H3 (Type-1: Resolution 2.35 angstrom) and in C222(1) (Type-2: Resolution 1.9 angstrom). Both are dimeric but Type-1 structure exhibited hexameric arrangement due to the presence of cadmium ions. A comparison has been made on the sequence and structures of all SAICAR synthetases to better understand the differences between mesophilic, thermophilic and hyperthermophilic SAICAR synthetases. These SAICAR synthetases are reasonably similar in sequence and three-dimensional structure; however, differences were visible only in the subtler details of percentage composition of the sequences, salt bridge interactions and non-polar contact areas. (c) 2012 Elsevier B.V. All rights reserved.