-
Using a state-of-the-art cosmological simulation of merging proto-galaxies at
high redshift from the FIRE project, with explicit treatments of star formation
and stellar feedback in the interstellar medium, we investigate the formation
of star clusters and examine one of the formation hypothesis of present-day
metal-poor globular clusters. We find that frequent mergers in high-redshift
proto-galaxies could provide a fertile environment to produce long-lasting
bound star clusters. The violent merger event disturbs the gravitational
potential and pushes a large gas mass of ~> 1e5-6 Msun collectively to high
density, at which point it rapidly turns into stars before stellar feedback can
stop star formation. The high dynamic range of the reported simulation is
critical in realizing such dense star-forming clouds with a small dynamical
timescale, t_ff <~ 3 Myr, shorter than most stellar feedback timescales. Our
simulation then allows us to trace how clusters could become virialized and
tightly-bound to survive for up to ~420 Myr till the end of the simulation.
Because the cluster's tightly-bound core was formed in one short burst, and the
nearby older stars originally grouped with the cluster tend to be
preferentially removed, at the end of the simulation the cluster has a small
age spread.
-
The CHIME Pathfinder is a new interferometric radio telescope that uses a
hybrid FPGA/GPU FX correlator. The GPU-based X-engine of this correlator
processes over 819 Gb/s of 4+4-bit complex astronomical data from N=256 inputs
across a 400 MHz radio band. A software framework is presented to manage this
real-time data flow, which allows each of 16 processing servers to handle 51.2
Gb/s of astronomical data, plus 8 Gb/s of ancillary data. Each server receives
data in the form of UDP packets from an FPGA F-engine over the eight 10 GbE
links, combines data from these packets into large (32MB-256MB) buffered
frames, and transfers them to multiple GPU co-processors for correlation. The
results from the GPUs are combined and normalized, then transmitted to a
collection server, where they are merged into a single file. Aggressive
optimizations enable each server to handle this high rate of data; allowing the
efficient correlation of 25 MHz of radio bandwidth per server. The solution
scales well to larger values of N by adding additional servers.
-
We present the radial-velocity follow-up of two Kepler planetary transiting
candidates (KOI-189 and KOI-686) carried out with the SOPHIE spectrograph at
the Observatoire de Haute Provence. These data promptly discard these objects
as viable planet candidates and show that the transiting objects are in the
regime of very low-mass stars, where a strong discrepancy between observations
and models persists for the mass and radius parameters. By combining the SOPHIE
spectra with the Kepler light curve and photometric measurements found in the
literature, we obtain a full characterization of the transiting companions,
their orbits, and their host stars. The two companions are in significantly
eccentric orbits with relatively long periods (30 days and 52.5 days), which
makes them suitable objects for a comparison with theoretical models, since the
effects invoked to understand the discrepancy with observations are weaker for
these orbital distances. KOI-189 B has a mass M = 0.0745 +/- 0.0033 Msun and a
radius R = 0.1025 +/- 0.0024 Rsun. The density of KOI-189 B is significantly
lower than expected from theoretical models for a system of its age. We explore
possible explanations for this difference. KOI-189 B is the smallest
hydrogen-burning star with such a precise determination of its fundamental
parameters. KOI-686 B is larger and more massive (M = 0.0915 +/- 0.0043 Msun; R
= 0.1201 +/- 0.0033 Rsun), and its position in the mass-radius diagram agrees
well with theoretical expectations.
-
Mapping Mg II resonance emission scattered by galactic winds offers a means
to determine the spatial extent and density of the warm outflow. Using
Keck/LRIS spectroscopy, we have resolved scattered Mg II emission to the east
of 32016857, a star-forming galaxy at z =0.9392 with an outflow. The Mg II
emission from this galaxy exhibits a P-Cygni profile, extends further than both
the continuum and [O II] emission along the eastern side of the slit, and has a
constant Doppler shift along the slit which does not follow the velocity
gradient of the nebular [O II] emission. Using the Sobolev approximation, we
derive the density of Mg+ ions at a radius of 12 to 18 kpc in the outflow. We
model the ionization correction and find that much of the outflowing Mg is in
Mg++. We estimate that the total mass flux could be as large as 330 - 500 solar
masses per year, with the largest uncertainties coming from the depletion of Mg
onto grains and the clumpiness of the warm outflow. We show that confining the
warm clouds with a hot wind reduces the estimated mass flux of the warm outflow
and indicates amass-loading factor near unity in the warm phase alone. Based on
the high blue luminosities that distinguish 32016857 and TKRS 4389, described
by Rubin et al. 2011, from other galaxies with P-Cygni emission, we suggest
that, as sensitivity to diffuse emission improves, scattering halos may prove
to be a generic property of star-forming galaxies at intermediate redshifts.
-
The whole of the Taurus region (a total area of 52 sq. deg.) has been
observed by the Herschel SPIRE and PACS instruments at wavelengths of 70, 160,
250, 350 and 500 {\mu}m as part of the Herschel Gould Belt Survey. In this
paper we present the first results from the part of the Taurus region that
includes the Barnard 18 and L1536 clouds. A new source-finding routine, the
Cardiff Source-finding AlgoRithm (CSAR), is introduced, which is loosely based
on CLUMPFIND, but that also generates a structure tree, or dendrogram, which
can be used to interpret hierarchical clump structure in a complex region.
Sources were extracted from the data using the hierarchical version of CSAR and
plotted on a mass-size diagram. We found a hierarchy of objects with sizes in
the range 0.024-2.7 pc. Previous studies showed that gravitationally bound
prestellar cores and unbound starless clumps appeared in different places on
the mass-size diagram. However, it was unclear whether this was due to a lack
of instrumental dynamic range or whether they were actually two distinct
populations. The excellent sensitivity of Herschel shows that our sources fill
the gap in the mass-size plane between starless and pre-stellar cores, and
gives the first clear supporting observational evidence for the theory that
unbound clumps and (gravitationally bound) prestellar cores are all part of the
same population, and hence presumably part of the same evolutionary sequence
(c.f. Simpson et al. 2011).
-
We present observational evidence for the inhibition of bar formation in
dispersion-dominated (dynamically hot) galaxies by studying the relationship
between galactic structure and host galaxy kinematics in a sample of 257
galaxies between 0.1 $<$ z $\leq$ 0.84 from the All-Wavelength Extended Groth
Strip International Survey (AEGIS) and the Deep Extragalactic Evolutionary
Probe 2 (DEEP2) survey. We find that bars are preferentially found in galaxies
that are massive and dynamically cold (rotation-dominated) and on the stellar
Tully-Fisher relationship, as is the case for barred spirals in the local
Universe. The data provide at least one explanation for the steep ($\times$3)
decline in the overall bar fraction from z=0 to z=0.84 in L$^*$ and brighter
disks seen in previous studies. The decline in the bar fraction at high
redshift is almost exclusively in the lower mass (10 $<$ log M$_{*}$(\Msun)$<$
11), later-type and bluer galaxies. A proposed explanation for this
"downsizing" of the bar formation / stellar structure formation is that the
lower mass galaxies may not form bars because they could be dynamically hotter
than more massive systems from the increased turbulence of accreting gas,
elevated star formation, and/or increased interaction/merger rate at higher
redshifts. The evidence presented here provides observational support for this
hypothesis. However, the data also show that not every disk galaxy that is
massive and cold has a stellar bar, suggesting that mass and dynamic coldness
of a disk are necessary but not sufficient conditions for bar formation -- a
secondary process, perhaps the interaction history between the dark matter halo
and the baryonic matter, may play an important role in bar formation.
-
Typical flows in stellar interiors are much slower than the speed of sound.
To follow the slow evolution of subsonic motions, various sound-proof equations
are in wide use, particularly in stellar astrophysical fluid dynamics. These
low-Mach number equations include the anelastic equations. Generally, these
equations are valid in nearly adiabatically stratified regions like stellar
convection zones, but may not be valid in the sub-adiabatic, stably stratified
stellar radiative interiors. Understanding the coupling between the convection
zone and the radiative interior is a problem of crucial interest and may have
strong implications for solar and stellar dynamo theories as the interface
between the two, called the tachocline in the Sun, plays a crucial role in many
solar dynamo theories. Here we study the properties of gravity waves in
stably-stratified atmospheres. In particular, we explore how gravity waves are
handled in various sound-proof equations. We find that some anelastic
treatments fail to conserve energy in stably-stratified atmospheres, instead
conserving pseudo-energies that depend on the stratification, and we
demonstrate this numerically. One anelastic equation set does conserve energy
in all atmospheres and we provide recommendations for converting low-Mach
number anelastic codes to this set of equations.
-
Young solar-type stars rotate rapidly and many are magnetically active; some
undergo magnetic cycles similar to the 22-year solar activity cycle. We conduct
simulations of dynamo action in rapidly rotating suns with the 3D MHD anelastic
spherical harmonic (ASH) code to explore dynamo action achieved in the
convective envelope of a solar-type star rotating at 5 times the current solar
rotation rate. Striking global-scale magnetic wreaths appear in the midst of
the turbulent convection zone and show rich time-dependence. The dynamo
exhibits cyclic activity and undergoes quasi-periodic polarity reversals where
both the global-scale poloidal and toroidal fields change in sense on a roughly
1500 day time scale. These magnetic activity patterns emerge spontaneously from
the turbulent flow and are more organized temporally and spatially than those
realized in our previous simulations of the solar dynamo. We assess in detail
the competing processes of magnetic field creation and destruction within our
simulations that contribute to the global-scale reversals. We find that the
mean toroidal fields are built primarily through an $\Omega$-effect, while the
mean poloidal fields are built by turbulent correlations which are not
necessarily well represented by a simple $\alpha$-effect. During a reversal the
magnetic wreaths propagate towards the polar regions, and this appears to arise
from a poleward propagating dynamo wave. The primary response in the convective
flows involves the axisymmetric differential rotation which shows variations
associated with the poleward propagating magnetic wreaths. In the Sun, similar
patterns are observed in the poleward branch of the torsional oscillations, and
these may represent poleward propagating magnetic fields deep below the solar
surface. [abridged]
-
Young solar-type stars rotate rapidly and are very magnetically active. The
magnetic fields at their surfaces likely originate in their convective
envelopes where convection and rotation can drive strong dynamo action. Here we
explore simulations of global-scale stellar convection in rapidly rotating suns
using the 3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields
built in these dynamos are organized on global-scales into wreath-like
structures that span the convection zone. We explore one case rotates five
times faster than the Sun in detail. This dynamo simulation, called case D5,
has repeated quasi-cyclic reversals of global-scale polarity. We compare this
case D5 to the broader family of simulations we have been able to explore and
discuss how future simulations and observations can advance our understanding
of stellar dynamos and magnetism.
-
When our Sun was young it rotated much more rapidly than now. Observations of
young, rapidly rotating stars indicate that many possess substantial magnetic
activity and strong axisymmetric magnetic fields. We conduct simulations of
dynamo action in rapidly rotating suns with the 3-D MHD anelastic spherical
harmonic (ASH) code to explore the complex coupling between rotation,
convection and magnetism. Here we study dynamo action realized in the bulk of
the convection zone for a system rotating at three times the current solar
rotation rate. We find that substantial organized global-scale magnetic fields
are achieved by dynamo action in this system. Striking wreaths of magnetism are
built in the midst of the convection zone, coexisting with the turbulent
convection. This is a surprise, for it has been widely believed that such
magnetic structures should be disrupted by magnetic buoyancy or turbulent
pumping. Thus, many solar dynamo theories have suggested that a tachocline of
penetration and shear at the base of the convection zone is a crucial
ingredient for organized dynamo action, whereas these simulations do not
include such tachoclines. We examine how these persistent magnetic wreaths are
maintained by dynamo processes and explore whether a classical mean-field
$\alpha$-effect explains the regeneration of poloidal field.
-
When stars like our Sun are young they rotate rapidly and are very
magnetically active. We explore dynamo action in rapidly rotating suns with the
3-D MHD anelastic spherical harmonic (ASH) code. The magnetic fields built in
these dynamos are organized on global-scales into wreath-like structures that
span the convection zone. Wreath-building dynamos can undergo quasi-cyclic
reversals of polarity and such behavior is common in the parameter space we
have been able to explore. These dynamos do not appear to require tachoclines
to achieve their spatial or temporal organization. Wreath-building dynamos are
present to some degree at all rotation rates, but are most evident in the more
rapidly rotating simulations.
-
[Abridged] Superluminous Supernovae (SN2006gy, SN2005gj, SN2005ap, SN2008fz,
SN2003ma) have been a challenge to explain by standard models. We present an
alternative scenario involving a quark-nova (QN), an explosive transition of
the newly born neutron star to a quark star in which a second explosion
(delayed) occurs inside the already expanding ejecta of a normal SN. The
reheated SN ejecta can radiate at higher levels for longer periods of time
primarily due to reduced adiabatic expansion losses, unlike the standard SN
case. Our model is successfully applied to SN2006gy, SN2005gj, SN2005ap,
SN2008fz, SN2003ma with encouraging fits to the lightcurves. There are four
predictions in our model: (i) superluminous SNe optical lightcurves should show
a double-hump with the SN hump at weaker magnitudes occurring days to weeks
before the QN; (ii) Two shock breakouts should be observed vis-a-vis one for a
normal SN. Depending on the time delay, this would manifest as two distinct
spikes in the X-ray region or a broadening of the first spike for extremely
short delays; (iii) The QN deposits heavy elements of mass number A> 130 at the
base of the preceeding SN ejecta. These QN r-processed elements should be
visible in the late spectrum (few days-weeks in case of strong ejecta mixing)
of the superluminous SN; (iv) The QN yield will also contain lighter elements
(Hydrogen and Helium). We expect the late spectra to include H_alpha emission
lines that should be distinct in their velocity signature from standard H_alpha
emission.
-
Primordial star formation appears to result in stars at least an order of
magnitude more massive than modern star formation. It has been proposed that
the transition from primordial to modern initial mass functions occurs due to
the onset of effective metal line cooling at a metallicity Z/Z_sun = 10^{-3.5}.
However, these simulations neglected molecular hydrogen cooling. We perform
simulations using the same initial conditions, but including molecular cooling,
using a complex network that follows molecular hydrogen formation and also
directly follows carbon monoxide and water. We find that molecular hydrogen
cooling allows roughly equivalent fragmentation to proceed even at zero
metallicity for these initial conditions. The apparent transition just
represents the point where metal line cooling becomes more important than
molecular cooling. In all cases, the fragments are massive enough to be
consistent with models of primordial stellar masses, suggesting that the
transition to the modern initial mass function may be determined by other
physics such as dust formation. We conclude that such additional cooling
mechanisms, combined with the exact initial conditions produced by cosmological
collapse are likely more important than metal line cooling in determining the
initial mass function, and thus that there is unlikely to be a sharp transition
in the initial mass function at Z/Z_sun = 10^{-3.5}.
-
Weak lensing measurements are starting to provide statistical maps of the
distribution of matter in the universe that are increasingly precise and
complementary to cosmic microwave background maps. The probability distribution
(PDF) provides a powerful tool to test non-Gausian features in the convergence
field and to discriminate the different cosmological models. In this letter, we
present a new PDF space Wiener filter approach to reconstruct the probability
density function of the convergence from the noisy convergence field. We find
that for parameters comparable to the CFHT legacy survey, the averaged PDF of
the convergence in a 3 degree field can be reconstructed with an uncertainty of
about 10%, even though the pointwise PDF is noise dominated.