Spline regression for zero-inflated models

We propose a regression model for count data when the classical generalized linear model approach is too rigid due to a high outcome of zero counts and a nonlinear influence of continuous covariates. Zero-Inflation is applied to take into account the presence of excess zeros with separate link functions for the zero and the nonzero component. Nonlinearity in covariates is captured by spline functions based on B-splines. Our algorithm relies on maximum-likelihood estimation and allows for adaptive box-constrained knots, thus improving the goodness of the spline fit and allowing for detection of sensitivity changepoints. A simulation study substantiates the numerical stability of the algorithm to infer such models. The AIC criterion is shown to serve well for model selection, in particular if nonlinearities are weak such that BIC tends to overly simplistic models. We fit the introduced models to real data of children's dental sanity, linking caries counts with the so-called Body-Mass-Index (BMI) and other socioeconomic factors. This reveals a puzzling nonmonotonic influence of BMI on caries counts which is yet to be explained by clinical experts

Bias-dependent Contact Resistance in Rubrene Single-Crystal Field-Effect Transistors

We report a systematic study of the bias-dependent contact resistance in rubrene single-crystal field-effect transistors with Ni, Co, Cu, Au, and Pt electrodes. We show that the reproducibility in the values of contact resistance strongly depends on the metal, ranging from a factor of two for Ni to more than three orders of magnitude for Au. Surprisingly, FETs with Ni, Co, and Cu contacts exhibits an unexpected reproducibility of the bias-dependent differential conductance of the contacts, once this has been normalized to the value measured at zero bias. This reproducibility may enable the study of microscopic carrier injection processes into organic semiconductors.Comment: 4 pages, 4 figure

A possible observational bias in the estimation of the virial parameter in virialized clumps

The dynamics of massive clumps, the environment where massive stars originate, is still unclear. Many theories predict that these regions are in a state of near-virial equilibrium, or near energy equi-partition, while others predict that clumps are in a sub-virial state. Observationally, the majority of the massive clumps are in a sub-virial state with a clear anti-correlation between the virial parameter $\alpha_{vir}$ and the mass of the clumps $M_{c}$, which suggests that the more massive objects are also the more gravitationally bound. Although this trend is observed at all scales, from massive clouds down to star-forming cores, theories do not predict it. In this work we show how, starting from virialized clumps, an observational bias is introduced in the specific case where the kinetic and the gravitational energies are estimated in different volumes within clumps and how it can contribute to the spurious $\alpha_{vir}-M_{c}$ anti-correlation in these data. As a result, the observed effective virial parameter $\tilde{\alpha}_{eff}<\alpha_{vir}$, and in some circumstances it might not be representative of the virial state of the observed clumps.Comment: A&A letter, accepte

Massive 70 micron quiet clumps I: evidence of embedded low/intermediate-mass star formation activity

Massive clumps, prior to the formation of any visible protostars, are the best candidates to search for the elusive massive starless cores. In this work we investigate the dust and gas properties of massive clumps selected to be 70 micron quiet, therefore good starless candidates. Our sample of 18 clumps has masses 300 < M < 3000 M_sun, radius 0.54 < R < 1.00 pc, surface densities Sigma > 0.05 g cm^-2 and luminosity/mass ratio L/M < 0.3. We show that half of these 70 micron quiet clumps embed faint 24 micron sources. Comparison with GLIMPSE counterparts shows that 5 clumps embed young stars of intermediate stellar mass up to ~5.5 M_sun. We study the clump dynamics with observations of N2H+ (1-0), HNC (1-0) and HCO+ (1-0) made with the IRAM 30m telescope. Seven clumps have blue-shifted spectra compatible with infall signatures, for which we estimate a mass accretion rate 0.04 < M_dot < 2.0 x 10^-3 M_sun yr^-1, comparable with values found in high-mass protostellar regions, and free-fall time of the order of t_ff = 3 x 10^5 yr. The only appreciable difference we find between objects with and without embedded 24 micron sources is that the infall rate appears to increase from 24 micron dark to 24 micron bright objects. We conclude that all 70 micron quiet objects have similar properties on clump scales, independently of the presence of an embedded protostar. Based on our data we speculate that the majority, if not all of these clumps may already embed faint, low-mass protostellar cores. If these clumps are to form massive stars, this must occur after the formation of these lower mass stars.Comment: 44 pages, 11 Figures. Accepted for publication in MNRA

Spitzer-IRAC GLIMPSE of high mass protostellar objects. I Infrared point sources and nebulae

The GLIMPSE archive was used to obtain 3.6--8.0micron, point source photometry and images for 381 massive protostellar candidates lying in the Galactic mid-plane. The colours, magnitudes and spectral indicies of sources in each of the 381 target fields were analysed and compared with the predictions of 2D radiative transfer model simulations. Although no discernable embedded clusters were found in any targets, multiple sources or associations of redenned young stellar objects were found in many sources indicating multiplicity at birth. The spectral index ($\alpha$) of these point sources in 3.6--8.0mum bands display large values of $\alpha$=2--5. A color-magnitude analog plot was used to identify 79 infrared counterparts to the HMPOs. Compact nebulae are found in 75% of the detected sources with morphologies that can be well described by core-halo, cometary, shell-like and bipolar geometries similar to those observed in ultra-compact HII regions. The IRAC band SEDs of the IR counterparts of HMPOs are best described to represent YSOs with a mass range of 8--20\msun in their Class I stages when compared with 2D radiative transfer models. They also suggest that the high $\alpha$ values represent reprocessed star/star+disk emission that is arising in the dense envelopes. Thus we are witnessing the luminous envelopes around the protostars rather than their photospheres or disks. We argue that the compact infrared nebulae likely reflect the underlying physical structure of the dense cores and are found to imitate the morphologies of known UCHII regions. Our results favour models of continuuing accretion involving both molecular and ionised accretion components to build the most massive stars rather than purely molecular rapid accretion flows.Comment: 13 pages, 7 figures, accepted by A&

The initial conditions of stellar protocluster formation. II. A catalogue of starless and protostellar clumps embedded in IRDCs in the Galactic longitude range 15<l<55

We present a catalogue of starless and protostellar clumps associated with infrared dark clouds (IRDCs) in a 40 degrees wide region of the inner Galactic Plane (b<1). We have extracted the far-infrared (FIR) counterparts of 3493 IRDCs with known distance in the Galactic longitude range 15<l<55 and searched for the young clumps using Hi-GAL, the survey of the Galactic Plane carried out with the Herschel satellite. Each clump is identified as a compact source detected at 160, 250 and 350 mum. The clumps have been classified as protostellar or starless, based on their emission (or lack of emission) at 70 mum. We identify 1723 clumps, 1056 (61%) of which are protostellar and 667 (39%) starless. These clumps are found within 764 different IRDCs, 375 (49%) of which are only associated with protostellar clumps, 178 (23%) only with starless clumps, and 211 (28%) with both categories of clumps. The clumps have a median mass of 250 M_sun and range up to >10^4$ M_sun in mass and up to 10^5 L_sun in luminosity. The mass-radius distribution shows that almost 30% of the starless clumps identified in this survey could form high-mass stars, however these massive clumps are confined in only ~4% of the IRDCs. Assuming a minimum mass surface density threshold for the formation of high-mass stars, the comparison of the numbers of massive starless clumps and those already containing embedded sources suggests an upper limit lifetime for the starless phase of 10^5 years for clumps with a mass M>500 M_sun.Comment: accepted for publication in MNRAS. Online catalogues available soon, please contact the authors if intereste

Effects of Chia Seed Supplementation on Biochemical Markers of Cardiometabolic Diseases in Spontaneously Hypertensive Rats

Some studies suggested a positive effect against cardiometabolic diseases of supplementation of alpha-linolenic acid (ALA, C18:3 n-3) rich foods in pathological subjects, even if the total literature is controversial. In order to clarify ALA-rich chia seed action in hypertensive model with the overt pathology and without drug interference, in the present study the biochemical markers of cardiometabolic diseases (endothelin-1, ET-1; nitric oxide, NO; and bradykinin, BK) in Spontaneously Hypertensive Rats (SHRs) were analysed after 5% chia seed dietary supplementation for five weeks, and compared with the staple raw material wheat and corn. At the end of the experimental period, also plasma antioxidant capacity and inflammatory condition were evaluated. Our results showed that the chia seed group was more oxidized. On the other hand, ET-1 significantly decreased in chia seed group, and there was no difference between groups for NO, BK, and the inflammatory C-reactive protein (CRP). In conclusion, some positive effects of chia seed consumption on cardiometabolic markers in SHRs were observed, despite this the association of chia seeds with antioxidants is suggested to avoid plasma oxidation increase

Enumeration of many-body skeleton diagrams

The many-body dynamics of interacting electrons in condensed matter and quantum chemistry is often studied at the quasiparticle level, where the perturbative diagrammatic series is partially resummed. Based on Hedin's equations for self-energy, polarization, propagator, effective potential, and vertex function in zero dimension of space-time, dressed Feynman (skeleton) diagrams are enumerated. Such diagram counts provide useful basic checks for extensions of the theory for future realistic simulations.Comment: 5 pages including 4 figure