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.

Reducing the parameter space for unparticle-inspired models using white dwarf masses

Based on astrophysical constraints derived from Chandrasekhar's mass limit for white dwarfs, we study the effects of the model on the parameters of unparticle-inspired gravity, on scales Lambda(U) > 1 TeV and d(U) approximate to 1.

Charcot foot: Skin temperature as a good clinical parameter for predicting disease outcome

Twenty-eight diabetics presenting with acute Charcot foot were immobilized and the temperature difference between limbs measured at each month. All patients had monthly follow-up visits for a year and the relapse rate was zero. We found that skin temperature is a good parameter to ensure safe immobilization withdrawal. (C) 2011 Elsevier Ireland Ltd. All rights reserved.

On the energy conversion efficiency in magnetic hyperthermia applications: A new perspective to analyze the departure from the linear regime

The success of magnetic hyperthermia cancer treatments rely strongly on the magnetic properties of the nanoparticles and their intricate dependence on the externally applied field. This is particularly more so as the response departs from the low field linear regime. In this paper we introduce a new parameter, referred to as the efficiency in converting electromagnetic energy into thermal energy, which is shown to be remarkably useful in the analysis of the system response, especially when the power loss is investigated as a function of the applied field amplitude. Using numerical simulations of dynamic hysteresis, through the stochastic Landau-Lifshitz model, we map in detail the efficiency as a function of all relevant parameters of the system and compare the results with simple-yet powerful-predictions based on heuristic arguments about the relaxation time. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4705392]

Transformer Operation at Deep Saturation: Model and Parameter Determination

This paper proposes a model that adequately describes the operation of the transformer at deep saturation, suitable for power-electronics applications, and a method for determining its parameters. Simulation and experimental results are presented to confirm the validity of the model and the method.

Linear parameter varying control design for rotating systems supported by journal bearings

Linear parameter varying (LPV) control is a model-based control technique that takes into account time-varying parameters of the plant. In the case of rotating systems supported by lubricated bearings, the dynamic characteristics of the bearings change in time as a function of the rotating speed. Hence, LPV control can tackle the problem of run-up and run-down operational conditions when dynamic characteristics of the rotating system change significantly in time due to the bearings and high vibration levels occur. In this work, the LPV control design for a flexible shaft supported by plain journal bearings is presented. The model used in the LPV control design is updated from unbalance response experimental results and dynamic coefficients for the entire range of rotating speeds are obtained by numerical optimization. Experimental implementation of the designed LPV control resulted in strong reduction of vibration amplitudes when crossing the critical speed, without affecting system behavior in sub- or supercritical speeds. (C) 2012 Elsevier Ltd. 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.

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.

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.

Mean-Field critical behavior and ergodicity break in a nonequilibrium one-dimensional RSOS growth model

We investigate the nonequilibrium roughening transition of a one-dimensional restricted solid-on-solid model by directly sampling the stationary probability density of a suitable order parameter as the surface adsorption rate varies. The shapes of the probability density histograms suggest a typical Ginzburg-Landau scenario for the phase transition of the model, and estimates of the "magnetic" exponent seem to confirm its mean-field critical behavior. We also found that the flipping times between the metastable phases of the model scale exponentially with the system size, signaling the breaking of ergodicity in the thermodynamic limit. Incidentally, we discovered that a closely related model not considered before also displays a phase transition with the same critical behavior as the original model. Our results support the usefulness of off-critical histogram techniques in the investigation of nonequilibrium phase transitions. We also briefly discuss in the appendix a good and simple pseudo-random number generator used in our simulations.

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.

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.

Deviation from quark number scaling of the anisotropy parameter v(2) of pions, kaons, and protons in Au+Au collisions at root s(NN)=200 GeV

Measurements of the anisotropy parameter v(2) of identified hadrons (pions, kaons, and protons) as a function of centrality, transverse momentum p(T), and transverse kinetic energy KET at midrapidity (vertical bar eta vertical bar < 0.35) in Au + Au collisions at root s(N N) = 200 GeV are presented. Pions and protons are identified up to p(T) = 6 GeV/c, and kaons up to p(T) = 4 GeV/c, by combining information from time-of-flight and aerogel Cerenkov detectors in the PHENIX Experiment. The scaling of v(2) with the number of valence quarks (n(q)) has been studied in different centrality bins as a function of transverse momentum and transverse kinetic energy. A deviation from previously observed quark-number scaling is observed at large values of KET/n(q) in noncentral Au + Au collisions (20-60%), but this scaling remains valid in central collisions (0-10%).

Constraining the dark energy and smoothness parameter with type Ia supernovae and gamma-ray bursts

The existence of inhomogeneities in the observed Universe modifies the distance-redshift relations thereby affecting the results of cosmological tests in comparison to the ones derived assuming spatially uniform models. By modeling the inhomogeneities through a Zeldovich-Kantowski-Dyer-Roeder approach which is phenomenologically characterized by a smoothness parameter alpha, we rediscuss the constraints on the cosmic parameters based on type Ia supernovae (SNe Ia) and gamma-ray bursts (GRBs) data. The present analysis is restricted to a flat Lambda CDM model with the reasonable assumption that Lambda does not clump. A chi(2) analysis using 557 SNe Ia data from the Union2 compilation data (R. Amanullah et al., Astrophys. J. 716, 712 (2010).) constrains the pair of parameters (Omega(m), alpha) to Omega(m) = 0.27(-0.03)(+0.08) (2 sigma) and alpha >= 0.25. A similar analysis based only on 59 Hymnium GRBs (H. Wei, J. Cosmol. Astropart. Phys. 08 (2010) 020.) constrains the matter density parameter to be Omega(m) = 0.35(-0.24)(+0.62) (2 sigma) while all values for the smoothness parameter are allowed. By performing a joint analysis, it is found that Omega(m) = 0.27(-0.06)(+0.06) and alpha >= 0.52. As a general result, although considering that current GRB data alone cannot constrain the smoothness alpha parameter, our analysis provides an interesting cosmological probe for dark energy even in the presence of inhomogeneities.

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.