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.

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

Low-lying singlet and triplet electronic states of RhB

The low-lying X-1 Sigma(+), a(3)Delta, A(1)Delta, b(3)Sigma(+), B-1 Pi, c(3)Pi, C-1 Phi, D-1 Sigma(+), E-1 Pi, d(1)Phi, and e(3)Pi electronic states of RhB have been investigated at the ab initio level, using the multistate multiconfigurational second-order perturbation (MS-CASPT2) theory, with extended atomic basis sets and inclusion of scalar relativistic effects. Among the eleven electronic states included in this work, only three (the X-1 Sigma(+), D-1 Sigma(+), and E-1 Pi states) have been investigated experimentally. Potential energy curves, spectroscopic constants, dipole moments, binding energies, and chemical bonding aspects are presented for all electronic states.

A theoretical study of the binding and electronic spectrum of the Mo-2 molecule

Multiconfigurational SCF and second-order perturbation theory have been employed to study seven low-lying singlet and triplet electronic states of the Mo-2 molecule. The bond order of the ground state has been analyzed based on the effective bond order (EBO), indicating that a fully developed sextuple bond is formed between the two Mo atoms. The experimentally observed excited states a(3)Sigma(+)(u) and A(1)Sigma(+)(u) have been determined and the so-called (3)Lambda excited state identified as the b(3)Sigma(+)(u) state, in agreement with experimental expectations. (C) 2007 Elsevier B.V. All rights reserved.

A CASSCF/MRCI Study of the Low-Lying Electronic States of the BeS Molecule

Accurate potential energy curves, dipole moment functions, dissociation energies, and molecular constants for several low-lying singlet and triplet electronic states of BeS were investigated using the CASSCF/MRCI methodology, and the cc-pV5Z basis set for beryllium, and the aug-cc-pV(5+d) Z set for sulfur. Besides presenting improved results for the three lowest lying states, this study presents the first theoretical characterization of another set of nine excited states so far unknown experimentally. Our results are sufficiently accurate to reliably guide the experimental search and characterization of these states, and also to confirm the experimental assignment of the B(1)Sigma(+) excited state. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 1694-1700, 2011

Dynamic portrait of the planetary 2/1 mean-motion resonance - I. Systems with a more massive outer planet

We analyse the global structure of the phase space of the planar planetary 2/1 mean-motion resonance in cases where the outer planet is more massive than its inner companion. Inside the resonant domain, we show the existence of two families of periodic orbits, one associated to the librational motion of resonant angle (sigma-family) and the other related to the circulatory motion of the difference in longitudes of pericentre (Delta pi-family). The well-known apsidal corotation resonances (ACR) appear as intersections between both families. A complex web of secondary resonances is also detected for low eccentricities, whose strengths and positions are dependent on the individual masses and spatial scale of the system. The construction of dynamical maps for various values of the total angular momentum shows the evolution of the families of stable motion with the eccentricities, identifying possible configurations suitable for exoplanetary systems. For low-moderate eccentricities, several different stable modes exist outside the ACR. For larger eccentricities, however, all stable solutions are associated to oscillations around the stationary solutions. Finally, we present a possible link between these stable families and the process of resonance capture, identifying the most probable routes from the secular region to the resonant domain, and discussing how the final resonant configuration may be affected by the extension of the chaotic layer around the resonance region.

DETECTABILITY AND ERROR ESTIMATION IN ORBITAL FITS OF RESONANT EXTRASOLAR PLANETS

We estimate the conditions for detectability of two planets in a 2/1 mean-motion resonance from radial velocity data, as a function of their masses, number of observations and the signal-to-noise ratio. Even for a data set of the order of 100 observations and standard deviations of the order of a few meters per second, we find that Jovian-size resonant planets are difficult to detect if the masses of the planets differ by a factor larger than similar to 4. This is consistent with the present population of real exosystems in the 2/1 commensurability, most of which have resonant pairs with similar minimum masses, and could indicate that many other resonant systems exist, but are currently beyond the detectability limit. Furthermore, we analyze the error distribution in masses and orbital elements of orbital fits from synthetic data sets for resonant planets in the 2/1 commensurability. For various mass ratios and number of data points we find that the eccentricity of the outer planet is systematically overestimated, although the inner planet`s eccentricity suffers a much smaller effect. If the initial conditions correspond to small-amplitude oscillations around stable apsidal corotation resonances, the amplitudes estimated from the orbital fits are biased toward larger amplitudes, in accordance to results found in real resonant extrasolar systems.

Alkaline phosphatase as a reporter of sigma(S) levels and rpoS polymorphisms in different E-coli strains

sigma(S) is responsible for the transcriptional regulation of genes related to protection against stresses and bacterial survival and it accumulates in the cell under conditions of stress, such as nutrient limitation. An increase in the levels of sigma(S) causes a reduction in the expression of genes that are transcribed by RNA polymerase associated with the principal sigma factor, sigma(70). phoA, that encodes alkaline phosphatase (AP) is expressed under phosphate shortage conditions, and is also repressed by sigma(S). Here we show that in a Pi-limited chemostat, accumulation of rpoS mutations is proportional to the intrinsic level of sigma(S) in the cells. Acquisition of mutations in rpoS relieves repression of the PHO genes. We also devised a non-destructive method based on the rpoS effect on AP that differentiates between rpo(S+) and rpoS mutants, as well as between high and low-sigma(S) producers. Using this method, we provide evidence that sigma(S) contributes to the repression of AP under conditions of Pi excess and that AP variation among different strains is at least partly due to intrinsic variation in sigma(S) levels. Consequently, a simple and non-destructive AP assay can be employed to differentiate between strains expressing different levels of sigma(S) on agar plates.

SIGMA THEORY AND TWISTED CONJUGACY CLASSES

Using Sigma theory we show that for large classes of groups G there is a subgroup H of finite index in Aut(G) such that for phi is an element of H the Reidemeister number R(phi) is infinite. This includes all finitely generated nonpolycyclic groups G that fall into one of the following classes: nilpotent-by-abelian groups of type FP(infinity); groups G/G `` of finite Prufer rank; groups G of type FP(2) without free nonabelian subgroups and with nonpolycyclic maximal metabelian quotient; some direct products of groups; or the pure symmetric automorphism group. Using a different argument we show that the result also holds for 1-ended nonabelian nonsurface limit groups. In some cases, such as with the generalized Thompson`s groups F(n,0) and their finite direct products, H = Aut(G).

(1)H NMR spectroscopy as a tool to determine accurate photoisomerization quantum yields of stilbene-like ligands coordinated to rhenium(I) polypyridyl complexes

In this work, the use of proton nuclear magnetic resonance, (1)H NMR, was fully described as a powerful tool to follow a photoreaction and to determine accurate quantum yields, so called true quantum yields (Phi(true)), when a reactant and photoproduct absorption overlap. For this, Phi(true) for the trans-cis photoisomerization process were determined for rhenium(I) polypyridyl complexes, fac-[Re(CO)(3)(NN)(trans-L)](+) (NN = 1,10-phenanthroline, phen, or 4,7-diphenyl-1,10-phenanthroline, ph(2)phen, and L = 1,2-bis(4-pyridyl) ethylene, bpe, or 4-styrylpyridine, stpy). The true values determined at 365 nm irradiation (e. g. Phi(NMR) = 0.80 for fac-[Re(CO)(3)(phen)(trans-bpe)](+)) were much higher than those determined by absorption spectral changes (Phi(UV-Vis) = 0.39 for fac-[Re(CO)(3)(phen)(trans-bpe)](+)). Phi(NMR) are more accurate in these cases due to the distinct proton signals of trans and cis-isomers, which allow the actual determination of each component concentration under given irradiation time. Nevertheless when the photoproduct or reactant contribution at the probe wavelength is negligible, one can determine Phi(true) by regular absorption spectral changes. For instance, Phi(313) nm for free ligand photoisomerization determined both by absorption and (1)H NMR variation are equal within the experimental error (bpe: Phi(UV-Vis) = 0.27, Phi(NMR) = 0.26; stpy: Phi(UV-Vis) = 0.49, Phi(NMR) = 0.49). Moreover, (1)H NMR data combined with electronic spectra allowed molar absorptivity determination of difficult to isolate cis-complexes. (C) 2009 Elsevier B. V. All rights reserved.