Image reconstruction in optical interferometry: Benchmarking the regularization

With the advent of infrared long-baseline interferometers with more than two telescopes, both the size and the completeness of interferometric data sets have significantly increased, allowing images based on models with no a priori assumptions to be reconstructed. Our main objective is to analyze the multiple parameters of the image reconstruction process with particular attention to the regularization term and the study of their behavior in different situations. The secondary goal is to derive practical rules for the users. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we performed multiple systematic tests, analyzing 11 regularization terms commonly used. The tests are made on different astrophysical objects, different (u,v) plane coverages and several signal-to-noise ratios to determine the minimal configuration needed to reconstruct an image. We establish a methodology and we introduce the mean-square errors (MSE) to discuss the results. From the ~24000 simulations performed for the benchmarking of image reconstruction with MiRA, we are able to classify the different regularizations in the context of the observations. We find typical values of the regularization weight. A minimal (u,v) coverage is required to reconstruct an acceptable image, whereas no limits are found for the studied values of the signal-to-noise ratio. We also show that super-resolution can be achieved with increasing performance with the (u,v) coverage filling. Using image reconstruction with a sufficient (u,v) coverage is shown to be reliable. The choice of the main parameters of the reconstruction is tightly constrained. We recommend that efforts to develop interferometric infrastructures should first concentrate on the number of telescopes to combine, and secondly on improving the accuracy and sensitivity of the arrays.Comment: 15 pages, 16 figures; accepted in A&

Testing Lorentz Invariance with Neutrinos from Ultrahigh Energy Cosmic Ray Interactions

We have previously shown that a very small amount of Lorentz invariance violation (LIV), which suppresses photomeson interactions of ultrahigh energy cosmic rays (UHECRs) with cosmic background radiation (CBR) photons, can produce a spectrum of cosmic rays that is consistent with that currently observed by the Pierre Auger Observatory (PAO) and HiRes experiments. Here, we calculate the corresponding flux of high energy neutrinos generated by the propagation of UHECR protons through the CBR in the presence of LIV. We find that LIV produces a reduction in the flux of the highest energy neutrinos and a reduction in the energy of the peak of the neutrino energy flux spectrum, both depending on the strength of the LIV. Thus, observations of the UHE neutrino spectrum provide a clear test for the existence and amount of LIV at the highest energies. We further discuss the ability of current and future proposed detectors make such observations.Comment: final version to appear in Astroparticle Physic

Complete Genome Sequences of Mycobacterium smegmatis Phages NihilNomen and Carlyle, Isolated in Las Vegas, Nevada

We present the complete genomes of the Mycobacterium smegmatis phages Carlyle and NihilNomen, isolated from soil in Las Vegas, Nevada. The phages were isolated and annotated by undergraduate students enrolled in the Phage Discovery course offered by the School of Life Sciences at the University of Nevada Las Vegas

The Dark Molecular Gas

The mass of molecular gas in an interstellar cloud is often measured using line emission from low rotational levels of CO, which are sensitive to the CO mass, and then scaling to the assumed molecular hydrogen H_2 mass. However, a significant H_2 mass may lie outside the CO region, in the outer regions of the molecular cloud where the gas phase carbon resides in C or C+. Here, H_2 self-shields or is shielded by dust from UV photodissociation, where as CO is photodissociated. This H_2 gas is "dark" in molecular transitions because of the absence of CO and other trace molecules, and because H_2 emits so weakly at temperatures 10 K < T < 100 K typical of this molecular component. This component has been indirectly observed through other tracers of mass such as gamma rays produced in cosmic ray collisions with the gas and far-infrared/submillimeter wavelength dust continuum radiation. In this paper we theoretically model this dark mass and find that the fraction of the molecular mass in this dark component is remarkably constant (~ 0.3 for average visual extinction through the cloud with mean A_V ~ 8) and insensitive to the incident ultraviolet radiation field strength, the internal density distribution, and the mass of the molecular cloud as long as mean A_V, or equivalently, the product of the average hydrogen nucleus column and the metallicity through the cloud, is constant. We also find that the dark mass fraction increases with decreasing mean A_V, since relatively more molecular H_2 material lies outside the CO region in this case.Comment: 38 page, 11 figures, Accepted for Publication in ApJ, corrected citation and typo in Appendix

A low optical depth region in the inner disk of the HerbigAe star HR5999

Circumstellar disks surrounding young stars are known to be the birthplaces of planets, and the innermost astronomical unit is of particular interest. We present new long-baseline spectro-interferometric observations of the HerbigAe star, HR5999, obtained in the H and K bands with the AMBER instrument at the VLTI, and aim to produce near-infrared images at the sub-AU spatial scale. We spatially resolve the circumstellar material and reconstruct images using the MiRA algorithm. In addition, we interpret the interferometric observations using models that assume that the near-infrared excess is dominated by the emission of a circumstellar disk. We compare the images reconstructed from the VLTI measurements to images obtained using simulated model data. The K-band image reveals three main elements: a ring-like feature located at ~0.65 AU, a low surface brightness region inside, and a central spot. At the maximum angular resolution of our observations (1.3 mas), the ring is resolved while the central spot is only marginally resolved, preventing us from revealing the exact morphology of the circumstellar environment. We suggest that the ring traces silicate condensation, i.e., an opacity change, in a circumstellar disk around HR 5999. We build a model that includes a ring at the silicate sublimation radius and an inner disk of low surface brightness responsible for a large amount of the near-infrared continuum emission. The model successfully fits the SED, visibilities, and closure phases, and provides evidence of a low surface brightness region inside the silicate sublimation radius. This study provides additional evidence that in HerbigAe stars, there is material in a low surface brightness region, probably a low optical depth region, located inside the silicate sublimation radius and of unknown nature.Comment: 11 pages, 10 figure

Observing a column-dependent zeta in dense interstellar sources: the case of the Horsehead Nebula

Context: Observations of small carbon-bearing molecules such as CCH, C4H, c-C3H2, and HCO in the Horsehead Nebula have shown these species to have higher abundances towards the edge of the source than towards the center. Aims: Given the determination of a wide range of values for zeta (s-1), the total ionization rate of hydrogen atoms, and the proposal of a column-dependent zeta(N_H), where N_H is the total column of hydrogen nuclei, we desire to determine if the effects of zeta(N_H) in a single object with spatial variation can be observable. We chose the Horsehead Nebula because of its geometry and high density. Method: We model the Horsehead Nebula as a near edge-on photon dominated region (PDR), using several choices for zeta, both constant and as a function of column. The column-dependent zeta functions are determined by a Monte Carlo model of cosmic ray penetration, using a steep power-law spectrum and accounting for ionization and magnetic field effects. We consider a case with low-metal elemental abundances as well as a sulfur-rich case. Results: We show that use of a column-dependent zeta(N_H) of 5(-15) s-1 at the surface and 7.5(-16) s-1 at Av = 10 on balance improves agreement between measured and theoretical molecular abundances, compared with constant values of zeta.Comment: 12 pages, 6 figures, 5 tables, accepted in A&

Milli-arcsecond images of the Herbig Ae star HD 163296

The very close environments of young stars are the hosts of fundamental physical processes, such as planet formation, star-disk interactions, mass accretion, and ejection. The complex morphological structure of these environments has been confirmed by the now quite rich data sets obtained for a few objects by near-infrared long-baseline interferometry. We gathered numerous interferometric measurements for the young star HD163296 with various interferometers (VLTI, IOTA, KeckI and CHARA), allowing for the first time an image independent of any a priori model to be reconstructed. Using the Multi-aperture image Reconstruction Algorithm (MiRA), we reconstruct images of HD 163296 in the H and K bands. We compare these images with reconstructed images obtained from simulated data using a physical model of the environment of HD 163296. We obtain model-independent $H$ and $K$-band images of the surroundings of HD 163296. The images present several significant features that we can relate to an inclined asymmetric flared disk around HD 163296 with the strongest intensity at about 4-5 mas. Because of the incomplete spatial frequency coverage, we cannot state whether each of them individually is peculiar in any way. For the first time, milli-arcsecond images of the environment of a young star are produced. These images confirm that the morphology of the close environment of young stars is more complex than the simple models used in the literature so far.Comment: 11 pages, 10 figures, accepted A&A pape

Mott-Superfluid transition in bosonic ladders

We show that in a commensurate bosonic ladder, a quantum phase transition occurs between a Mott insulator and a superfluid when interchain hopping increases. We analyse the properties of such a transition as well as the physical properties of the two phases. We discuss the physical consequences for experimental systems such as Josephson Junction arrays.Comment: 4 pages, 2 figures, revtex

On the Relationship Between Molecular Hydrogen and Carbon Monoxide Abundances in Molecular Clouds

The most usual tracer of molecular gas is line emission from CO. However, the reliability of that tracer has long been questioned in environments different from the Milky Way. We study the relationship between H2 and CO abundances using a fully dynamical model of magnetized turbulence coupled to a chemical network simplified to follow only the dominant pathways for H2 and CO formation and destruction, and including photodissociation using a six-ray approximation. We find that the abundance of H2 is primarily determined by the amount of time available for its formation, which is proportional to the product of the density and the metallicity, but insensitive to photodissociation. Photodissociation only becomes important at extinctions under a few tenths of a visual magnitude, in agreement with both observational and prior theoretical work. On the other hand, CO forms quickly, within a dynamical time, but its abundance depends primarily on photodissociation, with only a weak secondary dependence on H2 abundance. As a result, there is a sharp cutoff in CO abundance at mean visual extinctions A_V < 3. At lower values of A_V we find that the ratio of H2 column density to CO emissivity X_CO is proportional to A_V^(-3.5). This explains the discrepancy observed in low metallicity systems between cloud masses derived from CO observations and other techniques such as infrared emission. Our work predicts that CO-bright clouds in low metallicity systems should be systematically larger or denser than Milky Way clouds, or both. Our results further explain the narrow range of observed molecular cloud column densities as a threshold effect, without requiring the assumption of virial equilibrium.Comment: 16 pages, 11 figures. Updated to match version accepted by MNRA

Multi-timescale Solar Cycles and the Possible Implications

Based on analysis of the annual averaged relative sunspot number (ASN) during 1700 -- 2009, 3 kinds of solar cycles are confirmed: the well-known 11-yr cycle (Schwabe cycle), 103-yr secular cycle (numbered as G1, G2, G3, and G4, respectively since 1700); and 51.5-yr Cycle. From similarities, an extrapolation of forthcoming solar cycles is made, and found that the solar cycle 24 will be a relative long and weak Schwabe cycle, which may reach to its apex around 2012-2014 in the vale between G3 and G4. Additionally, most Schwabe cycles are asymmetric with rapidly rising-phases and slowly decay-phases. The comparisons between ASN and the annual flare numbers with different GOES classes (C-class, M-class, X-class, and super-flare, here super-flare is defined as $\geq$ X10.0) and the annal averaged radio flux at frequency of 2.84 GHz indicate that solar flares have a tendency: the more powerful of the flare, the later it takes place after the onset of the Schwabe cycle, and most powerful flares take place in the decay phase of Schwabe cycle. Some discussions on the origin of solar cycles are presented.Comment: 8 pages, 4 figure