Semiparametric Bayesian measurement error modeling

This work presents a Bayesian semiparametric approach for dealing with regression models where the covariate is measured with error. Given that (1) the error normality assumption is very restrictive, and (2) assuming a specific elliptical distribution for errors (Student-t for example), may be somewhat presumptuous; there is need for more flexible methods, in terms of assuming only symmetry of errors (admitting unknown kurtosis). In this sense, the main advantage of this extended Bayesian approach is the possibility of considering generalizations of the elliptical family of models by using Dirichlet process priors in dependent and independent situations. Conditional posterior distributions are implemented, allowing the use of Markov Chain Monte Carlo (MCMC), to generate the posterior distributions. An interesting result shown is that the Dirichlet process prior is not updated in the case of the dependent elliptical model. Furthermore, an analysis of a real data set is reported to illustrate the usefulness of our approach, in dealing with outliers. Finally, semiparametric proposed models and parametric normal model are compared, graphically with the posterior distribution density of the coefficients. (C) 2009 Elsevier Inc. All rights reserved.

CAN OLD GALAXIES AT HIGH REDSHIFTS AND BARYON ACOUSTIC OSCILLATIONS CONSTRAIN H(0)?

A new age-redshift test is proposed in order to constrain H(0) on the basis of the existence of old high-redshift galaxies (OHRGs). In the flat Lambda cold dark matter model, the value of H(0) is heavily dependent on the mass density parameter Omega(M) = 1- Omega(Lambda). Such a degeneracy can be broken through a joint analysis involving the OHRG and baryon acoustic oscillation signature. By assuming a galaxy incubation time, t(inc) = 0.8 +/- 0.4 Gyr, our joint analysis yields a value of H(0) = 71 +/- 4 km s(-1) Mpc(-1) (1 sigma) with the best-fit density parameter Omega(M) = 0.27 +/- 0.03. Such results are in good agreement with independent studies from the Hubble Space Telescope key project and recent estimates of the Wilkinson Microwave Anisotropy Probe, thereby suggesting that the combination of these two independent phenomena provides an interesting method to constrain the Hubble constant.

Cosmic distance duality relation and the shape of galaxy clusters

Context. Observations in the cosmological domain are heavily dependent on the validity of the cosmic distance-duality (DD) relation, eta = D(L)(z)(1+ z)(2)/D(A)(z) = 1, an exact result required by the Etherington reciprocity theorem where D(L)(z) and D(A)(z) are, respectively, the luminosity and angular diameter distances. In the limit of very small redshifts D(A)(z) = D(L)(z) and this ratio is trivially satisfied. Measurements of Sunyaev-Zeldovich effect (SZE) and X-rays combined with the DD relation have been used to determine D(A)(z) from galaxy clusters. This combination offers the possibility of testing the validity of the DD relation, as well as determining which physical processes occur in galaxy clusters via their shapes. Aims. We use WMAP (7 years) results by fixing the conventional Lambda CDM model to verify the consistence between the validity of DD relation and different assumptions about galaxy cluster geometries usually adopted in the literature. Methods. We assume that. is a function of the redshift parametrized by two different relations: eta(z) = 1+eta(0)z, and eta(z) = 1+eta(0)z/(1+z), where eta(0) is a constant parameter quantifying the possible departure from the strict validity of the DD relation. In order to determine the probability density function (PDF) of eta(0), we consider the angular diameter distances from galaxy clusters recently studied by two different groups by assuming elliptical (isothermal) and spherical (non-isothermal) beta models. The strict validity of the DD relation will occur only if the maximum value of eta(0) PDF is centered on eta(0) = 0. Results. It was found that the elliptical beta model is in good agreement with the data, showing no violation of the DD relation (PDF peaked close to eta(0) = 0 at 1 sigma), while the spherical (non-isothermal) one is only marginally compatible at 3 sigma. Conclusions. The present results derived by combining the SZE and X-ray surface brightness data from galaxy clusters with the latest WMAP results (7-years) favors the elliptical geometry for galaxy clusters. It is remarkable that a local property like the geometry of galaxy clusters might be constrained by a global argument provided by the cosmic DD relation.

Kinematics of galaxies in compact groups Studying the B-band Tully-Fischer relation

We obtained new Fabry-Perot data cubes and derived velocity fields, monochromatic, and velocity dispersion maps for 28 galaxies in the Hickson compact groups 37, 40, 47, 49, 54, 56, 68, 79, and 93. We also derived rotation curves for 9 of the studied galaxies, 6 of which are strongly asymmetric. Combining these new data with previously published 2D kinematic maps of compact group galaxies, we investigated the differences between the kinematic and morphological position angles for a sample of 46 galaxies. We find that one third of the unbarred compact group galaxies have position angle misalignments between the stellar and gaseous components. This and the asymmetric rotation curves are clear signatures of kinematic perturbations, probably because of interactions among compact group galaxies. A comparison between the B-band Tully-Fisher relation for compact group galaxies and for the GHASP field-galaxy sample shows that, despite the high fraction of compact group galaxies with asymmetric rotation curves, these lay on the TF relation defined by galaxies in less dense environments, although with more scatter. This agrees with previous results, but now confirmed for a larger sample of 41 galaxies. We confirm the tendency for compact group galaxies at the low-mass end of the Tully-Fisher relation (HCG 49b, 89d, 96c, 96d, and 100c) to have either a magnitude that is too bright for its mass (suggesting brightening by star formation) and/or a low maximum rotational velocity for its luminosity (suggesting tidal stripping). These galaxies are outside the Tully Fisher relation at the 1 sigma level, even when the minimum acceptable values of inclinations are used to compute their maximum velocities. Including such galaxies with nu < 100 km s(-1) in the determination of the zero point and slope of the compact group B-band Tully-Fisher relation would strongly change the fit, making it different from the relation for field galaxies, which has to be kept in mind when studying scaling relations of interacting galaxies, especially at high redshifts.

Constraints on cold dark matter accelerating cosmologies and cluster formation

We discuss the properties of homogeneous and isotropic flat cosmologies in which the present accelerating stage is powered only by the gravitationally induced creation of cold dark matter (CCDM) particles (Omega(m) = 1). For some matter creation rates proposed in the literature, we show that the main cosmological functions such as the scale factor of the universe, the Hubble expansion rate, the growth factor, and the cluster formation rate are analytically defined. The best CCDM scenario has only one free parameter and our joint analysis involving baryonic acoustic oscillations + cosmic microwave background (CMB) + SNe Ia data yields (Omega) over tilde = 0.28 +/- 0.01 (1 sigma), where (Omega) over tilde (m) is the observed matter density parameter. In particular, this implies that the model has no dark energy but the part of the matter that is effectively clustering is in good agreement with the latest determinations from the large- scale structure. The growth of perturbation and the formation of galaxy clusters in such scenarios are also investigated. Despite the fact that both scenarios may share the same Hubble expansion, we find that matter creation cosmologies predict stronger small scale dynamics which implies a faster growth rate of perturbations with respect to the usual Lambda CDM cosmology. Such results point to the possibility of a crucial observational test confronting CCDM with Lambda CDM scenarios through a more detailed analysis involving CMB, weak lensing, as well as the large-scale structure.

Kinematic constraints to the transition redshift from supernovae type Ia union data

The kinematic approach to cosmological tests provides direct evidence to the present accelerating stage of the Universe that does not depend on the validity of general relativity, as well as on the matter-energy content of the Universe. In this context, we consider here a linear two-parameter expansion for the decelerating parameter, q(z)=q(0)+q(1)z, where q(0) and q(1) are arbitrary constants to be constrained by the union supernovae data. By assuming a flat Universe we find that the best fit to the pair of free parameters is (q(0),q(1))=(-0.73,1.5) whereas the transition redshift is z(t)=0.49(-0.07)(+0.14)(1 sigma) +0.54-0.12(2 sigma). This kinematic result is in agreement with some independent analyses and more easily accommodates many dynamical flat models (like Lambda CDM).

Combining Monte Carlo simulation and density-functional theory to describe the spectral changes of Na(2) in liquid helium

Spectral changes of Na(2) in liquid helium were studied using the sequential Monte Carlo-quantum mechanics method. Configurations composed by Na(2) surrounded by explicit helium atoms sampled from the Monte Carlo simulation were submitted to time-dependent density-functional theory calculations of the electronic absorption spectrum using different functionals. Attention is given to both line shift and line broadening. The Perdew, Burke, and Ernzerhof (PBE1PBE, also known as PBE0) functional, with the PBE1PBE/6-311++G(2d,2p) basis set, gives the spectral shift, compared to gas phase, of 500 cm(-1) for the allowed X (1)Sigma(+)(g) -> B (1)Pi(u) transition, in very good agreement with the experimental value (700 cm(-1)). For comparison, cluster calculations were also performed and the first X (1)Sigma(+)(g) -> A (1)Sigma(+)(u) transition was also considered.

Effect of Femtosecond Laser Energy Level on Corneal Stromal Cell Death and Inflammation

PURPOSE: To analyze the effects of variations in femtosecond laser energy level on corneal stromal cell death. and inflammatory cell influx following flap creation in a rabbit model. METHODS: Eighteen rabbits were stratified in three different groups according to level of energy applied for flap creation (six animals per group). Three different energy levels were chosen for both the lamellar and side cut; 2.7 mu J (high energy), 1.6 mu J (intermediate energy), and 0.5 mu J (low energy) with a 60 kHz, model II, femtosecond laser (IntraLase). The opposite eye of each rabbit served as a control. At the 24-hour time point after surgery, all rabbits were euthanized and the comeoscleral rims were analyzed for the levels of cell death and inflammatory cell influx with the terminal uridine deoxynucleotidyl transferase dUTP-nick end labeling (TUNEL) assay and immunocytochemistry for monocyte marker CD11b, respectively. RESULTS: The high energy group (31.9 +/- 7.1 [standard error of mean (SEM) 2.9]) had significantly more TUNEL positive cells in the central flap compared to the intermediate (22.2 +/- 1.9 [SEM 0.8], P=.004), low (17.9 +/- 4.0 [SEM 1.6], P <= .001), and control eye (0.06 +/- 0.02 [SEM 0.009], P <= .001) groups. The intermediate and low energy groups also had significantly more TUNEL positive cells than the control groups (P <= .001). The difference between the intermediate and low energy levels was not significant (P=.56). The mean for CD11b-positive cells/400x field at the flap edge was 26.1 +/- 29.3 (SEM 11.9), 5.8 +/- 4.1 (SEM 1.6), 1.6 +/- 4.1 (SEM 1.6), and 0.005 +/- 0.01 (SEM 0.005) for high energy, intermediate energy, low energy, and control groups, respectively. Only the intermediate energy group showed statistically more inflammatory cells than control eyes (P = .015), most likely due to variability between eyes. CONCLUSIONS: Higher energy levels trigger greater cell death when the femtosecond laser is used to create corneal flaps: Greater corneal inflammatory cell infiltration is observed with higher femtosecond laser energy levels. [J Refract Surg. 2009;25:869-874.] doi:10.3928/1081597X-20090917-08

Reliability of electromyographic amplitude values of the upper limb muscles during closed kinetic chain exercises with stable and unstable surfaces

The purpose of the present study was to evaluate the intra and interday reliability of surface electromyographic amplitude values of the scapular girdle muscles and upper limbs during 3 isometric closed kinetic chain exercises, involving upper limbs with the fixed distal segment extremity on stable base of support and on a Swiss ball (relatively unstable). Twenty healthy adults performed the exercises push-up, bench-press and wall-press with different effort levels (80% and 100% maximal load). Subjects performed three maximal voluntary contractions (MVC) in muscular testing position of each muscle to obtain a reference value for root mean square (RMS) normalization. Individuals were instructed to randomly perform three isometric contraction series, in which each exercise lasted 6 s with a 2-min resting-period between series and exercises. Intra and interday reliabilities were calculated through the intraclass correlation coefficient (ICC 2.1), standard error of the measurement (SEM). Results indicated an excellent intraday reliability of electromyographic amplitude values (ICC >= 0.75). The interday reliability of normalized RMS values ranged between good and excellent (ICC 0.52-0.98). Finally, it is suggested that the reliability of normalized electromyographic amplitude values of the analyzed muscles present better values during exercises on a stable surface. However, load levels used during the exercises do not seem to have any influence on variability levels, possibly because the loads were quite similar. (C) 2007 Elsevier Ltd. All rights reserved.

Transition redshift: new kinematic constraints from supernovae

The transition redshift (deceleration/acceleration) is discussed by expanding the deceleration parameter to first order around its present value. A detailed study is carried out by considering two different parametrizations, q = q(0) + q(1)z and q = q(0) + q(1)z(1 + z)(-1), and the associated free parameters (q(0), q(1)) are constrained by three different supernovae (SNe) samples. A previous analysis by Riess et al. using the first expansion is slightly improved and confirmed in light of their recent data (Gold07 sample). However, by fitting the model with the Supernova Legacy Survey (SNLS) type Ia sample, we find that the best fit to the redshift transition is z(t) = 0.61, instead of z(t) = 0.46 as derived by the High-z Supernovae Search (HZSNS) team. This result based in the SNLS sample is also in good agreement with the sample of Davis et al., z(t) = 0.60(-0.11)(+0.28) (1 sigma). Such results are in line with some independent analyses and accommodate more easily the concordance flat model (Lambda CDM). For both parametrizations, the three SNe Ia samples considered favour recent acceleration and past deceleration with a high degree of statistical confidence level. All the kinematic results presented here depend neither on the validity of general relativity nor on the matter-energy contents of the Universe.

Bayesian analysis and constraints on kinematic models from union SNIa

The kinematic expansion history of the universe is investigated by using the 307 supernovae type Ia from the Union Compilation set. Three simple model parameterizations for the deceleration parameter ( constant, linear and abrupt transition) and two different models that are explicitly parametrized by the cosmic jerk parameter ( constant and variable) are considered. Likelihood and Bayesian analyses are employed to find best fit parameters and compare models among themselves and with the flat Lambda CDM model. Analytical expressions and estimates for the deceleration and cosmic jerk parameters today (q(0) and j(0)) and for the transition redshift (z(t)) between a past phase of cosmic deceleration to a current phase of acceleration are given. All models characterize an accelerated expansion for the universe today and largely indicate that it was decelerating in the past, having a transition redshift around 0.5. The cosmic jerk is not strongly constrained by the present supernovae data. For the most realistic kinematic models the 1 sigma confidence limits imply the following ranges of values: q(0) is an element of [-0.96, -0.46], j(0) is an element of [-3.2,-0.3] and z(t) is an element of [0.36, 0.84], which are compatible with the Lambda CDM predictions, q(0) = -0.57 +/- 0.04, j(0) = -1 and z(t) = 0.71 +/- 0.08. We find that even very simple kinematic models are equally good to describe the data compared to the concordance Lambda CDM model, and that the current observations are not powerful enough to discriminate among all of them.

SU(3) mass-splittings of heavy-baryons in QCD

We extract directly (for the first time) the charmed (C = 1) and bottom (B = -1) heavy-baryons (spin 1/2 and 3/2) mass-splittings due to SU(3) breaking using double ratios of QCD spectral sum rules (QSSR) in full QCD, which are less sensitive to the exact value and definition of the heavy quark mass, to the perturbative radiative corrections and to the QCD continuum contributions than the simple ratios commonly used for determining the heavy baryon masses. Noticing that most of the mass-splittings are mainly controlled by the ratio kappa <(S) over bars >/<(d) over bard > of the condensate, we extract this ratio, by allowing 1 sigma deviation from the observed masses of the Xi(c.b) and of the Omega(c). We obtain: kappa = 0.74(3), which improves the existing estimates: kappa = 0.70(10) from light hadrons. Using this value, we deduce M(Omega b) = 6078.5(27.4) MeV which agrees with the recent CDF data but disagrees by 2.4 sigma with the one from D0. Predictions of the Xi(Q)` and of the spectra of spin 3/2 baryons containing one or two strange quark are given in Table 2. Predictions of the hyperfine splittings Omega(Q)* - Omega(Q) and Xi(Q)* - Xi(Q) are also given in Table 3. Starting for a general choice of the interpolating currents for the spin 1/2 baryons, our analysis favours the optimal value of the mixing angle b similar or equal to (-1/5-0) found from light and non-strange heavy baryons. (C) 2010 Elsevier B.V. All rights reserved.

Electronic structure and chemical bonding in W(2) molecule

The electronic structure of the lowest-lying electronic states of W(2) were investigated at the CASPT2 level. The ground state is a X(1)Sigma(+)(g) state, followed by the a(3)Delta(u), b(3)Sigma(+)(u) and A(1)Delta(u) electronic states. Seven low-lying Omega-states were computed: (1)0(g)(+), (2)3(u), (3)2(u), (4)1(u), (5)0(u)(-), (6)1(u), and (7)2(u), with the ground state corresponding to the (1)0(g)(+)(X(1)Sigma(+)(g)) state. Comparison with the other VIB transition metal group dimers indicates a common pattern of electronic structure and spectroscopic properties. (C) 2010 Elsevier B.V. All rights reserved.

Structural and spectroscopic properties of the diazocarbene radical (CNN) and its ions CNN(+) and CNN(-): a high-level theoretical investigation

The diazocarbene radical, CNN, and the ions CNN(+) and CNN(-) were investigated at a high level of theory. Very accurate structural parameters for the states X(3)Sigma(-) and A(3)Pi of CNN, and X(2)Pi of both CNN(+) and CNN(-) were obtained with the UCCSD(T) method using correlated-consistent basis functions with extrapolations to the complete basis set limit, with valence only and also with all electrons correlated. Harmonic and anharmonic frequencies were obtained for all species and the Renner parameter and average frequencies evaluated for the Pi states. At the UCCSD(T)/CBS(T-5) level of theory, Delta(f)H(0 K) = 138.89 kcal/mol and Delta(f)H(298 K) = 139.65 kcal/mol were obtained for diazocarbene; for the ionization potential and the electron affinity of CNN, 10.969 eV (252.95 kcal/mol), and 1.743 eV (40.19 kcal/mol), respectively, are predicted. Geometry optimization was also carried out with the CASSCF/MRCI/CBS(T-5) approach for the states X(3)Sigma(-) A(3)Pi, and a(1)Delta of CNN, and with the CASSCF/MRSDCI/aug-cc-pVTZ approach for the states b(1)Sigma(+), c(1)Pi, d(1)Sigma(-), and B(3)Sigma(-), and excitation energies (T(e)) evaluated. Vertical energies were calculated for 15 electronic states, thus improving on the accuracy of the five transitions already described, and allowing for a reliable overview of a manifold of other states, which is expected to guide future spectroscopic experiments. This study corroborates the experimental assignment for the vertical transition X (3)Sigma(-) <- E (3)Pi.

Electronic Structure and Chemical Bonding in the Ground and Low-Lying Electronic States of Ta(2)

The electronic structure and chemical bonding of the ground and low-lying Lambda - S and Omega states of Ta(2) were investigated at the multiconfiguration second-order perturbation theory (CASSCF//CASPT2) level. The ground state of Ta(2) is computed to be a X(3)Sigma(-)(g) state (R(e) = 2.120 angstrom, omega(e) = 323 cm(-1), and D(e) = 4.65 eV), with two low-lying singlet states close to it (a(1) Sigma(+)(g) : T(e) = 409 cm(-1), R(e) = 2.131 angstrom, and omega(e) = 313 cm(-1); b(1) Gamma(g): T(e) = 1, 038 cm(-1), R(e) = 2.127 angstrom, and omega(e) = 316 cm(-1)). These electronic states are derived from the same electronic configuration: vertical bar 13 sigma(2)(g)14 sigma(2)(g)7 delta(2)(g)13 pi(4)(u)>. The effective bond order of the X(3) Sigma(-)(g) state is 4.52, which indicates that the Ta atoms are bound by a quintuple chemical bond. The a(1) Sigma(+)(g) state interacts strongly with the X(3)Sigma(-)(g) g ground state by a second-order spin-orbit interaction, giving rise to the (1)0(g)(+) (ground state) (dominated by the X(3)Sigma(-)(g) Lambda - S ground state) and (9)0(g)(+) (dominated by the a(1) Sigma(+)(g) Lambda - S state) Omega states. These results are in line with those reported for the group 5B homonuclear transition metal diatomics. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1306-1315, 2011