The Equilibrium Structure of Lithium Salt Solutions in Ether-Functionalized Ammonium Ionic Liquids

Molecular dynamics simulations have been performed for ionic liquids based on a ternary mixture of lithium and ammonium cations and a common anion, bis(trifluoromethylsulfonyl)imide, [Tf2N](-). We address structural changes resulting from adding Li+ in ionic liquids with increasing length of an ether-functionalized chain in the ammonium cation. The calculation of static structure factors reveals the lithium effect on charge ordering and intermediate range order in comparison with the neat ionic liquids. The charge ordering is modified in the lithium solution because the coordination of [Tf2N](-) toward Li+ is much stronger than ammonium cations. Intermediate range order is observed in neat ionic liquids based on ammonium cations with a long chain, but in the lithium solutions, there is also a nonhomogenous distribution of Li+ cations. The presence of Li+ enhances interactions between the ammonium cations due to correlations between the oxygen atom of the ether chain and the nitrogen atom of another ammonium cation.

Coupled Boltzmann calculation of mixed axion/neutralino cold dark matter production in the early universe

We calculate the relic abundance of mixed axion/neutralino cold dark matter which arises in R-parity conserving supersymmetric (SUSY) models wherein the strong CP problem is solved by the Peccei-Quinn (PQ) mechanism with a concommitant axion/saxion/axino supermultiplet. By numerically solving the coupled Boltzmann equations, we include the combined effects of 1. thermal axino production with cascade decays to a neutralino LSP, 2. thermal saxion production and production via coherent oscillations along with cascade decays and entropy injection, 3. thermal neutralino production and re-annihilation after both axino and saxion decays, 4. gravitino production and decay and 5. axion production both thermally and via oscillations. For SUSY models with too high a standard neutralino thermal abundance, we find the combined effect of SUSY PQ particles is not enough to lower the neutralino abundance down to its measured value, while at the same time respecting bounds on late-decaying neutral particles from BBN. However, models with a standard neutralino underabundance can now be allowed with either neutralino or axion domination of dark matter, and furthermore, these models can allow the PQ breaking scale f(a) to be pushed up into the 10(14) - 10(15) GeV range, which is where it is typically expected to be in string theory models.

A novel self consistent calculation approach for the capacitance-voltage characteristics of semiconductor quantum wire transistors based on a split-gate configuration

Purpose - The purpose of this paper is to develop an efficient numerical algorithm for the self-consistent solution of Schrodinger and Poisson equations in one-dimensional systems. The goal is to compute the charge-control and capacitance-voltage characteristics of quantum wire transistors. Design/methodology/approach - The paper presents a numerical formulation employing a non-uniform finite difference discretization scheme, in which the wavefunctions and electronic energy levels are obtained by solving the Schrodinger equation through the split-operator method while a relaxation method in the FTCS scheme ("Forward Time Centered Space") is used to solve the two-dimensional Poisson equation. Findings - The numerical model is validated by taking previously published results as a benchmark and then applying them to yield the charge-control characteristics and the capacitance-voltage relationship for a split-gate quantum wire device. Originality/value - The paper helps to fulfill the need for C-V models of quantum wire device. To do so, the authors implemented a straightforward calculation method for the two-dimensional electronic carrier density n(x,y). The formulation reduces the computational procedure to a much simpler problem, similar to the one-dimensional quantization case, significantly diminishing running time.

Supersymmetric Isospectral Formalism for the Calculation of Near-Zero Energy States: Application to the Very Weakly Bound He-4 Trimer Excited State

We propose a novel mathematical approach for the calculation of near-zero energy states by solving potentials which are isospectral with the original one. For any potential, families of strictly isospectral potentials (with very different shape) having desirable and adjustable features are generated by supersymmetric isospectral formalism. The near-zero energy Efimov state in the original potential is effectively trapped in the deep well of the isospectral family and facilitates more accurate calculation of the Efimov state. Application to the first excited state in He-4 trimer is presented.

Generalized Metropolis dynamics with a generalized master equation: An approach for time-independent and time-dependent Monte Carlo simulations of generalized spin systems

The extension of Boltzmann-Gibbs thermostatistics, proposed by Tsallis, introduces an additional parameter q to the inverse temperature beta. Here, we show that a previously introduced generalized Metropolis dynamics to evolve spin models is not local and does not obey the detailed energy balance. In this dynamics, locality is only retrieved for q = 1, which corresponds to the standard Metropolis algorithm. Nonlocality implies very time-consuming computer calculations, since the energy of the whole system must be reevaluated when a single spin is flipped. To circumvent this costly calculation, we propose a generalized master equation, which gives rise to a local generalized Metropolis dynamics that obeys the detailed energy balance. To compare the different critical values obtained with other generalized dynamics, we perform Monte Carlo simulations in equilibrium for the Ising model. By using short-time nonequilibrium numerical simulations, we also calculate for this model the critical temperature and the static and dynamical critical exponents as functions of q. Even for q not equal 1, we show that suitable time-evolving power laws can be found for each initial condition. Our numerical experiments corroborate the literature results when we use nonlocal dynamics, showing that short-time parameter determination works also in this case. However, the dynamics governed by the new master equation leads to different results for critical temperatures and also the critical exponents affecting universality classes. We further propose a simple algorithm to optimize modeling the time evolution with a power law, considering in a log-log plot two successive refinements.

STRUCTURE AND DYNAMICS: THE TRANSITION FROM NONEQUILIBRIUM TO EQUILIBRIUM IN INTEGRATE-AND-FIRE DYNAMICS

The transient and equilibrium properties of dynamics unfolding in complex systems can depend critically on specific topological features of the underlying interconnections. In this work, we investigate such a relationship with respect to the integrate-and-fire dynamics emanating from a source node and an extended network model that allows control of the small-world feature as well as the length of the long-range connections. A systematic approach to investigate the local and global correlations between structural and dynamical features of the networks was adopted that involved extensive simulations (one and a half million cases) so as to obtain two-dimensional correlation maps. Smooth, but diverse surfaces of correlation values were obtained in all cases. Regarding the global cases, it has been verified that the onset avalanche time (but not its intensity) can be accurately predicted from the structural features within specific regions of the map (i.e. networks with specific structural properties). The analysis at local level revealed that the dynamical features before the avalanches can also be accurately predicted from structural features. This is not possible for the dynamical features after the avalanches take place. This is so because the overall topology of the network predominates over the local topology around the source at the stationary state.

Liquid-Liquid Equilibrium Data for the System Lard plus Oleic Acid plus Ethanol + Water at 318.2 K: Cholesterol Distribution Coefficients

This research reports liquid liquid equilibrium data for the system lard (swine fat), cis-9-octadecenoic acid (oleic acid), ethanol, and water at 318.2 K, as well as their correlation with the nonrandom two-liquid (NRTL) and universal quasichemical activity coefficient (UNIQUAC) thermodynamic equations, which have provided global deviations of 0.41 % and 0.53 %, respectively. Additional equilibrium experiments were also performed to obtain cholesterol partition (or distribution) coefficients to verify the availability of the use of ethanol plus water to reduce the cholesterol content in lard. The partition experiments were performed with concentrations of free fatty acids (commercial oleic acid) that varied from (0 to 20) mass % and of water in the solvent that varied from (0 to 18) mass %. The percentage of free fatty acids initially present in lard had a slight effect on the distribution of cholesterol between the phases. Furthermore, the distribution coefficients decreased by adding water in the ethanol; specifically, it resulted in a diminution of the capability of the solvent to remove the cholesterol.

Fractionation of lemon essential oil by solvent extraction: Phase equilibrium for model systems at T=298.2 K

The fractioning of lemon essential oil can be performed by liquid-liquid extraction using hydrous ethanol as a solvent. A quaternary mixture composed of limonene, gamma-terpinene, beta-pinene, and citral was used to simulate lemon essential oil. In this paper, we present (liquid + liquid) equilibrium data that were experimentally determined for systems containing essential oil compounds, ethanol, and water at T = 298.2 K. The experimental data were correlated using the NRTL and UNIQUAC models, and the mean deviations between calculated and experimental data were less than 0.0053 in all systems, indicating the accuracy of these molecular models in describing our systems. The results show that as the water content in the solvent phase increased, the values of the distribution coefficients decreased, regardless of the type of compound studied. However, the oxygenated compound always showed the highest distribution coefficient among the components of the essential oil, thus making deterpenation of the lemon essential oil a feasible process. (C) 2012 Elsevier Ltd. All rights reserved.

Application of succinylated sugarcane bagasse as adsorbent to remove methylene blue and gentian violet from aqueous solutions - Kinetic and equilibrium studies

In a previous work, succinylated sugarcane bagasse (SCB 2) was prepared from sugarcane bagasse (B) using succinic anhydride as modifying agent. In this work the adsorption of cationic dyes onto SCB 2 from aqueous solutions was investigated. Methylene blue, MB, and gentian violet, GV, were selected as adsorbates. The capacity of SCB 2 to adsorb MB and GV from aqueous single dye solutions was evaluated at different contact times, pH, and initial adsorbent concentration. According to the obtained results, the adsorption processes could be described by the pseudo-second-order kinetic model. Adsorption isotherms were well fitted by Langmuir model. Maximum adsorption capacities for MB and GV onto SCB 2 were found to be 478.5 and 1273.2 mg/g, respectively. (C) 2011 Elsevier Ltd. All rights reserved.

Comparison of different dosimetric methods for red marrow absorbed dose calculation in thyroid cancer therapy

Several dosimetric methods have been proposed for estimating red marrow absorbed dose (RMAD) when radionuclide therapy is planned for differentiated thyroid cancer, although to date, there is no consensus as to whether dose calculation should be based on blood-activity concentration or not. Our purpose was to compare RMADs derived from methods that require collecting patients' blood samples versus those involving OLINDA/EXM software, thereby precluding this invasive procedure. This is a retrospective study that included 34 patients under treatment for metastatic thyroid disease. A deviation of 10 between RMADs was found, when comparing the doses from the most usual invasive dosimetric methods and those from OLINDA/EXM. No statistical difference between the methods was discovered, whereby the need for invasive procedures when calculating the dose is questioned. The use of OLINDA/EXM in clinical routine could possibly diminish data collection, thus giving rise to a simultaneous reduction in time and clinical costs, besides avoiding any kind of discomfort on the part of the patients involved.

Influence of Alloying Elements Upon the Solidification Interval of CA6NM Cast Martensitic Stainless Steel

The influence of different Cr and C contents upon the solidification interval of ASTM A352M-06 Grade CA6NM cast martensitic stainless steel has been investigated using computational thermodynamics, and checked against DTA measurements in samples taken from 13 large cast parts, in order to identify potential sources for improvement on the part castability. Calculation results suggest, indeed, that this would be the case for C: when its content increases from 0.018 to 0.044 wt.% C (within the allowed range in the alloy specification), the solidification intervals increases from 25 to 43 K, which suggests improved castability with decreasing C contents. DTA results, however, do not support this prediction, showing a fairly constant solidification interval around 23 K for all investigated samples. The results are discussed both regarding the impact in alloy processing and the fitness of the existing databases to reproduce experimental results in these limiting cases.

Influence of the initial condition in equilibrium last-passage percolation models

In this paper we consider an equilibrium last-passage percolation model on an environment given by a compound two-dimensional Poisson process. We prove an L-2-formula relating the initial measure with the last-passage percolation time. This formula turns out to be a useful tool to analyze the fluctuations of the last-passage times along non-characteristic directions.

On Equilibrium Distribution of a Reversible Growth Model

We study a probabilistic model of interacting spins indexed by elements of a finite subset of the d-dimensional integer lattice, da parts per thousand yen1. Conditions of time reversibility are examined. It is shown that the model equilibrium distribution converges to a limit distribution as the indexing set expands to the whole lattice. The occupied site percolation problem is solved for the limit distribution. Two models with similar dynamics are also discussed.

Calculation of parameter ranges for robust gain tuning of power system controllers

This work proposes a computational tool to assist power system engineers in the field tuning of power system stabilizers (PSSs) and Automatic Voltage Regulators (AVRs). The outcome of this tool is a range of gain values for theses controllers within which there is a theoretical guarantee of stability for the closed-loop system. This range is given as a set of limit values for the static gains of the controllers of interest, in such a way that the engineer responsible for the field tuning of PSSs and/or AVRs can be confident with respect to system stability when adjusting the corresponding static gains within this range. This feature of the proposed tool is highly desirable from a practical viewpoint, since the PSS and AVR commissioning stage always involve some readjustment of the controller gains to account for the differences between the nominal model and the actual behavior of the system. By capturing these differences as uncertainties in the model, this computational tool is able to guarantee stability for the whole uncertain model using an approach based on linear matrix inequalities. It is also important to remark that the tool proposed in this paper can also be applied to other types of parameters of either PSSs or Power Oscillation Dampers, as well as other types of controllers (such as speed governors, for example). To show its effectiveness, applications of the proposed tool to two benchmarks for small signal stability studies are presented at the end of this paper.

Shear force and torsion in reinforced concrete beam elements: theoretical analysis based on Brazilian Standard Code ABNT NBR 6118:2007

Reinforced concrete beam elements are submitted to applicable loads along their life cycle that cause shear and torsion. These elements may be subject to only shear, pure torsion or both, torsion and shear combined. The Brazilian Standard Code ABNT NBR 6118:2007 [1] fixes conditions to calculate the transverse reinforcement area in beam reinforced concrete elements, using two design models, based on the strut and tie analogy model, first studied by Mörsch [2]. The strut angle θ (theta) can be considered constant and equal to 45º (Model I), or varying between 30º and 45º (Model II). In the case of transversal ties (stirrups), the variation of angle α (alpha) is between 45º and 90º. When the equilibrium torsion is required, a resistant model based on space truss with hollow section is considered. The space truss admits an inclination angle θ between 30º and 45º, in accordance with beam elements subjected to shear. This paper presents a theoretical study of models I and II for combined shear and torsion, in which ranges the geometry and intensity of action in reinforced concrete beams, aimed to verify the consumption of transverse reinforcement in accordance with the calculation model adopted As the strut angle on model II ranges from 30º to 45º, transverse reinforcement area (Asw) decreases, and total reinforcement area, which includes longitudinal torsion reinforcement (Asℓ), increases. It appears that, when considering model II with strut angle above 40º, under shear only, transverse reinforcement area increases 22% compared to values obtained using model I.