Correlation between conductivity or diffusivity and activation energy in amorphous solids

There exist many investigations of ionic transport in a variety of glasses. These studies exhibit strong correlation between ionic conductivity and activation energy: Typically, it is found that higher conductivity is associated with lower activation energies and vice versa. Although there are explanations for this at a phenomenological level, there is no consistent physical picture to explain the correlation between conductivity and activation energy. We have carried out molecular dynamics simulation as a function of the size of the impurity atom or diffusant (both neutral and charged) in a host amorphous matrix. We find that there is a maximum in self-diffusivity as a function of the size of the impurity atom suggesting that there is an appropriate size for which the diffusivity is maximum. The activation energy is found to be the lowest for this size of the impurity. A similar maximum has been previously found in other condensed phases, such as confined fluids and dense liquids, and has its origin in the levitation effect. The implications of this result for understanding ionic conductivity in glasses are discussed. Our results suggest that there is a relation between microscopic structure of the amorphous solid, diffusivity or conductivity, and activation energy. The nature of this relationship is discussed in terms of the levitation parameter showing that diffusivity is maximum when the size of the neck or doorway radius is comparable with the size of the diffusant. Our computational results here are in excellent agreement with independent experimental results of Nascimento et al. [Braz. J. Phys. 35, 626 (2005)] that structural features of the glass are important in determining the ionic conductivity.

The mixed black hole partition function for the STU model

We evaluate the mixed partition function for dyonic BPS black holes using the recently proposed degeneracy formula for the STU model. The result factorizes into the OSV mixed partition function times a proportionality factor. The latter is in agreement with the measure factor that was recently conjectured for a class of N = 2 black holes that contains the STU model.

Time and Energy Complexity of Distributed Computation of a Class of Functions in Wireless Sensor Networks

We consider a scenario in which a wireless sensor network is formed by randomly deploying n sensors to measure some spatial function over a field, with the objective of computing a function of the measurements and communicating it to an operator station. We restrict ourselves to the class of type-threshold functions (as defined in the work of Giridhar and Kumar, 2005), of which max, min, and indicator functions are important examples: our discussions are couched in terms of the max function. We view the problem as one of message-passing distributed computation over a geometric random graph. The network is assumed to be synchronous, and the sensors synchronously measure values and then collaborate to compute and deliver the function computed with these values to the operator station. Computation algorithms differ in (1) the communication topology assumed and (2) the messages that the nodes need to exchange in order to carry out the computation. The focus of our paper is to establish (in probability) scaling laws for the time and energy complexity of the distributed function computation over random wireless networks, under the assumption of centralized contention-free scheduling of packet transmissions. First, without any constraint on the computation algorithm, we establish scaling laws for the computation time and energy expenditure for one-time maximum computation. We show that for an optimal algorithm, the computation time and energy expenditure scale, respectively, as Theta(radicn/log n) and Theta(n) asymptotically as the number of sensors n rarr infin. Second, we analyze the performance of three specific computation algorithms that may be used in specific practical situations, namely, the tree algorithm, multihop transmission, and the Ripple algorithm (a type of gossip algorithm), and obtain scaling laws for the computation time and energy expenditure as n rarr infin. In particular, we show that the computation time for these algorithms scales as Theta(radicn/lo- g n), Theta(n), and Theta(radicn log n), respectively, whereas the energy expended scales as , Theta(n), Theta(radicn/log n), and Theta(radicn log n), respectively. Finally, simulation results are provided to show that our analysis indeed captures the correct scaling. The simulations also yield estimates of the constant multipliers in the scaling laws. Our analyses throughout assume a centralized optimal scheduler, and hence, our results can be viewed as providing bounds for the performance with practical distributed schedulers.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Comparative analysis of thermophilic and mesophilic proteins using Protein Energy Networks

Background: Thermophilic proteins sustain themselves and function at higher temperatures. Despite their structural and functional similarities with their mesophilic homologues, they show enhanced stability. Various comparative studies at genomic, protein sequence and structure levels, and experimental works highlight the different factors and dominant interacting forces contributing to this increased stability. Methods: In this comparative structure based study, we have used interaction energies between amino acids, to generate structure networks called as Protein Energy Networks (PENs). These PENs are used to compute network, sub-graph, and node specific parameters. These parameters are then compared between the thermophile-mesophile homologues. Results: The results show an increased number of clusters and low energy cliques in thermophiles as the main contributing factors for their enhanced stability. Further more, we see an increase in the number of hubs in thermophiles. We also observe no community of electrostatic cliques forming in PENs. Conclusion: In this study we were able to take an energy based network approach, to identify the factors responsible for enhanced stability of thermophiles, by comparative analysis. We were able to point out that the sub-graph parameters are the prominent contributing factors. The thermophiles have a better-packed hydrophobic core. We have also discussed how thermophiles, although increasing stability through higher connectivity retains conformational flexibility, from a cliques and communities perspective.

Top Quark Mass Measurement in the lepton+jets Channel Using a Matrix Element Method and in situ Jet Energy Calibration

A precision measurement of the top quark mass m_t is obtained using a sample of ttbar events from ppbar collisions at the Fermilab Tevatron with the CDF II detector. Selected events require an electron or muon, large missing transverse energy, and exactly four high-energy jets, at least one of which is tagged as coming from a b quark. A likelihood is calculated using a matrix element method with quasi-Monte Carlo integration taking into account finite detector resolution and jet mass effects. The event likelihood is a function of m_t and a parameter DJES to calibrate the jet energy scale /in situ/. Using a total of 1087 events, a value of m_t = 173.0 +/- 1.2 GeV/c^2 is measured.

Asymtotic solution of the Frieman and Book equation for the pair-correlation function

Using Thomé's procedure, the asymptotic solutions of the Frieman and Book equation for the two-particle correlation in a plasma have been obtained in a complete form. The solution is interpreted in terms of the Lorentz distance. The exact expressions for the internal energy and pressure are evaluated and they are found to be a generalization of the result obtained earlier by others.

The function of NR3B and NR4A orphan nuclear receptors in osteoblasts

Nuclear receptors (NRs) comprise a large family of proteins that mediate the effects of small lipophilic molecules such as steroid hormones. In addition, there are a group of NRs which lack identified natural ligands and are referred as orphan NRs. In this thesis, the function of two such orphan NR families, the NR3B (ERRα, ERRβ and ERRγ) and the NR4A family (NGFI-B, Nurr1 and Nor1), was studied. NR3B and NR4A receptors regulate many biological processes such as energy metabolism and carcinogenesis. In addition, NR3B and NR4A receptors are expressed in bone. Therefore, the signaling and function of NR3B and NR4A orphan nuclear receptors was studied specifically in osteoblasts. NR4A receptors were found to be regulated by NR3B receptors and the Wnt/β-catenin signaling pathway as ERRα, ERRγ and β-catenin repressed the transcriptional activity of NR4A receptors in U2-OS cells. NGFI-B was found to repress the transcriptional activity of ERRγ in HeLa cells. The phytoestrogen equol was identified as a new agonist for ERRγ and ERRβ in PC-3, U2-OS, and SaOS-2 cells. Equol increased the transcriptional activity of ERRγ by increasing ERRγ co-activator binding and by inducing a conformational change in the ligand binding pocket of ERRγ. The growth inhibitory effect of equol on PC-3 prostate cancer cells was decreased by blocking ERRγ expression by siRNA. Therefore, ERRγ could mediate some of the beneficial health effects of equol. The Wnt/β-catenin signaling pathway is important for the differentiation and function of osteoblasts. NR3B and NR4A receptors were found to repress the transcriptional activity mediated by β-catenin in U2-OS cells. The mesenchymal stem cells (MSCs) isolated from ERRα knockout (KO) mice showed diminished proliferation and osteoblastic differentiation compared to the wild-type cells. The overexpression of ERRα in osteoblastic MC3T3-E1 cell line increased their mineralization. Bone sialoprotein (BSP) was shown to be a direct target gene for ERRα and ERRγ as the BSP promoter was activated by ERRα or ERRγ and PGC-1α in HeLa cells. The adipogenic differentiation of ERRα KO MSCs was also decreased and they expressed less adipogenic marker genes. In conclusion, the studies described in this thesis demonstrated that the transcriptional activity of NR3B and NR4A receptors can be regulated by other orphan NRs and signaling pathways in osteoblasts. NR3B receptors can also be regulated by ligands and a new agonist, equol, was identified for ERRβ and ERRγ. New roles for NR3B and NR4A were also identified as they were shown to converge with the Wnt signaling pathway in osteoblasts, ERRγ was shown to mediate the growth inhibitory effect of equol in prostate cancer cells, and ERRα was shown to regulate positively MSC proliferation, osteoblastic differentiation and adipogenesis.

Elucidation of the conformational free energy landscape in H.pylori LuxS and its implications to catalysis

Background: One of the major challenges in understanding enzyme catalysis is to identify the different conformations and their populations at detailed molecular level in response to ligand binding/environment. A detail description of the ligand induced conformational changes provides meaningful insights into the mechanism of action of enzymes and thus its function. Results: In this study, we have explored the ligand induced conformational changes in H. pylori LuxS and the associated mechanistic features. LuxS, a dimeric protein, produces the precursor (4,5-dihydroxy-2,3-pentanedione) for autoinducer-2 production which is a signalling molecule for bacterial quorum sensing. We have performed molecular dynamics simulations on H. pylori LuxS in its various ligand bound forms and analyzed the simulation trajectories using various techniques including the structure network analysis, free energy evaluation and water dynamics at the active site. The results bring out the mechanistic details such as co operativity and asymmetry between the two subunits, subtle changes in the conformation as a response to the binding of active and inactive forms of ligands and the population distribution of different conformations in equilibrium. These investigations have enabled us to probe the free energy landscape and identify the corresponding conformations in terms of network parameters. In addition, we have also elucidated the variations in the dynamics of water co-ordination to the Zn2+ ion in LuxS and its relation to the rigidity at the active sites. Conclusions: In this article, we provide details of a novel method for the identification of conformational changes in the different ligand bound states of the protein, evaluation of ligand-induced free energy changes and the biological relevance of our results in the context of LuxS structure-function. The methodology outlined here is highly generalized to illuminate the linkage between structure and function in any protein of known structure.

Gibbs energy of formation of orthorhombic CaZrO3

The Gibbs free energy of formation of the orthorhombic form of CaZrO3(o) from monoclinic ZrO2(m) and periclase CaO(p) has been determined as a function of temperature in the range 950-1225 K, using an electrochemical cell incorporating single-crystal CaF2 as the solid electrolyte. The results are corrected for the small solid solubility of CaO in ZrO2. For the reaction, ZrO2(m) + CaO(p) --> CaZrO3(o), DELTAG(phi) = -31590 -13.9T(+/- 180) J mol-1. The ''second-law'' enthalpy of formation of CaZrO3 obtained from the results of this study at a mean temperature of 1090 K is in excellent agreement with the high-temperature solution calorimetric measurements of Muromachi and Navrotsky at 1068 K (J. Solid State Chem., 72 (1988) 244), and the average value of the bomb and acid solution calorimetric studies of Lvova and Feodosev (Zh. Fiz. Khim., 38 (1964) 28), Korneev et al. (Izv. Akad. Nauk SSSR, Neorg. Mater., 7 (1971) 886) and Brown and Bennington (Thermochim. Acta, 106 (1986) 183). The standard entropy of CaZrO3(o) at 298.15 K from the free energy data is 96.4 (+/- 3.5) J K-1 mol-1. The results of this study are discussed in comparison with high-temperture e.m.f. measurements reported in the literature on cubic zirconia solid solutions.

Solvation dynamics, energy distribution and trapping of a light solute ion

In the theoretical treatments of the dynamics of solvation of a newly created ion in a dipolar solvent, the self-motion of the solute is usually ignored. Recently, it has been shown that for a light ion the translational motion of the ion can significantly enhance its own rate of solvation. Therefore, solvation itself may not be the rate determining step in the equilibration. Instead, the rate determining step is the search of the low energy configuration which serves to localize the light ion. In this article a microscopic calculation of the probability distribution of the interaction energy of the nascent charge with the dipolar solvent molecules is presented in order to address this problem of solute trapping. It is found that to a good approximation, this distribution is Gaussian and the second moment of this distribution is exactly equal to the half of its own solvation energy. It is shown that this is in excellent agreement with the simulation results that are available for the model Brownian dipolar lattice and for liquid acetonitrile. If the distortion of the solvent by the ion is negligible then the same relation gives the energy distribution for the solvated ion, with the average centered at the final equilibrium solvation energy. These results are expected to be useful in understanding various chemical processes in dipolar liquids. Another interesting outcome of the present study is a simple dynamic argument that supports Onsager's ''inverse snow-ball'' conjecture of solvation of a light ion. A simple derivation of the semi-phenomenological relation between the solvation time correlation function and the single particle orientation, reported recently by Maroncelli et al. (J. Phys. Chem. 97 (1993) 13), is also presented.

Effect of magnetic field on the binding energy of a hydrogenic donor in a vertical bar z vertical bar(2/3) quantum well

We have calculated the binding energy of a hydrogenic donor in a quantum well with potential shape proportional to \z\(2/3) as a function of the width of the quantum well and the barrier height under an applied uniform magnetic field along the a axis. As the well width decreases, the binding energy increases initially up to a critical well width (which is nearly the same for all magnetic fields) at which there is a turnover. The results are qualitatively similar to those of a hydrogenic donor in a rectangular well. We have also calculated [rho(2)](1/2) and [z(2)](1/2) for the donor electron. [rho(2)](1/2) is found to be strongly dependent on the magnetic field for a given well width and weakly dependent on the well width and the barrier height, for a given value of magnetic field [z(2)](1/2) is weakly dependent on the applied magnetic field. The probability of finding the donor electron inside the well shows a rapid decrease as the well width is reduced at nearly the well width at which the binding energy shows a maximum.

The structures of bacteriorhodopsin with different retinal-Schiff base orientations--computer modeling and energy minimization studies

Bacteriorhodopsin has been the subject of intense study in order to understand its photochemical function. The recent atomic model proposed by Henderson and coworkers based on electron cryo-microscopic studies has helped in understanding many of the structural and functional aspects of bacteriorhodopsin. However, the accuracy of the positions of the side chains is not very high since the model is based on low-resolution data. In this study, we have minimized the energy of this structure of bacteriorhodopsin and analyzed various types of interactions such as - intrahelical and interhelical hydrogen bonds and retinal environment. In order to understand the photochemical action, it is necessary to obtain information on the structures adopted at the intermediate states. In this direction, we have generated some intermediate structures taking into account certain experimental data, by computer modeling studies. Various isomers of retinal with 13-cis and/or 15-cis conformations and all possible staggered orientations of Lys-216 side chain were generated. The resultant structures were examined for the distance between Lys-216-schiff base nitrogen and the carboxylate oxygen atoms of Asp-96 - a residue which is known to reprotonate the schiff base at later stages of photocycle. Some of the structures were selected on the basis of suitable retinal orientation and the stability of these structures were tested by energy minimization studies. Further, the minimized structures are analyzed for the hydrogen bond interactions and retinal environment and the results are compared with those of the minimized rest state structure. The importance of functional groups in stabilizing the structure of bacteriorhodopsin and in participating dynamically during the photocycle have been discussed.

Escape of high mass ions due to initial thermal energy and its implications for axially symmetric RF ion trap mass analyzer design

This paper investigates the loss of high mass ions due to their initial thermal energy in ion trap mass analyzers. It provides an analytical expression for estimating the percentage loss of ions of a given mass at a particular temperature, in a trap operating under a predetermined set of conditions. The expression we developed can be used to study the loss of ions due to its initial thermal energy in traps which have nonlinear fields as well as those which have linear fields. The expression for the percentage of ions lost is shown to be a function of the temperature of the ensemble of ions, ion mass and ion escape velocity. An analytical expression for the escape velocity has also been derived in terms of the trapping field, drive frequency and ion mass. Because the trapping field is determined by trap design parameters and operating conditions, it has been possible to study the influence of these parameters on ion loss. The parameters investigated include ion temperature, magnitude of the initial potential applied to the ring electrode (which determines the low mass cut-off), trap size, dimensions of apertures in the endcap electrodes and RF drive frequency. Our studies demonstrate that ion loss due to initial thermal energy increases with increase in mass and that, in the traps investigated, ion escape occurs in the radial direction. Reduction in the loss of high mass ions is favoured by lower ion temperatures, increasing low mass cut-off, increasing trap size, and higher RF drive frequencies. However, dimensions of the apertures in the endcap electrodes do not influence ion loss in the range of aperture sizes considered. (C) 2010 Elsevier B.V. All rights reserved.

Absolute Entropy and Energy of Carbon Dioxide Using the Two-Phase Thermodynamic Model

The two-phase thermodynamic (2PT) model is used to determine the absolute entropy and energy of carbon dioxide over a wide range of conditions from molecular dynamics trajectories. The 2PT method determines the thermodynamic properties by applying the proper statistical mechanical partition function to the normal modes of a fluid. The vibrational density of state (DoS), obtained from the Fourier transform of the velocity autocorrelation function, converges quickly, allowing the free energy, entropy, and other thermodynamic properties to be determined from short 20-ps MD trajectories. The anharmonic effects in the vibrations are accounted for by the broadening of the normal modes into bands from sampling the velocities over the trajectory. The low frequency diffusive modes, which lead to finite DoS at zero frequency, are accounted for by considering the DoS as a superposition of gas-phase and solid-phase components (two phases). The analytical decomposition of the DoS allows for an evaluation of properties contributed by different types of molecular motions. We show that this 2PT analysis leads to accurate predictions of entropy and energy of CO2 over a wide range of conditions (from the triple point to the critical point of both the vapor and the liquid phases along the saturation line). This allows the equation of state of CO2 to be determined, which is limited only by the accuracy of the force field. We also validated that the 2PT entropy agrees with that determined from thermodynamic integration, but 2PT requires only a fraction of the time. A complication for CO2 is that its equilibrium configuration is linear, which would have only two rotational modes, but during the dynamics it is never exactly linear, so that there is a third mode from rotational about the axis. In this work, we show how to treat such linear molecules in the 2PT framework.