Evaluation of serial C-reactive protein measurements after surgical treatment of pleural empyema

OBJECTIVE: Serial C-reactive protein measurements have been used to diagnose and monitor the response to therapy in patients with pneumonia and other infectious diseases. Nonetheless, the role of C-reactive protein measurement after surgical treatment for pleural empyema is not well defined. The aim of this study is to describe the behavior of C-reactive protein levels after the surgical treatment of pleural empyema and to correlate this parameter with the patient's prognosis. METHODS: We retrospectively analyzed the records of patients with pleural empyema treated by either chest-tube drainage or surgery from January 2006 to December 2008. C-reactive protein levels were recorded preoperatively and 2 and 7 days postoperatively. The clinical outcome was binary: success or failure (mortality or the need for repeated pleural intervention). RESULTS: The study group comprised fifty-two patients. The median C-reactive protein values were as follows: 146 mg/L (pre-operative), 134 mg/L (post-operative day 2), and 116 mg/L (post-operative day 7). There was a trend toward a decrease in these values during the first week after surgery, but this difference was only statistically significant on day 7 after surgery. Over the first week after surgery, the C-reactive protein values decreased similarly in both groups (successful and failed treatment). No correlation between the preoperative C-reactive protein level and the clinical outcome was found. CONCLUSIONS: We observed that, in contrast to other medical conditions, C-reactive protein levels fall slowly during the first postoperative week in patients who have undergone surgical treatment for pleural empyema. No correlation between the perioperative C-reactive protein level and the clinical outcome was observed.

Quantum spherical model with competing interactions

We analyse the phase diagram of a quantum mean spherical model in terms of the temperature T, a quantum parameter g, and the ratio p = -J(2)/J(1) where J(1) > 0 refers to ferromagnetic interactions between first-neighbour sites along the d directions of a hypercubic lattice, and J(2) < 0 is associated with competing anti ferromagnetic interactions between second neighbours along m <= d directions. We regain a number of known results for the classical version of this model, including the topology of the critical line in the g = 0 space, with a Lifshitz point at p = 1/4, for d > 2, and closed-form expressions for the decay of the pair correlations in one dimension. In the T = 0 phase diagram, there is a critical border, g(c) = g(c) (p) for d >= 2, with a singularity at the Lifshitz point if d < (m + 4)/2. We also establish upper and lower critical dimensions, and analyse the quantum critical behavior in the neighborhood of p = 1/4. 2012 (C) Elsevier B.V. All rights reserved.

Parameter recovery for a skew-normal IRT model under a Bayesian approach: hierarchical framework, prior and kernel sensitivity and sample size

Item response theory (IRT) comprises a set of statistical models which are useful in many fields, especially when there is an interest in studying latent variables (or latent traits). Usually such latent traits are assumed to be random variables and a convenient distribution is assigned to them. A very common choice for such a distribution has been the standard normal. Recently, Azevedo et al. [Bayesian inference for a skew-normal IRT model under the centred parameterization, Comput. Stat. Data Anal. 55 (2011), pp. 353-365] proposed a skew-normal distribution under the centred parameterization (SNCP) as had been studied in [R. B. Arellano-Valle and A. Azzalini, The centred parametrization for the multivariate skew-normal distribution, J. Multivariate Anal. 99(7) (2008), pp. 1362-1382], to model the latent trait distribution. This approach allows one to represent any asymmetric behaviour concerning the latent trait distribution. Also, they developed a Metropolis-Hastings within the Gibbs sampling (MHWGS) algorithm based on the density of the SNCP. They showed that the algorithm recovers all parameters properly. Their results indicated that, in the presence of asymmetry, the proposed model and the estimation algorithm perform better than the usual model and estimation methods. Our main goal in this paper is to propose another type of MHWGS algorithm based on a stochastic representation (hierarchical structure) of the SNCP studied in [N. Henze, A probabilistic representation of the skew-normal distribution, Scand. J. Statist. 13 (1986), pp. 271-275]. Our algorithm has only one Metropolis-Hastings step, in opposition to the algorithm developed by Azevedo et al., which has two such steps. This not only makes the implementation easier but also reduces the number of proposal densities to be used, which can be a problem in the implementation of MHWGS algorithms, as can be seen in [R.J. Patz and B.W. Junker, A straightforward approach to Markov Chain Monte Carlo methods for item response models, J. Educ. Behav. Stat. 24(2) (1999), pp. 146-178; R. J. Patz and B. W. Junker, The applications and extensions of MCMC in IRT: Multiple item types, missing data, and rated responses, J. Educ. Behav. Stat. 24(4) (1999), pp. 342-366; A. Gelman, G.O. Roberts, and W.R. Gilks, Efficient Metropolis jumping rules, Bayesian Stat. 5 (1996), pp. 599-607]. Moreover, we consider a modified beta prior (which generalizes the one considered in [3]) and a Jeffreys prior for the asymmetry parameter. Furthermore, we study the sensitivity of such priors as well as the use of different kernel densities for this parameter. Finally, we assess the impact of the number of examinees, number of items and the asymmetry level on the parameter recovery. Results of the simulation study indicated that our approach performed equally as well as that in [3], in terms of parameter recovery, mainly using the Jeffreys prior. Also, they indicated that the asymmetry level has the highest impact on parameter recovery, even though it is relatively small. A real data analysis is considered jointly with the development of model fitting assessment tools. The results are compared with the ones obtained by Azevedo et al. The results indicate that using the hierarchical approach allows us to implement MCMC algorithms more easily, it facilitates diagnosis of the convergence and also it can be very useful to fit more complex skew IRT models.

Stripe-tetragonal phase transition in the two-dimensional Ising model with dipole interactions: Partition function zeros approach

We have performed multicanonical simulations to study the critical behavior of the two-dimensional Ising model with dipole interactions. This study concerns the thermodynamic phase transitions in the range of the interaction delta where the phase characterized by striped configurations of width h = 1 is observed. Controversial results obtained from local update algorithms have been reported for this region, including the claimed existence of a second-order phase transition line that becomes first order above a tricritical point located somewhere between delta = 0.85 and 1. Our analysis relies on the complex partition function zeros obtained with high statistics from multicanonical simulations. Finite size scaling relations for the leading partition function zeros yield critical exponents. that are clearly consistent with a single second-order phase transition line, thus excluding such a tricritical point in that region of the phase diagram. This conclusion is further supported by analysis of the specific heat and susceptibility of the orientational order parameter.

Solvation in aqueous binary mixtures: consequences of the hydrophobic character of the ionic liquids and the solvatochromic probes

Information on the solvation in mixtures of water, W, and the ionic liquids, ILs, 1-allyl-3-R-imidazolium chlorides; R = methyl, 1-butyl, and 1-hexyl, has been obtained from the responses of the following solvatochromic probes: 2,6-dibromo-4-[(E)-2-(1-R-pyridinium-4-yl)ethenyl] phenolate, R = methyl, MePMBr2; 1-octyl, OcPMBr(2), and the corresponding quinolinium derivative, MeQMBr(2). A model developed for solvation in binary mixtures of W and molecular solvents has been extended to the present mixtures. Our objective is to assess the relevance to solvation of hydrogen-bonding and the hydrophobic character of the IL and the solvatochromic probe. Plots of the medium empirical polarity, E-T(probe) versus its composition revealed non-ideal behavior, attributed to preferential solvation by the IL and, more efficiently, by the IL-W hydrogen-bonded complex. The deviation from linearity increases as a function of increasing number of carbon atoms in the alkyl group of the IL, and is larger than that observed for solvation by W plus molecular solvents (1-propanol and 2-(1-butoxy)ethanol) that are more hydrophobic than the ILs investigated. This enhanced deviation is attributed to the more organized structure of the ILs proper, which persists in their aqueous solutions. MeQMBr(2) is more susceptible to solvent lipophilicity than OcPMBr(2), although the former probe is less lipophilic. This enhanced susceptibility agrees with the important effect of annelation on the contributions of the quinonoid and zwitterionic limiting structures to the ground and excited states of the probe, hence on its response to both medium composition and lipophilicity of the IL.

Modeling the gluon propagator in Landau gauge: Lattice estimates of pole masses and dimension-two condensates

We present an analytic description of numerical results for the Landau-gauge SU(2) gluon propagator D(p(2)), obtained from lattice simulations (in the scaling region) for the largest lattice sizes to date, in d = 2, 3 and 4 space-time dimensions. Fits to the gluon data in 3d and in 4d show very good agreement with the tree-level prediction of the refined Gribov-Zwanziger (RGZ) framework, supporting a massive behavior for D(p(2)) in the infrared limit. In particular, we investigate the propagator's pole structure and provide estimates of the dynamical mass scales that can be associated with dimension-two condensates in the theory. In the 2d case, fitting the data requires a noninteger power of the momentum p in the numerator of the expression for D(p(2)). In this case, an infinite-volume-limit extrapolation gives D(0) = 0. Our analysis suggests that this result is related to a particular symmetry in the complex-pole structure of the propagator and not to purely imaginary poles, as would be expected in the original Gribov-Zwanziger scenario.

In-plane mapping of buried InGaAs quantum rings and hybridization effects on the electronic structure

This work reports the investigation on the structural differences between InAs quantum rings and their precursor quantum dots species as well as on the presence of piezoelectric fields and asymmetries in these nanostructures. The experimental results show significant reduction in the ring dimensions when the sizes of capped and uncapped ring and dot samples are compared. The iso-lattice parameter mapped by grazing-incidence x-ray diffraction has revealed the lateral extent of strained regions in the buried rings. A comparison between strain and composition of dot and ring structures allows inferring on how the ring formation and its final configuration may affect optical response parameters. Based on the experimental observations, a discussion has been introduced on the effective potential profile to emulate theoretically the ring-shape confinement. The effects of confinement and strain field modulation on electron and hole band structures are simulated by a multiband k.p calculation. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4733964]

Biaxial nematic phase in the Maier-Saupe model for a mixture of discs and cylinders

We analyze the global phase diagram of a Maier-Saupe lattice model with the inclusion of shape-disordered degrees of freedom to mimic a mixture of oblate and prolate molecules (discs and cylinders). In the neighborhood of a Landau multicritical point, solutions of the statistical problem can be written as a Landau-de Gennes expansion for the free energy. If the shape-disordered degrees of freedom are quenched, we confirm the existence of a biaxial nematic structure. If orientational and disorder degrees of freedom are allowed to thermalize, this biaxial solution becomes thermodynamically unstable. Also, we use a two-temperature formalism to mimic the presence of two distinct relaxation times, and show that a slight departure from complete thermalization is enough to stabilize a biaxial nematic phase.

Exploration of the parameter space in an agent-based model of tuberculosis spread: emergence of drug resistance in developing vs developed countries

In this work we present an agent-based model for the spread of tuberculosis where the individuals can be infected with either drug-susceptible or drug-resistant strains and can also receive a treatment. The dynamics of the model and the role of each one of the parameters are explained. The whole set of parameters is explored to check their importance in the numerical simulation results. The model captures the beneficial impact of the adequate treatment on the prevalence of tuberculosis. Nevertheless, depending on the treatment parameters range, it also captures the emergence of drug resistance. Drug resistance emergence is particularly likely to occur for parameter values corresponding to less efficacious treatment, as usually found in developing countries.

Structure and anomalous solubility for hard spheres in an associating lattice gas model

In this paper we investigate the solubility of a hard-sphere gas in a solvent modeled as an associating lattice gas. The solution phase diagram for solute at 5% is compared with the phase diagram of the original solute free model. Model properties are investigated both through Monte Carlo simulations and a cluster approximation. The model solubility is computed via simulations and is shown to exhibit a minimum as a function of temperature. The line of minimum solubility (TmS) coincides with the line of maximum density (TMD) for different solvent chemical potentials, in accordance with the literature on continuous realistic models and on the "cavity" picture. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4743635]

THE FIRST DETERMINATION OF THE VISCOSITY PARAMETER IN THE CIRCUMSTELLAR DISK OF A Be STAR

Be stars possess gaseous circumstellar decretion disks, which are well described using standard alpha-disk theory. The Be star 28 CMa recently underwent a long outburst followed by a long period of quiescence, during which the disk dissipated. Here we present the first time-dependent models of the dissipation of a viscous decretion disk. By modeling the rate of decline of the V-band excess, we determine that the viscosity parameter alpha = 1.0 +/- 0.2, corresponding to a mass injection rate (M) over dot = (3.5 +/- 1.3) x 10(-8) M-circle dot yr(-1). Such a large value of a suggests that the origin of the turbulent viscosity is an instability in the disk whose growth is limited by shock dissipation. The mass injection rate is more than an order of magnitude larger than the wind mass-loss rate inferred from UV observations, implying that the mass injection mechanism most likely is not the stellar wind, but some other mechanism.

The effects of the pressing step on the microstructure and aging of NdFeB bonded magnets

The effects of the compaction step on the (micro)structural features and aging behavior of polymer coated NdFeB-based bonded magnets is reported. Due to the fracture of the material during pressing, it is estimated an increase of at least 14% in the particles' area which is not coated. Such uncoated surfaces, when exposed to the environment, reduce the magnetic performance of the magnets aged/cured in air by 19% in the conditions evaluated in this investigation. Furthermore, XRD results interpreted by Rietveld analyses show a lattice parameter change in the tetragonal structure of the hard magnetic phase after pressing. Such change varies as a function of the height of the compacted part and it is ascribed to macro-elastic stress arising from the pressure distribution in the magnet. An aging/curing step during 24 h is able to relief such macro-elastic stress. (C) 2012 Elsevier B.V. All rights reserved.

Minimal Landau background gauge on the lattice

We present the first numerical implementation of the minimal Landau background gauge for Yang-Mills theory on the lattice. Our approach is a simple generalization of the usual minimal Landau gauge and is formulated for the general SU(N) gauge group. We also report on preliminary tests of the method in the four-dimensional SU(2) case, using different background fields. Our tests show that the convergence of the numerical minimization process is comparable to the case of a null background. The uniqueness of the minimizing functional employed is briefly discussed.

Binary and ternary inclusion complexes of dapsone in cyclodextrins and polymers: preparation, characterization and evaluation

Dapsone (DAP) is a synthetic sulfone drug with bacteriostatic activity, mainly against Mycobacterium leprae. In this study we have investigated the interactions of DAP with cyclodextrins, 2-hydroxypropyl-beta-cyclodextrin (HP beta CD) and beta-cyclodextrin (beta CD), in the presence and absence of water-soluble polymers, in order to improve its solubility and bioavailability. Solid systems DAP/HP beta CD and DAP/beta CD, in the presence or absence of polyvinylpyrrolidone (PVP K30) or hydroxypropyl methylcellulose (HPMC), were prepared. The binary and ternary systems were evaluated and characterized by SEM, DSC, XRD and NMR analysis as well as phase solubility assays, in order to investigate the interactions between DAP and the excipients in aqueous solution. This study revealed that inclusion complexes of DAP and cyclodextrins (HP beta CD and beta CD) can be produced in order to improve DAP solubility and bioavailability in the presence or absence of polymers (PVP K30 and HPMC). The more stable inclusion complex was obtained with HP beta CD, and consequently HP beta CD was more efficient in improving DAP solubility than beta CD, and the addition of polymers had no influence on DAP solubility or on the stability of the DAP/CDs complexes.

No-pole condition in Landau gauge: Properties of the Gribov ghost form factor and a constraint on the 2d gluon propagator

We study general properties of the Landau-gauge Gribov ghost form factor sigma(p(2)) for SU(N-c) Yang-Mills theories in the d-dimensional case. We find a qualitatively different behavior for d = 3, 4 with respect to the d = 2 case. In particular, considering any (sufficiently regular) gluon propagator D(p(2)) and the one-loop-corrected ghost propagator, we prove in the 2d case that the function sigma(p(2)) blows up in the infrared limit p -> 0 as -D(0) ln(p(2)). Thus, for d = 2, the no-pole condition sigma(p(2)) < 1 (for p(2) > 0) can be satisfied only if the gluon propagator vanishes at zero momentum, that is, D(0) = 0. On the contrary, in d = 3 and 4, sigma(p(2)) is finite also if D(0) > 0. The same results are obtained by evaluating the ghost propagator G(p(2)) explicitly at one loop, using fitting forms for D(p(2)) that describe well the numerical data of the gluon propagator in two, three and four space-time dimensions in the SU(2) case. These evaluations also show that, if one considers the coupling constant g(2) as a free parameter, the ghost propagator admits a one-parameter family of behaviors (labeled by g(2)), in agreement with previous works by Boucaud et al. In this case the condition sigma(0) <= 1 implies g(2) <= g(c)(2), where g(c)(2) is a "critical" value. Moreover, a freelike ghost propagator in the infrared limit is obtained for any value of g(2) smaller than g(c)(2), while for g(2) = g(c)(2) one finds an infrared-enhanced ghost propagator. Finally, we analyze the Dyson-Schwinger equation for sigma(p(2)) and show that, for infrared-finite ghost-gluon vertices, one can bound the ghost form factor sigma(p(2)). Using these bounds we find again that only in the d = 2 case does one need to impose D(0) = 0 in order to satisfy the no-pole condition. The d = 2 result is also supported by an analysis of the Dyson-Schwinger equation using a spectral representation for the ghost propagator. Thus, if the no-pole condition is imposed, solving the d = 2 Dyson-Schwinger equations cannot lead to a massive behavior for the gluon propagator. These results apply to any Gribov copy inside the so-called first Gribov horizon; i.e., the 2d result D(0) = 0 is not affected by Gribov noise. These findings are also in agreement with lattice data.