The research of exoplanets has entered an era in which we characterize extrasolar planets. This has become possible with measurements of radii and luminosities. Meanwhile, radial velocity surveys discover also very low-mass planets. Uniting all this observational data into one coherent picture to better understand planet formation is an important, but difficult undertaking. Our approach is to develop a model which can make testable predictions for all these observational methods. We continue to describe how we have extended our formation model into a self-consistently coupled formation and evolution model. We show how we calculate the internal structure of the solid core and radiogenic heating. We also improve the protoplanetary disk model. Finally, we conduct population synthesis calculations. We present how the planetary mass-radius relationship of planets with primordial H/He envelopes forms and evolves in time. The basic shape of the M-R relation can be understood from the core accretion model. Low-mass planets cannot bind massive envelopes, while super-critical cores necessarily trigger runway gas accretion, leading to "forbidden" zones in the M-R plane. For a given mass, there is a considerable diversity of radii. We compare the synthetic M-R relation with the observed one, finding good agreement for a>0.1 AU. The synthetic radius distribution is characterized by a strong increase towards small R, and a second, lower local maximum at ~1 Jovian radius. The increase towards small radii reflects the increase of the mass function towards low M. The second local maximum is due to the fact that radii are nearly independent of mass for giant planets. A comparison of the synthetic radius distribution with Kepler data shows agreement for R>2 Earth radii, but divergence for smaller radii. We predict that in the next few years, Kepler should find the second, local maximum at ~1 Jovian radius.
Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for structure formation is one of the most promising ways to test inflation and differentiate among different scenarios. The scale-dependent imprint of PNG on the large-scale clustering of galaxies and quasars has already been used to place significant constraints on the level of PNG in our observed Universe. Such measurements depend upon an accurate and robust theory for how PNG affects the bias of galactic halos relative to the underlying matter density field. We improve upon previous work by employing a more general analytical method - the path-integral extension of the excursion set formalism - which is able to account for the non-Markovianity caused by PNG in the random-walk model used to identify halos in the initial density field. This non-Markovianity encodes information about environmental effects on halo formation which have so far not been taken into account in analytical bias calculations. We compute both scale-dependent and -independent corrections to the halo bias, along the way presenting an expression for the conditional collapsed fraction for the first time, and a new expression for the conditional halo mass function. To leading order in our perturbative calculation, we recover the halo bias results of Desjacques et. al. (2011), including the new scale-dependent correction reported there. However, we show that the non-Markovian dynamics from PNG can lead to marked differences in halo bias when next-to-leading order terms are included. We quantify these differences here.
We study the impact of primordial non-Gaussianity generated during inflation on the bias of halos using excursion set theory. We recapture the familiar result that the bias scales as $k^{-2}$ on large scales but explicitly identify the approximations that go into this conclusion and the corrections to it. We solve the more complicated problem of non-spherical halos, for which the collapse threshold is scale dependent.
We present the discovery of a 0.2" companion to the young M dwarf GJ 3629 as part of our high contrast adaptive optics imaging search for giant planets around low-mass stars with the Keck-II and Subaru telescopes. Two epochs of imaging confirm the pair is co-moving and reveal signs of orbital motion. The primary exhibits saturated X-ray emission, which together with its UV photometry from GALEX point to an age younger than ~300 Myr. At these ages the companion lies below the hydrogen burning limit with a model-dependent mass of 46 +/- 16 Mjup based on the system's photometric distance of 22 +/- 3 pc. Resolved YJHK photometry of the pair indicates a spectral type of M7 +/- 2 for GJ 3629 B. With a projected separation of 4.4 +/- 0.6 AU and an estimated orbital period of 21 +/- 5 yr, GJ 3629 AB is likely to yield a dynamical mass in the next several years, making it one of only a handful of brown dwarfs to have a measured mass and an age constrained from the stellar primary.
We use stellar population synthesis modeling to analyze the host galaxy properties of a sample of 33 UV-selected, narrow-lined active galactic nuclei (AGNs) at z ~ 2 - 3. In order to quantify the contribution of AGN emission to host galaxy broadband spectral energy distributions (SEDs), we use the subsample of 11 AGNs with photometric coverage spanning from rest-frame UV through near-IR wavelengths. Modeling the SEDs of these objects with a linear combination of stellar population and AGN templates, we infer the effect of the AGN on derived stellar population parameters. We also estimate the typical bias in derived stellar populations for AGNs lacking rest-frame near-IR wavelength coverage, and develop a method for inferring the true host galaxy properties. We compare AGN host galaxy properties to those of a sample of UV-selected, star-forming non-AGNs in the same redshift range, including a subsample carefully matched in stellar mass. Although the AGNs have higher masses and SFRs than the full non-active sample, their stellar population properties are consistent with those of the mass-selected sample, suggesting that the presence of an AGN is not connected with the cessation of star-formation activity in star-forming galaxies at z ~ 2 - 3. We suggest that a correlation between M_BH and galaxy stellar mass is already in place at this epoch. Assuming a roughly constant Eddington ratio for AGNs at all stellar masses, we are unable to detect the AGNs in low-mass galaxies because they are simply too faint.
Rotation is expected to have an important influence on the structure and the evolution of stars. However, the mechanisms of angular momentum transport in stars remain theoretically uncertain and very complex to take into account in stellar models. To achieve a better understanding of these processes, we desperately need observational constraints on the internal rotation of stars, which until very recently were restricted to the Sun. In this paper, we report the detection of mixed modes - i.e. modes that behave both as g modes in the core and as p modes in the envelope - in the spectrum of the early red giant KIC7341231, which was observed during one year with the Kepler spacecraft. By performing an analysis of the oscillation spectrum of the star, we show that its non-radial modes are clearly split by stellar rotation and we are able to determine precisely the rotational splittings of 18 modes. We then find a stellar model that reproduces very well the observed atmospheric and seismic properties of the star. We use this model to perform inversions of the internal rotation profile of the star, which enables us to show that the core of the star is rotating at least five times faster than the envelope. This will shed new light on the processes of transport of angular momentum in stars. In particular, this result can be used to place constraints on the angular momentum coupling between the core and the envelope of early red giants, which could help us discriminate between the theories that have been proposed over the last decades.
We present a study of the Milky Way disk and halo magnetic field, determined from observations of Faraday rotation measure (RM) towards 641 polarized extragalactic radio sources in the Galactic longitude range 100-117 degs, within 30 degs of the Galactic plane. For |b| < 15 degs, we observe a symmetric RM distribution about the Galactic plane. This is consistent with a disk field in the Perseus arm of even parity across the Galactic mid-plane. In the range 15<|b|<30 degs, we find median rotation measures of -15+/-4 rad/m^2 and -62+/-5 rad/m^2 in the northern and southern Galactic hemispheres, respectively. If the RM distribution is a signature of the large-scale field parallel to the Galactic plane, this suggests that the halo magnetic field toward the outer Galaxy does not reverse direction across the mid-plane. The variation of RM as a function of Galactic latitude in this longitude range is such that RMs become more negative at larger |b|. This is consistent with an azimuthal magnetic field of strength 2 microGauss (7 microGauss) at a height 0.8-2 kpc above (below) the Galactic plane between the local and the Perseus spiral arm. We propose that the Milky Way could possess spiral-like halo magnetic fields similar to those observed in M51.
We present comprehensive seeing statistics for the San Pedro M\'artir site derived from the Thirty Meter Telescope site selection data. The observations were obtained between 2004 and 2008 with a Differential Image Motion Monitor (DIMM) and a Multi Aperture Scintillation Sensor (MASS) combined instrument (MASS--DIMM). The parameters that are statistically analised here are: whole atmosphere seeing -measured by the DIMM-; free atmosphere seeing --measured by the MASS--; and ground-layer seeing (GL) --difference between the total and free-atmosphere seeing--. We made a careful data coverage study along with statistical distributions of simultaneous MASS--DIMM seeing measurements, in order to investigate the nightly, monthly, seasonal, annual and global behaviour, as well as possible hourly seeing trends. Although this campaign covers five years, the sampling is uneven, being 2006 and 2007 the best sampled years in terms of seasonal coverage. The overall results yield a median seeing of 0.78 (DIMM), 0.37 (MASS) and 0.59 arcsec (GL). The strongest contribution to the whole atmosphere seeing comes, therefore, from a strong ground layer. We find that the best season is summer, while the worst one is winter, in accordance with previous studies. It is worth noting that the best yearly results are correlated with the best sampled years. The hourly analysis shows that there is no statistically significant tendency of seeing degradation towards dawn. The seeing values are slightly larger than those reported before. This may be caused by climate changes.
While the standard and most popular explanation for the flatness of galaxy rotation curves is dark matter, one cannot at this stage rule out an explanation based on a modified law of gravitation, which agrees with Newtonian gravitation on the scale of the solar system, but differs from it on larger length scales. Examples include Modfied Newtonian Dynamics [MOND] and Scalar-Tensor-Vector Gravity [STVG]. Here we report on a fourth order modification of the Poisson equation which yields the same Yukawa type modification of Newtonian gravity as STVG, and which can explain flat galaxy rotation curves for a large sample of galaxies, once specific values for two parameters have been chosen. We speculate on two possible origins for this modified Poisson equation: first, a possible fourth order modification of general relativity, and second, quadrupole gravitational polarization induced on a galaxy because of the pull of neighbouring galaxies.
Since its launch, the Fermi satellite has firmly identified 5 pulsar wind nebulae plus a large number of candidates, all powered by young and energetic pulsars. HESS J1857+026 is a spatially extended gamma-ray source detected by H.E.S.S. and classified as a possible pulsar wind nebula candidate powered by PSR J1856+0245. We search for gamma-ray pulsations from PSR J1856+0245 and explore the characteristics of its associated pulsar wind nebula. Using a rotational ephemeris obtained from the Lovell telescope at Jodrell Bank Observatory at 1.5 GHz, we phase-fold 36 months of gamma-ray data acquired by the Large Area Telescope (LAT) aboard Fermi. We also perform a complete gamma-ray spectral and morphological analysis. No gamma-ray pulsations were detected from PSR J1856+0245. However, significant emission is detected at a position coincident with the TeV source HESS J1857+026. The gamma-ray spectrum is well described by a simple power-law with a spectral index of 1.53 \pm 0.11_{\rm stat} \pm 0.55_{\rm syst}$ and an energy flux of $G(0.1$--100 GeV$)=(2.71 \pm 0.52_{\rm stat} \pm 1.51_{\rm syst}) \times 10^{-11}$ ergs cm$^{-2}$ s$^{-1}$. The $\gamma$-ray luminosity is $L_{PWN}^{\gamma} (0.1$--100 GeV$)=(2.5 \pm 0.5_{stat} \pm 1.5_{syst}) \times 10^{35} (\frac{d}{9 kpc})^2$ ergs s$^{-1}$, assuming a distance of 9 kpc. This implies a $\gamma-$ray efficiency of $\sim$ 5% for $\dot{E}=4.6 \times 10^{36}$ erg $s^{-1}$, in the range expected for pulsar wind nebulae. Detailed multi-wavelength modeling provides new constraints on its pulsar wind nebula nature.
The young, low-mass, triple system NTTS 155808-2219 (ScoPMS 20) was previously identified as a ~17-day period single-lined spectroscopic binary with a tertiary component at 0.21 arcseconds. Using high-resolution infrared spectra, acquired with NIRSPEC on Keck II, both with and without adaptive optics, we measured radial velocities of all three components. Reanalysis of the single-lined visible light observations, made from 1987 to 1993, also yielded radial velocity detections of the three stars. Combining visible light and infrared data to compute the orbital solution produces orbital parameters consistent with the single-lined solution and a mass ratio of q = 0.78 +/- 0.01 for the SB. We discuss the consistency between our results and previously published data on this system, our radial-velocity analysis with both observed and synthetic templates, and the possibility that this system is eclipsing, providing a potential method for the determination of the stars' absolute masses. Over the ~20 year baseline of our observations, we have measured the acceleration of the SB's center-of-mass in its orbit with the tertiary. Long-term, adaptive optics imaging of the tertiary will eventually yield dynamical data useful for component mass estimates.
We report on the analysis of two XMM-Newton observations of the recently discovered soft gamma repeater Swift J1834.9-0846, taken in September 2005 and one month after the source went into outburst on 2011 August 7. We performed timing and spectral analyses on the point source as well as on the extended emission. We find that the source period is consistent with an extrapolation of the Chandra ephemeris reported earlier and the spectral properties remained constant. The source luminosity decreased to a level of 1.6x10^34 erg s^-1 following a decay trend of $\propto t^{-0.5}$. Our spatial analysis of the source environment revealed the presence of two extended emission regions around the source. The first (Region A) is a symmetric ring around the point source, starting at 25arcsec and extending to ~50arcsec. We argue that Region A is a dust scattering halo. The second (Region B) has an asymmetrical shape extending between 50arcsec and 150arcsec, and is detected both in the pre- and post-outburst data. We argue that this region is a possible magnetar wind nebula (MWN). The X-ray efficiency of the MWN with respect to the rotation energy loss is substantially higher than those of rotation powered pulsars: $\eta_{\rm X}\equiv L_{\rm MWN,0.5-8 keV}/\dot{E}_{\rm rot}\approx0.7$. The higher efficiency points to a different energy source for the MWN of Swift J1834.9-0846, most likely bursting activity of the magnetar, powered by its high magnetic field, B=1.4x10^14 G.
The Core Mass Functions (CMFs) of low-mass star-forming regions are found to
resemble the shape of the Initial Mass Function (IMF). A similar result is
observed for the dust clumps in high-mass star forming regions, although at
spatial scales of clusters that do not resolve the substructure found in them.
The region IRAS 19410+2336 is one exception, having been observed at spatial
scales on the order of ~2500AU, resolving the clump substructure into
individual cores.
We mapped that region with the PdBI in the 1.4mm and 3mm continuum and
several transitions of H2CO and CH3CN. The H2CO transitions were also observed
with the IRAM 30m Telescope. We detected 26 continuum sources at 1.4mm with a
spatial resolution down to ~2200 AU, distributed in two protoclusters. With the
lines emission we derived the temperature structure of the region, ranging from
35 to 90K. With them we calculated the core masses of the detected sources,
ranging from ~0.7 to ~8 M_sun. These masses were strongly (~90%) affected by
the interferometer spatial filtering. Considering only the detected dense cores
we derived a CMF with a power-law index b=-2.3+-0.2. We resolve the Jeans
length of the protoclusters by one order of magnitude, and only find little
velocity dispersion between the different subsources.
Since we cannot unambiguously differentiate protostellar and prestellar
cores, the derived CMF is not prestellar. Also, because of the large missing
flux, we cannot establish a firm link between the CMF and the IMF. This implies
that future high-mass CMF studies will need to complement the interferometer
continuum data with the short spacing data, a task suitable for ALMA. We note
that the method of extracting temperatures using H2CO lines becomes less
applicable when reaching the dense core scales of the interferometric
observations because most of the H2CO appears to originate in the envelope
structure.
(Abridged) We study the polarisation properties, magnetic field strength, and synchrotron emission scale-height of Milky-Way-like galaxies in comparison with other spiral galaxies. We use our 3D-emission model of the Milky Way Galaxy for viewing the Milky Way from outside at various inclinations as spiral galaxies are observed. When seen edge-on the synchrotron emission from the Milky Way has an exponential scale-height of about 0.74 kpc, which is much smaller than the values obtained from previous models. We find that current analysis methods overestimate the scale-height of synchrotron emission of galaxies by about 10% at an inclination of 80 degree and about 40% at an inclination of 70 degree because of contamination from the disk. The observed RMs for face-on galaxies derived from high-frequency polarisation measurements approximate to the Faraday depths (FDs) when scaled by a factor of two. For edge-on galaxies, the observed RMs are indicative of the orientation of the large-scale magnetic field, but are not well related with the FDs. Assuming energy equipartition between the magnetic field and particles for the Milky Way results in an average magnetic-field strength, which is about two times larger than the intrinsic value for a K factor of 100. The number distribution of the integrated polarisation percentages of a large sample of unresolved Milky-Way-like galaxies peaks at about 4.2% at 4.8 GHz and at about 0.8% at 1.4GHz. Integrated polarisation angles rotated by 90 degree align very well with the position angles of the major axes, implying that unresolved galaxies do not have intrinsic RMs.
We present an analysis of over 150 ks of data on the planetary nebula WeBo 1 (PN G135.6+01.0) obtained with the Swift Ultraviolet Optical Telescope (UVOT). The central object of this nebula has previously been described as a late-type K giant barium star with a possible hot companion, most likely a young pre-white dwarf. UVOT photometry shows that while the optical photometry is consistent with a large cool object, the near-ultraviolet (UV) photometry shows far more UV flux than could be produced by any late-type object. Using model stellar atmospheres and a comparison to UVOT photometry for the pre-white dwarf PG 1159-035, we find that the companion has a temperature of at least 40,000 K and a radius of, at most, 0.056 R_sun. While the temperature and radius are consistent with a hot compact stellar remnant, they are lower and larger, respectively, than expected for a typical young pre-white dwarf. This likely indicates a deficiency in the assumed UV extinction curve. We find that higher temperatures more consistent with expectations for a pre-white dwarf can be derived if the foreground dust has a strong "blue bump" at 2175 AA and a lower R_V. Our results demonstrate the ability of Swift to both uncover and characterize hot hidden companion stars and to constrain the UV extinction properties of foreground dust based solely on UVOT photometry.
The atmospheres of hot Jupiters and other strongly-forced exoplanets are susceptible to a thermal instability in the presence of ohmic dissipation, weak magnetic drag and strong winds. The instability occurs in radiatively-dominated atmospheric regions when the ohmic dissipation rate increases with temperature faster than the radiative (cooling) rate. The instability domain covers a specific range of atmospheric pressures and temperatures, typically P ~ 3-300 mbar and T ~ 1500-2500K for hot Jupiters, which makes it a candidate mechanism to explain the dayside thermal "inversions" inferred for a number of such exoplanets. The instability is suppressed by high levels of non-thermal photoionization, in possible agreement with a recently established observational trend. We highlight several shortcomings of the instability treatment presented here. Understanding the emergence and outcome of the instability, which should result in locally hotter atmospheres with stronger levels of drag, will require global non-linear atmospheric models with adequate MHD prescriptions.
A complete understanding of Doppler shift in active region loops can help probe the basic physical mechanism involved into the heating of those loops. Here we present observations of upflows in coronal loops detected in a range of temperature temperatures (log T=5.8 - 6.2). The loop was not discernible above these temperatures. The speed of upflow was strongest at the footpoint and decreased with height. The upflow speed at the footpoint was about 20 km/s in Fe VIII which decreased with temperature being about 13 km/s in Fe X, about 8 km/s in Fe XII and about 4 km/s in FeXIII. To the best of our knowledge this is the first observation providing evidence of upflow of plasma in coronal loop structures at these temperatures. We interpret these observations as evidence of chromospheric evaporation in quasi-static coronal loops.
We describe the target selection algorithm for the low latitude disk portion of the LAMOST Pilot Survey, which aims to test systems in preparation for the LAMOST spectroscopic survey. We use the PPMXL (Roeser et al. 2010) astrometric catalog, which provides positions, proper motions, B/R/I magnitudes (mostly) from USNO-B (Monet et al. 2003) and J/H/Ks from The Two Micron All Sky Survey (2MASS, see Skrutskie et al. 2006) as well. We chose 8 plates along the Galactic plane, in the region $0^\circ<\alpha<67^\circ$ and $42^\circ<\delta<59^\circ$, that cover 22 known open clusters with a range of ages. Adjacent plates may have small overlapping. Each plate covers an area $2.5^\circ$ in radius,with central star (for Shark-Hartmann guider) brighter than $\sim8^{\rm th}$ magnitude. For each plate, we create an input catalog in the magnitude range $11.3<Imag<16.3$ and $Bmag$ available from PPMXL. The stars are selected to satisfy the requirements of the fiber positioning system and have a uniform distribution in the $I$ vs. $B-I$ color-magnitude diagram. Our final input catalog consists of 12,000 objects on each of 8 plates that are observable during the winter observing season in Xinglong Station of the National Astronomical Observatory of China.
Double radio relics in galaxy clusters are rare phenomena that trace shocks in the outskirts of merging galaxy clusters. We have carried out a spectral and polarization study of the spectacular double relics in the galaxy cluster A3376 using the Giant Metrewave Radio Telescope at 150 and 325 MHz and the Very Large Array at 1400 MHz. The polarization study at 1400 MHz reveals a high degree of polarization (~30%) and aligned magnetic field vectors (not corrected for Faraday rotation) in the eastern relic. A highly polarized (>60%) filamentary radio source of size ~300 kpc near the eastern relic and north of the bent-jet radio galaxy is detected for the first time. The western relic is less polarized and does not show aligned magnetic field vectors. The distribution of spectral indices between 325 and 1400 MHz over the radio relics show steepening from the outer to the inner edges of the relics. The spectral indices of the eastern and the western relics imply Mach numbers in the range 2.2 to 3.3. Remarkable features such as the inward filament extending from the eastern relic, the highly polarized filament, the complex polarization properties of the western relic and the separation of the BCG from the ICM by a distance >900 kpc are noticed in the cluster. A comparison with simulated cluster mergers is required to understand the complex properties of the double relics in the context of the merger in A3376. An upper limit (log(P(1.4GHz) W/Hz < 23.0) on the strength of a Mpc size radio halo in A3376 is estimated.
We present new [SII] images of the HH 30 jet and counterjet observed in 2006, 2007, and 2010 that allowed us to measure with improved accuracy the positions and proper motions of the jet and counterjet knots. Our results show that the motion of the knots is essentially ballistic, with the exception of the farthest knots, which trace the large scale C-shape bending of the jet. The observed bending of the jet can be produced by a relative motion of the HH 30 star with respect to its surrounding environment, caused either by a possible proper motion of the HH 30 star, or by the entrainment of environment gas by the red lobe of the nearby L1551-IRS 5 outflow. Alternatively, the bending can be produced by the stellar wind from a nearby CTTS, identified in the 2MASS catalog. The proper motion velocities of the knots of the counterjet show more variations than those of the jet. In particular, we identify two knots of the counterjet that have the same kinematic age but whose velocities differ by almost a factor of two. Thus, it appears that counterjet knots launched simultaneously can be ejected with very different velocities. We confirm that the observed wiggling of the jet and counterjet arises from the orbital motion of the jet source in a binary system. Precession is of secondary importance in shaping the jet. We derive an orbital period $\tau_o=114\pm2$ yr and a mass function $m\mu_c^3=0.014\pm0.006$ $M_\odot$. For a mass of the system of $m=0.45\pm0.04$ $M_\odot$ (the value inferred from the disk kinematics) we obtain a mass $m_j=0.31\pm0.04$ $M_\odot$ for the jet source, a mass $m_c=0.14\pm0.03$ $M_\odot$ for the companion, and a binary separation of $a=18.0\pm0.6$ AU. This binary separation coincides with the value required to account for the size of the inner hole observed in the disk, attributed to tidal truncation in a binary system.
We present the results of a visual search for galaxy-scale gravitational lenses in ~7 square degrees of Hubble Space Telescope (HST) images. The dataset comprises the whole imaging data ever taken with the Advanced Camera for Surveys (ACS) in the filter F814W (I-band) up to August 31st, 2011, i.e. 6.03 square degrees excluding the field of the Cosmic Evolution Survey (COSMOS). In addition, we have searched for lenses in the whole Wide Field Camera 3 (WFC3) near-IR imaging dataset in all filters 1.01 square degree up to the same date. Our primary goal is to provide a sample of lenses with a broad range of different morphologies and lens-source brightness contrast in order to design and train future automated lens finders in view of all-sky surveys. Our criteria to select lenses are purely morphological as we do not use any color or redshift information. The final candidate selection is very conservative hence leading to a nearly pure but incomplete sample. We find 49 new lens candidates: 40 in the ACS images and 9 in the WFC3 images. Out of these, 16 candidates are secure lenses owe to their highly recognizable morphology, 21 more are very good candidates, and 12 more have morphologies compatible with gravitational lensing. The imaging dataset is heterogeneous in depth and spans a broad range of galactic latitudes. It is therefore insensitive to cosmic variance and allows to estimate the number of galaxy-scale strong lenses on the sky for a putative survey depth. Because of the incompleteness of the sample, the estimated lensing rates should be taken as lower limits. Using these, we anticipate that a 15000 square degrees space survey such as Euclid will find at least 60000 galaxy-scale strong lenses down to a limiting AB magnitude of I=24.5 (10-sigma) or I=25.8 (3-sigma).
We observed a sample of 20 representative Herbig Ae/Be stars and five A-type debris discs with PACS onboard of Herschel. The observations were done in spectroscopic mode, and cover far-IR lines of [OI], [CII], CO, CH+, H2O and OH. We have a [OI]63 micron detection rate of 100% for the Herbig Ae/Be and 0% for the debris discs. [OI]145 micron is only detected in 25%, CO J=18-17 in 45% (and less for higher J transitions) of the Herbig Ae/Be stars and for [CII] 157 micron, we often found spatially variable background contamination. We show the first detection of water in a Herbig Ae disc, HD 163296, which has a settled disc. Hydroxyl is detected as well in this disc. CH+, first seen in HD 100546, is now detected for the second time in a Herbig Ae star, HD 97048. We report fluxes for each line and use the observations as line diagnostics of the gas properties. Furthermore, we look for correlations between the strength of the emission lines and stellar or disc parameters, such as stellar luminosity, UV and X-ray flux, accretion rate, PAH band strength, and flaring. We find that the stellar UV flux is the dominant excitation mechanism of [OI]63 micron, with the highest line fluxes found in those objects with a large amount of flaring and greatest PAH strength. Neither the amount of accretion nor the X-ray luminosity has an influence on the line strength. We find correlations between the line flux of [OI]63 micron and [OI]145 micron, CO J = 18-17 and [OI]6300 \AA, and between the continuum flux at 63 micron and at 1.3 mm, while we find weak correlations between the line flux of [OI]63 micron and the PAH luminosity, the line flux of CO J = 3-2, the continuum flux at 63 micron, the stellar effective temperature and the Brgamma luminosity. (Abbreviated version)
The R Coronae Borealis (RCB) stars are rare hydrogen-deficient, carbon-rich, supergiants, best known for their spectacular declines in brightness at irregular intervals. Efforts to discover more RCB stars have more than doubled the number known in the last few years and they appear to be members of an old, bulge population. Two evolutionary scenarios have been suggested for producing an RCB star, a double degenerate merger of two white dwarfs, or a final helium shell flash in a planetary nebula central star. The evidence pointing toward one or the other is somewhat contradictory, but the discovery that RCB stars have large amounts of 18O has tilted the scales towards the merger scenario. If the RCB stars are the product of white dwarf mergers, this would be a very exciting result since RCB stars would then be low-mass analogs of type Ia supernovae. The predicted number of RCB stars in the Galaxy is consistent with the predicted number of He/CO WD mergers. But, so far, only about 65 of the predicted 5000 RCB stars in the Galaxy have been discovered. The mystery has yet to be solved.
We exploit the SDSS galaxy groups catalogue of Yang et al. to study how the gas-phase metallicities of star-forming galaxies depend on environment. We find that satellite and central galaxies follow a qualitatively similar stellar mass (M_*) - gas-phase metallicity relation. Satellites, though, have higher gas-phase metallicities than equally massive centrals, and this difference increases with decreasing M_*. We also find that the gas-phase metallicity of satellites increases with halo mass at fixed stellar mass. This increment is more pronounced for less massive galaxies. We also show that low mass satellite galaxies have higher gas-phase metallicities than central galaxies of the same stellar metallicity. This difference becomes negligible for more massive galaxies of roughly solar metallicity. We demonstrate that the observed differences in gas-phase metallicity between centrals and satellites at fixed M_* are not a consequence of stellar mass stripping (advocated by Pasquali et al. in order to explain similar differences but in stellar metallicity), nor to the past star formation history of these galaxies as quantified by their surface mass density or gas mass fraction. Rather, we argue that these trends probably originate from a combination of three environmental effects: (i) strangulation, which prevents satellite galaxies from accreting new, low metallicity gas which would otherwise dilute their ISM, (ii) ram-pressure stripping of the outer gas disk, thereby inhibiting radial inflows of low-metallicity gas, and (iii) external pressure provided by the hot gas of the host halo which prevents metal-enriched outflows from escaping the galaxies. [abridged]
When searching for populations of rare and/or weak signals in noisy data, it is common to use a detection threshold to remove marginal events which are unlikely to be the signals of interest; or a detector might have limited sensitivity, causing it to not detect some of the population. In both cases a selection of data has occurred, which can potentially bias any inferences drawn from the remaining data, and this effect must be corrected for. We show how the selection bias is naturally avoided by using the full information from the search, considering both the selected data and our ignorance of the data that are thrown away, and considering all relevant signal and noise models. This approach produces unbiased estimates of parameters even in the presence of false alarms and incomplete data.
Highest posterior density (HPD) intervals are derived for the true eccentricities of spectroscopic binaries with measured values e near 0. This is achieved by adopting an informative prior representing orbital decay due to tidal dissipation. These HPD intervals yield upper limits when e is below the detection threshold e_th and seamlessly transform to upper and lower bounds when e > e_th. Sampling experiments confirm the validity of the HPD intervals.
Search for extraterrestrial life and intelligence constitutes one of the major endeavors in science, but has yet been quantitatively modeled only rarely and in a cursory and superficial fashion. We argue that probabilistic cellular automata (PCA) represent the best quantitative framework for modeling astrobiological history of the Milky Way and its Galactic Habitable Zone. The relevant astrobiological parameters are to be modeled as the elements of the input probability matrix for the PCA kernel. With the underlying simplicity of the cellular automata constructs, this approach enables a quick analysis of large and ambiguous input parameters' space. We perform a simple clustering analysis of typical astrobiological histories and discuss the relevant boundary conditions of practical importance for planning and guiding actual empirical astrobiological and SETI projects. In addition to showing how the present framework is adaptable to more complex situations and updated observational databases from current and near-future space missions, we demonstrate how numerical results could offer a cautious rationale for continuation of practical SETI searches.
The presence of two stellar populations in the Milky Way bulge has been reported recently. We aim at studying the abundances and kinematics of stars in the outer bulge, thereby providing additional constraints on models of its formation. Spectra of 401 red giant stars in a field at (l,b)=(0{\deg},-10{\deg}) were obtained with FLAMES at the VLT. Stars of luminosities down to below the two bulge red clumps (RCs) are included. From these spectra we measure general metallicities, abundances of Fe and the alpha-elements, and radial velocities (RV) of the stars. These measurements as well as photometric data are compared to simulations with the Besancon and TRILEGAL models of the Galaxy. We confirm the presence of two populations among our sample stars: i) a metal-rich one at [M/H] ~+0.3, comprising about 30% of the sample, with low RV dispersion and low alpha-abundance, and ii) a metal-poor population at [M/H] ~-0.6 with high RV dispersion and high alpha-abundance. The metal-rich population could be connected to the Galactic bar. We identify this population as the carrier of the double RC feature. We do not find a significant difference in metallicity or RV between the two RCs, a small difference in metallicity being probably due to a selection effect. The RV dispersion agrees well with predictions of the Besancon Galaxy model, but the metallicity of the "thick bulge" model component should be shifted to lower metallicity by 0.2 to 0.3dex to well reproduce the observations. We present evidence that the metallicity distribution function depends on the evolutionary state of the sample stars, suggesting that enhanced mass loss preferentially removes metal-rich stars. The enhancement in alpha-elements decreases with increasing metallicity, a plateau at low metallicity is however lacking in our data.
We study planetesimal evolution in circumbinary disks, focusing on the three systems Kepler 16, 34 and 35 where planets have been discovered recently. We show that for circumbinary planetesimals, in addition to secular forcing, eccentricities evolve on a dynamical timescale, which leads to orbital crossings even in the presence of gas drag. This makes the current locations of the circumbinary Kepler planets hostile to planetesimal accretion. We then present results from simulations including planetesimal formation and dust accretion, and show that even in the most favourable case of 100% efficient dust accretion, in situ growth starting from planetesimals smaller than ~10 km is difficult for Kepler 16b, Kepler 34b and Kepler 35b. These planets were likely assembled further out in the disk, and migrated inward to their current location.
SXP 1062, a newly discovered Be/X-ray binary in the Small Magellanic Cloud, provides the first example of a robust association with a supernova remnant (SNR). The short age estimated for the SNR qualifies SXP 1062 as the youngest known source in its class, ~ 10^4 yr. Here we discuss possible evolutionary scenarios for SXP 1062 in the attempt to reconcile its long spin period, P=1062 s, and short age. Our results indicate that SXP 1062 may host a neutron star born with a large initial magnetic field, typically in excess of ~ 10^14 G, which then decayed to ~ 10^13 G.
The new results of our observing campaign targeting the isolated neutron star 2XMM J104608.7-594306 in the Carina Nebula are used to understand how peculiar groups of isolated neutron stars relate to each other, as well as to the bulk of the normal radio pulsar population.
We describe a model of dark matter halo abundances and clustering which combines the two most widely used approaches to this problem: that based on peaks and the other based on excursion sets. Our approach can be thought of as addressing the cloud-in-cloud problem for peaks and/or modifying the excursion set approach so that it averages over a special subset, rather than all possible walks. In this respect, it seeks to account for correlations between steps in the walk as well as correlations between walks. We first show how the excursion set and peaks models can be written in the same formalism, and then use this correspondence to write our combined excursion set peaks model. We then give simple expressions for the mass function and bias, showing that even the linear halo bias factor is predicted to have a scale dependence as a consequence of the nonlocality associated with the peak constraint. At large masses, our model has little or no need to rescale the variable delta_c from the value associated with spherical collapse, and suggests a simple explanation for why the linear halo bias factor appears to lie above that based on the peak-background split at high masses when such a rescaling is assumed. Although we have concentrated on peaks, our analysis is more generally applicable to other traditionally single-scale analyses of large-scale structure.
We study the production of Hot Jupiters (HJs) in stellar binaries. We show that the "eccentric Kozai-Lidov" (EKL) mechanism can play a key role in the dynamical evolution of a star-planet-star triple system. We run a large set of Monte Carlo simulations including the secular evolution of the orbits, general relativistic precession, and tides, and we determine the semi-major axis, eccentricity, inclination and spin-orbit angle distributions of the HJs that are produced. We explore the effect of different tidal friction parameters on the results. We find that the efficiency of forming HJs when taking the EKL mechanism into account is higher then previously estimated. Accounting for the frequency of stellar binaries, we find that this production mechanism can account for about 30% of the observed HJ population. Current observations of spin-orbit angles are consistent with this mechanism producing \sim 30% of all HJs, and up to 100% of the misaligned systems. Based on the properties of binaries without a HJ in our simulations, we predict the existence of many Jupiter-like planets with moderately eccentric and inclined orbits and semi-major axes of several AU.
The physical relationship between the far-infrared and radio fluxes of star forming galaxies has yet to be definitively determined. The favored interpretation, the "calorimeter model," requires that supernova generated cosmic ray (CR) electrons cool rapidly via synchrotron radiation. However, this cooling should steepen their radio spectra beyond what is observed, and so enhanced ionization losses at low energies from high gas densities are also required. Further, evaluating the minimum energy magnetic field strength with the traditional scaling of the synchrotron flux may underestimate the true value in massive starbursts if their magnetic energy density is comparable to the hydrostatic pressure of their disks. Gamma-ray spectra of starburst galaxies, combined with radio data, provide a less ambiguous estimate of these physical properties in starburst nuclei. While the radio flux is most sensitive to the magnetic field, the GeV gamma-ray spectrum normalization depends primarily on gas density. To this end, spectra above 100 MeV were constructed for two nearby starburst galaxies, NGC 253 and M82, using Fermi data. Their nuclear radio and far-infrared spectra from the literature are compared to new models of the steady-state CR distributions expected from starburst galaxies. Models with high magnetic fields, favoring galaxy calorimetry, are overall better fits to the observations. These solutions also imply relatively high densities and CR ionization rates, consistent with molecular cloud studies.
Gamma Ray Burst prompt emission is believed to originate from electrons accelerated in a highly relativistic outflow. "Internal shocks" due to collisions between shells ejected by the central engine is a leading candidate for electron acceleration. While synchrotron radiation is generally invoked to interpret prompt gamma-ray emission within the internal shock model, synchrotron self-Compton (SSC) is also considered as a possible candidate of radiation mechanism. In this case, one would expect a synchrotron emission component at low energies, and the naked-eye GRB 080319B has been considered as such an example. In the view that the gamma-ray lightcurve of GRB 080319B is much more variable than its optical counterpart, in this paper we study the relative variability between the synchrotron and SSC components. We develop a "top-down" formalism by using observed quantities to infer physical parameters, and subsequently to study the temporal structure of synchrotron and SSC components of a GRB. We complement the formalism with a "bottom-up" approach where the synchrotron and SSC lightcurves are calculated through a Monte-Carlo simulations of the internal shock model. Both approaches lead to the same conclusion. Small variations in the synchrotron lightcurve can be only moderately amplified in the SSC lightcurve. The SSC model therefore cannot adequately interpret the gamma-ray emission properties of GRB 080319B.
The Earth is continuously showered by charged cosmic ray particles, naturally produced atomic nuclei moving with velocity close to the speed of light. Among these are ultra high energy cosmic ray particles with energy exceeding 5x10^19 eV, which is ten million times more energetic than the most energetic particles produced at the Large Hadron Collider. Astrophysical questions include: what phenomenon accelerates particles to such high energies, and what sort of nuclei are energized? Also, the magnetic deflection of the trajectories of the cosmic rays makes them potential probes of galactic and intergalactic magnetic fields. We develop a Bayesian hierarchical model that can be used to compare different association models between the cosmic rays and source population, using Bayes factors. A measurement model with directional uncertainties and accounting for non-uniform sky exposure is incoporated into the model. The methodology allows us to learn about astrophysical parameters, such as those governing the source luminosity function and the cosmic magnetic field.
We argue that from observations alone, only the transverse power spectrum $C_\ell(z_1,z_2)$ and the corresponding correlation function $\xi(\theta,z_1,z_2)$ can be measured and that these contain the full three dimensional information. We determine the two point galaxy correlation function at linear order in perturbation theory. Redshift space distortions are taken into account for arbitrary angular and redshift separations. We discuss the shape of the longitudinal and the transversal correlation functions which are very different from each other and from their real space counterpart. We then go on and suggest how to measure both, the Hubble parameter, $H(z)$, and the angular diameter distance, $D_A(z)$, separately from these correlation functions and perform an Alcock-Paczynski test.
We report the results of our numerical simulation of classical-dissipative dynamics of a charged particle subjected to a non-markovian stochastic forcing. We find that the system develops a steady-state orbital magnetic moment in the presence of a static magnetic field. Very significantly, the sign of the orbital magnetic moment turns out to be {\it paramagnetic} for our choice of parameters, varied over a wide range. This is shown specifically for the case of classical dynamics driven by a Kubo-Anderson type non-markovian noise. Natural spatial boundary condition was imposed through (1) a soft (harmonic) confining potential, and (2) a hard potential, approximating a reflecting wall. There was no noticeable qualitative difference. What appears to be crucial to the orbital magnetic effect noticed here is the non-markovian property of the driving noise chosen. Experimental realization of this effect on the laboratory scale, and its possible implications are briefly discussed. We would like to emphasize that the above steady-state classical orbital paramagnetic moment complements, rather than contradicts the Bohr-van Leeuwen (BvL) theorem on the absence of classical orbital diamagnetism in thermodynamic equilibrium.
We perform cosmological perturbation theory in Hassan-Rosen bimetric gravity for general homogeneous and isotropic backgrounds. In the de Sitter approximation, we obtain decoupled sets of massless and massive scalar gravitational fluctuations. Matter perturbations then evolve like in Einstein gravity. We perturb the future de Sitter regime by the ratio of matter to dark energy, producing quasi-de Sitter space. In this more general setting the massive and massless fluctuations mix. We argue that in the quasi-de Sitter regime, the growth of structure in bimetric gravity differs from that of Einstein gravity.
When a color superconductor of high density QCD is rotating, super- fluid vortices are inevitably created along the rotation axis. In the color-flavor locked phase realized at the asymptotically large chemical potential, there appear non-Abelian vortices carrying both circulations of superfluid and color magnetic fluxes. A family of solutions has a degeneracy characterized by the Nambu-Goldtone modes CP2, associ- ated with the color-flavor locked symmetry spontaneously broken in the vicinity of the vortex. In this paper, we study electromagnetic coupling of the non-Abelian vortices and find that the degeneracy is removed with the induced effective potential. We obtain one stable vortex solu- tion and a family of metastable vortex solutions, both of which carry ordinary magnetic fluxes in addition to color magnetic fluxes. We dis- cuss quantum mechanical decay of the metastable vortices by quantum tunneling, and compare the effective potential with the other known po- tentials, the quantum mechanically induced potential and the potential induced by the strange quark mass.
Links to: arXiv, form interface, find, astro-ph, recent, 1206, contact, help (Access key information)
(Abridged) We present measurements of the velocity dispersion profile (VDP) for different masses of galaxy groups, using ~16,000 groups selected from the final data release of the SDSS by Yang et al (2007). We divide the groups into subsamples according to the stellar mass of their central galaxy, and for each subsample we estimate the redshift-space cross-correlation function (CCF) with respect to a reference sample, xi(r_p,pi), as well as the projected CCF, w_p(r_p). An NFW profile plus a biased version of the linear mass autocorrelation function can well describe the observed w_p(r_p), providing an accurate determination of the real-space CCF xi(r). The one-dimensional VDP is extracted from the redshift distortion in xi(r_p,pi), by comparing xi(r_p,pi) with xi(r). We find that the velocity dispersion (VD) within the virial radius (R_200) shows a roughly flat profile, with a slight increase at <~ 0.3 R_200. The average VD within R_200, sigma_v, is a strongly increasing function of central galaxy mass, ranging from ~130km/s at M*~3x10^10 M_sun, up to ~650km/s at M*~5x10^11 M_sun. We extend our analysis to N-body simulations with the LCDM cosmology but different density fluctuation amplitudes, sigma_8. We apply the abundance matching model to assign a stellar mass to the central subhalos. We measure the VDP for the halos and compare the results to the data. We show that the sigma_v-M* relation provides stringent constraints on both sigma_8 and sigma_ms, the dispersion in log(M*) of central galaxies at fixed halo mass. Our data suggests sigma_8=0.86+/-0.03 and sigma_ms=0.16+/-0.03. The slightly higher sigma_8 compared to the WMAP7 result might be due to a smaller matter density parameter assumed in our simulations. Our sigma_v measurements also provide a direct measure of the halo masses, agreeing well with the results from the galaxy-galaxy lensing analysis by Mandelbaum et al. (2006).
We report the first three dimensional radiation magnetohydrodynamic (RMHD) simulations of protostellar collapse with and without the Ohmic dissipation. We take many physical processes required to study star formation processes including a realistic equation of state into account and we follow the evolution from molecular cloud cores until protostellar cores are formed with sufficiently high resolutions without introducing a sink particle. The physical processes involved in the simulations and adopted numerical methods are described in detail. We can calculate only about one year after the formation of the protostellar cores with our direct 3D RMHD simulations because of the extremely short timescale in the deep interior of the formed protostellar cores, but successfully reveal the early phase of star formation processes. The thermal evolution and the structure of the first and second (protostellar) cores are consistent with previous one dimensional simulations using full radiation transfer, and considerably differ from preceding multi-dimensional studies with the barotropic approximation. The protostellar cores evolve virtually spherically symmetric in the ideal MHD models due to efficient angular momentum transport by magnetic fields, but the Ohmic dissipation enables the formation of the circumstellar disks in the vicinity of the protostellar cores as in previous MHD studies with the barotropic approximation. We also confirm that two different types of outflows are naturally launched by magnetic fields from the first cores and protostellar cores in the resistive models.
We study the linear and nonlinear structure formation in the dilaton and symmetron modified gravity models using a generic parameterisation which describes a large class of scenarios using only a few parameters. For this we have modified the $N$-body code {\tt ECOSMOG} to perform a set of 110 simulations for different models and parameter values, including the default $\Lambda$CDM. These simulations enable us to explore a large portion of the parameter space. We have studied the effects of modified gravity on the matter power spectrum and mass function, and found rich and interesting phenomenology. The results reveal the qualitative behaviour of the models, and provide guidance for future more detailed studies.
LAMOST (Large sky Area Multi-Object fiber Spectroscopic Telescope) is a Chinese national scientific research facility operated by National Astronomical Observatories, Chinese Academy of Sciences (NAOC). After two years of commissioning beginning in 2009, the telescope, instruments, software systems and operations are nearly ready to begin the main science survey. Through a spectral survey of millions of objects in much of the northern sky, LAMOST will enable research in a number of contemporary cutting edge topics in astrophysics, such as: discovery of the first generation stars in the Galaxy, pinning down the formation and evolution history of galaxies especially theMilky Way and its central massive black hole, looking for signatures of dark matter distribution and possible sub-structures in the Milky Way halo. To maximize the scientific potential of the facility, wide national participation and international collaboration has been emphasized. The survey has two major components: the LAMOST ExtraGAlactic Survey (LEGAS), and the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE). Until LAMOST reaches its full capability, the LEGUE portion of the survey will use the available observing time, starting in 2012. An overview of the LAMOST project and the survey that will be carried out in next five to six years is presented in this paper. The science plan for the whole LEGUE survey, instrumental specifications, site conditions, the descriptions of the current on-going pilot survey, including its footprints and target selection algorithm, will be presented as separate papers in this volume.
Satellite galaxies in groups and clusters are more likely to have low star formation rates (SFR) and lie on the `red-sequence' than central (`field') galaxies. Using galaxy group/cluster catalogs from the Sloan Digital Sky Survey Data Release 7, together with a high-resolution, cosmological N-body simulation to track satellite orbits, we examine the star formation histories and quenching timescales of satellites of Mstar > 5 x 10^9 Msol at z ~ 0. We first explore satellite infall histories: group preprocessing and ejected orbits are critical aspects of satellite evolution, and properly accounting for these, satellite infall typically occurred at z ~ 0.5, or ~5 Gyr ago. To obtain accurate initial conditions for the SFRs of satellites at their time of first infall, we construct an empirical parametrization for the evolution of central galaxy SFRs. With this, we constrain empirically the amount and efficiency of satellite quenching as a function of satellite and host halo mass, finding that satellite quenching is the dominant process for building up all quiescent galaxies at Mstar < 10^10 Msol. We then constrain satellite star formation histories, finding a `delayed-then-rapid' quenching scenario: satellite SFRs evolve unaffected for 2 - 4 Gyr after infall, after which star formation quenches rapidly, with an e-folding time of < 0.8 Gyr. These quenching timescales are shorter for more massive satellites but do not depend on host halo mass: the observed increase in satellite quiescent fraction with halo mass arises simply because of group preprocessing, which is responsible for up to half of quenched satellites in massive clusters. Because of the long time delay before quenching starts, satellites experience significant stellar mass growth after infall, nearly identical to central galaxies. This fact is a key reason for the success of the subhalo abundance matching technique.
We demonstrate that dwarf galaxies (10^7 < M_stellar < 10^9 Msun) with no active star formation are extremely rare (<0.06%) in the field. Our sample is based on the NASA-Sloan Atlas which is a re-analysis of the Sloan Digital Sky Survey Data Release 8. We examine the relative number of quenched versus star forming dwarf galaxies, defining quenched galaxies as having no Halpha emission (EW_Halpha < 2 AA) and a strong 4000AA-break. The fraction of quenched dwarf galaxies decreases rapidly with increasing distance from a massive host, leveling off for distances beyond 1.5 Mpc. We define galaxies beyond 1.5 Mpc of a massive host galaxy to be in the field. We demonstrate that there is a stellar mass threshold of M_stellar < 1.0x10^9 Msun below which quenched galaxies do not exist in the field. Below this threshold, we find that none of the 2951 field dwarf galaxies are quenched; all field dwarf galaxies show evidence for recent star formation. Correcting for volume effects, this corresponds to a 1-sigma upper limit on the quenched fraction of 0.06%. In more dense environments, quenched galaxies account for 23% of the dwarf population over the same stellar mass range. The majority of quenched dwarf galaxies (often classified as dwarf elliptical galaxies) are within 2 virial radii of a massive galaxy, and only a few percent of quenched dwarf galaxies exist beyond 4 virial radii. Thus, for galaxies with stellar mass less than 1.0x10^9 Msun, ending star-formation requires the presence of a more massive neighbor, providing a stringent constraint on models of star formation feedback.
The weather at Xinglong Observing Station, where the Guo Shou Jing Telescope (GSJT) is located, is strongly affected by the monsoon climate in north-east China. The LAMOST survey strategy is constrained by these weather patterns. In this paper, we present a statistics on observing hours from 2004 to 2007, and the sky brightness, seeing, and sky transparency from 1995 to 2011 at the site. We investigate effects of the site conditions on the survey plan. Operable hours each month shows strong correlation with season: on average there are 8 operable hours per night available in December, but only 1-2 hours in July and August. The seeing and the sky transparency also vary with seasons. Although the seeing is worse in windy winters, and the atmospheric extinction is worse in the spring and summer, the site is adequate for the proposed scientific program of LAMOST survey. With a Monte Carlo simulation using historical data on the site condition, we find that the available observation hours constrain the survey footprint from 22h to 16h in right ascension; the sky brightness allows LAMOST to obtain the limit magnitude of V = 19.5mag with S/N = 10.
We describe the footprint and input catalog for bright nights in the LAMOST Pilot Survey, which began in October 2011. Targets are selected from two stripes in the north and south Galactic Cap regions, centered at $\alpha$= 29$^\circ$, with 10$^\circ$ width in declination, covering right ascension of 135$^\circ-290^\circ$ and -30$^\circ$ to 30$^\circ$ respectively. We selected spectroscopic targets from a combination of the SDSS and 2MASS point source catalogs. The catalog of stars defining the field centers (as required by the Shack-Hartmann wavefront sensor at the center of the LAMOST field) consists of all V < 8m stars from the Hipparcos catalog. We employ a statistical selection algorithm that assigns priorities to targets based on their positions in multidimensional color/magnitude space. This scheme overemphasizes rare objects and de-emphasizes more populated regions of magnitude and color phase space, while ensuring a smooth, well-understood selection function. A demonstration of plate design is presented based on the Shack-Hartmann star catalog and an input catalog that was generated by our target selection routines.
We outline the design of the dark nights portion of the LAMOST Pilot Survey, which began observations in October 2011. In particular, we focus on Milky Way stellar candidates that are targeted for the LEGUE (LAMOST Experiment for Galactic Understanding and Exploration) survey. We discuss the regions of sky in which spectroscopic candidates were selected, and the motivations for selecting each of these sky areas. Some limitations due to the unique design of the telescope are discussed, including the requirement that a bright (V < 8) star be placed at the center of each plate for wavefront sensing and active optics corrections. The target selection categories and scientific goals motivating them are briefly discussed, followed by a detailed overview of how these selection functions were realized. We illustrate the difference between the overall input catalog - Sloan Digital Sky Survey (SDSS) photometry - and the final targets selected for LAMOST observation.
We describe a general target selection algorithm that is applicable to any survey in which the number of available candidates is much larger than the number of objects to be observed. This routine aims to achieve a balance between a smoothly-varying, well-understood selection function and the desire to preferentially select certain types of targets. Some target-selection examples are shown that illustrate different possibilities of emphasis functions. Although it is generally applicable, the algorithm was developed specifically for the LAMOST Experiment for Galactic Understanding and Exploration (LEGUE) survey that will be carried out using the Chinese Guo Shou Jing Telescope. In particular, this algorithm was designed for the portion of LEGUE targeting the Galactic halo, in which we attempt to balance a variety of science goals that require stars at fainter magnitudes than can be completely sampled by LAMOST. This algorithm has been implemented for the halo portion of the LAMOST pilot survey, which began in October 2011.
We describe the current plans for a spectroscopic survey of millions of stars in the Milky Way galaxy using the Guo Shou Jing Telescope (GSJT, formerly the Large Area Multi-Object Spectroscopic Telescope - LAMOST). The survey will obtain spectra for 2.5 million stars brighter than $r<19$ during dark/grey time, and 5 million stars brighter than $r<17$ or $J<16$ on nights that are moonlit or have low transparency. The survey will begin in fall of 2012, and will run for at least four years. The telescope design constrains the optimal declination range for observations to $10^\circ<\delta<50^\circ$, and site conditions lead to an emphasis on stars in the direction of the Galactic anticenter. The survey is divided into three parts with different target selection strategies: disk, anticenter, and spheroid. The resulting dataset will be used to study the merger history of the Milky Way, the substructure and evolution of the disks, the nature of the first generation of stars through identification of the lowest metallicity stars, and star formation through study of open clusters and the OB associations. Detailed design of the LEGUE survey will be completed after a review of the results of the pilot survey in summer 2012.
Recent deep Hubble Space Telescope WFC3 imaging suggests that a majority of compact quiescent massive galaxies at z~2 may contain disks. To investigate this claim, we have compared the ellipticity distribution of 31 carefully selected high-redshift massive quiescent compact galaxies to a set of mass-selected ellipticity and Sersic index distributions obtained from 2D structural fits to ~40,000$ nearby galaxies from the Sloan Digital Sky Survey. A Kolmogorov-Smirnov test shows that the distribution of ellipticities for the high-redshift galaxies is consistent with the ellipticity distribution of a similarly chosen sample of massive early-type galaxies. However the distribution of Sersic indices for the high-redshift sample is inconsistent with that of local early-type galaxies, and instead resembles that of local disk-dominated populations. The mismatch between the properties of high-redshift compact galaxies and those of both local early-type and disk-dominated systems leads us to conclude that the basic structures of high-redshift compact galaxies probably do not closely resemble those of any single local galaxy population. Any galaxy population analog to the high-redshift compact galaxies that exists at the current epoch is either a mix of different types of galaxies, or possibly a unique class of objects on their own.
Kilometer-sized moonlets in Saturn's A ring create S-shaped wakes called "propellers" in surrounding material. The Cassini spacecraft has tracked the motions of propellers for several years and finds that they deviate from Keplerian orbits having constant semimajor axes. The inferred orbital migration is known to switch sign. We show using a statistical test that the time series of orbital longitudes of the propeller Bl\'eriot is consistent with that of a time-integrated Gaussian random walk. That is, Bl\'eriot's observed migration pattern is consistent with being stochastic. We further show, using a combination of analytic estimates and collisional N-body simulations, that stochastic migration of the right magnitude to explain the Cassini observations can be driven by encounters with ring particles 10-20 m in radius. That the local ring mass is concentrated in decameter-sized particles is supported on independent grounds by occultation analyses.
We use STEREO imagery to study the morphology of a shock driven by a fast coronal mass ejection (CME) launched from the Sun on 2011 March 7. The source region of the CME is located just to the east of a coronal hole. The CME ejecta is deflected away from the hole, in contrast with the shock, which readily expands into the fast outflow from the coronal hole. The result is a CME with ejecta not well centered within the shock surrounding it. The shock shape inferred from the imaging is compared with in situ data at 1 AU, where the shock is observed near Earth by the Wind spacecraft, and at STEREO-A. Shock normals computed from the in situ data are consistent with the shock morphology inferred from imaging.
We present the ground-based detection of the secondary eclipse of the transiting exoplanet WASP-19b. The observations were made in the Sloan z'-band using the ULTRACAM triple-beam CCD camera mounted on the NTT. The measurement shows a 0.088\pm0.019% eclipse depth, matching previous predictions based on H- and K-band measurements. We discuss in detail our approach to the removal of errors arising due to systematics in the data set, in addition to fitting a model transit to our data. This fit returns an eclipse centre, T0, of 2455578.7676 HJD, consistent with a circular orbit. Our measurement of the secondary eclipse depth is also compared to model atmospheres of WASP-19b, and is found to be consistent with previous measurements at longer wavelengths for the model atmospheres we investigated.
Observations of the thermal X-ray emission from old radio pulsars implicate
that the size of hot spots is much smaller then the size of the polar cap that
follows from the purely dipolar geometry of pulsar magnetic field. Plausible
explanation of this phenomena is an assumption that the magnetic field at the
stellar surface differs essentially from the purely dipolar field. Using the
conservation of the magnetic flux through the area bounded by open magnetic
field lines we can estimate the surface magnetic field as of the order of
$10^{14}$G. Based on observations that the hot spot temperature is about a few
million Kelvins the Partially Screened Gap (PSG) model was proposed which
assumes that the temperature of the actual polar cap equals to the so called
critical temperature. We discuss correlation between the temperature and
corresponding area of the thermal X-ray emission for a number of pulsars.
We have found that depending on the conditions in a polar cap region the gap
breakdown can be caused either by the Curvature Radiation (CR) or by the
Inverse Compton Scattering (ICS). When the gap is dominated by ICS the density
of secondary plasma with Lorentz factors $10^{2}-10^{3}$ is at least an order
of magnitude higher then in a CR scenario. We believe that two different gap
breakdown scenarios can explain the mode-changing phenomenon and in particular
the pulse nulling. Measurements of the characteristic spacing between
sub-pulses ($P_{2}$) and the period at which a pattern of pulses crosses the
pulse window ($P_{3}$) allowed us to determine more strict conditions for
avalanche pair production in the PSG.
In this {\it Letter}, we present sensitive millimeter SiO (J=5-4; $\nu$=0) line observations of the outflow arising from the enigmatic object Orion Source I made with the Atacama Large Millimeter/Submillimeter Array (ALMA). The observations reveal that at scales of a few thousand AU, the outflow has a marked "butterfly" morphology along a northeast-southwest axis. However, contrary to what is found in the SiO and H$_2$O maser observations at scales of tens of AU, the blueshifted radial velocities of the moving gas are found to the northwest, while the redshifted velocities are in the southeast. The ALMA observations are complemented with SiO (J=8-7; $\nu$=0) maps (with a similar spatial resolution) obtained with the Submillimeter Array (SMA). These observations also show a similar morphology and velocity structure in this outflow. We discuss some possibilities to explain these differences at small and large scales across the flow.
Quasars with redshifts greater than four are rare, and can provide us important helps in probing the nature of the early universe. Here we report the discovery of six new quasars with $i$-band magnitudes brighter than 19.5 and redshifts between 2.4 and 4.6 from the spectroscopic observations made with the YFOSC instrument of the Lijiang 2.4m telescope in February, 2012. These six quasars are in the list of high redshift ($z>3.6$) quasar candidates selected by using our proposed J-K/i-Y criterion and the photometric redshift estimations from the SDSS optical and UKIDSS near-IR photometric data. Nine candidates were observed by YFOSC, and five among six new quasars were identified as $z>3.6$ quasars. One of the other three objects was identified as a star and the rest two were unidentified due to the lower signal-to-noise ratio of their spectra. This is the first time to discover quasars with redshifts greater than four with a telescope in China. Thanks to the Chinese Telescope Access Program (TAP), the redshift of 4.6 for one of these new quasars was confirmed by the MMT Red Channel spectroscopy. The continuum and emission line properties of these six quasars, as well as their central black hole masses and Eddington ratios, were obtained.
In coming years, Australia may find the need to build new optical telescopes to continue local programmes, contribute to global survey projects, and form a local multi-wavelength connection for the new radio telescopes being built. In this study, we refine possible locations for a new optical telescope by studying remotely sensed meteorological infrared data to ascertain expected cloud coverage rates across Australia, and combine these data with a Digital Elevation Model using a Geographic Information System. We find that the best sites within Australia for building optical telescopes are likely to be on the highest mountains in the Hamersley Range in Northwest Western Australia, while the MacDonnell Ranges in the Northern Territory may also be appropriate. We believe that similar seeing values to Siding Spring should be obtainable and with significantly more observing time at the identified sites. We expect to find twice as many clear nights as at current telescope sites. These sites are thus prime locations for future on-site testing.
Characterizing atmospheres beyond the Solar System is an endeavor no longer confined to the realm of science fiction.
Using the Subaru 8.2m Telescope with an IRCS Echelle spectrograph, we obtained high-resolution (R=10,000) near-infrared (1.01-1.38 \mu m) spectra of images A and B of the gravitationally lensed QSO B1422+231 (z=3.628) consisting of four known lensed images. We detected MgII absorption lines at z=3.54, which show a large variance of column densities (~ 0.3 dex) and velocities (~ 10 km/s) between the sightlines A and B with a projected separation of only 8.4h_{70}^{-1} pc at the redshift. This is the smallest spatial structure of the high-z gas clouds ever detected after Rauch et al. found a 20-pc scale structure for the same z=3.54 absorption system using optical spectra of images A and C. The observed systematic variances imply that the system is an expanding shell as originally suggested by Rauch et al. By combining the data for three sightlines, we managed to constrain the radius and expansion velocity of the shell (~ 50-100 pc, 130 km/s), concluding that the shell is truly a supernova remnant (SNR) rather than other types of shell objects, such as a giant HII region. We also detected strong FeII absorption lines for this system, but with much broader Doppler width than that of \alpha-element lines. We suggest that this FeII absorption line originates in a localized FeII-rich gas cloud that is not completely mixed with plowed ambient interstellar gas clouds showing other \alpha-element low-ion absorption lines. Along with the Fe richness, we conclude that the SNR is produced by an SNIa explosion.
We present results of survey of interstellar absorptions towards supernova remnant (SNR) RX J0852.0-4622. The distribution of KI absorbers along the distance of the background stars is indicative of a local region (d<600pc) strongly depopulated by KI line-absorbing clouds. This fact is supported by the behavior of the interstellar extinction. We find four high-velocity CaII components with velocities of >100km/s towards three stars and identify them with shocked clouds of Vela SNR. We reveal and measure acceleration of two shocked clouds at the approaching and receding sides of Vela SNR along the same sight line. The clouds acceleration, velocity, and CaII column density are used to probe cloud parameters. The total hydrogen column density of both accelerating clouds is found to be similar (~6*10^{17} cm$^{-2}$) which indicates that possibly there is a significant amount of small-size clouds in the vicinity of Vela SNR.
In the present paper a self-consistent theory, interpreting the VERITAS observations of the very high energy pulsed emission from the Crab pulsar is considered. The photon spectrum between 10MeV and 400GeV can be described by two power-law functions with the spectral indexes equal to 2 and 3.8. The source of the pulsed emission above 10MeV is assumed to be the synchrotron radiation, which is generated near the light cylinder during the quasi-linear stage of the cyclotron instability. The emitting particles are the primary beam electrons with the Lorentz factors up to $10^{9}$. Such high energies by beam particles is supposed to be reached due to Landau damping of the centrifugally induced Langmuir waves. This mechanism provides simultaneous generation of low (radio) and high energy (10MeV-400GeV) emission on the light cylinder scales, in one location of the pulsar magnetosphere.
We discuss the feasibility of using solar-type main-sequence stars as probes of fundamental physics and unification. We use a simple polytropic stellar structure model and study its sensitivity to variations of the gravitational, strong and electroweak coupling constants in the context of unification scenarios. We quantify the sensitivity of the Sun's interior temperature to these variations, finding $|\Delta\alpha/\alpha|\lesssim1.3\times10^{-4}$ for a 'canonical' choice of unification scenario, and discuss prospects for future improvements.
Galaxy diversification proceeds by transforming events like accretion, interaction or mergers. These explain the formation and evolution of galaxies that can now be described with many observables. Multivariate analyses are the obvious tools to tackle the datasets and understand the differences between different kinds of objects. However, depending on the method used, redundancies, incompatibilities or subjective choices of the parameters can void the usefulness of such analyses. The behaviour of the available parameters should be analysed before an objective reduction of dimensionality and subsequent clustering analyses can be undertaken, especially in an evolutionary context. We study a sample of 424 early-type galaxies described by 25 parameters, ten of which are Lick indices, to identify the most structuring parameters and determine an evolutionary classification of these objects. Four independent statistical methods are used to investigate the discriminant properties of the observables and the partitioning of the 424 galaxies: Principal Component Analysis, K-means cluster analysis, Minimum Contradiction Analysis and Cladistics. (abridged)
We present new methods for radiative transfer on hierarchial grids. We develop a new method for calculating the scattered flux that employs the grid structure to speed up the computation. We describe a novel subiteration algorithm that can be used to accelerate calculations with strong dust temperature self-coupling. We compute two test models, a molecular cloud and a circumstellar disc, and compare the accuracy and speed of the new algorithms against existing methods. An adaptive model of the molecular cloud with less than 8 % of the cells in the uniform grid produced results in good agreement with the full resolution model. The relative RMS error of the surface brightness <4 % at all wavelengths, and in regions of high column density the relative RMS error was only 10^{-4}. Computation with the adaptive model was faster by a factor of ~5. The new method for calculating the scattered flux is faster by a factor of ~4 in large models with a deep hierarchy structure, when images of the scattered light are computed towards several observing directions. The efficiency of the subiteration algorithm is highly dependent on the details of the model. In the circumstellar disc test the speed-up was a factor of two, but much larger gains are possible. The algorithm is expected to be most beneficial in models where a large number of small, dense regions are embedded in an environment with a lower mean density.
The unipolar induction DC circuit model, originally developed by Goldreich & Lynden-Bell for the Jupiter-Io system, has been applied to different types of binary systems in recent years. We show that there exists an upper limit to the magnetic interaction torque and energy dissipation rate in such model. This arises because when the resistance of the circuit is too small, the large current flow severely twists the magnetic flux tube connecting the two binary components, leading to breakdown of the circuit. Applying this limit, we find that in coalescing neutron star binaries, magnetic interactions produce negligible correction to the phase evolution of the gravitational waveform, even for magnetar-like field strengths. However, energy dissipation in the binary magnetosphere may still give rise to electromagnetic radiation prior to the final merger. For ultra-compact white dwarf binaries, we find that DC circuit does not provide adequate energy dissipation to explain the observed X-ray luminosities of several sources. For exoplanetary systems containing close-in Jupiters or super-Earths, magnetic torque and dissipation are negligible, except possibly during the early T Tauri phase, when the stellar magnetic field is stronger than 10^3G.
Possibility of gravitational repulsion in General Relativity is discussed and astronomical data in favor of cosmological acceleration are described. The problem of vacuum energy is emphasized and possible ways of its solution are indicated. The main attention is payed to adjustment mechanism which in principle could compensate originally huge vacuum energy down to cosmologically acceptable value and to solve the coincidence problem of a close magnitudes of the non-compensated remnants of vacuum energy and the energy density of the universe at the present time. Finally possible modifications of gravity at large scales which could induce accelerated cosmological expansion are considered.
The ACME Spectra project provides absolutely calibrated, mostly empirical spectra of exoplanet host stars for use in analysis of the stars and their planets. Spectra are obtained from ground-based telescopes and are tied directly to calibrated ground- and space-based photometry. The spectra remain only "mostly" empirical because of telluric absorption, but interpolation of stellar models over the gaps in wavelength coverage provides continuous stellar spectra. Among other uses, the spectra are suitable for precisely converting observed secondary eclipses (occultations) into absolute flux units with minimal recourse to models. In this letter I introduce ACME's methods and present a calibrated spectrum of the nearby, super-Earth hosting star 55 Cancri that spans the range from 0.81-5.05 micron. With this spectrum I show that the brightness temperature of the transiting planet 55 Cnc e is 1950 (+260/-190) K at 4.5 micron (cooler than previously reported), which corresponds to a planetary flux of 0.44 (+0.12/-0.08) mJy. This result suggests the planet has some combination of a nonzero albedo, a moderately efficient redistribution of absorbed stellar irradiation, and/or an optically thick atmosphere.
The ultracompact minihalos would be formed during the earlier epoch of cosmology. The dark matter annihilation within them is very strong due to the high density. The electron and positron with high energy from the dark matter annihilation can inverse Compton scatter (ICS) the background photons, such as CMB photons, into the higher energy. On the other hand, the synchrotron radiation can also be produced when they meet the magnetic field. In this paper, we study the signals from the UCMHs due to the dark matter annihilation into the radio, x-ray and {\gamma}-ray bands. We found that for the radio emission the UCMHs can provide one kinds of sources of the radio excess observed by ARCADE 2. But for these dark matter models, the x-ray signals due to the ICS effect and the {\gamma}-ray signals mainly due to the prompt emission from dark matter would excess the present observations, such as Fermi, COMPTEL and CHANDRA.
In this study, we present the unpublished flare data collected from 222 flares detected in the B band observations of five stars and the results derived by statistical analysis and modeling of these data. Six basic properties have been found with a statistical analysis method applied to all models and analyses for the flares detected in the B band observation of UV Ceti type stars. We have also compared the U and B bands with the analysis results. This comparison allowed us to evaluate the methods used in the analyses. The analyses provided the following results. (1) The flares were separated into two types, fast and slow flares. (2) The mean values of the equivalent durations of the slow and the fast flares differ by a factor of 16.2 \pm 3.7. (3) Regardless of the total flare duration, the maximum flare energy can reach a different Plateau level for each star. (4) The Plateau values of EV Lac and EQ Peg are higher than the others. (5) The minimum values of the total flare duration increase toward the later spectral types. This value is called the Half-Life value in models. (6) Both the maximum flare rise times and the total flare duration obtained from the observed flares decrease toward the later spectral types.
LSI+61303 is one of the few GeV- and TeV-emitting X-ray binaries with a prominent, well-studied modulated radio and gamma-ray emission. Changes in its radio morphology suggested in the past the hypothesis of a precessing microquasar. In 2006, a set of VLBA observations performed all around the orbit were not used to study the precession because the souce was explained in the context of the pulsar model, the alternative model for this system. However, a successive radio spectral index analysis has confirmed the predictions of the microquasar scenario in LSI+61303. At the light of these results we reanalysed the set of VLBA observations that constitutes a unique tool to determine the precession period and render a better understanding of the physical mechanism behind the precession. We improved the dynamic range of the images by a factor of four using self-calibration, and the self-calibrated maps reveal, in six out of ten images, a double-sided structure. The double-sided structure has variable position angle and switches at some epochs to a one-sided structure. These variations indicate a scenario where the precessing jet, inducing variable Doppler boosting, points close to our line of sight - a microblazar, the galactic version of the extra-galactic blazars. High energy observations of LSI+61303 are consistent with the microblazar nature of this object. Moreover, we suggest in LSI+61303 the first case of core shift effect observed in a microquasar. Because of this effect, well known in AGN, the cm-core of the jet is rather displaced from the system center. In LSI+61303, the cm-core of the jet traces a large ellipse, 7 times larger than the orbit, in a period of about 28 d. Our hypothesis is that this ellipse is the cross-section of the precession cone of the jet at the distance of the 3.6 cm-core, and its period is the precession period.
A vertical X-shaped structure was recently reported in the Galactic bulge. Here we present evidence of a similar X-shaped structure in the Shen et al. (2010) bar/boxy bulge model that simultaneously matches the stellar kinematics successfully. The X-shaped structure is found in the central region of our bar/boxy bulge model, and is qualitatively consistent with the observed one in many aspects. End-to-end separations of the X-shaped structure in the radial and vertical directions are roughly 3 kpc and 2 kpc, respectively. The X-shaped structure contains about 7% of light in the boxy bulge region, but it is significant enough to be identified in observations. An X-shaped structure naturally arises in the formation of bar/boxy bulges, and is mainly associated with orbits trapped around the vertically-extended x_1 family. Like the bar in our model, the X-shaped structure tilts away from the Sun--Galactic center line by 20 degrees. The X-shaped structure becomes increasingly symmetric about the disk plane, so the observed symmetry may indicate that it formed at least a few billion years ago. The existence of the vertical X-shaped structure suggests that the formation of the Milky Way bulge is shaped mainly by internal disk dynamical instabilities.
A possible solution to the small scale problems of the cold dark matter (CDM) scenario is that the dark matter consists of two components, a cold and a warm one. We perform a set of high resolution simulations of the Milky Way halo varying the mass of the WDM particle ($m_{\rm WDM}$) and the cosmic dark matter mass fraction in the WDM component ($\bar{f}_{\rm W}$). The scaling ansatz introduced in combined analysis of LHC and astroparticle searches postulates that the relative contribution of each dark matter component is the same locally as on average in the Universe (e.g. $f_{\rm W,\odot} = \bar{f}_{\rm W}$). Here we find however, that the normalised local WDM fraction ($f_{\rm W,\odot}$ / $\bar{f}_{\rm W}$) depends strongly on $m_{\rm WDM}$ for $m_{\rm WDM} <$ 1 keV. Using the scaling ansatz can therefore introduce significant errors into the interpretation of dark matter searches. To correct this issue a simple formula that fits the local dark matter densities of each component is provided.
The cosmic-ray intensity record from Belgrade muon detectors in the period 2002-2004 is subjected to power spectra analysis. Statistically significant peaks are found in datasets from both ground and underground (25 m w.e.) detectors. The possible origin of these periodicities is discussed.
We present Hubble Space Telescope imaging and spectroscopy for the extended high-ionization cloud known as Hanny's Voorwerp, near the spiral galaxy IC 2497. WFC3 images show complex dust absorption near the nucleus of IC 2497. STIS spectra show a type 2 Seyfert AGN of rather low luminosity. The ionization parameter log U = -3.5 is in accord with its weak X-ray emission. We find no high-ionization gas near the nucleus, adding to evidence that the AGN is currently at low radiative output (perhaps now dominated by kinetic energy). The nucleus is accompanied by an expanding ring of ionized gas 500 pc in projected diameter on the side opposite Hanny's Voorwerp, with Doppler offset 300 km/s from the nucleus (kinematic age < 7 x10^5 years). [O III] and H-alpha + [N II] images show fine structure in Hanny's Voorwerp, with limb-brightened sections and small areas where H-alpha is strong. We identify these as regions ionized by recent star formation, in contrast to the AGN ionization of the entire cloud. These candidate "normal" H II regions contain blue continuum objects, whose colors are consistent with young stellar populations; they appear only in a 2-kpc region toward IC 2497 in projection. The ionization-sensitive ratio [O III]/H-alpha shows no discernible pattern near the prominent "hole" in the ionized gas. The independence of ionization and surface brightness suggests that substantial spatial structure remains unresolved, to such an extent that the surface brightness sample the number of denser filaments rather than the characteristic density in emission regions. These results fit with our picture of an ionization echo from an AGN whose ionizing luminosity has dropped by a factor > 100 (and possibly much more) within the last 1-2 x 10^5 years; we suggest a sequence of events and discuss implications of such rapid fluctuations for AGN demographics. (Abridged)
So far, no systematic long-term blazar monitoring programs and detailed variability studies exist at sub-mm wavelengths. Here, we present a new sub-mm blazar monitoring program using the APEX 12-m telescope. A sample of about 40 gamma-ray blazars has been monitored since 2007/2008 with the LABOCA bolometer camera at 345 GHz. First light curves, preliminary variability results and a first comparison with the longer cm/mm bands (F-GAMMA program) are presented, demonstrating the extreme variability characteristics of blazars at such short wavelengths.
Galaxy structure and morphology is nearly always studied using the light originating from stars, however ideally one is interested in measuring structure using the stellar mass distribution. Not only does stellar mass trace out the underlying distribution of matter, it also minimises the effects of star formation and dust on the appearance and structure of a galaxy. We present in this paper a study of the stellar mass distributions and structures of galaxies at z<1 as found within the GOODS fields. We use pixel by pixel K-corrections to construct stellar mass and mass-to-light ratio maps of 560 galaxies of known morphology at magnitudes z_{850}<24. We measure structural and size parameters using these stellar mass maps, as well as on ACS BViz band imaging. This includes investigating the structural CAS-Gini-M_{20} parameters and half-light radius for each galaxy. We compare structural parameters and half-light radii in the ACS z_{850}-band and stellar mass maps, finding no systematic bias introduced by measuring galaxy sizes in z_{850}. We furthermore investigate relations between structural parameters in the ACS BViz bands and stellar mass maps, and compare our result to previous morphological studies. Combinations of various parameters in stellar mass generally reveal clear separations between early and late type morphologies, but cannot easily distinguish between star formation and dynamically disturbed systems. We also show that while ellipticals and early-type spirals have fairly constant CAS values at z<1 we find a tendency for late-type spiral and peculiar morphological types to have a higher A(M_{*}) at higher redshift. We argue that this, and the large fraction of peculiars that appear spiral-like in stellar mass maps, are possible evidence for either an active bulge formation in some late-type disks at z<1 or the presence of minor merger events.
This brief historical review highlights the early work of Hayashi and his associates on the thermal physics of star-forming clouds, as summarized in the temperature-density diagrams first presented by this group. Some of the more recent developments in this subject, including its application to understanding stellar masses and to understanding the formation of the first stars, are also briefly reviewed.
The analysis of narrowband drifting of type III-like structures in radio bursts dynamic spectra allows to obtain unique information about primary energy release mechanisms in solar flares. The SSRT spatially resolved images and a high spectral and temporal resolution allow direct determination not only the positions of its sources but also the exciter velocities along the flare loop. Practically, such measurements are possible during some special time intervals when the SSRT (about 5.7 GHz) is observing the flare region in two high-order fringes; thus, two 1D scans are recorded simultaneously at two frequency bands. The analysis of type III-like bursts recorded during the flare 14 Apr 2002 is presented. Using-muliwavelength radio observations recorded by SSRT, SBRS, NoRP, RSTN we study an event with series of several tens of drifting microwave pulses with drift rates in the range from -7 to 13 GHz/s. The sources of the fast-drifting bursts were located near the top of the flare loop in a volume of a few Mm in size. The slow drift of the exciters along the flare loop suggests a high pitch-anisotropy of the emitting electrons.
The energy spectrum of cosmic rays between 10**16 eV and 10**18 eV, derived from measurements of the shower size (total number of charged particles) and the total muon number of extensive air showers by the KASCADE-Grande experiment, is described. The resulting all-particle energy spectrum exhibits strong hints for a hardening of the spectrum at approximately 2x10**16 eV and a significant steepening at c. 8x10**16 eV. These observations challenge the view that the spectrum is a single power law between knee and ankle. Possible scenarios generating such features are discussed in terms of astrophysical processes that may explain the transition region from galactic to extragalactic origin of cosmic rays.
We probe the possibility that Centaurus A (Cen A) is a point source of ultra-high energy cosmic rays (UHECR) observed by PAO, through the statistical analysis of the arrival direction distribution. For this purpose, we set up the Cen A dominance model for the UHECR sources, in which Cen A contributes the fraction $f_{\rm C}$ of the whole UHECR with energy above $5.5\times10^{19}\,{\rm eV}$ and the isotropic background contributes the remaining $1-f_{\rm C}$ fraction. The effect of the intergalactic magnetic fields on the bending of the trajectory of Cen A originated UHECR is parameterized by the gaussian smearing angle $\theta_s$. Using the correlational angular distance distribution (CADD), we identify the excess of UHECR in the Cen A direction and fit the CADD of the observed PAO data by varying two parameters $f_{\rm C}$ and $\theta_s$ of the Cen A dominance model. The best-fit parameter values are $f_{\rm C}\approx0.1$ (The corresponding Cen A fraction observed at PAO is $f_{\rm C,PAO}\approx0.15$, that is, about 10 out of 69 UHECR.) and $\theta_s=5^\circ$ with the maximum probability $P_{\rm max}=0.29$. Considering the uncertainty concerning the assumption of isotropic background in the Cen A dominance model, we extend the viable parameter ranges to the $2\sigma$ band, $0.09\lesssim f_{\rm C,PAO}\lesssim 0.25$ and $0^\circ\lesssim \theta_s\lesssim 20^\circ$. This result supports the existence of a point source extended by the intergalactic magnetic fields in the direction of Cen A. If Cen A is actually the source responsible for the observed excess of UHECR, the average deflection angle of the excess UHECR implies the order of $10\,{\rm nG}$ intergalactic magnetic field in the vicinity of Cen A.
We present initial results from the CO survey toward high redshift galaxies using the Nobeyama 45m telescope. Using the new wide bandwidth spectrometer equipped with a two-beam SIS receiver, we have robust new detections of three high redshift (z=1.6-3.4) submillimeter galaxies (SXDF 1100.001, SDP9, and SDP17), one tentative detection (SDSS J160705+533558), and one non-detection (COSMOS-AzTEC1). The galaxies observed during the commissioning phase are sources with known spectroscopic redshifts from previous optical or from wide-band submm spectroscopy. The derived molecular gas mass and line widths from Gaussian fits are ~10^11 Msun and 430-530 km/s, which are consistent with previous CO observations of distant submm galaxies and quasars. The spectrometer that allows a maximum of 32 GHz instantaneous bandwidth will provide new science capabilities at the Nobeyama 45m telescope, allowing us to determine redshifts of bright submm selected galaxies without any prior redshift information.
We present our analyses of 15 months of Kepler data on KIC 10139564. We detected 57 periodicities with a variety of properties not previously observed all together in one pulsating subdwarf B star. Ten of the periodicities were found in the low-frequency region, and we associate them with nonradial g-modes. The other periodicities were found in the high-frequency region, which are likely p-modes. We discovered that most of the periodicities are components of multiplets with a common spacing. Assuming that multiplets are caused by rotation, we derive a rotation period of 25.6(1.8) days. The multiplets also allow us to identify the pulsations to an unprecedented extent for this class of pulsator. We also detect l<=2 multiplets, which are sensitive to the pulsation inclination and can constrain limb darkening via geometric cancellation factors. While most periodicities are stable, we detected several regions that show complex patterns. Detailed analyses showed these regions are complicated by several factors. Two are combination frequencies that originate in the superNyquist region and were found to be reflected below the Nyquist frequency. The Fourier peaks are clear in the superNyquist region, but the orbital motion of Kepler smears the Nyquist frequency in the barycentric reference frame and this effect is passed on to the subNyquist reflections. Others are likely multiplets but unstable in amplitudes and/or frequencies. The density of periodicities also make KIC 10139564 challenging to explain using published models. This menagerie of properties should provide tight constraints on structural models, making this subdwarf B star the most promising for applying asteroseismology.
We confirm that there are at least three separate low-latitude over-densities of blue F turnoff stars near the Milky Way anti-center: the Monoceros Ring, the Anti-Center Stream (ACS), and the Eastern Banded Structure (EBS). There might also be a small number of normal thick disk stars at the same location. The ACS is a tilted component that extends to higher Galactic latitude at lower Galactic longitude, 10 kpc from the Sun towards the anti-center. It has a sharp cutoff on the high latitude side. Distance, velocity, and proper motion measurements are consistent with previous orbit fits. The mean metallicity is [Fe/H]$=-0.96 \pm 0.03$, which is lower than the thick disk and Monoceros Ring. The Monoceros Ring is a higher density substructure that is present at $15\arcdeg<b<22\arcdeg$ at all longitudes probed in this survey. The structure likely continues towards lower latitudes. The distances are consistent with a constant distance from the Galactic Center of 17.6 kpc. The mean line-of-sight velocity of the structure is consistent with a thick disk rotation. However, the velocity dispersion of these stars is $\sim 15$ km s$^{-1}$, and the metallicity is [Fe/H]$=-0.80 \pm 0.01$. Both of these quantities are lower than the canonical thick disk. We suggest that this ring structure is likely different from the thick disk, though its association with the disk cannot be definitively ruled out. The Eastern Banded Structure (EBS) is detected primarily photometrically, near $(l,b)=(225\arcdeg,30\arcdeg)$, at a distance of 10.9 kpc from the Sun.
We have conducted a systematic search of low-mass black holes (BHs) in active galactic nuclei (AGNs) with broad Halpha emission lines, aiming at building a homogeneous sample that is more complete than previous ones for fainter, less highly accreting sources. For this purpose, we developed a set of elaborate, automated selection procedures and applied it uniformly to the Fourth Data Release of the Sloan Digital Sky Survey. Special attention is given to AGN--galaxy spectral decomposition and emission-line deblending. We define a sample of 309 type 1 AGNs with BH masses in the range $8 \times 10^4$--$2 \times 10^6$ \msun (with a median of $1.2 \times 10^6$ solar mass), using the virial mass estimator based on the broad Halpha line. About half of our sample of low-mass BHs differs from that of Greene & Ho, with 61 of them discovered here for the first time. Our new sample picks up more AGNs with low accretion rates: the Eddington ratios of the present sample range from $<~0.01$ to ~1, with 30% below 0.1. This suggests that a significant fraction of low-mass BHs in the local Universe are accreting at low rates. The host galaxies of the low-mass BHs have luminosities similar to those of $L^*$ field galaxies, optical colors of Sbc spirals, and stellar spectral features consistent with a continuous star formation history with a mean stellar age of less than 1 Gyr.
We have modelled the inner surface brightness profiles of 39 alleged `core' galaxies with the core-Sersic model, and provide new physical parameters for the largest ever sample of `core' galaxies fit with this model. When present, additional nuclear components were simultaneously modelled and the typical rms scatter of the fits (out to ~10 arcsec) is 0.02 mag/arcsec^2. Model-independent estimates of each core's break radius are shown to agree with those from the core-Sersic model, and a comparison with the Nuker model is provided. We found an absence of cores in what amounts to 18% of the sample which are reclassified here as Sersic galaxies with low values of n (< ~ 4) and thus shallow inner profile slopes. In general, galaxies with n<3 and sigma < 183 km/s do not have depleted cores. We derive updated relations between core-Sersic break radii, their associated surface brightness, bulge luminosity, central velocity dispersion, and predicted black hole mass for galaxies with depleted cores. With the possible exception of NGC 584, we confirm that the inner negative logarithmic profile slopes gamma are < ~ 0.3 for the `core' galaxies, and 0 > gamma > -0.1 for six of these. Finally, the central stellar mass deficits are found to have values typically within a factor of 4 of the expected central black hole mass.
The 22Ne({\alpha},n)25Mg reaction is an important source of neutrons for the s-process. In massive stars responsible for the weak component of the s-process, 22Ne({\alpha},n)25Mg is the dominant source of neutrons, both during core helium burning and in shell carbon burning. For the main s-process component produced in Asymptotic Giant Branch (AGB) stars, the 13C({\alpha},n)16O reaction is the dominant source of neutrons operating during the interpulse period, with the 22Ne+{\alpha} source affecting mainly the s-process branchings during a thermal pulse. Rate uncertainties in the competing 22Ne({\alpha},n)25Mg and 22Ne({\alpha},{\gamma})26Mg reactions result in large variations of s-process nucleosynthesis. Here, we present up-to-date and statistically rigorous 22Ne+{\alpha} reaction rates using recent experimental results and Monte Carlo sampling. Our new rates are used in post-processing nucleosynthesis calculations both for massive stars and AGB stars. We demonstrate that the nucleosynthesis uncertainties arising from the new rates are dramatically reduced in comparison to previously published results, but several ambiguities in the present data must still be addressed. Recommendations for further study to resolve these issues are provided.
The recently discovered subdwarf B (sdB) pulsator KIC11558725 features a rich
g-mode frequency spectrum, with a few low-amplitude p-modes at short periods,
and is a promising target for a seismic study aiming to constrain the internal
structure of this star, and of sdB stars in general.
We have obtained ground-based spectroscopic Balmer-line radial-velocity
measurements of KIC11558725, spanning the 2010 and 2011 observing seasons. From
these data we have discovered that KIC11558725 is a binary with period P=10.05
d, and that the radial-velocity amplitude of the sdB star is 58 km/s.
Consequently the companion of the sdB star has a minimum mass of 0.63 M\odot,
and is therefore most likely an unseen white dwarf.
We analyse the near-continuous 2010-2011 Kepler light curve to reveal orbital
Doppler-beaming light variations at the 238 ppm level, which is consistent with
the observed spectroscopic orbital radial-velocity amplitude of the subdwarf.
We use the strongest 70 pulsation frequencies in the Kepler light curve of the
subdwarf as clocks to derive a third consistent measurement of the orbital
radial-velocity amplitude, from the orbital light-travel delay.
We use our high signal-to-noise average spectra to study the atmospheric
parameters of the sdB star, deriving Teff = 27 910K and log g = 5.41 dex, and
find that carbon, nitrogen and oxygen are underabundant relative to the solar
mixture.
Furthermore, we extract more than 160 significant frequencies from the Kepler
light curve. We investigate the pulsation frequencies for expected period
spacings and rotational splittings. We find period-spacing sequences of
spherical-harmonic degrees \ell=1 and \ell=2, and we associate a large fraction
of the g-modes in KIC11558725 with these sequences. From frequency splittings
we conclude that the subdwarf is rotating subsynchronously with respect to the
orbit.
We have searched and reviewed all multi- wavelength data available for the region towards the gamma-ray source 2FGL J2056.7+4939 in order to con- strain its possible counterpart at lower energies. As a result, only a point-like optical/infrared source with flat-spectrum radio emission is found to be consistent with all X-ray and gamma-ray error circles. Its struc- ture is marginally resolved at radio wavelengths at the sub-arcsecond level. An extragalactic scenario appears to be the most likely interpretation for this object.
We present the results of a photometric and astrometric study of the low mass stellar and substellar population of the young open cluster Blanco 1. We have exploited J band data, obtained recently with the Wide Field Camera (WFCAM) on the United Kingdom InfraRed Telescope (UKIRT), and 10 year old I and z band optical imaging from CFH12k and Canada France Hawaii Telescope (CFHT), to identify 44 candidate low mass stellar and substellar members, in an area of 2 sq. degrees, on the basis of their colours and proper motions. This sample includes five sources which are newly discovered. We also confirm the lowest mass candidate member of Blanco 1 unearthed so far (29MJup). We determine the cluster mass function to have a slope of alpha=+0.93, assuming it to have a power law form. This is high, but nearly consistent with previous studies of the cluster (to within the errors), and also that of its much better studied northern hemisphere analogue, the Pleiades.
We construct axisymmetric mass models for dwarf spheroidal (dSph) galaxies in the Milky Way to obtain plausible limits on the non-spherical structure of their dark halos. This is motivated by the fact that the observed luminous parts of the dSphs are actually non-spherical and Cold Dark Matter (CDM) models predict non-spherical virialized dark halos. Our models consider velocity anisotropy of stars $\bar{v^2_R} / \bar{v^2_{\phi}}$, which can vary with the adopted cylindrical coordinates under the assumption $\bar{v^2_z}=\bar{v^2_R}$ for simplicity, and also include an inclination of the system as a fitting parameter to explain the observed line-of-sight velocity dispersion profile. Applying these models to six of the bright dSphs in the Milky Way, we find that the best-fitting cases for most of the dSphs yield oblate and flattened dark halos, irrespective of assumed density profiles in their central parts. We also find that the total mass of the dSphs enclosed within a spheroid with major-axis length of 300 pc varies from $10^6M_{\odot}$ to $10^7M_{\odot}$, contrary to the conclusion from spherical models. This suggests the importance of considering shapes of dark halos in mass models of the dSphs. It is also found that dark halos of the Galactic dSphs may be more flattened than N-body predictions, thereby implying our yet incomplete understanding of baryonic and/or non-baryonic dark matter physics in dwarf galaxy scales.
The thermal stability of a weakly-magnetized, rotating, stratified, optically-thin plasma is studied by means of linear-perturbation analysis. We derive dispersion relations and criteria for stability against axisymmetric perturbations that generalize previous results on either non-rotating or unmagnetized fluids. The implications for the hot atmospheres of galaxies and galaxy clusters are discussed.
The relation between the ratio of infrared (IR) and ultraviolet (UV) flux densities (the infrared excess: IRX) and the slope of the UV spectrum (\beta) of galaxies plays a fundamental role in the evaluation of the dust attenuation of star forming galaxies especially at high redshifts. Many authors, however, pointed out that there is a significant dispersion and/or deviation from the originally proposed IRX-\beta relation depending on sample selection. We reexamined the IRX-\beta relation by measuring the far- and near-UV flux densities of the original sample galaxies with GALEX and AKARI imaging data, and constructed a revised formula. We found that the newly obtained IRX values were lower than the original relation because of the significant underestimation of the UV flux densities of the galaxies, caused by the small aperture of IUE, Further, since the original relation was based on IRAS data which covered a wavelength range of \lambda = 42--122\mum, using the data from AKARI which has wider wavelength coverage toward longer wavelengths, we obtained an appropriate IRX-\beta relation with total dust emission (TIR): \log(L_{\rm TIR}/L_{\rm FUV}) = \log [10^{0.4(3.06+1.58\beta)}-1] +0.22. This new relation is consistent with most of the preceding results for samples selected at optical and UV, though there is a significant scatter around it. We also found that even the quiescent class of IR galaxies follows this new relation, though luminous and ultraluminous IR galaxies distribute completely differently as well known before.
The phonon contribution to the shear viscosity $\eta$ in superfluid neutron stars is calculated by assuming neutron pairing in a $^1S_0$ channel. The shear viscosity is obtained by means of variational methods for the solution of the Boltzmann equation amended by a collision term which takes into account the binary collisions of phonons. Effective field theory techniques are used to extract the phonon scattering rates in terms of the equation of state (EoS) of the system. We find that $\eta \propto 1/T^5$, the proportionality factor depending on the EoS of the system. Our results indicate that the phonon contribution to $\eta$ might have important effects for the different oscillation modes of the star.
A posteriori anisotropy study of ultra-high energy cosmic rays (UHECRs) with the Pierre Auger Observatory (PAO) has shown evidence of excess of cosmic ray particles above 55 EeV within $18^{\circ}$ of the direction of the radio galaxy Centaurus A. However, the origin of the excess remains elusive. We simulate the propagation of different species of particles coming from the direction of Centaurus A in the Galactic magnetic fields, and find that only particles of nuclear charge $Z\la 10$ can avoid being deflected outside of the $18^{\circ}$ window of Centaurus A. On the other hand, considering the increasingly heavy composition of UHECRs at the highest energies measured by PAO, a plausible scenario for cosmic rays from the direction of Centaurus A can be found if they consist of intermediate-mass nuclei. The chemical composition of cosmic rays can be further constrained by lower-energy cosmic rays of the same rigidity. We find that cosmic ray acceleration in the lobes of Centaurus A is not favored, while acceleration in the stellar winds that are rich in intermediate-mass nuclei, could meet the requirement. This suggests that the observed excess may originate from cosmic ray accelerators induced by stellar explosions in the star-forming regions of Centaurus A and/or the Centaurus cluster located behind Centaurus A.
We analyzed the initial rising behaviors of X-ray outbursts from two transient low-mass X-ray binaries (LMXBs) containing a neutron-star (NS), Aql X-1 and 4U 1608-52, which are continuously being monitored by MAXI/GSC in 2--20 keV, RXTE/ASM in 2--10 keV, and Swift/BAT in 15--50 keV. We found that the observed ten outbursts are classified into two types by the patterns of the relative intensity evolutions in the two energy bands below/above 15 keV. One type behaves as the 15--50 keV intensity achieves the maximum during the initial hard-state period and drops greatly at the hard-to-soft state transition. On the other hand, the other type does as both the 2--15 keV and the 15--50 keV intensities achieve the maximums after the transition. The former have the longer initial hard-state ($\gtrsim$ 9 d) than the latter's ($\ltsim$5 d). Therefore, we named them as slow-type (S-type) and fast-type (F-type), respectively. These two types also show the differences in the luminosity at the hard-to-soft state transition as well as in the average luminosity before the outburst started, where the S-type are higher than the F-type in the both. These results suggest that the X-ray radiation during the pre-outburst period, which heats up the accretion disk and delays the disk transition (i.e., from a geometrically thick disk to a thin one), would determine whether the following outburst becomes S-type or F-type. The luminosity when the hard-to-soft state transition occurs is higher than $\sim 8 \times10^{36}$ erg s$^{-1}$ in the S-type, which corresponds to 4% of the Eddington luminosity for a 1.4 \Mo NS.
Gamma-ray bursts (GRBs) that emit photons at GeV energies form a small but significant population of GRBs. However, the number of GRBs whose GeV-emitting period is simultaneously observed in X-rays remains small. We report gamma-ray observations of GRB 110625A using Fermi's Large Area Telescope (LAT) in the energy range 100 MeV to 20 GeV. Gamma-ray emission at these energies was clearly detected using data taken between 180s and 580s after the burst, an epoch after the prompt emission phase. The GeV light curve differs from a simple power-law decay, and probably consists of two emission periods. Simultaneous Swift/XRT observations did not show flaring behaviors as in the case of GRB 100728A. We discuss the possibility that the GeV emission is the synchrotron self-Compton radiation of underlying ultraviolet flares.
Relativistic X-ray disk-lines have been found in multiple neutron star low-mass X-ray binaries, in close analogy with black holes across the mass-scale. These lines have tremendous diagnostic power and have been used to constrain stellar radii and magnetic fields, often finding values that are consistent with independent timing techniques. Here, we compare CCD-based data from Suzaku with Fe K line profiles from archival data taken with gas-based spectrometers. In general, we find good consistency between the gas-based line profiles from EXOSAT, BeppoSAX and RXTE and the CCD data from Suzaku, demonstrating that the broad profiles seen are intrinsic to the line and not broad due to instrumental issues. However, we do find that when fitting with a Gaussian line profile, the width of the Gaussian can depend on the continuum model in instruments with low spectral resolution, though when the different models fit equally well the line widths generally agree. We also demonstrate that three BeppoSAX observations show evidence for asymmetric lines, with a relativistic disk-line model providing a significantly better fit than a Gaussian. We test this by using the posterior predictive p-value method, and bootstrapping of the spectra to show that such deviations from a Gaussian are unlikely to be observed by chance.
"Propellers" are features in Saturn's A ring associated with moonlets that open partial gaps. They exhibit non-Keplerian motion (Tiscareno 2010); the longitude residuals of the best-observed propeller, "Bl\'eriot," appear consistent with a sinusoid of period ~4 years. Pan and Chiang (2010) proposed that propeller moonlets librate in "frog resonances" with co-orbiting ring material. By analogy with the restricted three-body problem, they treated the co-orbital material as stationary in the rotating frame and neglected non-co-orbital material. Here we use simple numerical experiments to extend the frog model, including feedback due to the gap's motion, and drag associated with the Lindblad disk torques that cause Type I migration. Because the moonlet creates the gap, we expect the gap centroid to track the moonlet, but only after a time delay t_diff, the time for a ring particle to travel from conjunction with the moonlet to the end of the gap. We find that frog librations can persist only if t_diff exceeds the frog libration period P_lib, and if damping from Lindblad torques balances driving from co-orbital torques. If t_diff << P_lib, then the libration amplitude damps to zero. In the case of Bl\'eriot, the frog resonance model can reproduce the observed libration period P_lib ~ 4 yr. However, our simple feedback prescription suggests that Bl\'eriot's t_diff ~ 0.01P_lib, which is inconsistent with the observed libration amplitude of 260 km. We urge more accurate treatments of feedback to test the assumptions of our toy models.
We present general relativistic magnetohydrodynamic (GRMHD) numerical
simulations of the accretion flow around the supermassive black hole in the
Galactic centre, Sagittarius A* (Sgr A*). The simulations include for the first
time radiative cooling processes (synchrotron, bremsstrahlung, and inverse
Compton) self-consistently in the dynamics, allowing us to test the common
simplification of ignoring all cooling losses in the modeling of Sgr A*. We
confirm that for Sgr A*, neglecting the cooling losses is a reasonable
approximation if the Galactic centre is accreting below ~10^{-8} Msun/yr i.e.
Mdot < 10^{-7} Mdot_Edd. But above this limit, we show that radiative losses
should be taken into account as significant differences appear in the dynamics
and the resulting spectra when comparing simulations with and without cooling.
This limit implies that most nearby low-luminosity active galactic nuclei are
in the regime where cooling should be taken into account.
We further make a parameter study of axisymmetric gas accretion around the
supermassive black hole at the Galactic centre. This approach allows us to
investigate the physics of gas accretion in general, while confronting our
results with the well studied and observed source, Sgr A*, as a test case. We
confirm that the nature of the accretion flow and outflow is strongly dependent
on the initial geometry of the magnetic field. For example, we find it
difficult, even with very high spins, to generate powerful outflows from discs
threaded with multiple, separate poloidal field loops.
We report on a detailed spectral analysis of all the available XMM-Newton data of RX J1856.5-3754, the brightest and most extensively observed nearby, thermally emitting neutron star. Very small variations (~1-2%) in the single-blackbody temperature are detected, but are probably due to an instrumental effect, since they correlate with the position of the source on the detector. Restricting the analysis to a homogeneous subset of observations, with the source at the same detector position, we place strong limits on possible spectral or flux variations from March 2005 to present-day. A slightly higher temperature (kT~61.5 eV, compared to the average value kT~61 eV) was instead measured in April 2002. If this difference is not of instrumental origin, it implies a rate of variation of about 0.15 eV/yr between April 2002 and March 2005. The high-statistics spectrum from the selected observations is well fit by the sum of two blackbody models, which extrapolate to an optical flux level in agreement with the observed value.
This summary reports on papers presented at the Cool Stars-16 meeting in the splinter session "Solar and Stellar flares." Although many topics were discussed, the main themes were the commonality of interests, and of physics, between the solar and stellar flare communities, and the opportunities for important new observations in the near future.
We consider models of accelerating Universe elaborated for Finsler like gravity theories constructed on tangent bundles to Lorentz manifolds. In the osculating approximation, certain locally anisotropic configurations are similar to those for f(R) gravity. This allows us to generalize a proposal (by Nojiri, Odintsov and Saez-Gomez, arXiv: 1108.0767) in order to reconstruct and compare two classes of Einstein-Finsler gravity, EFG, and f(R) gravity theories using modern cosmological data and realistic physical scenarios. We conclude that EFG provides inflation, acceleration and little rip evolution scenarios with realistic alternatives to standard Lambda CDM cosmology. The approach is based on a proof that there is a general decoupling property of gravitational field equations in EFG and modified theories which allows us to generate off-diagonal cosmological solutions.
We present Spitzer Space Telescope mid-infrared IRS spectra, supplemented by ground-based optical observations, of the classical novae V1974 Cyg, V382 Vel, and V1494 Aql more than 11, 8, and 4 years after outburst respectively. The spectra are dominated by forbidden emission from neon and oxygen, though in some cases, there are weak signatures of magnesium, sulfur, and argon. We investigate the geometry and distribution of the late time ejecta by examination of the emission line profiles. Using nebular analysis in the low density regime, we estimate lower limits on the abundances in these novae. In V1974 Cyg and V382 Vel, our observations confirm the abundance estimates presented by other authors and support the claims that these eruptions occurred on ONe white dwarfs. We report the first detection of neon emission in V1494 Aql and show that the system most likely contains a CO white dwarf.
We revisit the Polish doughnut model of accretion disks providing a comprehensive analytical description of the Polish doughnut structure. We describe a perfect fluid circularly orbiting around a Schwarzschild black hole, source of the gravitational field, by the effective potential approach for the exact gravitational and centrifugal effects. This analysis leads to a detailed, analytical description of the accretion disk, its toroidal surface, the thickness, the distance from the source. We determine the variation of these features with the effective potential and the fluid angular momentum. Many analytical formulas are given. In particular it turns out that the distance from the source of the inner surface of the torus increases with increasing fluid angular momentum but decreases with increasing energy function defined as the value of the effective potential for that momentum. The location of torus maximum thickness moves towards the external regions of the surface with increasing angular momentum, until it reaches a maximum an then decreases. Assuming a polytropic equation of state we investigate some specific cases.
We report the first two-dimensional mapping of the stellar population and non-stellar continua within the inner 180 pc (radius) of NGC 1068 at a spatial resolution of 8 pc, using integral field spectroscopy in the near-infrared. We have applied the technique of spectral synthesis to data obtained with the instrument NIFS and the adaptive optics module ALTAIR at the Gemini North Telescope. Two episodes of recent star formation are found to dominate the stellar population contribution: the first occurred 300 Myr ago, extending over most of the nuclear region; the second occurred just 30 Myr ago, in a ring-like structure at ~100 pc from the nucleus, where it is coincident with an expanding ring of H2 emission. Inside the ring, where a decrease in the stellar velocity dispersion is observed, the stellar population is dominated by the 300 Myr age component. In the inner 35 pc, the oldest age component (age > 2Gyr) dominates the mass, while the flux is dominated by black-body components with temperatures in the range 700 < T < 800 K which we attribute to the dusty torus. We also find some contribution from black-body and power-law components beyond the nucleus which we attribute to dust emission and scattered light.
We use numerical simulations to study the effects of the patchiness of a partly reionized intergalactic medium (IGM) on the observability of Ly-alpha emitters (LAEs) at high redshifts (z ~ 6). We present a new model that divides the Ly-alpha radiative transfer into a (circum-)galactic and an extragalactic (IGM) part, and investigate how the choice of intrinsic line model affects the IGM transmission results. We use our model to study the impact of neutral hydrogen on statistical observables such as the Ly-alpha restframe equivalent width (REW) distribution, the LAE luminosity function and the two-point correlation function. We find that if the observed changes in LAE luminosity functions and equivalent width distributions between z ~ 6 and z ~ 7 are to be explained by an increased IGM neutral fraction alone, we require an extremely late and rapid reionization scenario, where the Universe was ~ 40 % ionized at z = 7, ~ 50 % ionized at z = 6.5 and ~ 100 % ionized at z = 6. This is in conflict with other observations, suggesting that intrinsic LAE evolution at z > 6 cannot be completely neglected. We show how the two-point correlation function can provide more robust constraints once future observations obtain larger LAE samples, and provide predictions for the sample sizes needed to tell different reionization scenarios apart.
We present the first time-series study of the ultra-faint dwarf galaxy Hercules. Using a variety of telescope/instrument facilities we secured about 50 V and 80 B epochs. These data allowed us to detect and characterize 10 pulsating variable stars in Hercules. Our final sample includes 6 fundamental-mode (ab-type) and 3 first overtone (c-type) RR Lyrae stars, and one Anomalous Cepheid. The average period of the ab-type RR Lyrae stars, < Pab >= 0.68 d (sigma = 0.03 d), places Hercules in the Oosterhoff II group, as found for almost the totality of the ultra-faint dwarf galaxies investigated so far for variability. The RR Lyrae stars were used to obtain independent estimates of the metallicity, reddening and distance to Hercules, for which we find: [Fe/H] = -2.30+-0.15 dex, E(B -V) = 0.09+-0.02 mag, and (m-M)o = 20.6+-0.1 mag, in good agreement with the literature values. We have obtained a V, B - V color-magnitude diagram (CMD) of Hercules that reaches V ~ 25 mag and extends beyond the galaxy's half-light radius over a total area of 40' {\times} 36'. The CMD and the RR Lyrae stars indicate the presence of a population as old and metal-poor as (at least) the Galactic globular clusters M68.
Type IIP (Plateau) Supernovae are the most commonly observed variety of core collapse events. They have been detected in a wide range of wavelengths from radio, through optical to X-rays. The standard picture of a type IIP supernova has the blastwave interacting with the progenitor's circumstellar matter to produce a hot region bounded by a forward and a reverse shock. This region is thought to be responsible for most of the X-ray and radio emission from these objects. Yet the origin of X-rays from these supernovae is not well understood quantitatively. The relative contributions of particle acceleration and magnetic field amplification in generating the X-ray and radio emission need to be determined. In this work we analyze archival Chandra observations of SN 2004dj, the nearest supernova since SN 1987A, along with published radio and optical information. We determine the pre-explosion mass loss rate, blastwave velocity, electron acceleration and magnetic field amplification efficiencies. We find that a greater fraction of the thermal energy goes into accelerating electrons than into amplifying magnetic fields. We conclude that the X-ray emission arises out of a combination of inverse Compton scattering by non-thermal electrons accelerated in the forward shock and thermal emission from supernova ejecta heated by the reverse shock.
The time domain has been identified as one of the most important areas of astronomical research for the next decade. The Virtual Observatory is in the vanguard with dedicated tools and services that enable and facilitate the discovery, dissemination and analysis of time domain data. These range in scope from rapid notifications of time-critical astronomical transients to annotating long-term variables with the latest modeling results. In this paper, we will review the prior art in these areas and focus on the capabilities that the VAO is bringing to bear in support of time domain science. In particular, we will focus on the issues involved with the heterogeneous collections of (ancillary) data associated with astronomical transients, and the time series characterization and classification tools required by the next generation of sky surveys, such as LSST and SKA.
We report the detection of the Doppler shadow of the transiting hot Jupiter CoRoT-11b. Our analysis is based on line-profile tomography of time-series, Keck/HIRES high-resolution spectra acquired during the transit of the planet. We measured a sky-projected, spin-orbit angle of 0.1 +/- 2.6 degrees, which is consistent with a very low-inclined orbit with respect to the stellar rotation axis. We refined the physical parameters of the system using a Markov chain Monte Carlo simultaneous fitting of the available photometric and spectroscopic data. An analysis of the tidal evolution of the system shows how the currently measured obliquity and its uncertainty translate into an initial absolute value of less than about 10 degrees on the zero-age main sequence, for an expected average modified tidal quality factor of the star Q'* > 4 x 10^6. This is indicative of an inward migration scenario that would not have perturbed the primordial low obliquity of CoRoT-11b. Taking into account the effective temperature and mass of the planet host star (Teff=6440 K, M*=1.23 MSun), the system can be considered a new telling exception to the recently proposed trend, according to which relatively hot and massive stars (Teff>6250 K, M*>1.2 MSun) seem to be preferentially orbited by hot Jupiters with high obliquity.
The cosmological origin of both dark and baryonic matter can be explained through a unified mechanism called hylogenesis where baryon and antibaryon number are divided between the visible sector and a GeV-scale hidden sector, while the Universe remains net baryon symmetric. The "missing" antibaryons, in the form of exotic hidden states, are the dark matter. We study model-building, cosmological, and phenomenological aspects of this scenario within the framework of supersymmetry, which naturally stabilizes the light hidden sector and electroweak mass scales. Inelastic dark matter scattering on visible matter destroys nucleons, and nucleon decay searches offer a novel avenue for the direct detection of the hidden antibaryonic dark matter sea.
We suggest that galileon theories should have an additional self-coupling of the fields to the trace of their own energy-momentum tensor. We explore the classical features of one such model, in flat 4D spacetime, with emphasis on solutions that are scalar analogues of gravitational geons. We discuss the stability of these scalar geons, and some of their possible signatures, including shock fronts.
We study perturbations around some cosmological backgrounds in the dRGT theory of massive gravity. We develop a general formalism to calculate the perturbations around any background. We derive the Lagrangian for fluctuations in the small scale limit, and for the open FRW solution we repeat the analysis around the full background. We find that the perturbations display similar properties: the longitudinal modes of the massive graviton are instantaneous at quadratic level, but they acquire a time-kinetic term at cubic order.
We derive a phenomenological expression that predicts the final mass of the black-hole remnant resulting from the merger of a generic binary system of black holes on quasi-circular orbits. Besides recovering the correct test-particle limit for extreme mass-ratio binaries, our formula reproduces well the results of all the numerical-relativity simulations published so far, both when applied at separations of a few gravitational radii, and when applied at separations of tens of thousands of gravitational radii. These validations make our formula a useful tool in a variety of contexts ranging from gravitational-wave physics to cosmology. As representative examples, we first illustrate how it can be used to decrease the phase error of the effective-one-body waveforms during the ringdown phase. Secondly, we show that, when combined with the recently computed self-force correction to the binding energy of nonspinning black-hole binaries, it provides an estimate of the energy emitted during the merger and ringdown. Finally, we use it to calculate the energy radiated in gravitational waves by massive black-hole binaries as a function of redshift, using different models for the seeds of the black-hole population.
In de Sitter spacetime there exists an absolute minimum for the mass of a spin-2 field set by the Higuchi bound m^2 \geq 2H^2. We generalize this bound to arbitrary spatially flat FRW geometries in the context of the recently proposed ghostfree models of Massive Gravity by performing a Hamiltonian analysis for cosmological perturbations. We find that the bound generically indicates that spatially flat FRW solutions in massive gravity, which exhibit a Vainshtein mechanism in the background as required by consistency with observations, imply that the helicity zero mode is a ghost. In contradistinction to previous works, the tension between the Higuchi bound and the Vainshtein mechanism is equally strong regardless of the equation of state for matter.
Large astronomical objects such as stars or planets, produce approximately spherical shapes due to the large gravitational forces, and if the object is rotating rapidly, it becomes an oblate spheroid. In juxtaposition to this, we conduct a thought experiment regarding the properties of a planet being in the form of a perfect cube. We firstly calculate the gravitational potential and from the equipotentials, we deduce the shape of the lakes that would form on the surface of such an object. We then consider the formation of orbits around such objects both with a static and a rotating cube. A possible practical application of these results is that, because cuboid objects can be easily stacked together, we can calculate the field of more complicated shapes, using the principle of superposition, by simply adding the field from a set of component shapes.
We attempt an estimate for the distribution of the tensor mode fraction $r$ over the landscape of vacua in string theory. The dynamics of eternal inflation and quantum tunneling lead to a kind of democracy on the landscape, providing no bias towards large-field or small-field inflation regardless of the class of measure. The tensor mode fraction then follows the number frequency distributions of inflationary mechanisms of string theory over the landscape. We show that an estimate of the relative number frequencies for small-field vs large-field inflation, while unattainable on the whole landscape, may be within reach as a regional answer for warped Calabi-Yau flux compactifications of type IIB string theory.
Links to: arXiv, form interface, find, astro-ph, recent, 1206, contact, help (Access key information)
The zCOSMOS-bright 10k spectroscopic sample reveals a strong environmental dependence of close kinematic galaxy pair fractions in the redshift range 0.2 < z < 1. The fraction of close pairs is three times higher in the top density quartile than in the lowest one. This environmental variation in pair fractions will translate into merger fractions since merger timescales are shown, based on Millennium simulation catalogs, to be largely independent of environment. While galactic properties of close kinematic pairs (morphologies and star formation rates) may seem to be non-representative of an underlying galaxy population, they can be explained by taking into account well-known effects of environment, and changes caused by interactions. The latter is responsible for an increase of irregular galaxies in pairs by a factor of 50-75%, with a disproportionate increase in the number of irregular-irregular pairs (4-8 times), due to disturbance of about 15% of the disk galaxies in pairs. Another sign of interaction is an observed boost in specific star formation rate (factor 2-4) for the closest pairs. While significant for paired galaxies, this triggered star-formation due to interactions represents only about 5% of the integrated star-formation activity in our volume-limited sample. Although majority of close kinematic pairs are in dense environments, the effects of interactions appear to be strongest in the lower density environments. This may introduce strong biases into observational studies of mergers, especially those based on morphological criteria. Relative excess of post-starburst galaxies observed in paired galaxies (factor \sim2) as well as excess of AGNs (factor of over 2), linked with environmental dependence of the pair fractions could indicate that early phases of interactions and merging are plausible candidates for environmental quenching, observed in the global galaxy populations.
Using Keck laser guide star adaptive optics imaging, we have found that the T9 dwarf WISE J1217+1626 and T8 dwarf WISE J1711+3500 are exceptional binaries, with unusually wide separations (~0.8 arcsec, 8-15 AU), large near-IR flux ratios (~2-3 mags), and small mass ratios (~0.5). Keck/NIRSPEC H-band spectra give a spectral type of Y0 for WISE J1217+1626B, and photometric estimates suggest T9.5 for WISE J1711+3500B. The WISE J1217+1626AB system is very similar to the T9+Y0 binary CFBDSIR J1458+1013AB; these two systems are the coldest known substellar multiples, having secondary components of ~400 K and being planetary-mass binaries if their ages are <~1 Gyr. Both WISE J1217+1626B and CFBDSIR J1458+1013B have strikingly blue Y-J colors compared to previously known T dwarfs, including their T9 primaries. Combining all available data, we find that Y-J color drops precipitously between the very latest T dwarfs and the Y dwarfs. The fact that this is seen in (coeval, mono-metallicity) binaries demonstrates that the color drop arises from a change in temperature, not surface gravity or metallicity variations among the field population. Thus, the T/Y transition established by near-IR spectra coincides with a significant change in the ~1 micron fluxes of ultracool photospheres. One explanation is the depletion of potassium, whose broad absorption wings dominate the far-red optical spectra of T dwarfs. This large color change suggests that far-red data may be valuable for classifying objects of <~500 K.
As part of the Panchromatic Hubble Andromeda Treasury (PHAT) multi-cycle program, we observed a 12' \times 6.5' area of the bulge of M31 with the WFC3/UVIS filters F275W and F336W. From these data we have assembled a sample of \sim4000 UV-bright, old stars, vastly larger than previously available. We use updated Padova stellar evolutionary tracks to classify these hot stars into three classes: Post-AGB stars (P-AGB), Post-Early AGB (PE-AGB) stars and AGB-manqu\'e stars. P-AGB stars are the end result of the asymptotic giant branch (AGB) phase and are expected in a wide range of stellar populations, whereas PE-AGB and AGB-manqu\'e (together referred to as the hot post-horizontal branch; HP-HB) stars are the result of insufficient envelope masses to allow a full AGB phase, and are expected to be particularly prominent at high helium or {\alpha} abundances when the mass loss on the RGB is high. Our data support previous claims that most UV-bright sources in the bulge are likely hot (extreme) horizontal branch stars (EHB) and their progeny. We construct the first radial profiles of these stellar populations, and show that they are highly centrally concentrated, even more so than the integrated UV or optical light. However, we find that this UV-bright population does not dominate the total UV luminosity at any radius, as we are detecting only the progeny of the EHB stars that are the likely source of the UVX. We calculate that only a few percent of MS stars in the central bulge can have gone through the HP-HB phase and that this percentage decreases strongly with distance from the center. We also find that the surface density of hot UV-bright stars has the same radial variation as that of low-mass X-ray binaries. We discuss age, metallicity, and abundance variations as possible explanations for the observed radial variation in the UV-bright population.
We report new, sensitive observations of two z ~ 3-3.5 FIR-luminous radio galaxies, 6C1909+72 and B3J2330+3927, in 12CO J=1-0 with the Karl Jansky VLA and at 100-500um using Herschel, alongside new and archival 12CO J=4-3 observations from IRAM PdBI. We introduce a new colour-colour diagnostic plot to constrain the redshifts of several distant, dusty galaxies in our target fields. A bright SMG near 6C1909+72 likely shares the same node or filament as the signpost AGN, but it is not detected in CO despite ~20,000 km/s of velocity coverage. Also in the 6C1909+72 field, a large, red dust feature spanning ~500 kpc is aligned with the radio jet. We suggest several processes by which metal-rich material may have been transported, favouring a collimated outflow reminiscent of the jet-oriented metal enrichment seen in local cluster environments. Our interferometric imaging reveals a gas-rich companion to B3J2330+3927; indeed, all bar one of the eight z >~ 2 radio galaxies (or companions) detected in CO provide some evidence that starburst activity in radio-loud AGN at high redshift is driven by the interaction of two or more gas-rich systems in which a significant mass of stars has already formed, rather than via steady accretion of cold gas from the cosmic web. We find that the CO Tb ratios in radio-loud AGN host galaxies are significantly higher than those seen in similarly intense starbursts where AGN activity is less pronounced. Our most extreme example, where L'(CO4-3)/L'(CO1-0) > 2.7, provides evidence that significant energy is being deposited rapidly into the molecular gas via X-rays and/or mechanical (`quasar-mode') feedback from the AGN, leading to a high degree of turbulence globally and a low optical depth in 12CO - feedback that may lead to the cessation of star formation on a timescale commensurate with that of the jet activity, <~10 Myr.
We present spectroscopic observations in the rest-frame optical and near- to mid-infrared wavelengths of four gravitationally lensed infrared (IR) luminous star-forming galaxies at redshift 1 < z < 3 from the LUCIFER instrument on the Large Binocular Telescope and the Infrared Spectrograph on Spitzer. The sample was selected to represent pure, actively star-forming systems, absent of active galactic nuclei. The large lensing magnifications result in high signal-to-noise spectra that can probe faint IR recombination lines, including Pa-alpha and Br-alpha at high redshifts. The sample was augmented by three lensed galaxies with similar suites of unpublished data and observations from the literature, resulting in the final sample of seven galaxies. We use the IR recombination lines in conjunction with H-alpha observations to probe the extinction, Av, of these systems, as well as testing star formation rate (SFR) indicators against the SFR measured by fitting spectral energy distributions to far-IR photometry. Our galaxies occupy a range of Av from ~0 to 5.9 mag, larger than previously known for a similar range of IR luminosities at these redshifts. Thus, estimates of SFR even at z ~ 2 must take careful count of extinction in the most IR luminous galaxies. We also measure extinction by comparing SFR estimates from optical emission lines with those from far-IR measurements. The comparison of results from these two independent methods indicates a large variety of dust distribution scenarios at 1 < z < 3. Without correcting for dust extinction, the H-alpha SFR indicator underestimates the SFR; the size of the necessary correction depends on the IR luminosity and dust distribution scenario. Individual SFR estimates based on the 6.2 micron PAH emission line luminosity do not show a systematic discrepancy with extinction, although a considerable, ~0.2 dex scatter is observed.
Breakthrough direct detections of planetary companions orbiting A-type stars confirm the existence of massive planets at relatively large separations, but dedicated surveys are required to estimate the frequency of similar planetary systems. To measure the first estimation of the giant exoplanetary systems frequency at large orbital separation around A-stars, we have conducted a deep-imaging survey of young (8-400 Myr), nearby (19-84 pc) A- and F-stars to search for substellar companions in the 10-300 AU range. The sample of 42 stars combines all A-stars observed in previous AO planet search surveys reported in the literature with new AO observations from VLT/NaCo and Gemini/NIRI. It represents an initial subset of the International Deep Planet Survey (IDPS) sample of stars covering M- to B-stars. The data were obtained with diffraction-limited observations in H- and Ks-band combined with angular differential imaging to suppress the speckle noise of the central stars, resulting in typical 5-sigma detection limits in magnitude difference of 12 mag at 1", 14 mag at 2" and 16 mag at 5" which is sufficient to detect massive planets. A detailed statistical analysis of the survey results is performed using Monte Carlo simulations. Considering the planet detections, we estimate the fraction of A-stars having at least one massive planet (3-14 MJup) in the range 5-320 AU to be inside 5.9-18.8% at 68% confidence, assuming a flat distribution for the mass of the planets. By comparison, the brown dwarf (15-75 MJup) frequency for the sample is 2.0-8.9% at 68% confidence in the range 5-320 AU. Assuming power law distributions for the mass and semimajor axis of the planet population, the AO data are consistent with a declining number of massive planets with increasing orbital radius which is distinct from the rising slope inferred from radial velocity (RV) surveys around evolved A-stars.
We introduce a method for constructing end-to-end mock galaxy catalogues using a semi-analytical model of galaxy formation, applied to the halo merger trees extracted from a cosmological N-body simulation. The mocks that we construct are lightcone catalogues, in which a galaxy is placed according to the epoch at which it first enters the past lightcone of the observer, and incorporate the evolution of galaxy properties with cosmic time. We determine the position between the snapshot outputs at which a galaxy enters the observer's lightcone by interpolation. As an application, we consider the effectiveness of the BzK colour selection technique, which was designed to isolate galaxies in the redshift interval 1.4<z<2.5. The mock catalogue is in reasonable agreement with the observed number counts of all BzK galaxies, as well as with the observed counts of the subsample of BzKs that are star-forming galaxies. We predict that over 75 per cent of the model galaxies with K_{AB}<=23, and 1.4<z<2.5, are selected by the BzK technique. Interloper galaxies, outside the intended redshift range, are predicted to dominate bright samples of BzK galaxies (i.e. with K_{AB}<=21). Fainter K-band cuts are necessary to reduce the predicted interloper fraction. We also show that shallow B-band photometry can lead to confusion in classifying BzK galaxies as being star-forming or passively evolving. Overall, we conclude that the BzK colour selection technique is capable of providing a sample of galaxies that is representative of the 1.4<z<2.5 galaxy population.
We present the Pan-STARRS1 discovery and light curves, and follow-up MMT and Gemini spectroscopy of an ultra-luminous supernova (ULSN; dubbed PS1-11bam) at a redshift of z=1.566 with a peak brightness of M_UV=-22.3 mag. PS1-11bam is one of the highest redshift spectroscopically-confirmed SNe known to date. The spectrum is characterized by broad absorption features typical of previous ULSNe (e.g., CII, SiIII), and by strong and narrow MgII and FeII absorption lines from the interstellar medium (ISM) of the host galaxy, confirmed by an [OII]3727 emission line at the same redshift. The equivalent widths of the FeII2600 and MgII2803 lines are in the top quartile of the quasar intervening absorption system distribution, but are weaker than those of gamma-ray burst intrinsic absorbers (i.e., GRB host galaxies). We also detect the host galaxy in pre-explosion Pan-STARRS1 data and find that its UV spectral energy distribution is best fit with a young stellar population age of tau~15-45 Myr and a stellar mass of M \sim (1.1-2.6)x10^9 M_sun (for Z=0.05-1 Z_sun). The star formation rate inferred from the UV continuum and [OII]3727 emission line is ~10 M_sun/yr, higher than in any previous ULSN host. PS1-11bam provides the first direct demonstration that ULSNe can serve as probes of the interstellar medium in distant galaxies. At the present, the depth and red sensitivity of PS1 are uniquely suited to finding such events at cosmologically interesting redshifts (z~1-2); the future combination of LSST and 30-m class telescopes promises to extend this technique to z~4.
(Abridged) The 30 Dor region in the Large Magellanic Cloud (LMC) is the most vigorous star-forming region in the Local Group. Star formation in this region is taking place in low-metallicity molecular gas which is exposed to an extreme far--ultraviolet (FUV) radiation field powered by the massive compact star cluster R136. We used the NANTEN2 telescope to obtain high-angular resolution observations of the 12CO 4-3, 7-6, and 13CO 4-3 rotational lines and [CI] 3P1-3P0 and 3P2-3P1 fine-structure sub-millimeter transitions in 30Dor-10, the brightest CO and FIR-emitting cloud at the center of the 30Dor region. We derive the properties of the low-metallicity molecular gas using an excitation/radiative transfer code and find a self-consistent solution of the chemistry and thermal balance of the gas in the framework of a clumpy cloud PDR model. We compare the derived properties with those in the N159W region, which is exposed to a more moderate far-ultraviolet radiation field compared with 30Dor-10, but has similar metallicity. We also combine our CO detections with previously observed low-J CO transitions to derive the CO spectral-line energy distribution in 30Dor-10 and N159W. The separate excitation analysis of the submm CO lines and the neutral carbon fine structure lines shows that the mid-J CO and [CI]-emitting gas in the 30Dor-10 region has a temperature of about 160 K and a H2 density of about 10^4 cm^-3. We find that the molecular gas in 30Dor-10 is warmer and has a smaller beam filling factor compared to that of N159W, which might be a result of the effect of a strong FUV radiation field heating and disrupting the low--metallicity molecular gas. We use a clumpy PDR model (including the [CII] line intensity reported in the literature) to constrain the FUV intensity to about chi_0 ~ 3100 and an average total H density of the clump ensemble of about 10^5 cm^-3 in 30Dor-10.
(abridged) We present a clustering analysis of 370 high-confidence H-alpha emitters (HAEs) at z=2.23. The HAEs are detected in the Hi-Z Emission Line Survey (HiZELS), a large-area blank field 2.121um narrowband survey using the United Kingdom Infrared Telescope (UKIRT) Wide Field Camera (WFCAM). Averaging the two-point correlation function of HAEs in two ~1 degree scale fields (United Kingdom Infrared Deep Sky Survey/Ultra Deep Survey [UDS] and Cosmological Evolution Survey [COSMOS] fields) we find a clustering amplitude equivalent to a correlation length of r_0=3.7+/-0.3 Mpc/h for galaxies with star formation rates of >7 M_sun/yr. The data are well-fitted by the expected correlation function of Cold Dark Matter, scaled by a bias factor: omega_HAE=b^2 omega_DM where b=2.4^{+0.1}_{-0.2}. The corresponding 'characteristic' mass for the halos hosting HAEs is log(M_h/[M_sun/h])=11.7+/-0.1. Comparing to the latest semi-analytic GALFORM predictions for the evolution of HAEs in a LCDM cosmology, we find broad agreement with the observations, with GALFORM predicting a HAE correlation length of ~4 Mpc/h. Motivated by this agreement, we exploit the simulations to construct a parametric model of the halo occupation distribution (HOD) of HAEs, and use this to fit the observed clustering. Our best-fitting HOD can adequately reproduce the observed angular clustering of HAEs, yielding an effective halo mass and bias in agreement with that derived from the scaled omega_DM fit, but with the relatively small sample size the current data provide a poor constraint on the HOD. Our results support the broad picture that 'typical' (~L*) star-forming galaxies have been hosted by dark matter haloes with M_h<10^12 M_sun/h since z~2, but with a broad occupation distribution and clustering that is likely to be a strong function of luminosity.
Synchrotron radiation mechanism, when electrons are accelerated in a relativistic shock, is known to have serious problems to explain the observed gamma-ray spectrum below the peak for most Gamma-Ray Bursts (GRBs); the synchrotron spectrum below the peak is much softer than observed spectra. Recently, the possibility that electrons responsible for the radiation cool via Inverse Compton, but in the Klein-Nishina regime, has been proposed as a solution to this problem. We provide an analytical study of this effect and show that it leads to a hardening of the low energy spectrum but not by enough to make it consistent with the observed spectra for most GRBs (this is assuming that electrons are injected continuously over a time scale comparable to the dynamical time scale, as is expected for internal shocks of GRBs). In particular, we find that it is not possible to obtain a spectrum with \alpha>-0.1 (f_{\nu} \propto \nu^{\alpha}) whereas the typical observed value is \alpha\sim0. Moreover, extreme values for a number of parameters are required in order that \alpha\sim-0.1: the energy fraction in magnetic field needs to be less than about 10^{-4}, the thermal Lorentz factor of electrons should be larger than 10^6, and the radius where gamma-rays are produced should be not too far away from the deceleration radius. These difficulties suggest that the synchrotron radiation mechanism in internal shocks does not provide a self-consistent solution when \alpha>-0.2.
We present a detailed analysis (including redshift tomography) of the cosmic dipoles in the Keck+VLT quasar absorber and in the Union2 SnIa samples. We show that the fine structure constant cosmic dipole obtained through the Keck+VLT quasar absorber sample at $4.1\sigma$ level, is anomalously aligned with the corresponding dark energy dipole obtained through the Union2 sample at $2\sigma$ level. The angular separation between the two dipole directions is $11.3^\circ \pm 11.8^\circ$. We use Monte Carlo simulations to find the probability of obtaining the observed dipole magnitudes with the observed alignment, in the context of an isotropic cosmological model with no correlation between dark energy and fine structure constant. We find that this probability is less than one part in $10^6$. We propose a simple physical model (extended topological quintessence) which naturally predicts spherical inhomogeneous distribution for both dark energy density and fine structure constant values. The model is based on the existence of a recently formed giant global monopole with Hubble scale core which also couples non-minimally to electromagnetism. Aligned dipole anisotropies would naturally emerge for an off-center observer for both the fine structure constant and for dark energy density. This model smoothly reduces to \lcdm for proper limits of its parameters.
We obtain Keck HIRES spectroscopy of HVS5, one of the fastest unbound stars in the Milky Way halo. We show that HVS5 is a 3.62 +- 0.11 Msun main sequence B star at a distance of 50 +- 5 kpc. The difference between its age and its flight time from the Galactic center is 105 +-18(stat)+-30(sys) Myr; flight times from locations elsewhere in the Galactic disk are similar. This 10^8 yr `arrival time' between formation and ejection is difficult to reconcile with any ejection scenario involving massive stars that live for only 10^7 yr. For comparison, we derive arrival times of 10^7 yr for two unbound runaway B stars, consistent with their disk origin where ejection results from a supernova in a binary system or dynamical interactions between massive stars in a dense star cluster. For HVS5, ejection during the first 10^7 yr of its lifetime is ruled out at the 3-sigma level. Together with the 10^8 yr arrival times inferred for three other well-studied hypervelocity stars (HVSs), these results are consistent with a Galactic center origin for the HVSs. If the HVSs were indeed ejected by the central black hole, then the Galactic center was forming stars ~200 Myr ago, and the progenitors of the HVSs took ~100 Myr to enter the black hole's loss cone.
The evolution of multiple stellar systems can be driven by Kozai cycles and tidal friction (KCTF), which shrink the orbit of the inner binary. There is an interesting possibility that two close binaries on a common long-period orbit experience mutually-induced KCTF. We present the discovery of a possible new quadruple system composed of two unresolved eclipsing binaries (EBs), CzeV343 (V~13.5 mag). We obtained photometric observations of CzeV343 that completely cover the two orbital periods and we successfully model the light curves as the sum of two detached EBs. We provide confidence intervals for the model parameters and minima timings by bootstrap resampling of our data. One of the EBs shows a distinctly eccentric orbit with a total eccentricity of about 0.18. The two orbital periods, 1.20937 and 0.80693 days, are within 0.1% of a 3:2 ratio. We speculate that this might be the result of KCTF-driven evolution of a quadruple system and we discuss this hypothesis in the context of other quadruple systems composed of two EBs. We make our double EB fitting code publicly available to provide a tool for long-term monitoring of the mutual orbit in such systems.
Accretion onto a supermassive black hole of a rotating inflow is a particularly difficult problem to study because of the wide range of length scales involved. There has been some analytic and numerical treatment of the global properties of accretion flows, but detailed numerical simulations are required to address certain critical aspects. We use the ZEUS code to run hydrodynamical simulations of rotating, axisymmetric accretion flows with Bremsstrahlung cooling, considering solutions with and without viscous angular momentum transport, and also electron thermal conduction. Infalling gas is followed from well beyond R_Bondi down to the vicinity of the black hole. Absent viscous transport, when the centrifugal balance radius significantly exceeds R_Schwarzschild, the accretion rate is zero and the flow approaches a stationary solution in which pressure impedes inflow from large radii. With viscosity, we find two general classes of solutions: low inflow rate, hot, vertically extended disks with very low accretion and disk and conical wind outflows near R_Bondi, and strong inflow solutions which have cold, geometrically thin disks accreting at close to Mdot_Edd. We produce a continuum of solutions with respect to the Eddington ratio Mdot_Bondi/Mdot_Edd, and there is a sharp transition between the two general classes of solutions at Eddington ratio ~ few x 10^(-2). The low accretion inflow-outflow solutions are of two types. Equatorial outflow dominates when viscosity is larger than thermal conductivity, but polar outflow can be significant when the Prandtl number ~ 0.05. Our simulations have converged with respect to spatial resolution and temporal duration, and they do not depend strongly on our choice of boundary conditions. We also note the possibility that radiative feedback loops can cause the flow to switch between the hot and cold disk states, with potential applications to quasars.
We present the Spitzer Extragalactic Representative Volume Survey (SERVS), an 18 square degrees medium-deep survey at 3.6 and 4.5 microns with the post-cryogenic Spitzer Space Telescope to ~2 microJy (AB=23.1) depth of five highly observed astronomical fields (ELAIS-N1, ELAIS-S1, Lockman Hole, Chandra Deep Field South and XMM-LSS). SERVS is designed to enable the study of galaxy evolution as a function of environment from z~5 to the present day, and is the first extragalactic survey both large enough and deep enough to put rare objects such as luminous quasars and galaxy clusters at z>1 into their cosmological context. SERVS is designed to overlap with several key surveys at optical, near- through far-infrared, submillimeter and radio wavelengths to provide an unprecedented view of the formation and evolution of massive galaxies. In this paper, we discuss the SERVS survey design, the data processing flow from image reduction and mosaicing to catalogs, as well as coverage of ancillary data from other surveys in the SERVS fields. We also highlight a variety of early science results from the survey.
We explore the nature of heavily obscured quasar host galaxies at z~2 using deep Hubble Space Telescope WFC3/IR imaging of 28 Dust Obscured Galaxies (DOGs) to investigate the role of major mergers in driving black hole growth. The high levels of obscuration of the quasars selected for this study act as a natural coronagraph, blocking the quasar light and allowing a clear view of the underlying host galaxy. The sample of heavily obscured quasars represents a significant fraction of the cosmic mass accretion on supermassive black holes as the quasars have inferred bolometric luminosities around the break of the quasar luminosity function. We find that only a small fraction (4%, at most 11-25%) of the quasar host galaxies are major mergers. Fits to their surface brightness profiles indicate that 90% of the host galaxies are either disk dominated, or have a significant disk. This disk-like host morphology, and the corresponding weakness of bulges, is evidence against major mergers and suggests that secular processes are the predominant driver of massive black hole growth. Finally, we suggest that the co-incidence of mergers and AGN activity is luminosity dependent, with only the most luminous quasars being triggered mostly by major mergers.
It was recently shown that the power spectrum in redshift space can be written as a sum of cross-power spectra between number weighted velocity moments. We investigate the properties of these power spectra for simulated galaxies and dark matter halos and compare them to the dark matter power spectra, generalizing the concept of the bias. Because all of the quantities are number weighted this approach is well defined even for sparse systems such as massive halos, in contrasts to the previous approaches to RSD where velocity correlations have been explored. We find that the number density weighting leads to a strong scale dependence of the bias terms for momentum density auto-correlation and cross-correlation with density. This trend becomes more significant for the more biased halos and leads to an enhancement of RSD power relative to the linear theory. Fingers-of-god effects, which in this formalism come from the correlations of the higher order moments beyond the momentum density, lead to smoothing of the power spectrum and can reduce this enhancement of power, but are relatively small for halos with no small-scale velocity dispersion. In comparison, for a more realistic galaxy sample with satellites the velocity dispersion generated by satellite motions inside the halos leads to a larger power suppression on small scales, but this depends on the satellite fraction. We investigate several statistics such as the two-dimensional power spectrum, its multipole moments, its powers of mu^2, and configuration space statistics. Overall we find that the nonlinear effects in realistic galaxy samples such as luminous red galaxies affect the redshift space clustering on very large scales: for example, the quadrupole moment is affected by 10% for k<0.1h/mpc, which means that these effects need to be understood if we want to extract cosmological information from the redshift space distortions.
The Virtual Observatory (VO) is realizing global electronic integration of astronomy data. One of the long-term goals of the U.S. VO project, the Virtual Astronomical Observatory (VAO), is development of services and protocols that respond to the growing size and complexity of astronomy data sets. This paper describes how VAO staff are active in such development efforts, especially in innovative strategies and techniques that recognize the limited operating budgets likely available to astronomers even as demand increases. The project has a program of professional outreach whereby new services and protocols are evaluated.
The U.S. Virtual Astronomical Observatory (VAO; this http URL) has been in operation since May 2010. Its goal is to enable new science through efficient integration of distributed multi-wavelength data. This paper describes the management and organization of the VAO, and emphasizes the techniques used to ensure efficiency in a distributed organization. Management methods include using an annual program plan as the basis for establishing contracts with member organizations, regular communication, and monitoring of processes.
We present observations of the BL Lac object 1ES 0414+009 in the >200 GeV gamma-ray band by the VERITAS array of Cherenkov telescopes. 1ES 0414+009 was observed by VERITAS between January 2008 and February 2011, resulting in 56.2 hours of good quality pointed observations. These observations resulted in a detection of 822 events from the source corresponding to a statistical significance of 6.4 standard deviations (6.4 sigma) above the background. The source flux, showing no evidence for variability, is measured as 5.2 +/- 1.1_stat +/- 2.6_sys * 10^-12 photons cm^-2 s^-1 above 200 GeV, equivalent to approximately 2% of the Crab Nebula flux above this energy. The differential photon spectrum from 230 GeV to 850 GeV is well fit by a power law with an photon index of Gamma 3.4 +/- 0.5_stat +/- 0.3_sys and a flux normalization of 1.6 +/- 0.3_stat +/- 0.8_sys * 10^-11 photons cm^-2 s^-1 at 300 GeV. We also present multiwavelength results taken in the optical (MDM), X-ray (Swift-XRT), and GeV (Fermi-LAT) bands and use these results to construct a broadband spectral energy distribution (SED). Modeling of this SED indicates that homogenous one-zone leptonic scenarios are not adequate to describe emission from the system, with a lepto-hadronic model providing a better fit to the data.
Phase-induced amplitude apodization (PIAA) coronagraphs are a promising technology for imaging exoplanets, with the potential to detect Earth-like planets around Sun-like stars. A PIAA system nominally consists of a pair of mirrors which reshape incident light without attenuation, coupled with one or more apodizers to mitigate diffraction effects or provide additional beam-shaping to produce a desired output profile. We present a set of equations that allow apodizers to be chosen for any given pair of mirrors, or conversely mirror shapes chosen for given apodizers, to produce an arbitrary amplitude profile at the output of the system. We show how classical PIAA systems may be designed by this method, and present the design of a novel 4-mirror system with higher throughput than a standard 2-mirror system. We also discuss the limitations due to diffraction and the design steps that may be taken to mitigate them.
The review contains an analysis of the observed and model curves of the interstellar extinction and polarization. The observations mainly give information on dust in diffuse and translucent interstellar clouds. The features of various dust grain models including spherical/non-spherical, homogeneous/inhomogeneous particles are discussed. A special attention is devoted to the analysis of the grain size distributions, alignment mechanisms and magnetic field structure in interstellar clouds. It is concluded that the interpretation of interstellar extinction and polarization is not yet complete.
We present 2" - 10" imaging of eleven transitions from nine molecular species across the nuclear bar in Maffei 2. The data were obtained with the BIMA and OVRO interferometers. The ten detected transitions are compared with existing CO isotopologues, HCN, CS and millimeter continuum data. Dramatic spatial variations among the mapped species are observed across the nuclear bar. A principle component analysis is performed to characterize correlations between the transitions, star formation and molecular column density. The analysis reveals that HCN, HNC, HCO+ and 3 mm continuum are tightly correlated, indicating a direct connection to massive star formation. We find two main morphologically distinct chemical groups, CH3OH, SiO and HNCO comprising the grain chemistry molecules, versus HCN, HNC, HCO+ and C2H, molecules strong in the presence of star formation. The grain chemistry molecules, HNCO, CH3OH and SiO, trace hydrodynamical bar shocks. The near constancy of the HNCO/CH3OH, SiO/CH3OH and SiO/HNCO ratios argue that shock properties are uniform across the nucleus. HCN/HCO+, HCN/HNC, HCN/CS and HCN/CO ratios are explained primarily by variations in density. High HCO+/N2H+ ratios are correlated with the C2H line, suggesting that this ratio may be a powerful new dense photon-dominated region (PDR) probe in external galaxies. C2H reveals a molecular outflow along the minor axis. The morphology and kinematics of the outflow are consistent with an outflow age of 6-7 Myrs.
We present a detailed study of the Circinus Galaxy, investigating its star formation, dust and gas properties both in the inner and outer disk. To achieve this, we obtained high-resolution Spitzer mid-infrared images with the IRAC (3.6, 5.8, 4.5, 8.0 micron) and MIPS (24 and 70 micron) instruments and sensitive HI data from the Australia Telescope Compact Array (ATCA) and the 64-m Parkes telescope. These were supplemented by CO maps from the Swedish-ESO Submillimetre Telescope (SEST). Because Circinus is hidden behind the Galactic Plane, we demonstrate the careful removal of foreground stars as well as large- and small-scale Galactic emission from the Spitzer images. We derive a visual extinction of Av = 2.1 mag from the Spectral Energy Distribution of the Circinus Galaxy and total stellar and gas masses of 9.5 x 10^{10} Msun and 9 x 10^9 Msun, respectively. Using various wavelength calibrations, we find obscured global star formation rates between 3 and 8 Msun yr^{-1}. Star forming regions in the inner spiral arms of Circinus, which are rich in HI, are beautifully unveiled in the Spitzer 8 micron image. The latter is dominated by polycyclic aromatic hydrocarbon (PAH) emission from heated interstellar dust. We find a good correlation between the 8 micron emission in the arms and regions of dense HI gas. The (PAH 8 micron) / 24 micron surface brightness ratio shows significant variations across the disk of Circinus.
We present a detailed analysis of high resolution H I observations of the Magellanic spiral galaxies NGC 4618 and NGC 4625. While the H I disk of NGC 4625 is remarkably quiescent with a nearly uniform velocity dispersion and no evidence of H I holes, there is a dynamic interplay between star formation and the distribution of neutral hydrogen in NGC 4618. We calculate the critical density for widespread star formation in each galaxy and find that star formation proceeds even where the surface density of the atomic gas is well below the critical density necessary for global star formation. There are strong spatial correlations in NGC 4618 between UV emission, 1.4 GHz radio continuum emission, and peaks in the H I column density. Despite the apparent overlap of the outer disks of the two galaxies, we find that they are kinematically distinct, indicating that NGC 4618 and NGC 4625 are not interacting. The structure of NGC 4618 and, in particular, the nature of its outer ring, are highly suggestive of an interaction, but the timing and nature of such an interaction remain unclear.
We perform a systematic study of hybrid star configurations using several parametrizations of a relativistic mean-field hadronic EoS and the NJL model for three-flavor quark matter. For the hadronic phase we use the stiff GM1 and TM1 parametrizations, as well as the very stiff NL3 model. In the NJL Lagrangian we include scalar, vector and 't Hooft interactions. The vector coupling constant $g_v$ is treated as a free parameter. We also consider that there is a split between the deconfinement and the chiral phase transitions which is controlled by changing the conventional value of the vacuum pressure $- \Omega_0$ in the NJL thermodynamic potential by $- (\Omega_0 + \delta \Omega_0)$, being $\delta \Omega_0$ a free parameter. We find that, as we increase the value of $\delta \Omega_0$, hybrid stars have a larger maximum mass but are less stable, i.e. hybrid configurations are stable within a smaller range of central densities. For large enough $\delta \Omega_0$, stable hybrid configurations are not possible at all. The effect of increasing the coupling constant $g_v$ is very similar. We show that stable hybrid configurations with a maximum mass larger than the observed mass of the pulsar PSR J1614-2230 are possible for a large region of the parameter space of $g_v$ and $\delta \Omega_0$ provided the hadronic equation of state contains nucleons only. When the baryon octet is included in the hadronic phase, only a very small region of the parameter space allows to explain the mass of PSR J1614-2230. We compare our results with previous calculations of hybrid stars within the NJL model. We show that it is possible to obtain stable hybrid configurations also in the case $\delta \Omega_0=0$ that corresponds to the conventional NJL model for which the pressure and density vanish at zero temperature and chemical potential.
If a component of cosmological dark matter is made up of massive particles - such as sterile neutrinos - that decay with cosmological lifetime to emit photons, the reionization history of the universe would be affected, and cosmic microwave background anisotropies can be used to constrain such a decaying particle model of dark matter. The optical depth depends rather sensitively on the decaying dark matter particle mass m_{dm}, lifetime tau_{dm}, and the mass fraction of cold dark matter f that they account for in this model. Assuming that there are no other sources of reionization and using the WMAP 7-year data, we find that 250 eV < m_{dm} < 1 MeV, whereas 2.23*10^3 yr < tau_{dm} < 1.23*10^18 yr. The best fit values for m_{dm} and tau_{dm}/f are 17.3 keV and 2.03*10^16 yr respectively.
We update our prior work on the case B collisional-recombination spectrum of He I to incorporate \textit{ab initio} photoionisation cross-sections. This large set of accurate, self-consistent cross-sections represents a significant improvement in He I emissivity calculations because it largely obviates the piecemeal nature that has marked all modern works. A second, more recent set of \textit{ab initio} cross-sections is also available, but we show that those are less consistent with bound-bound transition probabilities than our adopted set. We compare our new effective recombination coefficients with our prior work and our new emissivities with those by other researchers, and we conclude with brief remarks on the effects of the present work on the He I error budget. Our calculations cover temperatures $5000 \le T_e \le 25000$ K and densities $10^1 \le n_e \le 10^{14}$ cm$^{-3}$. Full results are available online.
A new paradigm for active galactic jet kinematics has emerged through detailed investigations of BL Lac objects using very long baseline radio interferometry. In this new scheme, most, if not all, jet components appear to remain stationary with respect to the core but show significant non-radial motions. This paper presents results from our kinematic investigation of the jets of a statistically complete sample of radio-loud flat-spectrum active galaxies, focusing on the comparison between the jet kinematic properties of BL Lacs and flat-spectrum radio-quasars. It is shown that there is a statistically significant difference between the kinematics of the two AGN classes, with BL Lacs showing more bent jets, that are wider and show slower movement along the jet axis, compared to flat-spectrum radio-quasars. This is interpreted as evidence for helically structured jets.
We discuss the possible source of a highly-dispersed radio transient discovered in the Parkes Multi-beam Pulsar Survey (PMPS). The pulse has a dispersion measure of $746\mathrm{cm}^{-3}\mathrm{pc}$, a peak flux density of 400 mJy for the observed pulse width of 7.8 ms, and a flat spectrum across a 288-MHz band centred on 1374 MHz. The flat spectrum suggests that the pulse did not originate from a pulsar, but is consistent with radio-emitting magnetar spectra. The non-detection of subsequent bursts constrains any possible pulsar period to $\gtrsim1$ s, and the pulse energy distribution to being much flatter than typical giant pulse emitting pulsars. The burst is also consistent with the radio signal theorised from an annihilating mini black hole. Extrapolating the PMPS detection rate, provides a limit of $\Omega_{BH}\lesssim5\times10^{-14}$ on the density of these objects. We investigate the consistency of these two scenarios, plus several other possible solutions, as potential explanations to the origin of the pulse, as well as for another transient with similar properties: the Lorimer Burst.
Rotational mixing in massive main sequence stars is predicted to monotonically increase their surface nitrogen abundance with time. We use this effect to design a method for constraining the age and the inclination angle of massive main sequence stars, given their observed luminosity, effective temperature, projected rotational velocity and surface nitrogen abundance. This method relies on stellar evolution models for different metallicities, masses and rotation rates. We use the population synthesis code STARMAKER to show the range of applicability of our method. We apply this method to 79 early B-type main sequence stars near the LMC clusters NGC 2004 and N 11 and the SMC clusters NGC 330 and NGC 346. From all stars within the sample, 17 were found to be suitable for an age analysis. For ten of them, which are rapidly rotating stars without a strong nitrogen enhancement, it has been previously concluded that they did not evolve as rotationally mixed single stars. This is confirmed by our analysis, which flags the age of these objects as highly discrepant with their isochrone ages. For the other seven stars, their nitrogen and isochrone ages are found to agree within error bars, what validates our method. Constraints on the inclination angle have been derived for the other 62 stars,with the implication that the nitrogen abundances of the SMC stars, for which mostly only upper limits are known, fall on average significantly below those limits. Nitrogen chronology is found to be a new useful tool for testing stellar evolution and to constrain fundamental properties of massive main sequence stars. A web version of this tool is provided.
We show how the redshift and peak-flux distributions of gamma-ray bursts (GRBs) have an observation time dependence that can be used to discriminate between different burst populations. We demonstrate how observation time relations can be derived from the standard integral distributions and that they can differentiate between GRB populations detected by both the BATSE and \emph{Swift} satellites. Using \emph{Swift} data we show that a redshift--observation-time relation (log\,$Z$\,--\,log\,$T$) is consistent with both a peak-flux\,--\,observation time relation (log\,$P$\,--\,log\,$T$) and a standard log\,$N$\,--\,log\,$P$ brightness distribution. As the method depends only on rarer small-$z$ events, it is invariant to high-$z$ selection effects. We use the log\,$Z$\,--\,log\,$T$ relation to show that sub-luminous GRBs are a distinct population occurring at a higher rate of order $150^{+180}_{-90} \mathrm{Gpc}^{-3}\mathrm{yr}^{-1}$. Our analysis suggests that GRB 060505 -- a relatively nearby GRB observed without any associated supernova -- is consistent with a sub-luminous population of bursts. Finally, we suggest that our relations can be used as a consistency test for some of the proposed GRB spectral energy correlations.
Numerical simulations of hot accretion flow have shown that the mass accretion rate decreases with decreasing radius; consequently the density profile of accretion flow becomes flatter compared to the case of a constant accretion rate. This result has important theoretical and observational implications. However, because of technical difficulties, the radial dynamic range in almost all previous simulations usually spans at most two orders of magnitude. This small dynamical range, combined with the effects of boundary conditions, makes the simulation results suspectable. Especially, the radial profiles of density and accretion rate may not be precise enough to be used to compare with observations. In this paper we present a "two-zone" approach to expand the radial dynamical range from two to four orders of magnitude. We confirm previous results and find that from $r_s$ to $ 10^4r_s$ the radial profiles of accretion rate and density can be well described by $\dot{M}(r)\propto r^s$ and $\rho\propto r^{-p}$. The values of (s, p) are (0.48, 0.65) and (0.4, 0.85), for viscous parameter $\alpha=0.001$ and 0.01, respectively. We have looked up numerical simulation works in the literature and found that the values of $s$ and $p$ are all similar, no matter a magnetic field is included or not and what kind of initial conditions are adopted. The density profile we obtain is in good quantitative agreement with that obtained from the detailed observations and modeling to Sgr A* and NGC 3115. The origin of such a accretion rate profile will be investigated in a subsequent paper.
Previous hydrodynamical (HD) and magnetohydrodynamical (MHD) numerical simulations of hot accretion flows have shown that the mass accretion rate decreases with decreasing radius. Two models have been proposed to explain this result. In the ADIOS model, the inward decrease of accretion rate is because of the loss of gas in the outflow. In the CDAF model, the gas is assumed to be locked in convective eddies, which results in the inward decrease of the accretion rate. We investigate the nature of inward decrease of accretion rate using HD and MHD simulations. We calculate various properties of inflow and outflow, including the mass flux, radial and rotational velocities, temperature, and the Bernoulli parameter ($Be$). Systematic and significant differences between inflow and outflow are found. These results suggest that the inflow and outflow are not dominated by convective turbulence, but are systematic inward and outward motion. We have also analyzed the convective stability of MHD accretion flow and found that they are convectively stable. These results indicate that the ADIOS scenario is favored. The different properties of inflow and outflow also suggest that the mechanisms of producing outflow in HD and MHD flows are buoyancy and centrifugal force associated with the magnetic field, respectively. The latter mechanism is similar to the Blandford & Payne mechanism. We also study the effect of initial conditions in the simulations. We find that the value of $Be$ is mainly determined by the value of $Be$ of the initial condition. We discuss some possible observational applications of our outflow model. These observations include the Fermi bubble observed in the Galaxy center, and winds widely observed in AGNs and black hole X-ray binaries.
Microphysics of weakly magnetized relativistic collisionless shock waves, corroborated by recent high performance numerical simulations, indicate the presence of a microturbulent layer of large magnetic field strength behind the shock front, which must decay beyond some hundreds of skin depths. The present paper discusses the dynamics of such microturbulence, borrowing from these same numerical simulations, and calculates the synchrotron signature of a powerlaw of shock accelerated particles. The decaying microturbulent layer is found to leave distinct signatures in the spectro-temporal evolution of the spectrum $F_\nu \propto t^{-\alpha}\nu^{-\beta}$ of a decelerating blast wave, which are potentially visible in early multi-wavelength follow-up observations of gamma-ray bursts. This paper also discusses the influence of the evolving microturbulence on the acceleration process, with particular emphasis on the maximal energy of synchrotron afterglow photons, which falls in the GeV range for standard gamma-ray burst parameters. Finally, this paper argues that the evolving microturbulence plays a key role in shaping the spectra of recently observed gamma-ray bursts with extended GeV emission, such as GRB090510.
Primordial Black Holes (PBH's) can form in the early Universe from the collapse of large density fluctuations. Tight observational limits on their abundance constrain the amplitude of the primordial fluctuations on very small scales which can not otherwise be constrained, with PBH's only forming from the extremely rare large fluctuations. The number of PBH's formed is therefore sensitive to small changes in the shape of the tail of the fluctuation distribution, which itself depends on the amount of non-Gaussianity present. We study, for the first time, how quadratic and cubic local non-Gaussianity of arbitrary size (parameterised by f_nl and g_nl respectively) affects the PBH abundance and the resulting constraints on the amplitude of the fluctuations on very small scales. Intriguingly we find that even non-linearity parameters of order unity have a significant impact on the PBH abundance. The sign of the non-Gaussianity is particularly important, with the constraint on the allowed fluctuation amplitude tightening by an order of magnitude as f_nl changes from just -0.5 to 0.5. We find that if PBH's are observed in the future, then regardless of the amplitude of the fluctuations, non-negligible negative f_nl would be ruled out. Finally we show that g_nl can have an even larger effect on the number of PBH's formed than f_nl.
We study the evolution of debris created in the giant impacts expected during the final stages of terrestrial planet formation. The starting point is the debris created in a simulation of the Moon-forming impact. The dynamical evolution is followed for 10 Myr including the effects of Earth, Venus, Mars and Jupiter. The spatial distribution evolves from a clump in the first few months to an asymmetric ring for the first 10 kyr and finally becoming an axisymmetric ring by about 1 Myr after the impact. By 10 Myr after the impact 20% of the particles have been accreted onto Earth and 17% onto Venus, with 8% ejected by Jupiter and other bodies playing minor roles. However, the fate of the debris also depends strongly on how fast it is collisionally depleted, which depends on the poorly constrained size distribution of the impact debris. Assuming that the debris is made up of 30% by mass mm-cm-sized vapour condensates and 70% boulders up to 500 km, we find that the condensates deplete rapidly on ~1000 yr timescales, whereas the boulders deplete predominantly dynamically. By considering the luminosity of dust produced in collisions within the boulder-debris distribution we find that the Moon-forming impact would have been readily detectable around other stars in Spitzer 24 micron surveys for around 25 Myr after the impact, with levels of emission comparable to many known hot dust systems. The vapour condensates meanwhile produce a short-lived, optically thick, spike of emission. We use these surveys to make an estimate of the fraction of stars that form terrestrial planets, F_TPF. Since current terrestrial planet formation models invoke multiple giant impacts, the low fraction of 10-100 Myr stars found to have warm (~150 K) dust implies that F_TPF ~<10%.
A bouncing cosmology with an initial matter-dominated phase of contraction during which scales which are currently probed with cosmological observations exit the Hubble radius provides a mechanism alternative to inflation for producing a nearly scale-invariant spectrum of cosmological perturbations. In this review article I first discuss the evolution of cosmological fluctuations in the matter bounce scenario, and then discuss various attempts at realizing such a scenario. Observational signatures which will allow the matter bounce to be distinguished from the inflationary paradigm are also discussed.
We plan to install an infrared telescope at the new site of Tibet, China. The primary mirror diameter is 50cm, and the focal ratio F8. The Xenics 640\times512 near infrared camera is employed, equipped with a dedicated high speed InGaAs detector array, working up to 1.7{\mu}m. The new site is located on 5100m mountain, near Gar town, Ali, where is an excellent site for both infrared and submillimeter observations. The telescope will be remotely controlled through internet. The goal of IRT is to make site testing, detect variable stars, and search for extrasolar planets.
The large uncertainties associated with measuring the amount of high temperature emission in solar active regions represents a significant impediment to making progress on the coronal heating problem. Most current observations at temperatures of 3 MK and above are taken with broad band soft X-ray instruments. Such measurements have proven difficult to interpret unambiguously. Here we present the first spectroscopic observations of the Fe XVIII 974.86 AA emission line in an on-disk active region taken with then SUMER instrument on SOHO. Fe XVIII has a peak formation temperature of 7.1 MK and provides important constraints on the amount of impulsive heating in the corona. Detailed evaluation of the spectra and comparison of the SUMER data with soft X-ray images from the XRT on Hinode confirm that this line is unblended. We also compare the spectroscopic data with observations from the AIA 94 AA channel on SDO. The AIA 94 AA channel also contains Fe XVIII, but is blended with emission formed at lower temperatures. We find that is possible to remove the contaminating blends and form relatively pure Fe XVIII images that are consistent with the spectroscopic observations from SUMER. The observed spectra also contain the Ca XIV 943.63 AA line that, although a factor 2 to 6 weaker than the Fe XVIII 974.86 AA line, allows us to probe the plasma around 3.5 MK. The observed ratio between the two lines indicates (isothermal approximation) that most of the plasma in the brighter Fe XVIII active region loops is at temperatures between 3.5 and 4 MK.
We compare the observed mass functions and age distributions of star clusters in six well-studied galaxies: the Milky Way, Magellanic Clouds, M83, M51, and Antennae. In combination, these distributions span wide ranges of mass and age: $10^2\lea M/M_{\odot}\lea10^6$ and $10^6\lea\tau/yr \lea10^9$. We confirm that the distributions are well represented by power laws: $dN/dM\propto M^{\beta}$ with $\beta \approx-1.9$ and $dN/d\tau\propto\tau^{\gamma}$ with $\gamma\approx -0.8$. The mass and age distributions are approximately independent of each other, ruling out simple models of mass-dependent disruption. As expected, there are minor differences among the exponents, at a level close to the true uncertainties, $\epsilon_{\beta}\sim\epsilon_{\gamma}\sim$~0.1--0.2. However, the overwhelming impression is the similarity of the mass functions and age distributions of clusters in these different galaxies, including giant and dwarf, quiescent and interacting galaxies. This is an important empirical result, justifying terms such as "universal" or "quasi-universal." We provide a partial theoretical explanation for these observations in terms of physical processes operating during the formation and disruption of the clusters, including star formation and feedback, subsequent stellar mass loss, and tidal interactions with passing molecular clouds. A full explanation will require additional information about the molecular clumps and star clusters in galaxies beyond the Milky Way.
In this paper we describe the first data release of the the Visible and
Infrared Survey Telescope for Astronomy (VISTA) Deep Extragalactic Observations
(VIDEO) survey. VIDEO is a ~12degree^2 survey in the near-infrared Z,Y,J,H and
K_s bands, specifically designed to enable the evolution of galaxies and large
structures to be traced as a function of both epoch and environment from the
present day out to z=4, and active galactic nuclei (AGN) and the most massive
galaxies up to and into the epoch of reionization. With its depth and area,
VIDEO will be able to fully explore the period in the Universe where AGN and
starburst activity were at their peak and the first galaxy clusters were
beginning to virialize. VIDEO therefore offers a unique data set with which to
investigate the interplay between AGN, starbursts and environment, and the role
of feedback at a time when it was potentially most crucial.
We provide data over the VIDEO-XMM3 tile, which also covers the
Canada-France-Hawaii-Telescope Legacy Survey Deep-1 field (CFHTLS-D1). The
released VIDEO data reach a 5-sigma AB-magnitude depth of Z=25.7, Y=24.5,
J=24.4, H=24.1 and K_s=23.8 in 2 arcsec diameter apertures (the full depth of
Y=24.6 will be reached within the full integration time in future releases).
The data are compared to previous surveys over this field and we find good
astrometric agreement with the Two-Micron All Sky Survey, and source counts in
agreement with the recently released UltraVISTA survey data. The addition of
the VIDEO data to the CFHTLS-D1 optical data increases the accuracy of
photometric redshifts and significantly reduces the fraction of catastrophic
outliers over the redshift range 0<z<1 from 5.8 to 3.1 per cent in the absence
of an i-band luminosity prior. (Truncated Abstract)
This paper is an invited commentary on Tamas Budavari's presentation, "On statistical cross-identification in astronomy," for the Statistical Challenges in Modern Astronomy V conference held at Pennsylvania State University in June 2011. I begin with a brief review of previous work on probabilistic (Bayesian) assessment of directional and spatio-temporal coincidences in astronomy (e.g., cross-matching or cross-identification of objects across multiple catalogs). Then I discuss an open issue in the recent innovative work of Budavari and his colleagues on large-scale probabilistic cross-identification: how to assign prior probabilities that play an important role in the analysis. With a simple toy problem, I show how Bayesian multilevel modeling (hierarchical Bayes) provides a principled framework that justifies and generalizes pragmatic rules of thumb that have been successfully used by Budavari's team to assign priors.
To broaden the understanding of classical Cepheid structure, evolution and atmospheres, we have extended our continuing secret lives of Cepheids program by obtaining XMM/Chandra X-ray observations, and Hubble space telescope (HST) / cosmic origins spectrograph (COS) FUV-UV spectra of the bright, nearby Cepheids Polaris, {\delta} Cep and {\beta} Dor. Previous studies made with the international ultraviolet explorer (IUE) showed a limited number of UV emission lines in Cepheids. The well-known problem presented by scattered light contamination in IUE spectra for bright stars, along with the excellent sensitivity & resolution combination offered by HST/COS, motivated this study, and the spectra obtained were much more rich and complex than we had ever anticipated. Numerous emission lines, indicating 10^4 K up to ~3 x 10^5 K plasmas, have been observed, showing Cepheids to have complex, dynamic outer atmospheres that also vary with the photospheric pulsation period. The FUV line emissions peak in the phase range {\phi} ~ 0.8-1.0 and vary by factors as large as 10x. A more complete picture of Cepheid outer atmospheres is accomplished when the HST/COS results are combined with X-ray observations that we have obtained of the same stars with XMM-Newton & Chandra. The Cepheids detected to date have X-ray luminosities of log Lx ~ 28.5-29.1 ergs/sec, and plasma temperatures in the 2-8 x 10^6 K range. Given the phase-timing of the enhanced emissions, the most plausible explanation is the formation of a pulsation-induced shocks that excite (and heat) the atmospheric plasmas surrounding the photosphere. A pulsation-driven {\alpha}^2 equivalent dynamo mechanism is also a viable and interesting alternative. However, the tight phase-space of enhanced emission (peaking near 0.8-1.0 {\phi}) favor the shock heating mechanism hypothesis.
We report the first results of a multi-epoch search for wide (separations greater than a few tens of AU), low-mass tertiary companions of a volume-limited sample of 118 known spectroscopic binaries within 30 pc of the Sun, using the 2MASS Point Source Catalog and follow-up observations with the KPNO and CTIO 4m telescopes. Note that this sample is not volume-complete but volume-limited, and, thus, there is incompleteness in our reported companion rates. We are sensitive to common proper motion companions with separations from roughly 200 AU to 10,000 AU (~10" -> ~10'). From 77 sources followed-up to date, we recover 11 previously known tertiaries, three previously known candidate tertiaries, of which two are spectroscopically confirmed and one rejected, and three new candidates, of which two are confirmed and one rejected. This yields an estimated wide tertiary fraction of 19.5^+5.2%_-3.7%. This observed fraction is consistent with predictions set out in star formation simulations where the fraction of wide, low-mass companions to spectroscopic binaries is >10%, and is roughly twice the wide companion rate of single stars.
In this study, we analyze nine CMEs from the Sun to Earth as observed in both the remote sensing and in situ data sets. To date, this is the largest study of Earth impacting CMEs using the multi-view point remote sensing and in situ data. However, the remote sensing and in situ data of the same CME cannot be directly compared. Thus, we use several models to parameterize the two data sets. With the model results, we are able to compare the arrival time, Earth impact speed, internal magnetic field, size and orientation as derived from the remote sensing and in situ methods. From the derived kinematics, we compare the predicted arrival times and impact velocities with the in situ data. We find that even with nearly continuous observations and the best available model of the CME structure, there is still a significant error in the predicted values. We estimate the various forces acting on the CME as predicted by three theoretical models of CME propagation and expansion and compare these results with the observational results. We find that the flux rope model of Chen (1989) provides the best agreement with the observations. With the flux rope model, we are able to predict the internal magnetic field of the CME near Earth from the remote sensing data to an order of magnitude. Finally, we compare the size and orientation of the CMEs as predicted from the remote sensing and in situ data. We find very little agreement between the values derived from the two data sets.
We used WISE-derived geometric albedos (p_V) and diameters, as well as geometric albedos and diameters from the literature, to produce more accurate diurnal Yarkovsky drift predictions for 540 near-Earth asteroids (NEAs) out of the current sample of \sim 8,800 known objects. As ten of the twelve objects with the fastest predicted rates have observed arcs of less than a decade, we list upcoming apparitions of these NEAs to facilitate observations.
We study dynamical electroweak symmetry breaking in the setting of ultra-high--energy cosmic rays. The additional gauge interactions in these extension of the standard model can significantly modify the cosmic-ray--air cross section and thus lead to the direct detection or exclusion of a particular model as well as to explanations for features of the cosmic-ray spectrum.
In this letter, we study the effects of an interaction between dark matter and dark energy through a two scalar field model with a potential $V(\phi,\chi)=e^{-\lambda\phi}P(\phi,\chi)$, where $P(\phi,\chi)$ is a polynomial. We show that features of the present Universe are reproduced for a large range of the bare mass of the dark matter field. Simple modifications of the potential are studied, revealing important implications of the interaction, including the possibility of transient acceleration solutions.
Information field dynamics (IFD) is introduced here as a framework to derive numerical schemes for the simulation of physical and other fields. Any simulation scheme updates a discretized field representation, the data in a computer's memory, for the next time step according to a discretized, approximate representation of the underlying field dynamics. Assumptions about the continuum field behavior on sub-grid scales are reflected in these rules, e.g. the field might be assumed to be constant within a grid cell, or to be some weighted average of neighboring data points, and the like. In contrast to such parametrized sub-grid field structures, IFD constructs non-parametric sub-grid field configurations from the combination of the data, representing constraints on possible field configurations, and prior assumptions on the sub-grid field statistics. Each of these field configurations can formally be evolved to a later moment since any differential operator of the dynamics can act on fields living in continuous space. However, these virtually evolved fields need again a representation by data in computer memory. The maximum entropy principle guides the construction of updated datasets via entropic matching, optimally representing these field configurations at the later time. The field dynamics thereby become represented by a finite set of evolution equations for the data that can be solved numerically. These should provide a more accurate description of the physical field dynamics, due to the more rigorous accounting of sub-grid physics and the space discretization process. The IFD approach is illustrated using the example of a coarsely discretized representation of a thermally excited classical Klein-Gordon field. The next steps towards the construction of IFD simulation schemes for more complex systems, e.g. for turbulent hydrodynamics, are also briefly discussed here.
Links to: arXiv, form interface, find, astro-ph, recent, 1206, contact, help (Access key information)
Dwarf spheroidal galaxies are often conjectured to be the sites of the first stars. The best current contenders for finding the chemical imprints from the enrichment by those massive objects are the "ultrafaint dwarfs" (UFDs). Here we present evidence for remarkably low heavy element abundances in the metal poor Hercules UFD. Combined with other peculiar abundance patterns this indicates that Hercules was likely only influenced by very few, massive explosive events - thus bearing the traces of an early, localized chemical enrichment with only very little other contributions from other sources at later times.
We present the first detailed study of the stellar populations of star-forming galaxies at z~1.5, which are selected by their [O II] emission line, detected in narrow-band surveys. We identified ~1,300 [O II] emitters at z=1.47 and z=1.62 in the Subaru Deep Field with rest-frame EWs above 13\AA. Optical and near-infrared spectroscopic observations for ~10% of our samples show that our two-color identification of [O II] emission-line galaxies is 99% successful. We analyze the multi-wavelength properties of a subset of ~1,200 galaxies with the best photometry. They have average rest-frame EW of 45\AA, stellar mass of 3 x 10^9 M_sun, and stellar age of 100 Myr. In addition, our SED fitting and broad-band colors indicate that [O II] emitters span the full range of galaxy populations at z~1.5. We also find that 80% of [O II] emitters are also photometrically classified as "BX/BM" (UV) galaxies and/or the star-forming "BzK" (near-IR) galaxies. Our [O II] emission line survey produces a far more complete, and somewhat deeper sample of z~1.5 galaxies than either the BX/BM or sBzK selection alone. We constructed average SEDs and find that higher [O II] EW galaxies have somewhat bluer continua. SED model-fitting shows that they have on average half the stellar mass of galaxies with lower [O II] EW. The observed [O II] luminosity is well-correlated with the far-UV continuum with a slope slightly less than one (0.89). However, this offset is only marginally significant (< 1\sigma). The fit also suggests that [O II] can be used as a SFR indicator, but with an additional systematic uncertainty of ~0.2 dex.
We study the shapes of galaxy dark matter haloes by measuring the anisotropy of the weak gravitational lensing signal around galaxies in the second Red-sequence Cluster Survey (RCS2). We determine the average shear anisotropy within the virial radius for three lens samples: all galaxies with 19<m_r'<21.5, and the `red' and `blue' samples, whose lensing signals are dominated by massive low-redshift early-type and late-type galaxies, respectively. To study the environmental dependence of the lensing signal, we separate each lens sample into an isolated and clustered part and analyse them separately. We also measure the azimuthal dependence of the distribution of physically associated galaxies around the lens samples. We find that these satellites preferentially reside near the major axis of the lenses, and constrain the angle between the major axis of the lens and the average location of the satellites to <theta>=43.7 deg +/- 0.3 deg for the `all' lenses, <theta>=41.7 deg +/- 0.5 deg for the `red' lenses and <theta>=42.0 deg +/- 1.4 deg for the `blue' lenses. For the `all' sample, we find that the anisotropy of the galaxy-mass cross-correlation function <f-f_45>=0.23 +/- 0.12, providing weak support for the view that the average galaxy is embedded in, and preferentially aligned with, a triaxial dark matter halo. Assuming an elliptical Navarro-Frenk-White (NFW) profile, we find that the ratio of the dark matter halo ellipticity and the galaxy ellipticity f_h=e_h/e_g=1.50+1.03-1.01, which for a mean lens ellipticity of 0.25 corresponds to a projected halo ellipticity of e_h=0.38+0.26-0.25 if the halo and the lens are perfectly aligned. For isolated galaxies of the `all' sample, the average shear anisotropy increases to <f-f_45>=0.51+0.26-0.25 and f_h=4.73+2.17-2.05, whilst for clustered galaxies the signal is consistent with zero. (abridged)
Near-IR observations are important for the detection and characterization of exoplanets using the transit technique, either in surveys of large numbers of stars or for follow-up spectroscopic observations of individual planets. In a controlled laboratory experiment, we imaged $\sim 10^4$ critically sampled spots onto an Teledyne Hawaii-2RG (H2RG) detector to emulate an idealized star-field. We obtained time-series photometry of up to $\simeq 24$ hr duration for ensembles of $\sim 10^3$ pseudo-stars. After rejecting correlated temporal noise caused by various disturbances, we measured a photometric performance of $<$50 ppm-hr$^{-1/2}$ limited only by the incident photon rate. After several hours we achieve a photon-noise limited precision level of $10\sim20$ ppm after averaging many independent measurements. We conclude that IR detectors such as the H2RG can make the precision measurements needed to detect the transits of terrestrial planets or detect faint atomic or molecular spectral features in the atmospheres of transiting extrasolar planets.
Ground-based optical surveys such as PanSTARRS, DES, and LSST, will produce large catalogs to limiting magnitudes of r > 24. Star-galaxy separation will pose a major challenge to such surveys because galaxies---even very compact galaxies---outnumber halo stars at these depths. Here we investigate photometric classification techniques on stars and galaxies with intrinsic FWHM < 0.2 arcsec. We consider unsupervised SED template fitting and supervised, data-driven Support Vector Machines (SVM). For template fitting, we use a Maximum Likelihood (ML) method and a new Hierarchical Bayesian (HB) method, in which we learn the prior distribution of template probabilities by optimizing the likelihood for the entire dataset. SVM requires training data to classify unknown sources; ML and HB don't. We consider both i.) a best-case scenario (SVM_best) in which the training data is (unrealistically) a random sampling of the data in both signal-to-noise and demographics, and ii.) a more realistic scenario in which the SVM is trained on higher signal-to-noise data (SVM_real) at brighter apparent magnitudes. Testing with COSMOS ugriz data we find that HB outperforms ML, delivering ~80% completeness in both star and galaxy samples, with purity of ~40-90% and ~70-90% for stars and galaxies, respectively. We find no algorithm delivers perfect performance, and that studies of M-giant and metal-poor main-sequence turnoff stars may be most affected by poor star-galaxy separation. We measure the area under the ROC curve to assess the relative performance of the approaches and find a best-to-worst ranking of SVM_best, HB, ML, and SVM_real. We conclude, therefore, that a well trained SVM will outperform template-fitting methods. However, a normally trained SVM performs worse. Thus, Hierarchical Bayesian template fitting may prove to be the optimal method for source classification in future surveys.
Advancements in infrared IR open up the possibility to spatially resolve AGN on the parsec-scale level and study the circumnuclear dust distribution, commonly referred to as the "dust torus", that is held responsible for the type 1/type 2 dichotomy of AGN. We used the mid-IR beam combiner MIDI together with the 8m telescopes at the VLTI to observe the nucleus of the Seyfert 2 galaxy NGC 424, achieving an almost complete coverage of the uv-plane accessible by the available telescope configurations. We detect extended mid-IR emission with a relatively baseline- and model-independent mid-IR half-light radius of (2.0 \pm 0.2) pc \times (1.5 \pm 0.3) pc (averaged over the 8-13 {\mu}m wavelength range). The extended mid-IR source shows an increasing size with wavelength. The orientation of the major axis in position angle -27deg is closely aligned with the system axis as set by optical polarization observations. Torus models typically favor extension along the mid-plane at mid-IR wavelengths instead. Therefore, we conclude that the majority of the pc-scale mid-IR emission (>~60%) in this type 2 AGN originates from optically-thin dust in the polar region of the AGN, a scenario consistent with the near- to far-IR SED. We suggest that a radiatively-driven dusty wind, possibly launched in a puffed-up region of the inner hot part of the torus, is responsible for the polar dust. In this picture, the torus dominates the near-IR emission up to about 5 {\mu}m, while the polar dust is the main contributor to the mid-IR flux. Our results of NGC 424 are consistent with recent observations of the AGN in the Circinus galaxy and resemble large-scale characteristics of other objects. If our results reflect a general property of the AGN population, the current paradigm for interpreting and modeling the IR emission of AGN have to be revised. (abridged)
Observations from multiple gamma-ray telescopes have uncovered a high energy gamma-ray source spatially coincident with the Galactic center. Recently, a compelling model for the broad-band gamma-ray emission has been formulated which posits that high energy protons emanating from Sgr A* could produce gamma-rays through pi0 decays resulting from inelastic collisions with the traversed interstellar gas in the region. Models of the gas distribution in the Galactic center region imply that the resulting gamma-ray morphology would be observed as a point source with all current telescopes, but that the upcoming Cherenkov Telescope Array (CTA) may be able to detect an extended emission profile with an unmistakable morphology. Here, we critically evaluate this claim, employing a three dimensional gas distribution model and a detailed Monte Carlo simulation, and using the anticipated effective area and angular resolution of CTA. We find that the impressive angular resolution of CTA will be key to test hadronic emission models conclusively against, for example, point source or dark matter annihilation scenarios. We comment on the relevance of this result for searches for dark matter annihilation in the Galactic center region.
Recent work has suggested that the stellar initial mass function (IMF) is not universal, but rather is correlated with galaxy stellar mass, stellar velocity dispersion, or morphological type. In this paper, we investigate variations of the IMF within individual galaxies. For this purpose, we use strong lensing and gas kinematics to measure independently the normalization of the IMF of the bulge and disk components of a sample of 5 massive spiral galaxies with substantial bulge components taken from the SWELLS survey. We find that the stellar mass of the bulges is tightly constrained by the lensing and kinematic data and consistent with that obtained by stellar population synthesis models fits to their colours only by assuming a Salpeter-like normalisation of the IMF. Conversely, the disk masses are less well constrained due to degeneracies with the dark matter halo, but are consistent with Milky Way type IMFs in agreement with previous studies. The disks are submaximal at 2.2 disk scale lengths, but due to the contribution of the bulges, the galaxies are baryon dominated at 2.2 disk scale lengths. Globally, our inferred IMF normalisation is consistent with that found for early-type galaxies of comparable stellar mass (>10^11 M_sun). Our results provide evidence for a non-universal IMF within the different components of spiral galaxies, adding to the well-known differences in stellar populations between disks and bulges.
To investigate possible variations in the stellar initial mass function (IMF) in red-sequence galaxies, we have obtained infrared spectroscopy with Subaru/FMOS for a sample of 92 red-sequence galaxies in the Coma cluster. Velocity dispersions, ages and element abundances for these galaxies have been previously determined from optical data. By stacking the FMOS spectra in the rest frame, removing sky-subtraction residuals and other artefacts fixed in the observed frame, we derive composite spectra in the 9600-10500 Angstrom range for galaxies grouped according to their velocity dispersion or Mg/Fe ratio. We measure the Wing-Ford band of FeH and a new index centred on a CaI line at 10345 Angstrom; these features are strong in cool dwarf stars, and hence reflect the form of the IMF at low mass (<0.5M_sun) The CaI line, unlike the Wing-Ford band and other `classical' IMF indicators (NaI doublet, CaII triplet), is unaffected by the abundance of sodium. We compare the measured indices against predictions from spectral synthesis models matched to the element abundances estimated from the optical data. Binning galaxies by velocity dispersion, we find that both the Wing-Ford and CaI index measurements are best reproduced by models with the Salpeter IMF. There is no clear evidence for an increase in dwarf-star content with velocity dispersion. Binning the observed galaxies instead by Mg/Fe ratio, the behaviour of both indices implies a trend of IMF from Chabrier-like, at abundance ratios close to solar, to Salpeter or heavier for highly alpha-enhanced populations. At face value, this suggests that the IMF depends on the mode of star formation, with intense rapid star-bursts generating a larger population of low-mass stars.
As brown dwarfs cool, a variety of species condense in their atmospheres, forming clouds. Iron and silicate clouds shape the emergent spectra of L dwarfs, but these clouds dissipate at the L/T transition. A variety of other condensates are expected to form in cooler T dwarf atmospheres. These include Cr, MnS, Na2S, ZnS, and KCl, but the opacity of these optically thinner clouds has not been included in previous atmosphere models. Here, we examine their effect on model T and Y dwarf atmospheres. The cloud structures and opacities are calculated using the Ackerman & Marley (2001) cloud model, which is coupled to an atmosphere model to produce atmospheric pressure-temperature profiles in radiative-convective equilibrium. We generate a suite of models between Teff = 400 and 1300 K, log g=4.0 and 5.5, and condensate sedimentation efficiencies from fsed=2 to 5. Model spectra are compared to two red T dwarfs, Ross 458C and UGPS 0722-05; models that include clouds are found to match observed spectra significantly better than cloudless models. The emergence of sulfide clouds in cool atmospheres, particularly Na2S, may be a more natural explanation for the "cloudy" spectra of these objects, rather than the re-emergence of silicate clouds that wane at the L-to-T transition. We find that sulfide clouds provide a mechanism to match the near- and mid-infrared colors of observed T dwarfs. Our results indicate that including the opacity of condensates in T dwarf atmospheres is necessary to accurately determine the physical characteristics of many of the observed objects.
We present one-dimensional models of the hot gas in dark-matter halos, which both predict the existence of cool cores and explain their structure. Our models are directly applicable to semi-analytic models (SAMs) of galaxy formation. We have previously argued that filaments of cold (~10^4 K) gas condense out of the intracluster medium (ICM) in hydrostatic and thermal equilibrium when the ratio of the thermal instability timescale to the free-fall time $t_{TI}/t_{ff}$ falls below 5-10. This criterion corresponds to an upper limit on the density of the ICM and motivates a model in which a density core forms wherever $t_{TI}/t_{ff} \lesssim 10$. Consistent with observations and numerical simulations, this model predicts larger and more tenuous cores for lower-mass halos---while the core density in a cluster may be as large as ~ 0.1 cm^{-3}, the core density in the Galactic halo should not exceed ~ 10^{-4} cm^{-3}. Our models produce a favorable match to the observational X-ray luminosity-temperature (Lx-Tx) relation. For halo masses $\lesssim 10^{13}$ Msun the core size approaches the virial radius. Thus, most of the baryons in such halos cannot be in the hot ICM, but either in the form of stars or in the form of hot gas beyond the virial radius. Because of the smaller mass in the ICM and much larger mass available for star formation, the majority of the baryons in low mass halos ($\lesssim 10^{13}$ Msun) can be expelled beyond the virial radius due to supernova feedback. This can account for the baryons `missing' from low mass halos, such as the Galactic halo.
It is well known that current spectroscopic determinations of the chemical composition of the Sun are starkly at odds with the metallicity implied by helioseismology. We propose that the discrepancy may be due to conversion of photons to a new light boson in the solar photosphere. We investigate the impact of particles with axion-like interactions with the photon on the inferred photospheric abundances. Conversion of photons into new light bosons can in principle easily produce the +0.2 dex change in derived abundances required to reconcile spectroscopic and helioseismological determinations of the solar metallicity. We show that this is however not possible for any of the presently-allowed parameter space of either the QCD axion or a standard axion-like particle. We speculate that other models involving photon-boson mixing, such as hidden photons, might prove more successful.
We study the secular, hierarchical three-body problem to first-order in a post-Newtonian expansion of General Relativity. We expand the first-order post-Newtonian Hamiltonian to leading-order in the ratio of the semi-major axis of the two orbits. In addition to the well-known terms that correspond to the GR precession of the inner and outer orbits, we find a new secular post-Newtonian interaction term that can affect the long-term evolution of the triple. We explore the parameter space for highly inclined and eccentric systems, where the Kozai-Lidov mechanism can produce large-amplitude oscillations in the eccentricities. The standard lore, i.e., that General Relativity effects suppress eccentricity, is only consistent with the parts of phase space where the General Relativity timescales are several orders of magnitude shorter than the secular Newtonian one. In other parts of phase space, however, post-Newtonian corrections combined with the three body ones, can excite eccentricities. In particular, for systems where the General Relativity timescale is comparable to the secular Newtonian timescales, the three-body interactions give rise to a resonant-like eccentricity excitation. Furthermore, for triples with a comparable-mass inner binary, where the eccentric Kozai-Lidov mechanism is suppressed, post-Newtonian corrections can further increase the eccentricity and lead to orbital flips even when the timescale of the former is much longer than the timescale of the secular Kozai-Lidov quadrupole perturbations.
We study the matter and velocity divergence power spectra in a f(R) gravity theory and their time evolution measured from several large-volume N-body simulations with varying box sizes and resolution. We find that accurate prediction of the matter power spectrum in f(R) gravity places stronger requirements on the simulation than is the case with LCDM, because of the nonlinear nature of the fifth force. Linear perturbation theory is shown to be a poor approximation for the f(R) models, except when the chameleon effect is very weak. We show that the relative differences from the fiducial LCDM model are much more pronounced in the nonlinear tail of the velocity divergence power spectrum than in the matter power spectrum, which suggests that future surveys which target the collection of peculiar velocity data will open new opportunities to constrain modified gravity theories. A close investigation of the time evolution of the power spectra shows that there is a pattern in the evolution history, which can be explained by the properties of the chameleon-type fifth force in f(R) gravity. Varying the model parameter |f_R0|, which quantifies the strength of the departure from standard gravity, mainly varies the epoch marking the onset of the fifth force, as a result of which the different f(R) models are in different stages of the same evolutionary path at any given time
Vista Variables in The Via Lactea (VVV) is an ESO variability survey that is performing observations in near infrared bands (ZYJHKs) towards the Galactic bulge and part of the disk with the completeness limits at least 3 mag deeper than 2MASS. In the present work, we searched in the VVV survey data for background galaxies near the Galactic plane using ZYJHKs photometry that covers 1.636 square degrees. We identified 204 new galaxy candidates by analyzing colors, sizes, and visual inspection of multi-band (ZYJHKs) images. The galaxy candidates colors were also compared with the predicted ones by star counts models considering a more realistic extinction model at the same completeness limits observed by VVV. A comparison of the galaxy candidates with the expected one by Milennium simulations is also presented. Our results increase the number density of known galaxies behind the Milky Way by more than one order of magnitude. A catalog with galaxy properties including ellipticity, Petrosian radii and ZYJHKs magnitudes is provided, as well as comparisons of the results with other surveys of galaxies towards Galactic plane.
The present observational understanding of the evolution of the mass accretion rates (Macc) in pre-main sequence stars is limited by the lack of accurate measurements of Macc over homogeneous and large statistical samples of young stars. Such observational effort is needed to properly constrain the theory of star formation and disk evolution. Based on HST/WFPC2 observations, we present a study of Macc for a sample of \sim 700 sources in the Orion Nebula Cluster, ranging from the Hydrogen-burning limit to M\ast \sim 2M\odot. We derive Macc from both the U-band excess and the H{\alpha} luminosity (LH{\alpha}), after determining empirically both the shape of the typical accretion spectrum across the Balmer jump and the relation between the accretion luminosity (Lacc) and LH{\alpha}, that is Lacc/L\odot = (1.31\pm0.03)\cdotLH{\alpha}/L\odot + (2.63\pm 0.13). Given our large statistical sample, we are able to accurately investigate relations between Macc and the parameters of the central star such as mass and age. We clearly find Macc to increase with stellar mass, and decrease over evolutionary time, but we also find strong evidence that the decay of Macc with stellar age occurs over longer timescales for more massive PMS stars. Our best fit relation between these parameters is given by: log(Macc/M\odot\cdotyr)=(-5.12 \pm 0.86) -(0.46 \pm 0.13) \cdot log(t/yr) -(5.75 \pm 1.47)\cdot log(M\ast/M\odot) + (1.17 \pm 0.23)\cdot log(t/yr) \cdot log(M\ast/M\odot). These results also suggest that the similarity solution model could be revised for sources with M\ast > 0.5M\odot. Finally, we do not find a clear trend indicating environmental effects on the accretion properties of the sources.
We have used deep, HST, near-IR imaging to study the morphological properties of the most massive galaxies at high z, modelling the WFC3/IR H-band images of the ~200 galaxies in the CANDELS-UDS field with 1 < z_phot < 3, and stellar masses M_star > 10^11 M_sun. We have used both single-Sersic and bulge+disk models, have investigated the errors/biases introduced by uncertainties in the background and the PSF, and have obtained formally-acceptable model fits to >90% of the galaxies. Our results indicate that these massive galaxies at 1 < z < 3 lie both on and below the local size-mass relation, with a median R_e~2.6 kpc, a factor of ~2.3 smaller than comparably-massive local galaxies. Moreover, we find that bulge-dominated objects in particular show evidence for a growing bimodality in the size-mass relation with increasing z, and by z > 2 the compact bulges display effective radii a factor ~4 smaller than local ellipticals of comparable mass. These trends appear to extend to the bulge components of disk-dominated galaxies, and vice versa. We also find that, while such massive galaxies at low z are bulge-dominated, at 1 < z < 2 they are predominantly mixed bulge+disk systems, and by z > 2 they are mostly disk-dominated. The majority of the disk-dominated galaxies are actively forming stars, but this is also true for many of the bulge-dominated systems. Interestingly, however, while most of the quiescent galaxies are bulge-dominated, we find that a significant fraction (25-40%) of the most quiescent galaxies have disk-dominated morphologies. Thus, while our results show that the massive galaxy population is undergoing dramatic changes at this crucial epoch, they also suggest that the physical mechanisms which quench star-formation activity are not simply connected to those responsible for the morphological transformation of massive galaxies into present-day giant ellipticals.
It has been posited that lunar eclipse observations may help predict the in-transit signature of Earth-like extrasolar planets. However, a comparative analysis of the two phenomena addressing in detail the transport of stellar light through the planet's atmosphere has not yet been presented. Here, we proceed with the investigation of both phenomena by making use of a common formulation. Our starting point is a set of previously unpublished near-infrared spectra collected at various phases during the August 2008 lunar eclipse. We then take the formulation to the limit of an infinitely distant observer in order to investigate the in-transit signature of the Earth-Sun system as being observed from outside our Solar System. The refraction-bending of sunlight rays that pass through the Earth's atmosphere is a critical factor in the illumination of the eclipsed Moon. Likewise, refraction will have an impact on the in-transit transmission spectrum for specific planet-star systems depending on the refractive properties of the planet's atmosphere, the stellar size and the planet's orbital distance. For the Earth-Sun system, at mid-transit, refraction prevents the remote observer's access to the lower ~12-14 km of the atmosphere and, thus, also to the bulk of the spectroscopically-active atmospheric gases. We demonstrate that the effective optical radius of the Earth in transit is modulated by refraction and varies by ~12 km from mid-transit to 2nd contact. The refractive nature of atmospheres, a property which is rarely accounted for in published investigations, will pose additional challenges to the characterization of Earth-like extrasolar planets. Refraction may have a lesser impact for Earth-like extrasolar planets within the habitable zone of some M-type stars.
Recent high resolution spectropolarimetric observations from Hinode detected the presence of supersonic downflows in a sunspot light bridge (Louis et al. 2009). These downflows occurred in localized patches, close to regions where the field azimuth changed by a large value. This apparent discontinuity in the field azimuth was seen along a thin ridge running along the western edge of the light bridge. Some, but not all, of these downflowing patches were co-spatial with chromospheric brightness enhancements seen in Ca II H filtergrams. The presence of magnetic inhomogeneities at scales of 0.3 arcsec could facilitate the reconnection of field lines in the lower chromosphere whose signatures might be the supersonic downflows and the brightness enhancements that have been observed.
We report radio interferometric observations of the 12C16O 1.3 mm J = 2-1 emission line in the circumstellar envelope of the M supergiant Alpha Ori and have detected and separated both the S1 and S2 flow components for the first time. Observations were made with the Combined Array for Research in Millimeter-wave Astronomy (CARMA) interferometer in the C, D, and E antenna configurations. We obtain good u-v coverage (5-280 klambda) by combining data from all three configurations allowing us to trace spatial scales as small as 0.9\arcsec over a 32\arcsec field of view. The high spectral and spatial resolution C configuration line profile shows that the inner S1 flow has slightly asymmetric outflow velocities ranging from -9.0 km s-1 to +10.6 km s-1 with respect to the stellar rest frame. We find little evidence for the outer S2 flow in this configuration because the majority of this emission has been spatially-filtered (resolved out) by the array. We also report a SOFIA-GREAT CO(J= 12-11) emission line profile which we associate with this inner higher excitation S1 flow. The outer S2 flow appears in the D and E configuration maps and its outflow velocity is found to be in good agreement with high resolution optical spectroscopy of K I obtained at the McDonald Observatory. We image both S1 and S2 in the multi-configuration maps and see a gradual change in the angular size of the emission in the high absolute velocity maps. We assign an outer radius of 4\arcsec to S1 and propose that S2 extends beyond CARMA's field of view (32\arcsec at 1.3 mm) out to a radius of 17\arcsec which is larger than recent single-dish observations have indicated. When azimuthally averaged, the intensity fall-off for both flows is found to be proportional to R^{-1}, where R is the projected radius, indicating optically thin winds with \rho \propto R^{-2}.
We study Hinode/SOT-FG observations of intensity fluctuations in Ca II H-line and G-band image sequences and their relation to simultaneous and co-spatial magnetic field measurements. We explore the G-band and H-line intensity oscillation spectra both separately and comparatively via their relative phase differences, time delays and cross-coherences. In the non-magnetic situations, both sets of fluctuations show strong oscillatory power in the 3 - 7 mHz band centered at 4.5 mHz, but this is suppressed as magnetic field increases. A relative phase analysis gives a time delay of H-line after G-band of 20\pm1 s in non-magnetic situations implying a mean effective height difference of 140 km. The maximum coherence is at 4 - 7 mHz. Under strong magnetic influence the measured delay time shrinks to 11 s with the peak coherence near 4 mHz. A second coherence maximum appears between 7.5 - 10 mHz. Investigation of the locations of this doubled-frequency coherence locates it in diffuse rings outside photospheric magnetic structures. Some possible interpretations of these results are offered.
We report on a discovery of a diffuse nebula containing a pointlike source in the southern blowout region of the Cygnus Loop supernova remnant, based on Suzaku and XMM-Newton observations. The X-ray spectra from the nebula and the pointlike source are well represented by an absorbed power-law model with photon indices of 2.2+/-0.1 and 1.6+/-0.2, respectively. The photon indices as well as the flux ratio of F_nebula/F_pointlike ~ 4 lead us to propose that the system is a pulsar wind nebula, although pulsations have not yet been detected. If we attribute its origin to the Cygnus Loop supernova, then the 0.5-8 keV luminosity of the nebula is computed to be 2.1e31 (d/540pc)^2 ergs/s, where d is the distance to the Loop. This implies a spin-down loss-energy E_dot ~ 2.6e35 (d/540pc)^2 ergs/s. The location of the neutron star candidate, ~2 degrees away from the geometric center of the Loop, implies a high transverse velocity of ~1850 (d/540pc)(t/10kyr)^{-1} km/s, assuming the currently accepted age of the Cygnus Loop.
Jupiter is expected to pulsate in a spectrum of acoustic modes and recent re-analysis of a spectroscopic time series has identified a regular pattern in the spacing of the frequencies \citep{gaulme2011}. This exciting result can provide constraints on gross Jovian properties and warrants a more in-depth theoretical study of the seismic structure of Jupiter. With current instrumentation, such as the SYMPA instrument \citep{schmider2007} used for the \citet{gaulme2011} analysis, we assume that, at minimum, a set of global frequencies extending up to angular degree $\ell=25$ could be observed. In order to identify which modes would best constrain models of Jupiter's interior and thus help motivate the next generation of observations, we explore the sensitivity of derived parameters to this mode set. Three different models of the Jovian interior are computed and the theoretical pulsation spectrum from these models for $\ell\leq 25$ is obtained. We compute sensitivity kernels and perform linear inversions to infer details of the expected discontinuities in the profiles in the Jovian interior. We find that the amplitude of the sound-speed jump of a few percent in the inner/outer envelope boundary seen in two of the applied models should be reasonably inferred with these particular modes. Near the core boundary where models predict large density discontinuities, the location of such features can be accurately measured, while their amplitudes have more uncertainty. These results suggest that this mode set would be sufficient to infer the radial location and strength of expected discontinuities in Jupiter's interior, and place strong constraints on the core size and mass. We encourage new observations to detect these Jovian oscillations.
Despite expanding research activity in gravitational lens modeling, there is no particular software which is considered a standard. Much of the gravitational lens modeling software is written by individual investigators for their own use. Some gravitational lens modeling software is freely available for download but is widely variable with regard to ease of use and quality of documentation. This review of 12 software packages was undertaken to provide a single source of information. Gravitational lens models are classified as parametric models or non-parametric models, and can be further divided into research and educational software. Software used in research includes GravLens / LensModel, Lenstool, Lensperfect, Glafic, PixeLens, SimpLens, Lensview, and GRALE. In this review, GravLensHD, G-Lens, Gravitational Lensing and Lens are categorized as educational programs that are useful for demonstrating various aspects of lensing. Each of the 12 software packages is reviewed with regard to software features (installation, documentation, files provided, etc.) and lensing features (type of model, input data, output data, etc.) as well as a brief review of studies where they have been used. Recent studies have demonstrated the utility of strong gravitational lensing data for mass mapping, and suggest increased use of these techniques in the future. Coupled with the advent of greatly improved imaging, new approaches to modeling of strong gravitational lens systems are needed. This is the first systematic review of strong gravitational lens modeling software, providing investigators with a starting point for future software development to further advance gravitational lens modeling research.
Observations with Subaru telescope have detected surprisingly strong Lyman continuum (LyC; $\sim900$ \AA\ in the rest-frame) from some Lyman $\alpha$ emitters (LAEs) at $z=3.1$. We have examined the stellar population which simultaneously accounts for the strength of the LyC and the spectral slope of non-ionizing ultraviolet of the LAEs. As a result, we have found that stellar populations with metallicity $Z\geq1/50 Z_\odot$ can explain the observed LyC strength only with a very top-heavy initial mass function (IMF; $< m > \sim 50 M_\odot$). However, the critical metallicity for such an IMF is expected to be much lower. A very young ($\sim1$ Myr) and massive ($\sim100$ $M_\odot$) extremely metal-poor ($Z\leq5\times10^{-4}Z_\odot$) or metal-free (so-called Population III) stellar population can also reproduce the observed LyC strength if the mass fraction of such 'primordial' stellar population is $\sim1$% in total stellar mass of the LAEs.
Models of spontaneous wave function collapse modify the linear Schr\"{o}dinger equation of standard Quantum Mechanics by adding stochastic non-linear terms to it. The aim of such models is to describe the quantum (linear) nature of microsystems along with the classical nature (violation of superposition principle) of macroscopic ones. The addition of such non-linear terms in the Schr\"{o}dinger equation leads to non-conservation of energy of the system under consideration. Thus, a striking feature of collapse models is to heat non-relativistic particles with a constant rate. If such a process is physical, then it has the ability to perturb the well-understood thermal history of the universe. In this article we will try to investigate the impacts of such heating terms on standard evolution of non-relativistic matter and on the formation of CMBR. We will use the CSL model, the most widely used collapse model. We will also put constraints on the CSL collapse rate $\lambda$ by considering that the standard evolution of non-relativistic matter is not hampered and the observed precise blackbody spectrum of CMBR would not get distorted (in the form of $\mu-$type and $y-$type distortions) so as to violate the observed bounds.
I outline the prospects for performing weak lensing studies with the new generation of radio telescopes that are coming online now and in the future. I include a description of a proposed technique to use polarization observations in radio weak lensing analyses which could prove extremely useful for removing a contaminating signal from intrinsic alignments. Ultimately, the Square Kilometre Array promises to be an exceptional instrument for performing weak lensing studies due to the high resolution, large area surveys which it will perform. In the nearer term, the e-MERLIN instrument in the UK offers the high sensitivity and sub-arcsec resolution required to prove weak lensing techniques in the radio band. I describe the SuperCLASS survey -- a recently accepted e-MERLIN legacy programme which will perform a pioneering radio weak lensing analysis of a supercluster of galaxies.
In 2008 AGILE and Fermi detected gamma-ray flaring activity from the unidentified EGRET source 3EG J1236+0457, recently associated with a flat spectrum radio quasar GB6 J1239+0443 at z=1.762. The optical counterpart of the gamma-ray source underwent a flux enhancement of a factor 15-30 in 6 years, and of ~10 in six months. We interpret this flare-up in terms of a transition from an accretion-disk dominated emission to a synchrotron-jet dominated one. We analysed a Sloan Digital Sky Survey (SDSS) archival optical spectrum taken during a period of low radio and optical activity of the source. We estimated the mass of the central black hole using the width of the CIV emission line. In our work, we have also investigated SDSS archival optical photometric data and UV GALEX observations to estimate the thermal-disk emission contribution of GB6 J1239+0443. Our analysis of the gamma-ray data taken during the flaring episodes indicates a flat gamma-ray spectrum, with an extension of up to 15 GeV, with no statistically-relevant sign of absorption from the broad line region, suggesting that the blazar-zone is located beyond the broad line region. This result is confirmed by the modeling of the broad-band spectral energy distribution (well constrained by the available multiwavelength data) of the flaring activity periods and by the accretion disk luminosity and black hole mass estimated by us using archival data.
We present CCD photometric observations of the slow irregular type variable KU Her. We have also used observations from the ASAS Database and the AAVSO International Database. The V-band observations span eight years from August of 2002 to September of 2010. The length and density of the data enable us to look for variations on time scales ranging from days to years. The V-band observations of KU Her are analysed to derive periodicity. This analysis was an approach to find the fundamental periods of the variability of KU Her. Analysis of the KU Her light variation indicates a complex combination of different periods. Five main and three harmonic periods are identified. In addition to the period analysis, we present the colour variations of the star.
Spectral distortion of the cosmic microwave background provides a unique opportunity to probe primordial perturbations on very small scales by performing large-scale measurements. We discuss in a systematic and pedagogic way all the relevant physical phenomena involved in the production and evolution of the mu-type spectral distortion. Our main results agree with previous estimates (in particular we show that a recently found factor of 3/4 arises from relativistic corrections to the wave energy). We also discuss several subleading corrections such as adiabatic cooling and the effects of bulk viscosity, baryon loading and photon heat conduction. Finally we calculate the transfer function for mu-distortions between the end of the mu-era and now.
We present results of broad band photometric reverberation mapping (RM) to measure the radius of the broad line region, and subsequently the black hole mass (M$_{\rm BH}$), in the nearby, low luminosity active galactic nuclei (AGN) NGC 4395. Using the Wise Observatory's 1m telescope equipped with the SDSS g$'$, r$'$ and i$'$ broad band filters, we monitored NGC 4395 for 9 consecutive nights and obtained 3 light curves each with over 250 data points. The g$'$ and r$'$ bands include time variable contributions from H$\beta$ and H$\alpha$ (respectively) plus continuum. The i$'$ band is free of broad lines and covers exclusively continuum. We show that by looking for a peak in the difference between the cross-correlation and the auto-correlation functions for all combinations of filters, we can get a reliable estimate of the time lag necessary to compute M$_{\rm BH}$. We measure the time lag for H$\alpha$ to be $3.6 \pm 0.8 $ hours, comparable to previous studies using the line resolved spectroscopic RM method. We argue that this lag implies a black hole mass of M$_{\rm BH} = (4.9 \pm 2.6) \times 10^{4}$ \Msun .
Operation of the US Virtual Astronomical Observatory shares some issues with modern physical observatories, e.g., intimidating data volumes and rapid technological change, and must also address unique concerns like the lack of direct control of the underlying and scattered data resources, and the distributed nature of the observatory itself. In this paper we discuss how the VAO has addressed these challenges to provide the astronomical community with a coherent set of science-enabling tools and services. The distributed nature of our virtual observatory-with data and personnel spanning geographic, institutional and regime boundaries-is simultaneously a major operational headache and the primary science motivation for the VAO. Most astronomy today uses data from many resources. Facilitation of matching heterogeneous datasets is a fundamental reason for the virtual observatory. Key aspects of our approach include continuous monitoring and validation of VAO and VO services and the datasets provided by the community, monitoring of user requests to optimize access, caching for large datasets, and providing distributed storage services that allow user to collect results near large data repositories. Some elements are now fully implemented, while others are planned for subsequent years. The distributed nature of the VAO requires careful attention to what can be a straightforward operation at a conventional observatory, e.g., the organization of the web site or the collection and combined analysis of logs. Many of these strategies use and extend protocols developed by the international virtual observatory community.
The probability density function (PDF) of the gas density in subsonic and supersonic, isothermal, driven turbulence is analyzed with a systematic set of hydrodynamical grid simulations with resolutions up to 1024^3 cells. We performed a series of numerical experiments with root mean square (r.m.s.) Mach number M ranging from the nearly incompressible, subsonic (M=0.1) to the highly compressible, supersonic (M=15) regime. We study the influence of two extreme cases for the driving mechanism by applying a purely solenoidal (divergence-free) and a purely compressive (curl-free) forcing field to drive the turbulence. We find that our measurements fit the linear relation between the r.m.s. Mach number and the standard deviation of the density distribution in a wide range of Mach numbers, where the proportionality constant depends on the type of the forcing. In addition, we propose a new linear relation between the standard deviation of the density distribution and the standard deviation of the velocity in compressible modes, i.e. the compressible component of the r.m.s. Mach number. In this relation the influence of the forcing is significantly reduced, suggesting a linear relation between the standard deviation of the density distribution and the standard deviation of the velocity in compressible modes, independent of the forcing, ranging from the subsonic to the supersonic regime.
We propose two detection techniques that take advantage of a small sky area approximation and are based on modifications of the "internal linear combination" (ILC) method, an approach widely used in Cosmology for the separation of the various components that contribute to the microwave background. The main advantage of the proposed approach, especially in handling multi-frequency maps of the same region, is that it does not require the "a priori" knowledge of the spatial power-spectrum of either the CMB and/or the Galactic foreground. Hence, it is more robust, easier and more intuitive to use. The performance of the proposed algorithms is tested with numerical experiments that mimic the physical scenario expected for high Galactic latitude observations with the Atacama Large Millimeter/submillimeter Array (ALMA).
We analyze two regions of the quiet Sun (35.6 x 35.6 Mm^2) observed at high spatial resolution (~100 km) in polarized light by the IMaX spectropolarimeter onboard the Sunrise balloon. We identify 497 small-scale (~400 km) magnetic loops, appearing at an effective rate of 0.25 loop h^{-1} arcsec^{-2}; further, we argue that this number and rate are underestimated by ~30%. However, we find that these small dipoles do not appear uniformly on the solar surface: their spatial distribution is rather filamentary and clumpy, creating dead calm areas, characterized by a very low magnetic signal and a lack of organized loop-like structures at the detection level of our instruments, that cannot be explained as just statistical fluctuations of a Poisson spatial process. We argue that this is an intrinsic characteristic of the mechanism that generates the magnetic fields in the very quiet Sun. The spatio-temporal coherences and the clumpy structure of the phenomenon suggest a recurrent, intermittent mechanism for the generation of magnetic fields in the quietest areas of the Sun.
We present APEX SABOCA 350micron and LABOCA 870micron observations of 11 representative examples of the rare, extremely bright (S_1.4mm > 15mJy), dust-dominated millimeter-selected galaxies recently discovered by the South Pole Telescope (SPT). All 11 sources are robustly detected with LABOCA with 40 < S_870micron < 130mJy, approximately an order of magnitude higher than the canonical submillimeter galaxy (SMG) population. Six of the sources are also detected by SABOCA at >3sigma, with the detections or upper limits providing a key constraint on the shape of the spectral energy distribution (SED) near its peak. We model the SEDs of these galaxies using a simple modified blackbody and perform the same analysis on samples of SMGs of known redshift from the literature. These calibration samples inform the distribution of dust temperature for similar SMG populations, and this dust temperature prior allows us to derive photometric redshift estimates and far infrared luminosities for the sources. We find a median redshift of <z> = 3.0, higher than the <z> = 2.2 inferred for the normal SMG population. We also derive the apparent size of the sources from the temperature and apparent luminosity, finding them to appear larger than our unlensed calibration sample, which supports the idea that these sources are gravitationally magnified by massive structures along the line of sight.
We present a general method for identifying the pre-main-sequence population of any star-forming region, unbiased with respect to the presence or absence of disks, in contrast to samples selected primarily via their mid-infrared emission from Spitzer surveys. We have applied this technique to a new, deep, wide-field, near-infrared imaging survey of the Rho Ophiuchi cloud core to search for candidate low mass members. In conjunction with published Spitzer IRAC photometry, and least squares fits of model spectra (COND, DUSTY, NextGen, and blackbody) to the observed spectral energy distributions, we have identified 948 candidate cloud members within our 90% completeness limits of J=20.0, H=20.0, and K_S=18.50. This population represents a factor of ~3 increase in the number of known young stellar objects in the Rho Ophiuchi cloud. A large fraction of the candidate cluster members (81% +/- 3%) exhibit infrared excess emission consistent with the presence of disks, thus strengthening the possibility of their being bona fide cloud members. Spectroscopic follow-up will confirm the nature of individual objects, better constrain their parameters, and allow an initial mass function to be derived.
We present a new search for variable stars in the Galactic globular cluster M28 (NGC 6626). The search is based on a series of BVI images obtained with the SMARTS Consortium's 1.3m telescope at Cerro Tololo Inter-American Observatory, Chile. The search was carried out using the ISIS v2.2 image subtraction package. We find a total of 25 variable stars in the field of the cluster, 9 being new discoveries. Of the newly found variables, 1 is an ab-type RR Lyrae star, 6 are c-type RR Lyrae, and 2 are long-period/semi-regular variables. V22, previously classified as a type II Cepheid, appears as a bona-fide RRc in our data. In turn, V20, previously classified as an ab-type RR Lyrae, could not be properly phased with any reasonable period. The properties of the ab-type RR Lyrae stars in M28 appear most consistent with an Oosterhoff-intermediate classification, which is unusual for bona-fide Galactic globulars clusters. However, the cluster's c-type variables do not clearly support such an Oosterhoff type, and a hybrid Oosterhoff I/II system is accordingly another possibility, thus raising the intriguing possibility of multiple populations being present in M28. Coordinates, periods, and light curves in differential fluxes are provided for all the detected variables.
The sources of ultra-high energy cosmic rays (UHECRs) are unknown but are likely nearby galaxies. To assess association of UHECRs and candidate sources, we developed a multilevel Bayesian framework. We demonstrate this framework using simple models similar to those in previous studies, but our results suggest a need for more complex models; these also can be fit within our framework. We use MCMC methods to implement the approach to model data on 69 UHECRs observed by the Pierre Auger Observatory (PAO) during 2004-2009, using a volume-complete catalog of 17 nearby active galactic nuclei (AGN) as candidate sources. The reported PAO data are incomplete; an early portion (Period 1) was used to set an energy cut maximizing anisotropy, and only cosmic rays above that cut are reported. This data-tuning proves problematic for independent analyses. Assuming a common, isotropic UHECR emission rate ("standard candles") and using the untuned data after Period 1, there is no significant evidence for association of UHECRs with nearby AGNs versus an isotropic population of sources. If the association model is adopted, the fraction of UHECRs associated with these AGN is likely nonzero but well below 50%. Relatively small magnetic deflections are favored; models that assign a large fraction of UHECRs to a single nearby source are ruled out unless very large deflections are specified a priori. Including Period 1 data alters the conclusions significantly, and a simulation study suggests that the Period 1 data are anomalous, presumably due to the tuning. Accurate and optimal analysis of future data will likely require more complete disclosure of the data.
Infrared 3.6 to 8 micron images of the Extended Groth Strip yield plausible counterpart identifications for all but one of 510 radio sources in the AEGIS20 S(1.4 GHz) > 50 micro-Jy sample. This is the first such deep sample that has been effectively 100% identified. Achieving the same identification rate at R-band would require observations reaching R_AB > 27. Spectroscopic redshifts are available for 46% of the sample and photometric redshifts for an additional 47%. Almost all of the sources with 3.6 micron AB magnitudes brighter than 19 have spectroscopic redshifts z < 1.1, while fainter objects predominantly have photometric redshifts with 1 \lapprox z \lapprox 3. Unlike more powerful radio sources that are hosted by galaxies having large stellar masses within a relatively narrow range, the AEGIS20 counterparts have stellar masses spanning more than a factor of 10 at z \sim 1. The sources are roughly 10--15% starbursts at z \lapprox 0.5 and 20--25% AGNs mostly at z > 1 with the remainder of uncertain nature.
Steady accretion of a radiating gas onto a central mass point is described and compared to classic Bondi accretion. Radiation losses are essential for accretion flows to be observed. Unlike Bondi flows, radiating Bondi flows pass through a sonic point at a finite radius and become supersonic near the center. The morphology of all radiating Bondi flows is described by a single dimensionless parameter. In radiating Bondi flows the mass accretion rate varies approximately as the first power of the central mass -- this differs significantly from the quadratic dependence on the central mass in classical Bondi flows. Mass accretion rates onto galaxy or cluster-centered black holes estimated from traditional and radiating Bondi flows are significantly different. In radiating Bondi flows the gas temperature increases at large radii, as in the cores of many galaxy groups and clusters, allowing radiating Bondi flows to merge naturally with gas arriving from their cluster environments. Some radiating flows cool completely before reaching the center of the flow, and this also occurs in cooling site flows in which there is no central gravitating mass.
We calculate the baryon asymmetry produced at the electroweak phase transition by quasi-degenerate third generation sfermions in the minimal supersymmetric extension of the Standard Model. We evaluate constraints from Higgs searches, from collider searches for supersymmetric particles, and from null searches for the permanent electric dipole moment (EDM) of the electron, of the neutron and of atoms. We find that resonant sfermion sources can in principle provide a large enough baryon asymmetry in various corners of the sfermion parameter space, and we focus, in particular, on the case of large $\tan\beta$, where third-generation down-type (s)fermions become relevant. We show that in the case of stop and sbottom sources, the viable parameter space is ruled out by constraints from the non-observation of the Mercury EDM. We introduce a new class of CP violating sources, quasi-degenerate staus, that escapes current EDM constraints while providing large enough net chiral currents to achieve successful "slepton-mediated" electroweak baryogenesis.
We study the spectrum of cosmological fluctuations in scenarios such as Galilean Genesis in which a spectator scalar field acquires a scale-invariant spectrum of perturbations during an early phase which asymptotes in the far past to Minkowski space-time. In the case of minimal coupling to gravity and standard scalar field Lagrangian, the induced curvature fluctuations depend quadratically on the spectator field and are hence non-scale-invariant and highly non-Gaussian. We show that if higher dimensional operators (the same operators that lead to the {\eta}-problem for inflation) are considered, a linear coupling between background and spectator field fluctuations is induced which leads to scale-invariant and Gaussian curvature fluctuations.
In previous works, it was shown that the Lagrangians and Hamiltonians of cosmological linear scalar, vector and tensor perturbations of homogeneous and isotropic space-times with flat spatial sections containing a perfect fluid can be put in a simple form through the implementation of canonical transformations and redefinitions of the lapse function, without ever using the background classical equations of motion. In this paper, we generalize this result to general fluids, which includes entropy perturbations, and to arbitrary spacelike hyper-surfaces through a new method together with the Faddeev-Jackiw procedure for the constraint reduction. A simple second order Hamiltonian involving the Mukhanov-Sasaki variable is obtained, again without ever using the background equations of motion.
We present large scale electronic structure based molecular dynamics simulations of liquid methane at planetary conditions. In particular, we address the controversy of whether or not the interior of Uranus and Neptune consists of diamond. In our simulations we find no evidence for the formation of diamond, but rather sp2-bonded polymeric carbon. Furthermore, we predict that at high tem- perature hydrogen may exist in its monoatomic and metallic state. The implications of our finding for the planetary models of Uranus and Neptune are in detail discussed.
Links to: arXiv, form interface, find, astro-ph, recent, 1206, contact, help (Access key information)
We present Keck/NIRC2 $K_{s}$ band high-contrast coronagraphic imaging of the luminous debris disk around the nearby, young A star HD 32297 resolved at a projected separation of $r$ = 0.3-2.5\arcsec{} ($\approx$ 35-280 AU). The disk is highly warped to the north and exhibits a complex, "wavy" surface brightness profile interior to $r$ $\approx$ 110 AU, where the peaks/plateaus in the profiles are shifted between the NE and SW disk lobes. The SW side of the disk is 50--100% brighter at $r$ = 35-80 AU, and the location of its peak brightness roughly coincides with the disk's mm emission peak. Spectral energy distribution modeling suggests that HD 32297 has at least two dust populations that may originate from two separate belts likely at different locations, possibly at distances coinciding with the surface brightness peaks. A disk model for a single dust belt including a phase function with two components and a 5-10 AU pericenter offset explains the disk's warped structure and reproduces some of the surface brightness profile's shape (e.g. the overall "wavy" profile, the SB peak/plateau shifts) but more poorly reproduces the disk's brightness asymmetry. Although there may be alternate explanations, agreement between the SW disk brightness peak and disk's peak mm emission is consistent with an overdensity of very small, sub-blowout-sized dust and large, 0.1-1 mm-sized grains at $\approx$ 45 AU tracing the same parent population of planetesimals. New near-IR and submm observations may be able to clarify whether even more complex grain scattering properties or dynamical sculpting by an unseen planet are required to explain HD 32297's disk structure.
We present IRAM Plateau de Bure Interferometer observations of the 12CO(3-2) emission from two far-infrared luminous QSOs at z ~ 2.5 selected from the Herschel-ATLAS survey. These far-infrared bright QSOs were selected to have supermassive black holes (SMBHs) with masses similar to those thought to reside in sub-millimetre galaxies (SMGs) at z ~ 2.5; making them ideal candidates as systems in transition from an ultraluminous infrared galaxy phase to a sub-mm faint, unobscured, QSO. We detect 12CO(3-2) emission from both QSOs and we compare their baryonic, dynamical and SMBH masses to those of SMGs at the same epoch. We find that these far-infrared bright QSOs have similar dynamical but lower gas masses than SMGs. In particular we find that far-infrared bright QSOs have ~50+-23% less warm/dense gas than SMGs, which combined with previous results showing the QSOs lack the extended, cool reservoir of gas seen in SMGs, suggests that they are at a different evolutionary stage. This is consistent with the hypothesis that far-infrared bright QSOs represent a short (~1Myr) but ubiquitous phase in the transformation of dust obscured, gas-rich, starburst-dominated SMGs into unobscured, gas-poor, QSOs.
We predict the space density of molecular gas reservoirs in the Universe, and place a lower limit on the number counts of carbon monoxide (CO), hydrogen cyanide (HCN) molecular and [CII] atomic emission lines in blind redshift surveys in the submillimeter-centimeter spectral regime. Our model uses: (a) recently available HCN Spectral Line Energy Distributions (SLEDs) of local Luminous Infrared Galaxies (LIRGs, L_IR>10^11 L_sun), (b) a value for epsilon=SFR/M_dense(H_2) provided by new developments in the study of star formation feedback on the interstellar medium and (c) a model for the evolution of the infrared luminosity density. Minimal 'emergent' CO SLEDs from the dense gas reservoirs expected in all star-forming systems in the Universe are then computed from the HCN SLEDs since warm, HCN-bright gas will necessarily be CO-bright, with the dense star-forming gas phase setting an obvious minimum to the total molecular gas mass of any star-forming galaxy. We include [CII] as the most important of the far-infrared cooling lines. Optimal blind surveys with the Atacama Large Millimeter Array (ALMA) could potentially detect very distant (z~10-12) [CII] emitters in the >ULIRG galaxy class at a rate of ~0.1-1 per hour (although this prediction is strongly dependent on the star formation and enrichment history at this early epoch), whereas the (high-frequency) Square Kilometer Array (SKA) will be capable of blindly detecting z>3 low-J CO emitters at a rate of ~40-70 per hour. The [CII] line holds special promise for the detection of metal-poor systems with extensive reservoirs of CO-dark molecular gas where detection rates with ALMA can reach up to 2-7 per hour in Bands 4-6.
During their main sequence lifetimes, the majority of all Galactic Disc field stars must endure at least one stellar intruder passing within a few hundred AU. Mounting observations of planet-star separations near or beyond this distance suggest that these close encounters may fundamentally shape currently-observed orbital architectures and hence obscure primordial orbital features. We consider the commonly-occurring fast close encounters of two single-planet systems in the Galactic Disc, and investigate the resulting change in the planetary eccentricity and semimajor axis. We derive explicit 4-body analytical limits for these variations and present numerical cross-sections which can be applied to localized regions of the Galaxy. We find that each wide-orbit planet has a few percent chance of escape and an eccentricity that will typically change by at least 0.1 due to these encounters. The orbital properties established at formation of millions of tight-orbit Milky Way exoplanets are likely to be disrupted.
We report on the long-term dynamical evolution of the two-planet Kepler-36 system, which we studied through numerical integrations of initial conditions that are consistent with observations of the system. The orbits are chaotic with a Lyapunov time of only ~10 years. The chaos is a consequence of a particular set of orbital resonances, with the inner planet orbiting 34 times for every 29 orbits of the outer planet. The rapidity of the chaos is due to the interaction of the 29:34 resonance with the nearby first order 6:7 resonance, in contrast to the usual case in which secular terms in the Hamiltonian play a dominant role. Only one contiguous region of phase space, accounting for ~4.5% of the sample of initial conditions studied, corresponds to planetary orbits that do not show large scale orbital instabilities on the timescale of our integrations (~200 million years). The long-lived subset of the allowed initial conditions are those that satisfy the Hill stability criterion by the largest margin. Any successful theory for the formation of this system will need to account for why its current state is so close to unstable regions of phase space.
We study the properties of low-velocity material in the line of sight towards
nearby Type Ia Supernovae (SNe Ia) that have measured late phase nebular
velocity shifts (v_neb), thought to be an environment-independent observable.
We have found that the distribution of equivalent widths of narrow blended Na I
D1 & D2 and Ca II H & K absorption lines differs significantly between those
SNe Ia with negative and positive v_neb, with generally stronger absorption for
SNe Ia with v_neb > 0. A similar result had been found previously for the
distribution of colors of SNe Ia, which was interpreted as a dependence of the
temperature of the ejecta with viewing angle. Our work suggests that: 1) a
significant part of these differences in color should be attributed to
extinction, 2) this extinction is caused by an asymmetric distribution of
circumstellar material (CSM) and 3) the CSM absorption is generally stronger on
the side of the ejecta opposite to where the ignition occurs.
Since it is difficult to explain 3) via any known physical processes that
occur before explosion, we argue that the asymmetry of the CSM is originated
after explosion by a stronger ionizing flux on the side of the ejecta where
ignition occurs, probably due to a stronger shock breakout and/or more exposed
radioactive material on one side of the ejecta. This result has important
implications for both progenitor and explosion models.
Recent observational results for magnetic fields in molecular clouds reviewed by Crutcher (2012) seem to be inconsistent with the predictions of the ambipolar diffusion theory of star formation. These include the measured decrease in mass to flux ratio between envelopes and cores, the failure to detect any self-gravitating magnetically subcritical clouds, the determination of the flat PDF of the total magnetic field strengths implying that there are many clouds with very weak magnetic fields, and the observed scaling $B \propto \rho^{2/3}$ that implies gravitational contraction with weak magnetic fields. We consider the problem of magnetic field evolution in turbulent molecular clouds and discuss the process of magnetic field diffusion mediated by magnetic reconnection. For this process that we termed "reconnection diffusion" we provide a simple physical model and explain that this process is inevitable in view of the present day understanding of MHD turbulence. We address the issue of the expected magnetization of cores and envelopes in the process of star formation and show that reconnection diffusion provides an efficient removal of magnetic flux that depends only on the properties of MHD turbulence in the core and the envelope. As a result, the magnetic flux trapped during the collapse in the envelope is being released faster than the flux trapped in the core, resulting in much weaker fields in envelopes than in cores, as observed. We provide simple semi-analytical model calculations which support this conclusion and qualitatively agree with the observational results. We argue that magnetic reconnection provides a solution to the magnetic flux problem of star formation that agrees better with observations than the long-standing ambipolar diffusion paradigm.
We make use of our 'minimal' cold interstellar medium (ISM) emission line model that predicts the molecular and atomic line emission per unit dense, star-forming gas mass (Geach & Papadopoulos 2012; Paper I) to examine the utility of key line ratios in surveys of the so-called star formation 'mode' as traced by xi_SF = M_dense(H_2)/M_total(H_2). We argue that xi_SF and its proxies provide very sensitive, extinction-free discriminators of rapid starburst/merger-driven versus secular quiescent/disk-like stellar mass assembly, with the most promising diagnostic to be applied in the near-future being CO(4-3)/[CI](1-0). These lines are accessible across nearly the full range 0<z<2 (thus covering the bulk of galaxy evolution) with the Atacama Large Millimeter Array. In addition to their diagnostic power, another advantage of this combination is the similar observed frequencies (Delta nu_0 ~ 30 GHz) of the lines, resulting in nearly spatially-matched beams for a fixed aperture, thus mitigating the effects of resolution/morphology bias in the interpretation of galaxy-averaged line ratios. Finally we discuss the capability of deep blind redshift surveys with the high frequency component of the Square Kilometer Arrray (SKA) in discovering H_2-rich galaxies with very low xi_SF values. These could be the progenitors of starburst galaxies seen prior to the onset of star formation; such galaxies could be a class of extreme outliers from local (gas surface density)-(star formation rate) scaling laws, which would exclude them from current star formatation or stellar mass selected samples. Our conservative model suggests that SKA could detect such systems residing at z~3 at a rate of 20-200 per hour.
We have developed and characterized an imaging instrument to measure the spatial properties of the diffuse near-infrared extragalactic background light in a search for fluctuations from z > 6 galaxies during the epoch of reionization. The instrument is part of the Cosmic Infrared Background Experiment (CIBER), designed to observe the extragalactic background light above the Earth's atmosphere during a suborbital sounding rocket flight. The imaging instrument incorporates a 2x2 degree field of view, to measure fluctuations over the predicted peak of the spatial power spectrum at 10 arcminutes, and 7"x7" pixels, to remove lower redshift galaxies to a depth sufficient to reduce the low-redshift galaxy clustering foreground below instrumental sensitivity. The imaging instrument employs two cameras with \Delta \lambda / \lambda ~0.5 bandpasses centered at 1.1 and 1.6 microns to spectrally discriminate reionization extragalactic background fluctuations from local foreground fluctuations. CIBER operates at wavelengths where the electromagnetic spectrum of the reionization extragalactic background is thought to peak, and complements fluctuations measurements by AKARI and Spitzer at longer wavelengths. We have characterized the instrument in the laboratory, including measurements of the sensitivity, flat-field response, stray light performance, and noise properties. Several modifications were made to the instrument following a first flight in 2009 February. The instrument performed to specifications in subsequent flights in 2010 July and 2012 March, and the scientific data are now being analyzed.
Magnetohydrodynamic turbulence is able to create hierarchical structures in the interstellar medium that are correlated on a wide range of scales via the energy cascade. We use hierarchical tree diagrams known as dendrograms to characterize structures in synthetic Position-Position-Velocity (PPV) emission cubes of optically thin isothermal magnetohydrodynamic turbulence. We show that the structures and degree of hierarchy observed in PPV space are related to the physics of the gas, i.e. self-gravity and the global sonic and Alfvenic Mach number. Simulations with higher Alfvenic Mach number, self-gravity and supersonic flows display enhanced hierarchical structure. We observed a strong sonic and Alfvenic dependency when we apply the the statistical moments (i.e. mean, variance, skewness, kurtosis) to the dendrogram distribution. Larger magnetic field and sonic Mach number correspond to larger values of the moments. Application of the dendrogram to 3D density cubes, also known as Position-Position-Position cubes (PPP), reveals that the dominant emission contours in PPP and PPV are related for supersonic gas but not for subsonic. We also explore the effects of smoothing, thermal broadening and velocity resolution on the dendrograms in order to make our study more applicable to observational data. These results all point to hierarchical tree diagrams as being a promising additional tool for studying ISM turbulence and star forming regions in the direction of obtaining information on the degree of self-gravity, the Mach numbers and the complicated relationship between PPV and PPP.
A fundamental problem of cosmic ray (CR) physics is the determination of the average properties of Galactic CRs outside the Solar system. Starting from COS-B data in the 1980's, gamma-ray observations of molecular clouds in the Gloud Belt above the Galactic plane have been used to deduce the Galactic CR energy spectrum. We reconsider this problem in view of the improved precision of observational data which in turn require a more precise treatment of photon production in proton-proton scatterings. We show that the spectral shape $dN/dp\propto p^{-2.85}$ of CR protons as determined by the PAMELA collaboration in the energy range 80 GeV<pc<230 GeV is consistent with the photon spectra from molecular clouds observed with Fermi-LAT down to photon energies E\sim 1-2 GeV. Adding a break of the CR flux at 3 GeV, caused by a corresponding change of the diffusion coefficient, improves further the agreement in the energy range 0.2-3 GeV.
In the effort to characterize the masses, radii, and atmospheres of potentially habitable exoplanets, there is an urgent need to find examples of such planets transiting nearby M dwarfs. The MEarth Project is an ongoing effort to do so, as a ground-based photometric survey designed to detect exoplanets as small as 2 Earth radii transiting mid-to-late M dwarfs within 33 pc of the Sun. Unfortunately, identifying transits of such planets in photometric monitoring is complicated both by the intrinsic stellar variability that is common among these stars and by the nocturnal cadence, atmospheric variations, and instrumental systematics that often plague Earth-bound observatories. Here we summarize the challenges MEarth faces, and address them with a new framework to detect shallow exoplanet transits in wiggly and irregularly-spaced light curves. In contrast to previous methods that clean trends from light curves before searching for transits, this framework assesses the significance of individual transits simultaneously while modeling variability, systematics, and the photometric quality of individual nights. Our Method for Including Starspots and Systematics in the Marginalized Probability of a Lone Eclipse (MISS MarPLE) uses a computationally efficient semi-Bayesian approach to explore the vast probability space spanned by the many parameters of this model, naturally incorporating the uncertainties in these parameters into its evaluation of candidate events. We show how to combine individual transits processed by MISS MarPLE into periodic transiting planet candidates and compare our results to the popular Box-fitting Least Squares (BLS) method with simulations. By applying MISS MarPLE to observations from the MEarth Project, we demonstrate the utility of this framework for robustly assessing the false alarm probability of transit signals in real data.
In the Solar system the planets' compositions vary with orbital distance, with rocky planets in close orbits and lower-density gas giants in wider orbits. The detection of close-in giant planets around other stars was the first clue that this pattern is not universal, and that planets' orbits can change substantially after their formation. Here we report another violation of the orbit-composition pattern: two planets orbiting the same star with orbital distances differing by only 10%, and densities differing by a factor of 8. One planet is likely a rocky `super-Earth', whereas the other is more akin to Neptune. These planets are thirty times more closely spaced--and have a larger density contrast--than any adjacent pair of planets in the Solar system.
Some Galactic models predict a significant population of radio pulsars close to the our galactic center. Beams from these pulsars could get strongly deflected by the supermassive black hole (SMBH) believed to reside at the galactic center and reach the Earth. Earlier work assuming a Schwarzschild SMBH gave marginal chances of observing this exotic phenomenon with current telescopes and good chances with future telescopes. Here we calculate the odds of observability for a rotating SMBH. We find that the estimates of observation are not affected by the SMBH spin, but a pulsar timing analysis of deflected pulses might be able to provide an estimate of the spin of the central black hole.
We study the dependence of the delayed neutrino-heating mechanism for core-collapse supernovae on the equation of state. Using a simplified treatment of the neutrino physics with a parameterized neutrino luminosity, we explore the relationship between explosion time, mass accretion rate, and neutrino luminosity for a 15 Msun progenitor in 1D and 2D. We test three different equations of state commonly used in core-collapse simulations: the models of Lattimer & Swesty (1991) with incompressibility of 180 MeV and 220 MeV, and the model of Shen et al. (1998), in order of increasing stiffness. We find that for a given neutrino luminosity the time after bounce until explosion increases with the stiffness of the equation of state: the Lattimer & Swesty EOS explode more easily than that of Shen et al. We find this holds in both 1D and 2D, while for all models explosions are obtained more easily in 2D than in 1D. We also discuss the relevance of approximate instability criteria to realistic simulations.
We present 230 GHz Submillimeter Array continuum and molecular line observations of the newly discovered FUor candidate HBC722. We report the detection of seven 1.3 mm continuum sources in the vicinity of HBC722, none of which correspond to HBC722 itself. We compile infrared and submillimeter continuum photometry of each source from previous studies and conclude that three are Class 0 embedded protostars, one is a Class I embedded protostar, one is a Class I/II transition object, and two are either starless cores or very young, very low luminosity protostars or first hydrostatic cores. We detect a northwest-southeast outflow, consistent with the previous detection of such an outflow in low-resolution, single-dish observations, and note that its axis may be precessing. We show that this outflow is centered on and driven by one of the nearby Class 0 sources rather than HBC722, and find no conclusive evidence that HBC722 itself is driving an outflow. The non-detection of HBC722 in the 1.3 mm continuum observations suggests an upper limit of 0.02 solar masses for the mass of the circumstellar disk. This limit is consistent with typical T Tauri disks and with a disk that provides sufficient mass to power the burst.
We study the relation between the measured anisotropies in the extragalactic diffuse gamma-ray background (DGRB) and the DGRB spectral intensity, and their potential origin from the unresolved blazar population. Using a physical-evolution model for blazars with a luminosity dependent density evolution (LDDE) and an observationally-determined luminosity-dependent blazar spectral energy distribution (SED), we find that blazars can account for the observed anisotropy of the DGRB consistent with their observed source-count distribution, but are in turn constrained in contributing significantly to the observed DGRB intensity. For the best-fit LDDE model accounting for the DGRB anisotropy and source-count distribution, blazars only contribute 4.3^{+4.1}_{-1.1}% (68% CL) of the DGRB intensity above 1 GeV. Requiring a higher fraction of the DGRB intensity contribution by blazars overproduces the DGRB anisotropy, and therefore blazars in the LDDE+SED-sequence model cannot simultaneously account for the DGRB intensity as well as anisotropy. We discuss the limitations of LDDE models. However, these models do not require the many unjustified and observationally-inconsistent simplifying assumptions---including a single power law for all blazar spectra and a simple broken power-law model for their source-count distribution---that are present in much previous work.
A significant fraction of the total photospheric light in nearby galaxy clusters is thought to be contained within the diffuse intracluster light (ICL), which extends 100s of kpc from cluster cores. The study of the ICL can reveal details of the evolutionary histories and processes occurring within galaxy clusters, however since it has a very low surface brightness it is often difficult to detect. We present here the first measurements of the ICL as a fraction of total cluster light at z \sim 1 using deep J-band (1.2 {\mu}m) imaging from HAWK-I on the VLT. We investigate the ICL in 6 X-ray selected galaxy clusters at 0.8< z <1.2 and find that the ICL below isophotes {\mu}(J) = 22 mag/arcsec2 constitutes 1-4% of the total cluster light within a radius R500. This is broadly consistent with simulations of the ICL at a similar redshift and when compared to nearby observations suggests that the fraction of the total cluster light that is in the ICL has increased by a factor 2 - 4 since z\sim1. We also find the fraction of the total cluster light contained within the Brightest Cluster Galaxy (BCG) to be 2.0-6.3% at these redshifts, which in 5 out of 6 cases is larger than the fraction of the ICL component, in contrast to results from nearby clusters. This suggests that the evolution in cluster cores involves substantial stripping activity at late times, in addition to the early build up of the BCG stellar mass through merging. The presence of significant amounts of stellar light at large radii from these BCGs may help towards solving the recent disagreement between the semi-analytic model predictions of BCG mass growth (e.g. De Lucia & Blaziot, 2007) and the observed large masses and scale sizes reported for BCGs at high redshift.
In a recent computational campaign [Ng et al., Astrophys. J. 747, 109, 2012] to investigate a three-dimensional model of coronal heating using reduced magnetohydrodynamics (RMHD), we have obtained scaling results of heating rate versus Lundquist number based on a series of runs in which random photospheric motions are imposed for hundreds to thousands of Alfv\'en time in order to obtain converged statistical values. Using this collection of numerical data, we have performed additional statistical analysis related to the formation of current sheets and heating events, or nanoflares [Parker, Astrophys. J. 330, 474, 1988]. While there have been many observations of the energy distribution of solar flares, there have not been many results based on large-scale three-dimensional direct simulations due to obvious numerical difficulties. We will present energy distributions and other statistics based on our simulations, calculated using a method employed in [Dmitruk & G\'omez, Astrophys. J., 484, L83, 1997]. We will also make comparisons of our results with observations.
Recently, a tentative 130 GeV $\gamma$-ray line signal was identified by quite a few groups. If correct it would constitute a "smoking gun" for dark matter annihilations. Interestingly, the spectra of the sum of cosmic ray electrons and positrons detected by ATIC and PAMELA both show small wiggle-like structure at $\sim 100$ GeV, which could be the result of the annihilation of $\sim 140$ GeV dark matter particles into electrons/positrons with a velocity-weighted cross section $<\sigma v>_{\rm \chi\chi\rightarrow e^{+}e^{-}}\sim 10^{-26}-10^{-25} {\rm cm^{3} s^{-1}}$. Accurate measurements of the total spectrum of electron and positron cosmic rays by AMS-2 and the upcoming missions such as DAMPE and CALET are highly needed to pin down the profile of the wiggle-like structure and then its physical origin.
We present the analysis of supernova remnants (SNRs) and candidates in M31 identified in the XMM-Newton large programme survey of M31. SNRs are among the bright X-ray sources in a galaxy. They are good indicators of recent star formation activities of a galaxy and of the interstellar environment in which they evolve. By combining the X-ray data of sources in M31 with optical data as well as with optical and radio catalogues, we aim to compile a complete, revised list of SNRs emitting X-rays in M31 detected with XMM-Newton, study their luminosity and spatial distribution, and understand the X-ray spectrum of the brightest SNRs. We analysed the X-ray spectra of the twelve brightest SNRs and candidates using XMM-Newton data. The four brightest sources allowed us to perform a more detailed spectral analysis and the comparison of different models to describe their spectrum. For all M31 large programme sources we searched for optical counterparts on the Ha, [Sii], and [Oiii] images of the Local Group Galaxy Survey. We confirm 21 X-ray sources as counterparts of known SNRs. In addition, we identify five new X-ray sources as X-ray and optically emitting SNRs. Seventeen sources are no longer considered as SNR candidates. We have thus created a list of 26 X-ray SNRs and 20 candidates in M31 based on their X-ray, optical, and radio emission, which is the most recent complete list of X-ray SNRs in M31. The brightest SNRs have X-ray luminosities of up to 8 x 10^36 erg/s in the 0.35 - 2.0 keV band.
We present the results of a comprehensive analysis of the structure and kinematics of six Galactic globular clusters. By comparing the results of the most extensive photometric and kinematical surveys available to date with suited dynamical models, we determine the stellar and dynamical masses of these stellar systems taking into account for the effect of mass segregation, anisotropy and unresolved binaries. We show that the stellar masses of these clusters are on average smaller than those predicted by canonical integrated stellar evolution models because of the shallower slope of their mass functions. The derived stellar masses are found to be also systematically smaller than the dynamical masses by ~40%, although the presence of systematics affecting our estimates cannot be excluded. If confirmed, this evidence can be linked to an increased fraction of retained dark remnants or to the presence of a modest amount of dark matter.
Interstellar filaments are an important part of star formation. To understand the structure of filaments, cross-section profiles are often fitted with Plummer profiles. This profiling is often done with submm studies, such as Herschel. It would be convenient if filament properties could also be studied using groundbased NIR data. We compare the filament profiles obtained by NIR extinction and submm observations to find out if reliable profiles can be derived using NIR data. We use J-, H-, and K-band data of a filament north of TMC-1 to derive an extinction map from colour excesses of background stars. We compare the Plummer profiles obtained from extinction maps with Herschel dust emission maps. We present 2 methods to estimate profiles from NIR: Plummer profile fits to median Av of stars or directly to the Av of individual stars. We compare the methods by simulations. In simulations extinction maps and the new methods give correct results to within ~10-20 for modest densities. Direct fit to data on individual stars gives more accurate results than extinction map, and can work in higher density. In profile fits to real observations, values of Plummer parameters are generally similar to within a factor of ~2. Although parameter values can vary significantly, estimates of filament mass usually remain accurate to within some tens of per cent. Our results for TMC-1 are in agreement with earlier results. High resolution NIR data give more details, but 2MASS data can be used to estimate profiles. NIR extinction can be used as an alternative to submm observations to profile filaments. Direct fits of stars can also be a valuable tool. Plummer profile parameters are not always well constrained, and caution should be taken when making fits. In the evaluation of Plummer parameters, one can use the independence of dust emission and NIR data and the difference in the shapes of the confidence regions.
J1211+743 is a giant radio galaxy with a one-sided jet and two asymmetric lobes, one of which is of Fanaroff-Riley (FR) type II with a hotspot and the other is a diffuse relic devoid of a hotspot. The jet points towards the latter lobe, which is difficult to explain in a standard way within the double-lobed radio source paradigm. Here, I propose to assume that the nucleus of J1211+743 has undergone a re-ignition of activity and its lobes, presumably both originally of FR II type, represent an earlier active phase, while the jet represents the current one. The asymmetry of the lobes is a consequence of the orientation of the source combined with an activity switch-off that occurred between two active periods. The relic lobe is on the near side with regard to the observer, whereas the radiation from the far-side lobe arrives significantly later owing to its longer distance to the observer. The far-side lobe is thus perceived to have not yet decayed. On the other hand, the jet behaves in a standard way, i.e. its projected orientation reflects the near side of the source. Hence, we are able to explain why the location of the relic lobe correlates with the direction of the jet.
Trapping of two-armed ($m=2$) vertical p-mode oscillations in relativistic disks are examined. The disks are assumed to be isothermal in the vertical direction, but are truncated at a certain height by the presence of corona. The same issues have been examined in a previous paper (Kato 2012a). In this paper, unlike the previous paper, however, we do not use the approximation that the oscillations are nearly vertical, but limit to a simpler case of no magnetic field. As in the previous paper, the results suggest that the two basic oscillation modes [both are the fundamental ($n=1$) in the vertical direction but in the horizontal direction one is the fundamental ($n_{\rm r}=0$) and the other the first overtone ($n_{\rm r}=1$)] correspond to the twin kHz QPOs. Second, the oscillation mode which is the first overtone $(n=2)$ in the vertical direction and the fundamental in the horizontal direction ($n_{\rm r}=0$) will correspond to the horizontal branch oscillation (HBO) of Z-sources. The latter suggests that the horizontal branch of Z-sources is a sequence of temperature change in disks whose vertical thickness is strongly terminated. The temperature increases leftward along the sequence from the apex between normal and horizontal branches.
We use high (0."65 x 0."52,(65x52pc)) resolution SubMillimeter Array (SMA)
observations to image the CO and 13CO 2-1 line emission of the extreme
FIR-excess galaxy NGC 1377. We find bright, complex CO 2-1 line emission in the
inner 400 pc of the galaxy. The CO 2-1 line has wings that are tracing a
kinematical component which appears perpendicular to that of the line core.
Together with an intriguing X-shape of the integrated intensity and dispersion
maps, this suggests that the molecular emission of NGC 1377 consists of a
disk-outflow system. Lower limits to the molecular mass and outflow rate are
M_out(H2)>1e7 Msun and dM/dt>8 Msun/yr. The age of the proposed outflow is
estimated to 1.4Myrs, the extent to 200pc and the outflow speed to 140 km/s.
The total molecular mass in the SMA map is estimated to M_tot(H2)=1.5e8 Msun
(on a scale of 400 pc) while in the inner r=29 pc the molecular mass is
M_core(H2)=1.7e7 Msun with a corresponding H2 column density of N(H2)=3.4e23
cm-2 and an average CO 2-1 brightness temperature of 19K.
Observing the molecular properties of the FIR-excess galaxy NGC 1377 allows
us to probe the early stages of nuclear activity and the onset of feedback in
active galaxies. The age of the outflow supports the notion that the current
nuclear activity is young - a few Myrs. The outflow may be powered by radiation
pressure from a compact, dust enshrouded nucleus, but other driving mechanisms
are possible. The buried source may be an AGN or an extremely young (1Myr)
compact starburst. Limitations on size and mass lead us to favour the AGN
scenario, but further studies are required to settle the issue. In either case,
the wind with its implied mass outflow rate will quench the nuclear power
source within a very short time of 5-25 Myrs. It is however possible that the
gas is unable to escape the galaxy and may eventually fall back onto NGC 1377
again.
In the last years several exoplanets have been discovered that orbit one component of a compact binary system (separation < 50 astronomical units), the probably best-known case is gamma-Cephei. So far, all attempts to explain the in-situ formation of these planets has been unsuccessful, in part because of the strong gravitational perturbations of the secondary star on any initial planetesimal swarm. Here we test whether planetary bodies in compact binaries, in particular gamma-Cephei, could have originated from a close encounter with a passing star, assuming initial configurations for the stellar system suitable for planetary formation. In other words, we analyze whether the orbital configuration of the current binary system might have been generated after the formation of the planet, and as a consequence of a close encounter with a third star in hyperbolic orbit. We performed a series of time-reverse N-body simulations of stellar scattering events in which the present-day configuration of gamma-Cephei was used as the initial condition plus a hypothetical third star as an impactor. We analyzed which configurations and system parameters could have given birth to the current system. Depending on the maximum impact velocity allowed for accretional collisions, we find that between 1% and 5% of stellar encounters correspond to an "original" system in which planetary formation around the primary star is not inhibited by the secondary, but is acceptable within the classical core-accretion scenario. Thus, although not highly probable, it is plausible that stellar encounters may have played a significant role in shaping these types of exoplanetary systems.
On 2010 October 13, the X-ray astronomical satellite Rossi XTE, during the observation of the newly discovered accretion powered X-ray pulsar IGR J17480--2446, detected a lunar occultation of the source. From knowledge of lunar topography and Earth, Moon, and spacecraft ephemeris at the epoch of the event, we determined the source position with an accuracy of 40 mas (1{\sigma} c.l.), which is interesting, given the very poor imaging capabilities of RXTE (\sim 1\circ). For the first time, using a non-imaging X-ray observatory, the position of an X-ray source with a subarcsecond accuracy is derived, demonstrating the neat capabilities of a technique that can be fruitfully applied to current and future X-ray missions.
Backreaction effects of the large scale structure on the background dynamics have been claimed to lead to a renormalization of the background dynamics that may account for the late time acceleration of the cosmic expansion. This article emphasizes that generically the averaged flow is locally anisotropic, a property that can be related to observation. Focusing on perturbation theory, the spatially averaged shear, that characterizes the anisotropy of the flow, is computed. It is shown that this shear arising from backreaction differs from a homogeneous shear: its time evolution is different and its amplitude is completely determined by the cosmological parameters and the matter power spectrum. It ranges within (2-37)% at a redshift of order 0.5 so that the isotropy of the Hubble flow may allow to constrain the backreaction approach to dark energy.
The most powerful explosions on the Sun [...] drive the most severe space-weather storms. Proxy records of flare energies based on SEPs in principle may offer the longest time base to study infrequent large events. We conclude that one suggested proxy, nitrate concentrations in polar ice cores, does not map reliably to SEP events. Concentrations of select radionuclides measured in natural archives may prove useful in extending the time interval of direct observations up to ten millennia, but as their calibration to solar flare fluences depends on multiple poorly known properties and processes, these proxies cannot presently be used to help determine the flare energy frequency distribution. Being thus limited to the use of direct flare observations, we evaluate the probabilities of large-energy solar explosions by combining solar flare observations with an ensemble of stellar flare observations. We conclude that solar flare energies form a relatively smooth distribution from small events to large flares, while flares on magnetically-active, young Sun-like stars have energies and frequencies markedly in excess of strong solar flares, even after an empirical scaling with the mean activity level of these stars. In order to empirically quantify the frequency of uncommonly large solar flares extensive surveys of stars of near-solar age need to be obtained, such as is feasible with the Kepler satellite. Because the likelihood of flares larger than approximately X30 remains empirically unconstrained, we present indirect arguments, based on records of sunspots and on statistical arguments, that solar flares in the past four centuries have likely not substantially exceeded the level of the largest flares observed in the space era, and that there is at most about a 10% chance of a flare larger than about X30 in the next 30 years.
The presence of a dark matter core in the central kiloparsec of many dwarf galaxies has been a long standing problem in galaxy formation theories based on the standard cold dark matter paradigm. Recent cosmological simulations, based on Smooth Particle Hydrodynamics and rather strong feedback recipes have shown that it was indeed possible to form extended dark matter cores using baryonic processes related to a more realistic treatment of the interstellar medium. Using adaptive mesh refinement, together with a new, stronger supernovae feedback scheme that we have recently implemented in the RAMSES code, we show that it is also possible to form a prominent dark matter core within the well-controlled framework of an isolated, initially cuspy, 10 billion solar masses dark matter halo. Although our numerical experiment is idealized, it allows a clean and unambiguous identification of the dark matter core formation process. Our dark matter inner profile is well fitted by a pseudo-isothermal profile with a core radius of 800 pc. The core formation mechanism is consistent with the one proposed by Pontzen & Governato (2012). We highlight two key observational predictions of all simulations that find cusp-core transformations: (i) a bursty star formation history with peak to trough ratio of 5 to 10 and a duty cycle comparable to the local dynamical time; and (ii) a stellar distribution that is hot with v/sigma=1. We compare the observational properties of our model galaxy with recent measurements of the isolated dwarf WLM. We show that the spatial and kinematical distribution of stars and HI gas are in striking agreement with observations, supporting the fundamental role played by stellar feedback in shaping both the stellar and dark matter distribution.
We analyse the results of four simulations of isolated galaxies: two with a rigid spiral potential of fixed pattern speed, but with different degrees of star-formation induced feedback, one with an axisymmetric galactic potential and one with a `live' self-gravitating stellar component. Since we use a Lagrangian method we are able to select gas that lies within giant molecular clouds (GMCs) at a particular timeframe, and to then study the properties of this gas at earlier and later times. We find that gas which forms GMCs is not typical of the interstellar medium at least 50 Myr before the clouds form and reaches mean densities within an order of magnitude of mean cloud densities by around 10 Myr before. The gas in GMCs takes at least 50 Myr to return to typical ISM gas after dispersal by stellar feedback, and in some cases the gas is never fully recycled. We also present a study of the two-dimensional, vertically-averaged velocity fields within the ISM. We show that the velocity fields corresponding to the shortest timescales (that is, those timescales closest to the immediate formation and dissipation of the clouds) can be readily understood in terms of the various cloud formation and dissipation mechanisms. Properties of the flow patterns can be used to distinguish the processes which drive converging flows (e.g.\ spiral shocks, supernovae) and thus molecular cloud formation, and we note that such properties may be detectable with future observations of nearby galaxies.
Many chemical models of dense interstellar clouds predict that the majority of gas-phase elemental nitrogen should be present as N2, with an abundance approximately five orders of magnitude less than that of hydrogen. As a homonuclear diatomic molecule, N2 is difficult to detect spectroscopically through infrared or millimetre-wavelength transitions so its abundance is often inferred indirectly through its reaction product N2H+. Two main formation mechanisms each involving two radical-radical reactions are the source of N2 in such environments. Here we report measurements of the low temperature rate constants for one of these processes, the N + CN reaction down to 56 K. The effect of the measured rate constants for this reaction and those recently determined for two other reactions implicated in N2 formation are tested using a gas-grain model employing a critically evaluated chemical network. We show that the amount of interstellar nitrogen present as N2 depends on the competition between its gas-phase formation and the depletion of atomic nitrogen onto grains. As the reactions controlling N2 formation are inefficient, we argue that N2 does not represent the main reservoir species for interstellar nitrogen. Instead, elevated abundances of more labile forms of nitrogen such as NH3 should be present on interstellar ices, promoting the eventual formation of nitrogen-bearing organic molecules.
In early October 2008, the Soft Gamma Repeater SGRJ1550 - 5418 (1E 1547.0 - 5408, AXJ155052 - 5418, PSR J1550 - 5418) became active, emitting a series of bursts which triggered the Fermi Gamma-ray Burst Monitor (GBM) after which a second especially intense activity period commenced in 2009 January and a third, less active period was detected in 2009 March-April. Here we analyze the GBM data all the bursts from the first and last active episodes. We performed temporal and spectral analysis for all events and found that their temporal characteristics are very similar to the ones of other SGR bursts, as well the ones reported for the bursts of the main episode (average burst durations \sim 170 ms). In addition, we used our sample of bursts to quantify the systematic uncertainties of the GBM location algorithm for soft gamma-ray transients to < 8 deg. Our spectral analysis indicates significant spectral evolution between the first and last set of events. Although the 2008 October events are best fit with a single blackbody function, for the 2009 bursts an Optically Thin Thermal Bremsstrahlung (OTTB) is clearly preferred. We attribute this evolution to changes in the magnetic field topology of the source, possibly due to effects following the very energetic main bursting episode.
The identification of strong gravitational lenses in large surveys has historically been a rather time consuming exercise. Early data from the Herschel Astrophysical Terahertz Large Area Survey (Herschel-ATLAS) demonstrate that lenses can be identified efficiently at submillimetre wavelengths using a simple flux criteria. Motivated by that development, this work considers the statistical properties of strong gravitational lens systems which have been, and will be, found by the Herschel-ATLAS. Analytical models of lens statistics are tested with the current best estimates for the various model ingredients. These include the cosmological parameters, the mass function and the lens density profile, for which we consider the singular isothermal sphere (SIS) and the Navarro, Frenk & White (NFW) approximations. The five lenses identified in the Herschel-ATLAS Science Demonstration Phase suggest a SIS density profile is preferred, but cannot yet constrain \Omega_\Lambda to an accuracy comparable with other methods. The complete Herschel-ATLAS data set should be sufficient to yield competitive constraints on \Omega_\Lambda. Whilst this huge number of lenses has great potential for constraining cosmological parameters, they will be most powerful in constraining uncertainty in astrophysical processes. Further investigation is needed to fully exploit this unprecedented data set.
We present an empirical relations for correcting the estimated masses, effective temperatures, and radii of chromospherically active low-mass stars and brown dwarfs. We base our corrections on a large set of low-mass stars in the field with Halpha activity measurements, and on a set of low-mass eclipsing binaries with X-ray activity measurements from which we indirectly infer the Halpha activity. Both samples yield consistent relations linking the amount by which an active object's temperature is suppressed, and the amount by which its radius is inflated, to the strength of its Halpha emission. Bolometric luminosity is found to be approximately preserved by these temperature and radius corrections. We apply these relations to the peculiar brown-dwarf eclipsing binary 2M0535-05, in which the active, higher-mass brown dwarf has a cooler temperature than its inactive, lower-mass companion. We find that the Halpha-corrected temperatures bring the inferred masses of the brown dwarfs into agreement with theoretical isochrones. These empirical relations are applicable to brown dwarfs and low-mass stars with masses below 0.8 Msun and for which the activity, as measured by the fractional Halpha luminosity, is in the range -4.6 < log LHa/Lbol < -3.3. We expect these corrections to be most useful for improving temperatures and radii of low-mass stars and brown dwarfs over their active lifetimes (few Gyr) and thereby also the inferred masses of objects with unknown ages or distances (and therefore unknown luminosities). We also discuss the implications of this work for improved determinations of young cluster initial mass functions.
The Cosmic Infrared Background (CIB) provides an opportunity to constrain many properties of the high redshift (z>6) stellar population as a whole. This background, specifically, from 1 to 200 microns, will contain any information about the era of reionization and the stars responsible for producing these ionizing photons. In this paper, we look at the fractional anisotropy delta I/I of this high redshift population, which is the ratio of the magnitude of the fluctuations (delta I) and the mean intensity (I). We show that this can be used to constrain the escape fraction of the population as a whole. The magnitude of the fluctuations of the CIB depend on the escape fraction, while the mean intensity does not. This results in lower values of the escape fraction producing higher values of the fractional anisotropy. This difference is predicted to be larger at the longer wavelengths bands (above 10 microns), albeit it is also much harder to observe in that range. We show that the fractional anisotropy can also be used to separate a dusty from a dust-free population. Finally, we discuss the constraints provided by current observations on the CIB fractional anisotropy.
We report the discovery of a unique object, BD+48 740, a lithium overabundant giant with A(Li)=2.33 +/- 0.04 (where A(Li) = log(n_Li/n_H) + 12), that exhibits radial velocity (RV) variations consistent with a 1.6 M_J companion in a highly eccentric, e = 0.67 +/- 0.17 and extended, a=1.89 AU (P=771 d), orbit. The high eccentricity of the planet is uncommon among planetary systems orbiting evolved stars and so is the high lithium abundance in a giant star. The ingestion by the star of a putative second planet in the system originally in a closer orbit, could possibly allow for a single explanation to these two exceptional facts. If the planet candidate is confirmed by future RV observations, it might represent the first example of the remnant of a multiple planetary system possibly affected by stellar evolution.
We propose a mechanism for the observed non-keplerian motion (Tiscareno et al. 2010, ApJL) of "propeller" moons embedded in Saturn's rings. Our mechanism, in which radial variations in surface density -- external to, and unaffected by, the embedded moon -- result in an equilibrium semimajor axis for the moon due to "Type I" angular momentum exchange (Crida et al. 2010, AJ), provides a good fit to the observations. Future observations should distinguish between our model and others recently proposed.
We study the redshift distribution of two samples of early-type gravitational lenses, extracted from a larger collection of 122 systems, to constrain the cosmological constant in the LCDM model and the parameters of a set of alternative dark energy models (XCDM, Dvali-Gabadadze-Porrati and Ricci dark energy models), under a spatially flat universe. The likelihood is maximized for $\Omega_\Lambda= 0.70 \pm 0.09$ when considering the sample excluding the SLACS systems (known to be biased towards large image-separation lenses) and no-evolution, and $\Omega_\Lambda= 0.81\pm 0.05$ when limiting to gravitational lenses with image separation larger than 2" and no-evolution. In both cases, results accounting for galaxy evolution are consistent within 1$\sigma$. The present test supports the accelerated expansion, by excluding the null-hypothesis (i.e., $\Omega_\Lambda = 0 $) at more than 4$\sigma$, regardless of the chosen sample and assumptions on the galaxy evolution. A comparison between competitive world models is performed by means of the Bayesian information criterion. This shows that the simplest cosmological constant model - that has only one free parameter - is still preferred by the available data on the redshift distribution of gravitational lenses. We perform an analysis of the possible systematic effects, finding that the systematic errors due to sample incompleteness, galaxy evolution and model uncertainties approximately equal the statistical errors, with present-day data. We find that the largest sources of systemic errors are the dynamical normalization and the high-velocity cut-off factor, followed by the faint-end slope of the velocity dispersion function.
The frequency and properties of multiple star systems offer powerful tests of star formation models. Multiplicity surveys over the past decade have shown that binary properties vary strongly with mass, but the functional forms and the interplay between frequency and semimajor axis remain largely unconstrained. We present the results of a large-scale survey of multiplicity at the bottom of the IMF in several nearby young associations, encompassing 78 very low mass members observed with Keck laser guide star adaptive optics. Our survey confirms the overall trend observed in the field for lower-mass binary systems to be less frequent and more compact, including a null detection for any substellar binary systems with separations wider than ~7 AU. Combined with a Bayesian re-analysis of existing surveys, our results demonstrate that the binary frequency and binary separations decline smoothly between masses of 0.5 Msun and 0.02 Msun, though we can not distinguish the functional form of this decline due to a degeneracy between the total binary frequency and the mean binary separation. We also show that the mass ratio distribution becomes progressively more concentrated at q~1 for declining masses, though a small number of systems appear to have unusually wide separations and low mass ratios for their mass. Finally, we compare our results to synthetic binary populations generated by smoothed particle hydrodynamic simulations, noting the similarities and discussing possible explanations for the differences.
Misaligned blazars have been the subject of some of the most successful radio and gamma-ray multiwavelength campaigns. These campaigns have included many of the major ground and space based gamma-ray telescopes and span decades of energy. Even though misaligned blazars account for only a small number of the total AGN detected at VHE, they provide a unique view on the AGN population. By viewing blazars at larger angles to our line of sight, they become a unique laboratory for the study of AGN jet substructure and the morphology of non-thermal emission processes. This contribution will discuss our understanding of three VHE misaligned blazars.
We describe approximate axisymmetric computations of the dynamical evolution of material inside radio lobes and X-ray cluster gas cavities in Fanaroff-Riley II sources such as Cygnus A. All energy is delivered by a jet to the lobe/cavity via a moving hotspot where jet energy dissipates in a reverse shock. Our calculations describe the evolution of hot plasma, cosmic rays (CRs) and toroidal magnetic fields flowing from the hotspot into the cavity. Many observed features are explained. Gas, CRs and field flow back along the cavity surface in a "boundary backflow" consistent with detailed FRII observations. Computed ages of backflowing CRs are consistent with observed radio-synchrotron age variations only if shear instabilities in the boundary backflow are damped and we assume this is done with viscosity of unknown origin. Magnetic fields estimated from synchrotron self-Compton (SSC) X-radiation observed near the hotspot evolve into radio lobe fields. Computed profiles of radio synchrotron lobe emission perpendicular to the jet are dramatically limb-brightened in excellent agreement with FRII observations although computed lobe fields exceed those observed. Strong winds flowing from hotspots naturally create kpc-sized spatial offsets between hotspot inverse Compton (IC-CMB) X-ray emission and radio synchrotron emission that peaks 1-2 kpc ahead where the field increases due to wind compression. In our computed version of Cygnus A, nonthermal X-ray emission increases from the hotspot (some IC-CMB, mostly SSC) toward the offset radio synchrotron peak (mostly SSC). A faint thermal jet along the symmetry axis may be responsible for redirecting the Cygnus A non-thermal jet.
Latent force models (LFMs) are flexible models that combine mechanistic modelling principles (i.e., physical models) with non-parametric data-driven components. Several key applications of LFMs need non-linearities, which results in analytically intractable inference. In this work we show how non-linear LFMs can be represented as non-linear white noise driven state-space models and present an efficient non-linear Kalman filtering and smoothing based method for approximate state and parameter inference. We illustrate the performance of the proposed methodology via two simulated examples, and apply it to a real-world problem of long-term prediction of GPS satellite orbits.
We recently proposed a chameleonic solution to the cosmological constant problem - Phys. Rev. D82 (2010) 044006. We further elaborate on our proposal discussing also some new results. The basic ingredients of the model are A) a chameleonic string dilaton parametrizing the amount of scale symmetry of the problem in the Einstein frame; B) a dual nature of the concept of particle where a splitting in local and global components is introduced (for every local particle a corresponding global particle is present in the theory). In this paper we proceed stepwise: 1) we point out that the value of the chameleonic dilaton in the Einstein frame parametrizes also the string length and, in particular, the string mass is an increasing function of the matter density; 2) the concept of global particle is clarified (indeed global particles are simply cosmic strings); 3) in the last quantization step, in the Feynman diagrams of our model, the UV cut-off (which is chosen to be the string mass) is a function of the chameleonic dilaton (naturally this stringy regularization is fully compatible with the chameleonic scale invariance of the model); 4) we show that a large number of particle species might be useful to keep locally under control the dangerous variations of the fundamental constants, however, more research efforts are necessary to make this point clear; 5) the chameleonic behaviour of matter fields and of the de Broglie wavelength of local matter particles is pointed out; 6) we briefly touch upon some ideas (which are still waiting for a full confirmation) regarding a potential connection between the shortest length scale of nature and the cosmological constant. A detailed phenomenological analysis of the entire model is required to test these ideas.
In this paper we discuss the restrictions of the spacetime for the standard model of cosmology by using results of the differential topology of 3- and 4-manifolds. The smoothness of the cosmic evolution is the strongest restriction. The Poincare model (dodecaeder model), the Picard horn and the 3-torus are ruled out by the restrictions but a sum of two Poincare spheres is allowed.
We discuss the Hamiltonian dynamics for cosmologies coming from Extended Theories of Gravity. In particular, minisuperspace models are taken into account searching for Noether symmetries. The existence of conserved quantities gives selection rule to recover classical behaviors in cosmic evolution according to the so called Hartle criterion, that allows to select correlated regions in the configuration space of dynamical variables. We show that such a statement works for general classes of Extended Theories of Gravity and is conformally preserved. Furthermore, the presence of Noether symmetries allows a straightforward classification of singularities that represent the points where the symmetry is broken. Examples of nonminimally coupled and higher-order models are discussed.
One possible interpretation of the holographic principle is the equality of the number of degrees of freedom in a bulk region of space and the number of degrees of freedom on the boundary surface. It is known that such an equality is maintained on equipotential surfaces in any static spacetime in the form of an equipartition law N_{bulk}= N_{sur}. In the cosmological context, the de Sitter universe obeys the same holographic equipartition. I argue that the difference between the surface degrees of freedom and the bulk degrees of freedom in a region of space (which has already emerged) drives the accelerated expansion of the universe through a simple equation dV/dt = (N_{sur} - N_{bulk}) where V is the Hubble volume in Planck units and t is the cosmic time in Planck units. This equation reproduces the standard evolution of the universe. This approach provides a novel paradigm to study the emergence of space and cosmology and has far reaching implications.
A variation of sneutrino inflation based on chi^2 potential is considered where the inflaton and the late-decaying field are sneutrinos of different generations. The lighter, late-decaying sneutrino dilutes the gravitinos over-produced after inflaton decay and generates the matter asymmetry. It can also significantly contribute to the curvature perturbation, realizing the mixed inflaton-curvaton case. The cosmic microwave background (CMB) observables can distinguish this case from inflation with chi^2 potential provided that the initial value of the late-decaying sneutrino is either an order of magnitude smaller or larger than the reduced Planck scale.
Links to: arXiv, form interface, find, astro-ph, recent, 1206, contact, help (Access key information)