PO and PN in the wind of the oxygen-rich AGB star IK Tau

Phosphorus-bearing compounds have only been studied in the circumstellar environments (CSEs) of the asymptotic giant branch (AGB) star IRC +10216 and the protoplanetary nebula CRL 2688, both C-rich objects, and the O-rich red supergiant VY CMa. The current chemical models cannot reproduce the high abundances of PO and PN derived from observations of VY CMa. No observations have been reported of phosphorus in the CSEs of O-rich AGB stars. We aim to set observational constraints on the phosphorous chemistry in the CSEs of O-rich AGB stars, by focussing on the Mira-type variable star IK Tau. Using the IRAM 30m telescope and the Submillimeter Array (SMA), we observed four rotational transitions of PN (J=2-1,3-2,6-5,7-6) and four of PO (J=5/2-3/2,7/2-5/2,13/2-11/2,15/2-13/2). The IRAM 30m observations were dedicated line observations, while the SMA data come from an unbiased spectral survey in the frequency range 279-355 GHz. We present the first detections of PN and PO in an O-rich AGB star and estimate abundances X(PN/H2) of about 3x10^-7 and X(PO/H2) in the range 0.5-6.0x10^-7. This is several orders of magnitude higher than what is found for the C-rich AGB star IRC +10216. The diameter (<=0.7") of the PN and PO emission distributions measured in the interferometric data corresponds to a maximum radial extent of about 40 stellar radii. The abundances and the spatial occurrence of the molecules are in very good agreement with the results reported for VY CMa. We did not detect PS or PH3 in the survey. We suggest that PN and PO are the main carriers of phosphorus in the gas phase, with abundances possibly up to several 10^-7. The current chemical models cannot account for this, underlining the strong need for updated chemical models that include phosphorous compounds.Comment: Accepted for publication in Astronomy & Astrophysics, 10 pages, 8 figure

The Spitzer c2d Survey of Nearby Dense Cores: VI. The Protostars of Lynds Dark Nebula 1221

Observations of Lynds Dark Nebula 1221 from the Spitzer Space Telescope are presented. These data show three candidate protostars towards L1221, only two of which were previously known. The infrared observations also show signatures of outflowing material, an interpretation which is also supported by radio observations with the Very Large Array. In addition, molecular line maps from the Five College Radio Astronomy Observatory are shown. One-dimensional dust continuum modelling of two of these protostars, IRS1 and IRS3, is described. These models show two distinctly different protostars forming in very similar environments. IRS1 shows a higher luminosity and larger inner radius of the envelope than IRS3. The disparity could be caused by a difference in age or mass, orientation of outflow cavities, or the impact of a binary in the IRS1 core.Comment: accepted for publication in Ap

Location of the Multicritical Point for the Ising Spin Glass on the Triangular and Hexagonal Lattices

A conjecture is given for the exact location of the multicritical point in the phase diagram of the +/- J Ising model on the triangular lattice. The result p_c=0.8358058 agrees well with a recent numerical estimate. From this value, it is possible to derive a comparable conjecture for the exact location of the multicritical point for the hexagonal lattice, p_c=0.9327041, again in excellent agreement with a numerical study. The method is a variant of duality transformation to relate the triangular lattice directly with its dual triangular lattice without recourse to the hexagonal lattice, in conjunction with the replica method.Comment: 9 pages, 1 figure; Minor corrections in notatio

Dynamics of the Wang-Landau algorithm and complexity of rare events for the three-dimensional bimodal Ising spin glass

We investigate the performance of flat-histogram methods based on a multicanonical ensemble and the Wang-Landau algorithm for the three-dimensional +/- J spin glass by measuring round-trip times in the energy range between the zero-temperature ground state and the state of highest energy. Strong sample-to-sample variations are found for fixed system size and the distribution of round-trip times follows a fat-tailed Frechet extremal value distribution. Rare events in the fat tails of these distributions corresponding to extremely slowly equilibrating spin glass realizations dominate the calculations of statistical averages. While the typical round-trip time scales exponential as expected for this NP-hard problem, we find that the average round-trip time is no longer well-defined for systems with N >= 8^3 spins. We relate the round-trip times for multicanonical sampling to intrinsic properties of the energy landscape and compare with the numerical effort needed by the genetic Cluster-Exact Approximation to calculate the exact ground state energies. For systems with N >= 8^3 spins the simulation of these rare events becomes increasingly hard. For N >= 14^3 there are samples where the Wang-Landau algorithm fails to find the true ground state within reasonable simulation times. We expect similar behavior for other algorithms based on multicanonical sampling.Comment: 9 pages, 12 figure

Properties of solar polar coronal plumes constrained by Ultraviolet Coronagraph Spectrometer data

We investigate the plasma dynamics (outflow speed and turbulence) inside polar plumes. We compare line profiles (mainly of \ion{O}{6}) observed by the UVCS instrument on SOHO at the minimum of solar cycle 22-23 with model calculations. We consider Maxwellian velocity distributions with different widths in plume and inter-plume regions. Electron densities are assumed to be enhanced in plumes and to approach inter-plume values with increasing height. Different combinations of the outflow and turbulence velocity in the plume regions are considered. We compute line profiles and total intensities of the \ion{H}{1} Ly$\alpha$ and the \ion{O}{6} doublets. The observed profile shapes and intensities are reproduced best by a small solar wind speed at low altitudes in plumes that increases with height to reach ambient inter-plume values above roughly 3-4 R_\sun combined with a similar variation of the width of the velocity distribution of the scattering atoms/ions. We also find that plumes very close to the pole give narrow profiles at heights above 2.5 R_\sun, which are not observed. This suggests a tendency for plumes to be located away from the pole. We find that the inclusion of plumes in the model computations provides an improved correspondence with the observations and confirms previous results showing that published UVCS observations in polar coronal holes can be roughly reproduced without the need for large temperature anisotropy. The latitude distributions of plumes and magnetic flux distributions are studied by analyzing data from different instruments on SOHO and with SOLIS.Comment: 11 figure

SCUBA Mapping of Spitzer c2d Small Clouds and Cores

We present submillimeter observations of dark clouds that are part of the Spitzer Legacy Program, From Molecular Cores to Planet-Forming Disks (c2d). We used the Submillimetre Common User's Bolometer Array to map the regions observed by Spitzer by the c2d program to create a census of dense molecular cores including data from the infrared to the submillimeter. In this paper, we present the basic data from these observations: maps, fluxes, and source attributes. We also show data for an object just outside the Perseus cloud that was serendipitously observed in our program. We propose that this object is a newly discovered, evolved protostar.Comment: 37 pages, accepted to The Astronomical Journa

Derivatives and inequalities for order parameters in the Ising spin glass

Identities and inequalities are proved for the order parameters, correlation functions and their derivatives of the Ising spin glass. The results serve as additional evidence that the ferromagnetic phase is composed of two regions, one with strong ferromagnetic ordering and the other with the effects of disorder dominant. The Nishimori line marks a crossover between these two regions.Comment: 10 pages; 3 figures; new inequalities added, title slightly change

Imaging of Ultraluminous Infrared Galaxies in the Near-UV

We present the first ground-based U' (3410 angstroms) images of Ultraluminous Infrared Galaxies (ULIGs). Strong U' emission (median total M_U' = -20.8) is seen in all systems and in some cases the extended tidal features (both the smooth stellar distribution and compact star-forming features) contribute up to 60-80% of the total flux. The star-forming regions in both samples are found to have ages based on spectral synthesis models in the range 10-100 Myrs, and most differences in color between them can be attributed to the effects of dust reddening. Additionally, it is found that star-formation in compact knots in the tidal tails is most prominent in those ULIGs which have double nuclei, suggesting that the star-formation rate in the tails peaks prior to the actual coalescence of the galaxy nuclei and diminishes quickly thereafter. Similar to results at other wavelengths, the observed star formation at U' can only account for a small fraction of the known bolometric luminosity of the ULIGs. Azimuthally averaged radial light profiles at U' are characterized by a sersic law with index n=2, which is intermediate between an exponential disk and an r^(-1/4) law and closely resembles the latter at large radii. The implications of this near-ultraviolet imaging for optical/near-infrared observations of high redshift counterparts of ULIGs are discussed.Comment: 30 pages, 4 tables, and 9 figures, 2 of which are JPEGs. To appear in the August, 2000 edition of the Astronomical Journa

Existence of a phase transition under finite magnetic field in the long-range RKKY Ising spin glass Dyx_{x}Y1−x_{1-x}Ru2_{2}Si2_{2}

A phase transition of a model compound of the long-range Ising spin glass (SG) Dy$_{x}$Y$_{1-x}$Ru$_{2}$Si$_{2}$, where spins interact via the RKKY interaction, has been investigated. The static and the dynamic scaling analyses reveal that the SG phase transition in the model magnet belongs to the mean-field universality class. Moreover, the characteristic relaxation time in finite magnetic fields exhibits a critical divergent behavior as well as in zero field, indicating a stability of the SG phase in finite fields. The presence of the SG phase transition in field in the model magnet strongly syggests that the replica symmetry is broken in the long-range Ising SG.Comment: 4 pages, 4 figures, to be published in JPSJ (2010