Table of contents

Both observations and numerical simulations are discordant with predictions of conventional stellar evolution codes for the latest stages of a massive star's life before core collapse. The most dramatic example of this disconnect is in the eruptive mass loss occurring in the decade preceding Type IIn supernovae. We outline the key empirical evidence that indicates severe pre-supernova instability in massive stars, and we suggest that the chief reason that these outbursts are absent in stellar evolution models may lie in the treatment of turbulent convection in these codes. The mixing length theory that is used ignores (1) finite amplitude fluctuations in velocity and temperature and (2) their nonlinear interaction with nuclear burning. Including these fluctuations is likely to give rise to hydrodynamic instabilities in the latest burning sequences, which prompts us to discuss a number of far-reaching implications for the fates of massive stars. In particular, we explore connections to enhanced pre-supernova mass loss, unsteady nuclear burning and consequent eruptions, swelling of the stellar radius that may trigger violent interactions with a companion star, and potential modifications to the core structure that could dramatically alter calculations of the core-collapse explosion mechanism itself. These modifications may also impact detailed nucleosynthesis and measured isotopic anomalies in meteorites, as well as the interpretation of young core-collapse supernova remnants. Understanding these critical instabilities in the final stages of evolution may make possible the development of an early warning system for impending core collapse, if we can identify their asteroseismological or eruptive signatures.

As part of the Binaries Are Not Always Neatly Aligned project (BANANA), we have found that the eclipsing binary CV Velorum has misaligned rotation axes. Based on our analysis of the Rossiter–McLaughlin effect, we find sky-projected spin-orbit angles of β_p = −52° ± 6° and β_s = 3° ± 7° for the primary and secondary stars (B2.5V + B2.5V, P = 6.9 days). We combine this information with several measurements of changing projected stellar rotation speeds (vsin i_⋆) over the last 30 yr, leading to a model in which the primary star's obliquity is ≈65°, and its spin axis precesses around the total angular momentum vector with a period of about 140 yr. The geometry of the secondary star is less clear, although a significant obliquity is also implicated by the observed time variations in the vsin i_⋆. By integrating the secular tidal evolution equations backward in time, we find that the system could have evolved from a state of even stronger misalignment similar to DI  Herculis, a younger but otherwise comparable binary.

We use a parent sample of 118 gamma-ray burst (GRB) afterglows, with known redshift and host galaxy extinction, to separate afterglows with and without signatures of dominant reverse-shock (RS) emission and to determine which physical conditions lead to a prominent reverse-shock emission. We identify 10 GRBs with reverse-shock signatures: 990123, 021004, 021211, 060908, 061126, 080319B, 081007, 090102, 090424, and 130427A. By modeling their optical afterglows with reverse- and forward-shock analytic light curves and using Monte Carlo simulations, we estimate the parameter space of the physical quantities describing the ejecta and circumburst medium. We find that physical properties cover a wide parameter space and do not seem to cluster around any preferential values. Comparing the rest-frame optical, X-ray, and high-energy properties of the larger sample of non-RS-dominated GRBs, we show that the early-time (<1 ks) optical spectral luminosity, X-ray afterglow luminosity, and γ-ray energy output of our reverse-shock dominated sample do not differ significantly from the general population at early times. However, the GRBs with dominant reverse-shock emission have fainter than average optical forward-shock emission at late times (>10 ks). We find that GRBs with an identifiable reverse-shock component show a high magnetization parameter R_B = ε_{B, r}/ε_{B, f} ∼ 2–10⁴. Our results are in agreement with the mildly magnetized baryonic jet model of GRBs.

We present a detailed study of the interaction process of two coronal mass ejections (CMEs) successively launched on 2011 February 14 (CME1) and 2011 February 15 (CME2). Reconstructing the three-dimensional shape and evolution of the flux ropes, we verify that the two CMEs interact. The frontal structure of both CMEs, measured along different position angles (PAs) over the entire latitudinal extent, reveals differences in the kinematics for the interacting flanks and the apexes. The interaction process is strongly PA-dependent in terms of timing as well as kinematical evolution. The central interaction occurs along PA-100°, which shows the strongest changes in kinematics. During interaction, CME1 accelerates from ∼400 km s⁻¹ to ∼700 km s⁻¹ and CME2 decelerates from ∼1300 km s⁻¹ to ∼600 km s⁻¹. Our results indicate that a simplified scenario such as inelastic collision may not be sufficient to describe the CME–CME interaction. The magnetic field structures of the intertwining flux ropes and the momentum transfer due to shocks each play an important role in the interaction process.

To determine the observable radio signatures of the fast sausage standing wave, we examine gyrosynchrotron (GS) emission modulation using a linear three-dimensional magnetohydrodynamic model of a plasma cylinder. Effects of the line-of-sight angle and instrumental resolution on perturbations of the GS intensity are analyzed for two models: a base model with strong Razin suppression and a low-density model in which the Razin effect was unimportant. Our finding contradicts previous predictions made with simpler models: an in-phase variation of intensity between low (f < f_peak) and high (f > f_peak) frequencies is found for the low-density model and an anti-phase variation for the base model in the case of a viewing angle of 45°. The spatially inhomogeneous character of the oscillating emission source and the spatial resolution of the model are found to have a significant effect on the resulting intensity.

A magnetohydrodynamic model is used to determine conditions under which the Lorentz force accelerates plasma to type 2 spicule speeds in the chromosphere. The model generalizes a previous model to include a more realistic pre-spicule state, and the vertical viscous force. Two cases of acceleration under upper chromospheric conditions are considered. The magnetic field strength for these cases is ⩽12.5 and 25 G. Plasma is accelerated to terminal vertical speeds of 66 and 78 km s⁻¹ in 100 s, compared with 124 and 397 km s⁻¹ for the case of zero viscosity. The flows are localized within horizontal diameters ∼80 and 50 km. The total thermal energy generated by viscous dissipation is ∼10 times larger than that due to Joule dissipation, but the magnitude of the total cooling due to rarefaction is ≳ this energy. Compressive heating dominates during the early phase of acceleration. The maximum energy injected into the corona by type 2 spicules, defined as the energy flux in the upper chromosphere, may largely balance total coronal energy losses in quiet regions, possibly also in coronal holes, but not in active regions. It is proposed that magnetic flux emergence in intergranular regions drives type 2 spicules.

Solar filaments are magnetic structures often observed in the solar atmosphere and consist of plasma that is cooler and denser than their surroundings. They are visible for days—even weeks—which suggests that they are often in equilibrium with their environment before disappearing or erupting. Several eruption models have been proposed that aim to reveal what mechanism causes (or triggers) these solar eruptions. Validating these models through observations represents a fundamental step in our understanding of solar eruptions. We present an analysis of the observation of a filament eruption that agrees with the torus instability model. This model predicts that a magnetic flux rope embedded in an ambient field undergoes an eruption when the axis of the flux rope reaches a critical height that depends on the topology of the ambient field. We use the two vantage points of the Solar Dynamics Observatory (SDO) and the Solar TErrestrial RElations Observatory to reconstruct the three-dimensional shape of the filament, to follow its morphological evolution, and to determine its height just before eruption. The magnetograms acquired by SDO/Helioseismic and Magnetic Imager are used to infer the topology of the ambient field and to derive the critical height for the onset of the torus instability. Our analysis shows that the torus instability is the trigger of the eruption. We also find that some pre-eruptive processes, such as magnetic reconnection during the observed flares and flux cancellation at the neutral line, facilitated the eruption by bringing the filament to a region where the magnetic field was more vulnerable to the torus instability.

The PICARD spacecraft was launched on 2010 June 15 with the scientific objective of studying the geometry of the Sun. It is difficult to measure solar oblateness because images are affected by optical distortion. Rolling the satellite, as done in previous space missions, determines the contribution of the telescope by assuming that the geometry of the Sun is constant during the observations. The optical response of the telescope is considered to be time-invariant during the roll operations. This is not the case for PICARD because an orbital signature is clearly observed in the solar radius computed from its images. We take this effect into account and provide the new space value of solar oblateness from PICARD images recorded in the solar continuum at 535.7 nm on 2011 July 4–5. The equator–pole radius difference is 8.4 ± 0.5 mas, which corresponds to an absolute radius difference of 6.1 km. This coincides with the mean value of all solar oblateness measurements obtained during the last two decades from the ground, balloons, and space. It is also consistent with values determined from models using helioseismology data.

We present numerical simulations of active region scale flux emergence covering a time span of up to 6 days. Flux emergence is driven by a bottom boundary condition that advects a semi-torus of magnetic field with 1.7 × 10²² Mx flux into the computational domain. The simulations show that, even in the absence of twist, the magnetic flux is able the rise through the upper 15.5 Mm of the convection zone and emerge into the photosphere to form spots. We find that spot formation is sensitive to the persistence of upflows at the bottom boundary footpoints, i.e., a continuing upflow would prevent spot formation. In addition, the presence of a torus-aligned flow (such flow into the retrograde direction is expected from angular momentum conservation during the rise of flux ropes through the convection zone) leads to a significant asymmetry between the pair of spots, with the spot corresponding to the leading spot on the Sun being more axisymmetric and coherent, but also forming with a delay relative to the following spot. The spot formation phase transitions directly into a decay phase. Subsurface flows fragment the magnetic field and lead to intrusions of almost field free plasma underneath the photosphere. When such intrusions reach photospheric layers, the spot fragments. The timescale for spot decay is comparable to the longest convective timescales present in the simulation domain. We find that the dispersal of flux from a simulated spot in the first two days of the decay phase is consistent with self-similar decay by turbulent diffusion.

The classic optical nebular diagnostics [N ii], [O ii], [O iii], [S ii], [S iii], and [Ar iii] are employed to search for evidence of non-Maxwellian electron distributions, namely κ distributions, in a sample of well-observed Galactic H ii regions. By computing new effective collision strengths for all these systems and A-values when necessary (e.g., S ii), and by comparing with previous collisional and radiative data sets, we have been able to obtain realistic estimates of the electron-temperature dispersion caused by the atomic data, which in most cases are not larger than ∼10%. If the uncertainties due to both observation and atomic data are then taken into account, it is plausible to determine for some nebulae a representative average temperature while in others there are at least two plasma excitation regions. For the latter, it is found that the diagnostic temperature differences in the high-excitation region, e.g., T_e(O iii), T_e(S iii), and T_e(Ar iii), cannot be conciliated by invoking κ distributions. For the low-excitation region, it is possible in some, but not all, cases to arrive at a common, lower temperature for [N ii], [O ii], and [S ii] with κ ≈ 10, which would then lead to significant abundance enhancements for these ions. An analytic formula is proposed to generate accurate κ-averaged excitation rate coefficients (better than 10% for κ ⩾ 5) from temperature tabulations of the Maxwell–Boltzmann effective collision strengths.

We present three-dimensional atmospheric circulation models of GJ 1214b, a 2.7 Earth-radius, 6.5 Earth-mass super Earth detected by the MEarth survey. Here we explore the planet's circulation as a function of atmospheric metallicity and atmospheric composition, modeling atmospheres with a low mean molecular weight (MMW; i.e., H₂-dominated) and a high MMW (i.e., water- and CO₂-dominated). We find that atmospheres with a low MMW have strong day–night temperature variations at pressures above the infrared photosphere that lead to equatorial superrotation. For these atmospheres, the enhancement of atmospheric opacities with increasing metallicity lead to shallower atmospheric heating, larger day–night temperature variations, and hence stronger superrotation. In comparison, atmospheres with a high MMW have larger day–night and equator-to-pole temperature variations than low MMW atmospheres, but differences in opacity structure and energy budget lead to differences in jet structure. The circulation of a water-dominated atmosphere is dominated by equatorial superrotation, while the circulation of a CO₂-dominated atmosphere is instead dominated by high-latitude jets. By comparing emergent flux spectra and light curves for 50× solar and water-dominated compositions, we show that observations in emission can break the degeneracy in determining the atmospheric composition of GJ 1214b. The variation in opacity with wavelength for the water-dominated atmosphere leads to large phase variations within water bands and small phase variations outside of water bands. The 50× solar atmosphere, however, yields small variations within water bands and large phase variations at other characteristic wavelengths. These observations would be much less sensitive to clouds, condensates, and hazes than transit observations.

Though there are now many hundreds of confirmed exoplanets known, the binarity of exoplanet host stars is not well understood. This is particularly true of host stars that harbor a giant planet in a highly eccentric orbit since these are more likely to have had a dramatic dynamical history that transferred angular momentum to the planet. Here we present observations of four exoplanet host stars that utilize the excellent resolving power of the Differential Speckle Survey Instrument on the Gemini North telescope. Two of the stars are giants and two are dwarfs. Each star is host to a giant planet with an orbital eccentricity >0.5 and whose radial velocity (RV) data contain a trend in the residuals to the Keplerian orbit fit. These observations rule out stellar companions 4–8 mag fainter than the host star at passbands of 692 nm and 880 nm. The resolution and field of view of the instrument result in exclusion radii of 005–14, which excludes stellar companions within several AU of the host star in most cases. We further provide new RVs for the HD 4203 system that confirm that the linear trend previously observed in the residuals is due to an additional planet. These results place dynamical constraints on the source of the planet's eccentricities, place constraints on additional planetary companions, and inform the known distribution of multiplicity amongst exoplanet host stars.

The lithium abundances for 378 G/K giants are derived with non-local thermodynamic equilibrium correction considered. Among these are 23 stars that host planetary systems. The lithium abundance is investigated, as a function of metallicity, effective temperature, and rotational velocity, as well as the impact of a giant planet on G/K giants. The results show that the lithium abundance is a function of metallicity and effective temperature. The lithium abundance has no correlation with rotational velocity at v sin  i < 10 km s⁻¹. Giants with planets present lower lithium abundance and slow rotational velocity (v sin  i < 4 km s⁻¹). Our sample includes three Li-rich G/K giants, 36 Li-normal stars, and 339 Li-depleted stars. The fraction of Li-rich stars in this sample agrees with the general rate of less than 1% in the literature, and the stars that show normal amounts of Li are supposed to possess the same abundance at the current interstellar medium. For the Li-depleted giants, Li-deficiency may have already taken place at the main sequence stage for many intermediate mass (1.5–5 M_☉) G/K giants. Finally, we present the lithium abundance and kinematic parameters for an enlarged sample of 565 giants using a compilation of the literature, and confirm that the lithium abundance is a function of metallicity and effective temperature. With the enlarged sample, we investigate the differences between the lithium abundance in thin-/thick-disk giants, which indicate that the lithium abundance in thick-disk giants is more depleted than that in thin-disk giants.

Type IIb supernova (SN) 2011dh, with conclusive detection of an unprecedented yellow supergiant (YSG) progenitor, provides an excellent opportunity to deepen our understanding on the massive star evolution in the final centuries toward the SN explosion. In this paper, we report on detection and analyses of thermal X-ray emission from SN IIb 2011dh at ∼500 days after the explosion on Chandra archival data, providing a solidly derived mass-loss rate of a YSG progenitor for the first time. We find that the circumstellar media should be dense, more than that expected from a Wolf–Rayet (W-R) star by one order of magnitude. The emission is powered by a reverse shock penetrating into an outer envelope, fully consistent with the YSG progenitor but not with a W-R progenitor. The density distribution at the outermost ejecta is much steeper than that expected from a compact W-R star, and this finding must be taken into account in modeling the early UV/optical emission from SNe IIb. The derived mass-loss rate is ∼3 × 10⁻⁶ M_☉ yr⁻¹ for the mass-loss velocity of ∼20 km s⁻¹ in the final ∼1300 yr before the explosion. The derived mass-loss properties are largely consistent with the standard wind mass-loss expected for a giant star. This is not sufficient to be a main driver to expel nearly all the hydrogen envelope. Therefore, the binary interaction, with a huge mass transfer having taken place at ≳ 1300 yr before the explosion, is a likely scenario to produce the YSG progenitor.

The cross sections of the ²H(d,p)³H and ²H(d,n)³He reactions have been measured via the Trojan Horse method applied to the quasi-free ²H(³He,p³H)¹H and ²H(³He,n³He)¹H processes at 18 MeV off the proton in ³He. For the first time, the bare nucleus S(E) factors have been determined from 1.5 MeV, across the relevant region for standard Big Bang nucleosynthesis, down to the thermal energies of deuterium burning in the pre-main-sequence (PMS) phase of stellar evolution, as well as of future fusion reactors. Both the energy dependence and the absolute value of the S(E) factors deviate by more than 15% from the available direct data and existing fitting curves, with substantial variations in the electron screening by more than 50%. As a consequence, the reaction rates for astrophysics experience relevant changes, with a maximum increase of up to 20% at the temperatures of the PMS phase. From a recent primordial abundance sensitivity study, it turns out that the ²H(d,n)³He reaction is quite influential on ⁷Li, and the present change in the reaction rate leads to a decrease in its abundance by up to 10%. The present reaction rates have also been included in an updated version of the FRANEC evolutionary code to analyze their influence on the central deuterium abundance in PMS stars with different masses. The largest variation of about 10%–15% pertains to young stars (⩽1 Myr) with masses ⩾1 M_☉.

We identified a general course of classical nova outbursts in the B − V versus U − B color–color diagram. It is reported that novae show spectra similar to those of A–F supergiants near optical light maximum. However, they do not follow the supergiant sequence in the color–color diagram, neither the blackbody nor the main-sequence sequence. Instead, we found that novae evolve along a new sequence in the pre-maximum and near-maximum phases, which we call "the nova-giant sequence." This sequence is parallel to but Δ(U − B) ≈ −0.2 mag bluer than the supergiant sequence. This is because the mass of a nova envelope is much (∼10⁻⁴ times) less than that of a normal supergiant. After optical maximum, its color quickly evolves back blueward along the same nova-giant sequence and reaches the point of free–free emission (B − V = −0.03, U − B = −0.97), which coincides with the intersection of the blackbody sequence and the nova-giant sequence, and remains there for a while. Then the color evolves leftward (blueward in B − V but almost constant in U − B), owing mainly to the development of strong emission lines. This is the general course of nova outbursts in the color–color diagram, which was deduced from eight well-observed novae in various speed classes. For a nova with unknown extinction, we can determine a reliable value of the color excess by matching the observed track of the target nova with this general course. This is a new and convenient method for obtaining the color excesses of classical novae. Using this method, we redetermined the color excesses of 20 well-observed novae. The obtained color excesses are in reasonable agreement with the previous results, which in turn support the idea of our general track of nova outbursts. Additionally, we estimated the absolute V magnitudes of about 30 novae using a method for time-stretching nova light curves to analyze the distance-reddening relations of the novae.

After the big bang nucleosynthesis, the first heavy element enrichment in the universe was made by a supernova (SN) explosion of a population (Pop) III star (Pop III SN). The abundance ratios of elements produced from Pop III SNe are recorded in abundance patterns of extremely metal-poor (EMP) stars. The observations of the increasing number of EMP stars have made it possible to statistically constrain the explosion properties of Pop III SNe. We present Pop III SN models whose nucleosynthesis yields well reproduce, individually, the abundance patterns of 48 such metal-poor stars as [Fe/H]  ≲   − 3.5. We then derive relations between the abundance ratios of EMP stars and certain explosion properties of Pop III SNe: the higher [(C + N)/Fe] and [(C + N)/Mg] ratios correspond to the smaller ejected Fe mass and the larger compact remnant mass, respectively. Using these relations, the distributions of the abundance ratios of EMP stars are converted to those of the explosion properties of Pop III SNe. Such distributions are compared with those of the explosion properties of present day SNe: the distribution of the ejected Fe mass of Pop III SNe has the same peak as that of the present day SNe but shows an extended tail down to ∼10⁻²–10⁻⁵ M_☉, and the distribution of the mass of the compact remnant of Pop III SNe is as wide as that of the present-day, stellar-mass black holes. Our results demonstrate the importance of large samples of EMP stars obtained by ongoing and future EMP star surveys and subsequent high-dispersion spectroscopic observations in clarifying the nature of Pop III SNe in the early universe.

Modeling the spectral emission of low-charge iron group ions enables the diagnostic determination of the local physical conditions of many cool plasma environments such as those found in H ii regions, planetary nebulae, active galactic nuclei, etc. Electron-impact excitation drives the population of the emitting levels and, hence, their emissivities. By carrying-out Breit-Pauli and intermediate coupling frame transformation (ICFT) R-matrix calculations for the electron-impact excitation of Fe^{2 +}, which both use the exact same atomic structure and the same close-coupling expansion, we demonstrate the validity of the application of the powerful ICFT method to low-charge iron group ions. This is in contradiction to the finding of Bautista et al., who carried-out ICFT and Dirac R-matrix calculations for the same ion. We discuss possible reasons.

We present new results on the chemical composition of the Galactic ring nebula NGC 6888 surrounding the WN6(h) star WR136. The data are based on deep spectroscopical observations taken with the High Dispersion Spectrograph at the 8.2 m Subaru Telescope. The spectra cover the optical range from 3700 to 7400 Å. The effect of the CNO cycle is well-identified in the abundances of He, N, and O, while elements not involved in the synthesis such as Ar, S, and Fe present values consistent with the solar vicinity and the ambient gas. The major achievement of this work is the first detection of the faint C ii λ4267 recombination line in a Wolf–Rayet nebula. This allows us to estimate the C abundance in NGC 6888 and therefore investigate for the first time the trace of the CNO cycle in a ring nebula around a Wolf–Rayet star. Although the detection of the C ii line has a low signal-to-noise ratio, the C abundance seems to be higher than the predictions of recent stellar evolution models of massive stars. The Ne abundance also shows a puzzling pattern with an abundance of about 0.5 dex lower than the solar vicinity, which may be related to the action of the NeNa cycle. Attending to the constraints imposed by the dynamical timescale and the He/H and N/O ratios of the nebula, the comparison with stellar evolution models indicates that the initial mass of the stellar progenitor of NGC 6888 is between 25 M_☉ and 40 M_☉.

We examine whether the magnetorotational instability (MRI) can serve as a mechanism of angular momentum transport in circumplanetary disks. For the MRI to operate the ionization degree must be sufficiently high and the magnetic pressure must be sufficiently lower than the gas pressure. We calculate the spatial distribution of the ionization degree and search for the MRI-active region where the two criteria are met. We find that there can be thin active layers at the disk surface depending on the model parameters, however, we find hardly any region which can sustain well-developed MRI turbulence; when the magnetic field is enhanced by MRI turbulence at the disk surface layer, a magnetically dominated atmosphere encroaches on a lower altitude and a region of well-developed MRI turbulence becomes smaller. We conclude that if there are no angular momentum transfer mechanisms other than MRI in gravitationally stable circumplanetary disks, gas is likely to pile up until disks become gravitationally unstable, and massive disks may survive for a long time.

We present α-element abundances of Mg, Si, and Ti for a large sample of field stars in two outer fields of the Fornax dwarf spheroidal (dSph) galaxy, obtained with Very Large Telescope/GIRAFFE (R ∼ 16, 000). Due to the large fraction of metal-poor (MP) stars in our sample, we are able to follow the α-element evolution from [Fe/H] ≈ −2.5 continuously to [Fe/H] ≈ −0.7. For the first time we are able to resolve the turnover from the Type II supernovae (SNe) dominated, α-enhanced plateau down to subsolar [α/Fe] values, due to the onset of SNe Ia, and thus to trace the chemical enrichment efficiency of the galaxy. Our data support the general concept of an α-enhanced plateau at early epochs, followed by a well-defined "knee" caused by the onset of SNe Ia, and finally a second plateau with sub-solar [α/Fe] values. We find the position of this knee to be at [Fe/H] ≈ −1.9 and therefore significantly more MP than expected from comparison with other dSphs and standard evolutionary models. Surprisingly, this value is rather comparable to the knee in Sculptor, a dSph ∼10 times less luminous than Fornax. Using chemical evolution models, we find that the position of the knee and the subsequent plateau at the sub-solar level can hardly be explained unless the galaxy experienced several discrete star formation (SF) events with a drastic variation in SF efficiency, while a uniform SF can be ruled out. One possible evolutionary scenario is that Fornax experienced one or several major accretion events from gas-rich systems in the past, so that its current stellar mass is not indicative of the chemical evolution environment at ancient times. If Fornax is the product of several smaller buildings blocks, this may also have implications for the understanding of the formation process of dSphs in general.

Aside from the grand-design stellar spirals appearing in the disk of M81, a pair of stellar spiral arms situated well inside the bright bulge of M81 has been recently discovered by Kendall et al. The seemingly unrelated pairs of spirals pose a challenge to the theory of spiral density waves. To address this problem, we have constructed a three-component model for M81, including the contributions from a stellar disk, a bulge, and a dark matter halo subject to observational constraints. Given this basic state for M81, a modal approach is applied to search for the discrete unstable spiral modes that may provide an understanding for the existence of both spiral arms. It is found that the apparently separated inner and outer spirals can be interpreted as a single trailing spiral mode. In particular, these spirals share the same pattern speed 25.5 km s⁻¹ kpc⁻¹ with a corotation radius of 9.03 kpc. In addition to the good agreement between the calculated and the observed spiral pattern, the variation of the spiral amplitude can also be naturally reproduced.

We describe redMaPPer, a new red sequence cluster finder specifically designed to make optimal use of ongoing and near-future large photometric surveys. The algorithm has multiple attractive features: (1) it can iteratively self-train the red sequence model based on a minimal spectroscopic training sample, an important feature for high-redshift surveys. (2) It can handle complex masks with varying depth. (3) It produces cluster-appropriate random points to enable large-scale structure studies. (4) All clusters are assigned a full redshift probability distribution P(z). (5) Similarly, clusters can have multiple candidate central galaxies, each with corresponding centering probabilities. (6) The algorithm is parallel and numerically efficient: it can run a Dark Energy Survey-like catalog in ∼500 CPU hours. (7) The algorithm exhibits excellent photometric redshift performance, the richness estimates are tightly correlated with external mass proxies, and the completeness and purity of the corresponding catalogs are superb. We apply the redMaPPer algorithm to ∼10, 000 deg² of SDSS DR8 data and present the resulting catalog of ∼25,000 clusters over the redshift range z ∈ [0.08, 0.55]. The redMaPPer photometric redshifts are nearly Gaussian, with a scatter σ_z ≈ 0.006 at z ≈ 0.1, increasing to σ_z ≈ 0.02 at z ≈ 0.5 due to increased photometric noise near the survey limit. The median value for |Δz|/(1 + z) for the full sample is 0.006. The incidence of projection effects is low (⩽5%). Detailed performance comparisons of the redMaPPer DR8 cluster catalog to X-ray and Sunyaev–Zel'dovich catalogs are presented in a companion paper.

Merging white dwarfs are a possible progenitor of Type Ia supernovae (SNe Ia). Numerical models suggest that a detonation might be initiated before the stars have coalesced to form a single compact object. Here we study such prompt detonations by means of numerical simulations, modeling the disruption and nucleosynthesis of the stars until the ejecta reach the coasting phase, and generating synthetic light curves and spectra. Three models are considered with primary masses 0.96 M_☉, 1.06 M_☉, and 1.20 M_☉. Of these, the 0.96 M_☉ dwarf merging with a 0.81 M_☉ companion, with an ⁵⁶Ni yield of 0.58 M_☉, is the most promising candidate for reproducing common SNe Ia. The more massive mergers produce unusually luminous SNe Ia with peak luminosities approaching those attributed to "super-Chandrasekhar" mass SNe Ia. While the synthetic light curves and spectra of some of the models resemble observed SNe Ia, the significant asymmetry of the ejecta leads to large orientation effects. The peak bolometric luminosity varies by more than a factor of two with the viewing angle, and the velocities of the spectral absorption features are lower when observed from angles where the light curve is brightest. The largest orientation effects are seen in the ultraviolet, where the flux varies by more than an order of magnitude. The set of three models roughly obeys a width–luminosity relation, with the brighter light curves declining more slowly in the B band. Spectral features due to unburned carbon from the secondary star are also seen in some cases.

We have investigated a supra-arcade structure associated with an M1.6 flare, which occurred on the south-east limb on 2010 November 4. It is observed in EUV with the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory, microwaves at 17 and 34 GHz with the Nobeyama Radioheliograph (NoRH), and soft X-rays of 8–20 keV with RHESSI. Interestingly, we found exceptional properties of the supra-arcade thermal plasma from the AIA 131 Å and the NoRH: (1) plasma upflows along large coronal loops and (2) enhancing microwave emission. RHESSI detected two soft X-ray sources, a broad one in the middle of the supra-arcade structure and a bright one just above the flare-arcade. We estimated the number density and thermal energy for these two source regions during the decay phase of the flare. In the supra-arcade source, we found that there were increases of the thermal energy and the density at the early and last stages, respectively. On the contrary, the density and thermal energy of the source on the top of the flare-arcade decreases throughout. The observed upflows imply that there is continuous energy supply into the supra-arcade structure from below during the decay phase of the flare. It is hard to explain by the standard flare model in which the energy release site is located high in the corona. Thus, we suggest that a potential candidate of the energy source for the hot supra-arcade structure is the flare-arcade, which has exhibited a predominant emission throughout.

In an effort to provide an accurate spectroscopic characterization of oxirane, state-of-the-art computational methods and approaches have been employed to determine highly accurate fundamental vibrational frequencies and rotational parameters. Available experimental data were used to assess the reliability of our computations, and an accuracy on average of 10 cm⁻¹ for fundamental transitions as well as overtones and combination bands has been pointed out. Moving to rotational spectroscopy, relative discrepancies of 0.1%, 2%–3%, and 3%–4% were observed for rotational, quartic, and sextic centrifugal-distortion constants, respectively. We are therefore confident that the highly accurate spectroscopic data provided herein can be useful for identification of oxirane in Titan's atmosphere and the assignment of unidentified infrared bands. Since oxirane was already observed in the interstellar medium and some astronomical objects are characterized by very high D/H ratios, we also considered the accurate determination of the spectroscopic parameters for the mono-deuterated species, oxirane-d1. For the latter, an empirical scaling procedure allowed us to improve our computed data and to provide predictions for rotational transitions with a relative accuracy of about 0.02% (i.e., an uncertainty of about 40 MHz for a transition lying at 200 GHz).

Spitzer UltRa Faint SUrvey Program is a joint Spitzer and Hubble Space Telescope Exploration Science program using 10 galaxy clusters as cosmic telescopes to study z ≳ 7 galaxies at intrinsically lower luminosities, enabled by gravitational lensing, than blank field surveys of the same exposure time. Our main goal is to measure stellar masses and ages of these galaxies, which are the most likely sources of the ionizing photons that drive reionization. Accurate knowledge of the star formation density and star formation history at this epoch is necessary to determine whether these galaxies indeed reionized the universe. Determination of the stellar masses and ages requires measuring rest-frame optical light, which only Spitzer can probe for sources at z ≳ 7, for a large enough sample of typical galaxies. Our program consists of 550 hr of Spitzer/IRAC imaging covering 10 galaxy clusters with very well-known mass distributions, making them extremely precise cosmic telescopes. We combine our data with archival observations to obtain mosaics with ∼30 hr exposure time in both 3.6 μm and 4.5 μm in the central 4' × 4' field and ∼15 hr in the flanking fields. This results in 3σ sensitivity limits of ∼26.6 and ∼26.2 AB magnitudes for the central field in the IRAC 3.6 and 4.5 μm bands, respectively. To illustrate the survey strategy and characteristics we introduce the sample, present the details of the data reduction and demonstrate that these data are sufficient for in-depth studies of z ≳ 7 sources (using a z = 9.5 galaxy behind MACS J1149.5+2223 as an example). For the first cluster of the survey (the Bullet Cluster) we have released all high-level data mosaics and IRAC empirical point-spread function models. In the future we plan to release these data products for the entire survey.

We recently reported extremely rapid changes in chromospheric fine structure observed using the IBIS instrument in the red wing of Hα. Here, we examine data obtained during the same observing run (2010 August 7), of a mature active region NOAA 11094. We analyze more IBIS data including wavelength scans and data from the Solar Dynamics Observatory, all from within a 30 minute interval. Using a slab radiative transfer model, we investigate the physical nature of fibrils in terms of tube-like versus sheet-like structures. Principal Component Analysis shows that the very rapid Hα variations in the line wings depend mostly on changes of line width and line shift, but for Ca ii 854.2 the variations are dominated by changes in column densities. The tube model must be rejected for a small but significant class of fibrils undergoing very rapid changes. If our wing data arise from the same structures leading to "type II spicules," our analysis calls into question much recent work. Instead, the data do not reject the hypothesis that some fibrils are optical superpositions of plasma collected into sheets. We review how Parker's theory of tangential discontinuities naturally leads to plasma collecting into sheets, and show that the sheet picture is falsifiable. Chromospheric fine structures seem to be populated by both tubes and sheets. We assess the merits of spectral imaging versus slit spectroscopy for future studies.

Magnetic Rayleigh–Taylor (MRT) instabilities may play a relevant role in many astrophysical problems. In this work the effect of magnetic shear on the growth rate of the MRT instability is investigated. The eigenmodes of an interface and a slab model under the presence of gravity are analytically calculated assuming that the orientation of the magnetic field changes in the equilibrium, i.e., there is magnetic shear. We solve the linearized magnetohydrodynamic equations in the incompressible regime. We find that the growth rate is bounded under the presence of magnetic shear. We have derived simple analytical expressions for the maximum growth rate, corresponding to the most unstable mode of the system. These expressions provide the explicit dependence of the growth rate on the various equilibrium parameters. For small angles the growth time is linearly proportional to the shear angle, and in this regime the single interface problem and the slab problem tend to the same result. On the contrary, in the limit of large angles and for the interface problem the growth time is essentially independent of the shear angle. In this regime we have also been able to calculate an approximate expression for the growth time for the slab configuration. Magnetic shear can have a strong effect on the growth rates of the instability. As an application of the results found in this paper we have indirectly determined the shear angle in solar prominence threads using their lifetimes and the estimation of the Alfvén speed of the structure.

We derive physical properties of 10 submillimeter galaxies located in the CANDELS coverage of the GOODS-S field. The galaxies were first identified as submillimeter sources with the LABOCA bolometer and subsequently targeted for 870 μm continuum observation with ALMA. The high angular resolution of the ALMA imaging allows secure counterparts to be identified in the CANDELS multiband data set. The CANDELS data provide deep photometric data from UV through near-infrared wavelengths. Using synthetic spectral energy distributions, we derive photometric redshifts, stellar masses, extinction, ages, and the star formation history. The redshift range is z = 1.65–4.76, with two of the galaxies located at z > 4. Two submillimeter galaxy (SMG) counterparts have stellar masses 2–3 orders of magnitude lower than the rest. The remaining SMG counterparts have stellar masses around 1 × 10¹¹ M_☉. The stellar population in the SMGs is typically older than the expected duration of the submillimeter phase, suggesting that the star formation history of SMGs is more complex than a single burst. Non-parametric morphology indices suggest that the SMG counterparts are among the most asymmetric systems compared with galaxies of the same stellar mass and redshift. The Hubble Space Telescope images show that three of the SMGs are associated with ongoing mergers. The remaining counterparts are isolated. Estimating the dust and molecular gas mass from the submillimeter fluxes, and comparing with our stellar masses shows that the gas mass fraction of SMGs is ∼28% and that the final stellar mass is likely to be ∼(1–2) × 10¹¹ M_☉.

Recent observations by the Fermi satellite suggest that a photosphere emission component is contributing to the observed spectrum of many gamma-ray bursts (GRBs). One important question is whether the photosphere component can interpret the typical "Band" function of GRBs with a typical low energy photon spectral index α ∼ −1. We perform a detailed study of the photosphere emission spectrum by progressively introducing several physical ingredients previously not fully incorporated, including the probability distribution of the location of a dynamically evolving photosphere, superposition of emission from an equal arrival time "volume" in a continuous wind, the evolution of optical depth of a wind with finite but evolving outer boundary, as well as the effect of different top-hat wind luminosity (L_w) profiles. By assuming a comoving blackbody spectrum emerging from the photosphere, we find that for an outflow with a constant or increasing L_w, the low-energy spectrum below the peak energy (E_p), can be modified to F_ν ∼ ν^1.5 (α ∼ +0.5). A softer (−1 < α < +0.5) or flat (α = −1) spectrum can be obtained during the L_w decreasing phase or high-latitude-emission-dominated phase. We also study the evolution of E_p as a function of wind and photosphere luminosity in this photosphere model. An E_p − L tracking pattern can be reproduced if a certain positive dependence between the dimensionless entropy η and L_w is introduced. However, the hard-to-soft evolution pattern cannot be reproduced unless a contrived condition is invoked. In order to interpret the Band spectrum, a more complicated photosphere model or a different energy dissipation and radiation mechanism is needed.

GRB 120326A is an unusual gamma-ray burst (GRB) that has a long plateau and a very late rebrightening in both X-ray and optical bands. The similar behavior of the optical and X-ray light curves suggests that they may share a common origin. The long plateau starts at several hundred seconds and ends at tens of thousands of seconds, and the peak time of the late rebrightening is about 30,000 s. We analyze the energy injection model by means of numerical and analytical solutions, considering both the wind environment and the interstellar medium environment for GRB afterglows. We particularly study the influence of the injection starting time, ending time, stellar wind density (or density of the circumburst environment), and injection luminosity on the shape of the afterglow light curves, respectively. In the wind model, we find that the light curve is largely affected by the parameters and that there is a "bump" in the late stage. In the wind environment, we found that the longer the energy is injected, the more obvious the rebrightening will be. We also find that the peak time of the bump is determined by the stellar wind density. We use the late continuous injection model to interpret the unusual afterglow of GRB 120326A. The model fits the observational data well; however, we find that the timescale of the injection must be higher than 10,000 s, which implies that the timescale of the central engine activity must also be more than 10,000 s. This information can give useful constraints on the central engines of GRBs—we consider a newborn millisecond pulsar with a strong magnetic field to be the central engine. On the other hand, our results suggest that the circumburst environment of GRB 120326A is very likely a stellar wind.

We report on PTF1 J191905.19+481506.2, a newly discovered, partially eclipsing, outbursting AM CVn system found in the Palomar Transient Factory synoptic survey. This is only the second known eclipsing AM CVn system. We use high-speed photometric observations and phase-resolved spectroscopy to establish an orbital period of 22.4559(3) minutes. We also present a long-term light curve and report on the normal and super-outbursts regularly seen in this system, including a super-outburst recurrence time of 36.8(4) days. We use the presence of the eclipse to place upper and lower limits on the inclination of the system and discuss the number of known eclipsing AM CVn systems versus what would be expected.

Observations indicate that most massive galaxies contain a supermassive black hole, and theoretical studies suggest that when such galaxies have a major merger, the central black holes will form a binary and eventually coalesce. Here we discuss two spectral signatures of such binaries that may help distinguish them from ordinary active galactic nuclei. These signatures are expected when the mass ratio between the holes is not extreme and the system is fed by a circumbinary disk. One such signature is a notch in the thermal continuum that has been predicted by other authors; we point out that it should be accompanied by a spectral revival at shorter wavelengths and also discuss its dependence on binary properties such as mass, mass ratio, and separation. In particular, we note that the wavelength λ_n at which the notch occurs depends on these three parameters in such a way as to make the number of systems displaying these notches $\propto \lambda _n^{16/3}$; longer wavelength searches are therefore strongly favored. A second signature, first discussed here, is hard X-ray emission with a Wien-like spectrum at a characteristic temperature ∼100 keV produced by Compton cooling of the shock generated when streams from the circumbinary disk hit the accretion disks around the individual black holes. We investigate the observability of both signatures. The hard X-ray signal may be particularly valuable as it can provide an indicator of black hole merger a few decades in advance of the event.

The secular dynamical evolution of a hierarchical three-body system in which a distant third object orbits around a binary has been studied extensively, demonstrating that the inner orbit can undergo large eccentricity and inclination oscillations. It was shown before that starting with a circular inner orbit, large mutual inclination (40°–140°) can produce long timescale modulations that drive the eccentricity to extremely large values and can flip the orbit. Here, we demonstrate that starting with an almost coplanar configuration, for eccentric inner and outer orbits, the eccentricity of the inner orbit can still be excited to high values, and the orbit can flip by ∼180°, rolling over its major axis. The ∼180° flip criterion and the flip timescale are described by simple analytic expressions that depend on the initial orbital parameters. With tidal dissipation, this mechanism can produce counter-orbiting exoplanetary systems. In addition, we also show that this mechanism has the potential to enhance the tidal disruption or collision rates for different systems. Furthermore, we explore the entire e₁ and i₀ parameter space that can produce flips.

Tentative evidence for statistical anisotropy in the Wilkinson Microwave Anisotropy Probe data was alleged to be due to "insufficient handling of beam asymmetries." In this paper, we investigate this issue and develop a method to estimate the shape of the inflight effective beam, particularly the asymmetry and azimuthal orientation. We divide the whole map into square patches and exploit the information in the Fourier space. For patches containing bright extragalactic point sources, we can directly estimate their shapes, from which the inflight effective beam can be estimated. For those without, we estimate the pattern from iso-power contours in two-dimensional Fourier space. We show that the inflight effective beam convolving the signal is indeed non-symmetric for most of the sky, and it is not randomly oriented. Around the ecliptic poles, however, the asymmetry is smaller due to the averaging effect from different orientations of the beam from the scan strategy. The orientations of the effective beam with significant asymmetry are parallel to the lines of ecliptic longitude. In the foreground-cleaned Internal Linear Combination map, however, the systematics caused by beam effect is significantly lessened.

We report on XMM-Newton and Chandra observations of the Galactic supernova remnant candidate G38.7-1.4, together with complementary radio, infrared, and γ-ray data. An approximately elliptical X-ray structure is found to be well correlated with a radio shell as seen by the Very Large Array. The X-ray spectrum of G38.7-1.4 can be well described by an absorbed collisional ionization equilibrium plasma model, which suggests the plasma is shock heated. Based on the morphology and the spectral behavior, we suggest that G38.7-1.4 is indeed a supernova remnant belonging to a mix-morphology category.

We present the results of searches for gravitational waves from a large selection of pulsars using data from the most recent science runs (S6, VSR2 and VSR4) of the initial generation of interferometric gravitational wave detectors LIGO (Laser Interferometric Gravitational-wave Observatory) and Virgo. We do not see evidence for gravitational wave emission from any of the targeted sources but produce upper limits on the emission amplitude. We highlight the results from seven young pulsars with large spin-down luminosities. We reach within a factor of five of the canonical spin-down limit for all seven of these, whilst for the Crab and Vela pulsars we further surpass their spin-down limits. We present new or updated limits for 172 other pulsars (including both young and millisecond pulsars). Now that the detectors are undergoing major upgrades, and, for completeness, we bring together all of the most up-to-date results from all pulsars searched for during the operations of the first-generation LIGO, Virgo and GEO600 detectors. This gives a total of 195 pulsars including the most recent results described in this paper.

Three luminous blue variables (LBVs) are located in and near the Quintuplet Cluster at the Galactic center: the Pistol Star, G0.120-0.048, and qF362. We present imaging at 19, 25, 31, and 37 μm of the region containing these three LBVs, obtained with SOFIA using FORCAST. We argue that Pistol and G0.120-0.048 are identical "twins" that exhibit contrasting nebulae due to the external influence of their different environments. Our images reveal the asymmetric, compressed shell of hot dust surrounding the Pistol Star and provide the first detection of the thermal emission from the symmetric, hot dust envelope surrounding G0.120-0.048. However, no detection of hot dust associated with qF362 is made. Dust and gas composing the Pistol nebula are primarily heated and ionized by the nearby Quintuplet Cluster stars. The northern region of the Pistol nebula is decelerated due to the interaction with the high-velocity (2000 km s⁻¹) winds from adjacent Wolf–Rayet Carbon (WC) stars. From fits to the spectral energy distribution (SED) of the Pistol nebula with the DustEM code we determine that the Pistol nebula is composed of a distribution of very small, transiently heated grains (10 to  ∼ 35 Å) having a total dust mass of 0.03 M_☉, and that it exhibits a gradient of decreasing grain size from south to north due to differential sputtering by the winds from the WC stars. The total IR luminosity of the Pistol nebula is 5.2 × 10⁵ L_☉. Dust in the G0.120-0.048 nebula is primarily heated by the central star; however, the nebular gas is ionized externally by the Arches Cluster. Unlike the Pistol nebula, the G0.120-0.048 nebula is freely expanding into the surrounding medium. A grain size distribution identical to that of the non-sputtered region of the Pistol nebula satisfies the constraints placed on the G0.120-0.048 nebula from DustEM model fits to its SED and implies a total dust mass of 0.021 M_☉. The total IR luminosity of the G0.120-0.048 nebula is ∼10⁵ L_☉. From Paschen-α and 6 cm observations we determine a total gas mass of 9.3 M_☉ and 6.2 M_☉ for the Pistol and G0.120-0.048 nebulae, respectively. Given the independent dust and gas mass estimates we find that the Pistol and G0.120-0.048 nebulae exhibit similar gas-to-dust mass ratios of $310^{+77}_{-52}$ and $293^{+73}_{-101}$, respectively. Both nebulae share identical size scales (∼0.7 pc) which suggests that they have similar dynamical timescales of ∼10⁴ yr, assuming a shell expansion velocity of v_exp = 60 km s⁻¹.

We present the first Chandra and Swift X-ray study of the spiral galaxy NGC 4088 and its ultraluminous X-ray source (ULX N4088–X1). We also report very long baseline interferometry (VLBI) observations at 1.6 and 5 GHz performed quasi-simultaneously with the Swift and Chandra observations, respectively. Fifteen X-ray sources are detected by Chandra within the D25 ellipse of NGC 4088, from which we derive the X-ray luminosity function (XLF) of this galaxy. We find the XLF is very similar to those of star-forming galaxies and estimate a star-formation rate of 4.5 M_☉ yr⁻¹. The Chandra detection of the ULX yields its most accurate X-ray position, which is spatially coincident with compact radio emission at 1.6 GHz. The ULX Chandra X-ray luminosity, L_{0.2–10.0 keV} = 3.4 × 10³⁹ erg s⁻¹, indicates that N4088–X1 could be located at the high-luminosity end of the high-mass X-ray binary (HMXB) population of NGC 4088. The estimates of the black hole (BH) mass and ratio of radio to X-ray luminosity of N4088–X1 rule out a supermassive BH nature. The Swift X-ray spectrum of N4088–X1 is best described by a thermal Comptonization model and presents a statistically significant high-energy cutoff. We conclude that N4088–X1 is most likely a stellar remnant BH in an HMXB, probably fed by Roche lobe overflow, residing in a super-Eddington ultraluminous state. The 1.6 GHz VLBI source is consistent with radio emission from possible ballistic jet ejections in this state.

We perform a systematic study of the dynamics of dust particles in protoplanetary disks with embedded planets using global two-dimensional and three-dimensional inviscid hydrodynamic simulations. Lagrangian particles have been implemented into the magnetohydrodynamic code Athena with cylindrical coordinates. We find two distinct outcomes depending on the mass of the embedded planet. In the presence of a low-mass planet (8 M_⊕), two narrow gaps start to open in the gas on each side of the planet where the density waves are shocked. These shallow gaps can dramatically affect particle drift speed and cause significant, roughly axisymmetric dust depletion. On the other hand, a more massive planet (>0.1 M_J) carves out a deeper gap with sharp edges, which are subject to Rossby wave instability leading to vortex formation. Particles with a wide range of sizes (0.02 < Ωt_s < 20) are trapped and settle to the midplane in the vortex, with the strongest concentration for particles with Ωt_s ∼ 1. The dust concentration is highly elongated in the ϕ direction, and can be as wide as four disk scale heights in the radial direction. Dust surface density inside the vortex can be increased by more than a factor of 10² in a very non-axisymmetric fashion. For very big particles (Ωt_s ≫ 1) we find strong eccentricity excitation, in particular around the planet and in the vicinity of the mean motion resonances, facilitating gap openings there. Our results imply that in weakly turbulent protoplanetary disk regions (e.g., the "dead zone") dust particles with a very wide range of sizes can be trapped at gap edges and inside vortices induced by planets with M_p < M_J, potentially accelerating planetesimal and planet formation there, and giving rise to distinctive features that can be probed by ALMA and the Extended Very Large Array.

Three-dimensional (3D) simulations of core-collapse supernovae (CCSNe) are granting new insight into the as-yet-uncertain mechanism that drives successful explosions. While there is still debate about whether explosions are obtained more easily in 3D than in 2D, it is undeniable that there exist qualitative and quantitative differences between the results of 3D and 2D simulations. We present an extensive set of high-resolution 1D, 2D, and 3D CCSN simulations with multispecies neutrino leakage carried out in two different progenitors. Our simulations confirm the results of Couch indicating that 2D explodes more readily than 3D. We argue that this is due to the inadequacies of 2D to accurately capture important aspects of the 3D dynamics. We find that without artificially enhancing the neutrino heating rate, we do not obtain explosions in 3D. We examine the development of neutrino-driven convection and the standing accretion shock instability (SASI) and find that, in separate regimes, either instability can dominate. We find evidence for growth of the SASI for both 15 M_☉ and 27 M_☉ progenitors; however, it is weaker in 3D exploding models. The growth rate of both instabilities is artificially enhanced along the symmetry axis in 2D as compared with our axis-free 3D Cartesian simulations. Our work highlights the growing consensus that CCSNe must be studied in 3D if we hope to solve the mystery of how the explosions are powered.

The early phase of the hydrodynamic evolution following the collision of two stars is analyzed. Two strong shocks propagate from the contact surface and move toward the center of each star at a velocity that is a small fraction of the velocity of the approaching stars. The shocked region near the contact surface has a planar symmetry and a uniform pressure. The density vanishes at the (Lagrangian) surface of contact, and the speed of sound diverges there. The temperature, however, reaches a finite value, since as the density vanishes, the finite pressure is radiation dominated. For carbon–oxygen white dwarf (CO WD) collisions, this temperature is too low for any appreciable nuclear burning shortly after the collision, which allows for a significant fraction of the mass to be highly compressed to the density required for efficient ⁵⁶Ni production in the detonation wave that follows. This property is crucial for the viability of collisions of typical CO WD as progenitors of type Ia supernovae, since otherwise only massive (>0.9 M_☉) CO WDs would have led to such explosions (as required by all other progenitor models). The divergence of the speed of sound limits numerical studies of stellar collisions, as it makes convergence tests exceedingly expensive unless dedicated schemes are used. We provide a new one-dimensional Lagrangian numerical scheme to achieve this. A self-similar planar solution is derived for zero-impact parameter collisions between two identical stars, under some simplifying assumptions (including a power-law density profile), which is the planar version of previous piston problems that were studied in cylindrical and spherical symmetries.

The C to O ratio is a crucial determinant of the chemical properties of planets. The recent observation of WASP 12b, a giant planet with a C/O value larger than that estimated for its host star, poses a conundrum for understanding the origin of this elemental ratio in any given planetary system. In this paper, we propose a mechanism for enhancing the value of C/O in the disk through the transport and distribution of volatiles. We construct a model that computes the abundances of major C- and O-bearing volatiles under the influence of gas drag, sublimation, vapor diffusion, condensation, and coagulation in a multi-iceline 1+1D protoplanetary disk. We find a gradual depletion in water and carbon monoxide vapors inside the water's iceline, with carbon monoxide depleting slower than water. This effect increases the gaseous C/O and decreases the C/H ratio in this region to values similar to those found in WASP 12b's day side atmosphere. Giant planets whose envelopes were accreted inside the water's iceline should then display C/O values larger than those of their parent stars, making them members of the class of so-called carbon-rich planets.

In this paper we search for distant massive companions to known transiting gas giant planets that may have influenced the dynamical evolution of these systems. We present new radial velocity observations for a sample of 51 planets obtained using the Keck HIRES instrument, and find statistically significant accelerations in fifteen systems. Six of these systems have no previously reported accelerations in the published literature: HAT-P-10, HAT-P-22, HAT-P-29, HAT-P-32, WASP-10, and XO-2. We combine our radial velocity fits with Keck NIRC2 adaptive optics (AO) imaging data to place constraints on the allowed masses and orbital periods of the companions responsible for the detected accelerations. The estimated masses of the companions range between 1–500 M_Jup, with orbital semi-major axes typically between 1–75 AU. A significant majority of the companions detected by our survey are constrained to have minimum masses comparable to or larger than those of the transiting planets in these systems, making them candidates for influencing the orbital evolution of the inner gas giant. We estimate a total occurrence rate of 51%  ±  10% for companions with masses between 1–13 M_Jup and orbital semi-major axes between 1–20 AU in our sample. We find no statistically significant difference between the frequency of companions to transiting planets with misaligned or eccentric orbits and those with well-aligned, circular orbits. We combine our expanded sample of radial velocity measurements with constraints from transit and secondary eclipse observations to provide improved measurements of the physical and orbital characteristics of all of the planets included in our survey.

Large-scale magnetic fields are key ingredients of magnetically driven disk accretion. We study how large-scale poloidal fields evolve in accretion disks, with the primary aim of quantifying the viability of magnetic accretion mechanisms in protoplanetary disks. We employ a kinematic mean-field model for poloidal field transport and focus on steady states where inward advection of a field balances with outward diffusion due to effective resistivities. We analytically derive the steady-state radial distribution of poloidal fields in highly conducting accretion disks. The analytic solution reveals an upper limit on the strength of large-scale vertical fields attainable in steady states. Any excess poloidal field will diffuse away within a finite time, and we demonstrate this with time-dependent numerical calculations of the mean-field equations. We apply this upper limit to large-scale vertical fields threading protoplanetary disks. We find that the maximum attainable strength is about 0.1 G at 1 AU, and about 1 mG at 10 AU from the central star. When combined with recent magnetic accretion models, the maximum field strength translates into the maximum steady-state accretion rate of ∼10⁻⁷ M_☉ yr⁻¹, in agreement with observations. We also find that the maximum field strength is ∼1 kG at the surface of the central star provided that the disk extends down to the stellar surface. This implies that any excess stellar poloidal field of strength ≳ kG can be transported to the surrounding disk. This might in part resolve the magnetic flux problem in star formation.

Even though the recently discovered high-magnification event MOA-2010-BLG-311 had complete coverage over its peak, confident planet detection did not happen due to extremely weak central perturbations (EWCPs, fractional deviations of ≲ 2%). For confident detection of planets in EWCP events, it is necessary to have both high cadence monitoring and high photometric accuracy better than those of current follow-up observation systems. The next-generation ground-based observation project, Korea Microlensing Telescope Network (KMTNet), satisfies these conditions. We estimate the probability of occurrence of EWCP events with fractional deviations of ⩽2% in high-magnification events and the efficiency of detecting planets in the EWCP events using the KMTNet. From this study, we find that the EWCP events occur with a frequency of >50% in the case of ≲ 100 M_E planets with separations of 0.2 AU ≲ d ≲ 20 AU. We find that for main-sequence and sub-giant source stars, ≳ 1 M_E planets in EWCP events with deviations ⩽2% can be detected with frequency >50% in a certain range that changes with the planet mass. However, it is difficult to detect planets in EWCP events of bright stars like giant stars because it is easy for KMTNet to be saturated around the peak of the events because of its constant exposure time. EWCP events are caused by close, intermediate, and wide planetary systems with low-mass planets and close and wide planetary systems with massive planets. Therefore, we expect that a much greater variety of planetary systems than those already detected, which are mostly intermediate planetary systems, regardless of the planet mass, will be significantly detected in the near future.

The expected anisotropy in the 1 to 10⁴ TeV energy range is calculated for Galactic cosmic rays with both anisotropy in the diffusion tensor and source discreteness taken into account. We find that if the sources are distributed radially (but with azimuthal symmetry) in proportion to Galactic pulsars, the expected anisotropy almost always exceeds the observational limits by one order of magnitude in the case of isotropic diffusion. If the radial diffusion is more than an order of magnitude smaller than the azimuthal diffusion rate, the radial gradient of the sources can be accommodated about 5% of the time. If the sources are concentrated in the spiral arms, then the anisotropy depends on our location between them, but in some spatial window, roughly equidistant from adjacent spiral arms, the observational constraints on anisotropy are obeyed roughly 20%–30% of the time for extremely anisotropic diffusion. The solar system is in that window less than 10% of the time, but it may be there now. Under the assumption of isotropic diffusion, nearby supernovae are found to produce a discreteness anisotropy that is nearly two orders of magnitude in excess of the observational limit if all supernovae are assumed to contribute equally with a source rate 1 in every 100 years.

The processes responsible for the broadband radiation of the young supernova remnant Cas A are explored by using a new code that is designed for a detailed treatment of the diffusive shock acceleration of particles in the nonlinear regime. The model is based on spherically symmetric hydrodynamic equations complemented with transport equations for relativistic particles. Electrons, protons, and the oxygen ions accelerated by forward and reverse shocks are included in the numerical calculations. We show that the available multi-wavelength observations in the radio, X-ray, and gamma-ray bands can be best explained by invoking particle acceleration by both forward and reversed shocks. Although the TeV gamma-ray observations can be interpreted by interactions of both accelerated electrons and protons/ions, the measurements by Fermi Large Area Telescope at energies below 1 GeV give a tentative preference to the hadronic origin of gamma-rays. Then, the acceleration efficiency in this source, despite the previous claims, should be very high; 25% of the explosion energy (or approximately 3 × 10⁵⁰ erg) should already be converted to cosmic rays, mainly by the forward shock. At the same time, the model calculations do not provide extension of the maximum energy of accelerated protons beyond 100 TeV. In this model, the acceleration of electrons is dominated by the reverse shock; the required 10⁴⁸ erg can be achieved under the assumption that the injection of electrons (positrons) is supported by the radioactive decay of ⁴⁴Ti.

Recent observations strongly suggest that the millisecond pulsar binary PSR J1023+0038 has developed an accretion disk since 2013 June. We present a multi-wavelength analysis of PSR J1023+0038, which reveals that (1) its gamma-rays suddenly brightened within a few days in 2013 June/July and has remained at a high gamma-ray state for several months; (2) both UV and X-ray fluxes have increased by roughly an order of magnitude; and (3) the spectral energy distribution has changed significantly after the gamma-ray sudden flux change. Time variabilities associated with UV and X-rays are on the order of 100–500 s and 50–100 s, respectively. Our model suggests that a newly formed accretion disk, due to the sudden increase of the stellar wind, could explain the changes of all these observed features. The increase of UV is emitted from the disk, and a new component in gamma-rays is produced by inverse Compton scattering between the new UV component and pulsar wind. The increase of X-rays results from the enhancement of injection pulsar wind energy into the intra-binary shock due to the increase of the stellar wind. We also predict that the radio pulses may be blocked by the evaporated winds from the disk, and the pulsar is still powered by rotation.

We extend our approach of modeling spectral energy distribution (SED) and light curves of blazars to include external Compton (EC) emission due to inverse Compton scattering of an external anisotropic target radiation field. We describe the time-dependent impact of such seed photon fields on the evolution of multifrequency emission and spectral variability of blazars using a multi-zone time-dependent leptonic jet model, with radiation feedback, in the internal shock model scenario. We calculate accurate EC-scattered high-energy (HE) spectra produced by relativistic electrons throughout the Thomson and Klein–Nishina regimes. We explore the effects of varying the contribution of (1) a thermal Shakura–Sunyaev accretion disk, (2) a spherically symmetric shell of broad-line clouds, the broad-line region (BLR), and (3) a hot infrared emitting dusty torus (DT), on the resultant seed photon fields. We let the system evolve to beyond the BLR and within the DT and study the manifestation of the varying target photon fields on the simulated SED and light curves of a typical blazar. The calculations of broadband spectra include effects of γ–γ absorption as γ-rays propagate through the photon pool present inside the jet due to synchrotron and inverse Compton processes, but neglect γ–γ absorption by the BLR and DT photon fields outside the jet. Thus, our account of γ–γ absorption is a lower limit to this effect. Here, we focus on studying the impact of parameters relevant for EC processes on HE emission of blazars.

The outskirts of galaxy clusters are continuously disturbed by mergers and gas infall along filaments, which in turn induce turbulent flow motions and shock waves. We examine the properties of shocks that form within r₂₀₀ in sample galaxy clusters from structure formation simulations. While most of these shocks are weak and inefficient accelerators of cosmic rays (CRs), there are a number of strong, energetic shocks which can produce large amounts of CR protons via diffusive shock acceleration. We show that the energetic shocks reside mostly in the outskirts and a substantial fraction of them are induced by infall of the warm-hot intergalactic medium from filaments. As a result, the radial profile of the CR pressure in the intracluster medium is expected to be broad, dropping off more slowly than that of the gas pressure, and might be even temporarily inverted, peaking in the outskirts. The volume-integrated momentum spectrum of CR protons inside r₂₀₀ has the power-law slope of 4.25–4.5, indicating that the average Mach number of the shocks of main CR production is in the range of <M_s > _CR ≈ 3–4. We suggest that some radio relics with relatively flat radio spectrum could be explained by primary electrons accelerated by energetic infall shocks with M_s ≳ 3 induced in the cluster outskirts.

We investigate the connection between the reionization era and the present-day universe by examining the mass reionization histories of z = 0 dark matter halos. In a 600³ Mpc³ volume, we combine a dark matter N-body simulation with a three-dimensional seminumerical reionization model. This tags each particle with a reionization redshift, so that individual present-day halos can be connected to their reionization histories and environments. We find that the vast majority of present-day halos with masses larger than ∼ few × 10¹¹ M_☉ reionize earlier than the rest of the universe. We also find significant halo-to-halo diversity in mass reionization histories, and find that in realistic inhomogeneous models, the material within a given halo is not expected to reionize at the same time. In particular, the scatter in reionization times within individual halos is typically larger than the scatter among halos. From our fiducial reionization model, we find that the typical 68% scatter in reionization times within halos is ∼115 Myr for 10^{12 ± 0.25} M_☉ halos, decreasing slightly to ∼95 Myr for 10^{15 ± 0.25} M_☉ halos. We find a mild correlation between reionization history and environment: halos with shorter reionization histories are typically in more clustered environments, with the strongest trend on a scale of ∼20 Mpc. Material in Milky Way mass halos with short reionization histories is preferentially reionized in relatively large H ii regions, implying reionization mostly by sources external to the progenitors of the present-day halo. We investigate the impact on our results of varying the reionization model parameters, which span a range of reionization scenarios with varying timing and morphology.

We present new Herschel observations of the (6,6) and (9,9) inversion transitions of the hydronium ion toward Sagittarius B2(N) and W31C. Sensitive observations toward Sagittarius B2(N) show that the high, ∼500 K, rotational temperatures characterizing the population of the highly excited metastable H₃O⁺ rotational levels are present over a wide range of velocities corresponding to the Sagittarius B2 envelope, as well as the foreground gas clouds between the Sun and the source. Observations of the same lines toward W31C, a line of sight that does not intersect the Central Molecular Zone but instead traces quiescent gas in the Galactic disk, also imply a high rotational temperature of ∼380 K, well in excess of the kinetic temperature of the diffuse Galactic interstellar medium. While it is plausible that some fraction of the molecular gas may be heated to such high temperatures in the active environment of the Galactic center, characterized by high X-ray and cosmic-ray fluxes, shocks, and high degree of turbulence, this is unlikely in the largely quiescent environment of the Galactic disk clouds. We suggest instead that the highly excited states of the hydronium ion are populated mainly by exoergic chemical formation processes and the temperature describing the rotational level population does not represent the physical temperature of the medium. The same arguments may be applicable to other symmetric top rotors, such as ammonia. This offers a simple explanation of the long-standing puzzle of the presence of a pervasive, hot molecular gas component in the central region of the Milky Way. Moreover, our observations suggest that this is a universal process not limited to the active environments associated with galactic nuclei.

We measure the extinction curve in the central 200 pc of M31 at mid-ultraviolet to near-infrared wavelengths (from 1928 Å to 1.5 μm), using Swift/UVOT and Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3)/Advanced Camera for Surveys (ACS) observations in 13 bands. Taking advantage of the high angular resolution of the HST/WFC3 and ACS detectors, we develop a method to simultaneously determine the relative extinction and the fraction of obscured starlight for five dusty complexes located in the circumnuclear region. The extinction curves of these clumps (R_V = 2.4–2.5) are steeper than the average Galactic one (R_V = 3.1), but are similar to optical and near-infrared curves recently measured toward the Galactic bulge (R_V ∼ 2.5). This similarity suggests that steep extinction curves may be common in the inner bulge of galaxies. In the ultraviolet, the extinction curves of these clumps are also unusual. We find that one dusty clump (size < 2 pc) exhibits a strong UV bump (extinction at 2175 Å), more than three standard deviation higher than that predicted by common models. Although the high stellar metallicity of the M31 bulge indicates that there are sufficient carbon and silicon to produce large dust grains, the grains may have been destroyed by supernova explosions or past activity of the central supermassive black hole, resulting in the observed steepened extinction curve.

We present evidence that recurrent spiral activity, long manifested in simulations of disk galaxies, results from the superposition of a few transient spiral modes. Each mode lasts between 5 and 10 rotations at its corotation radius where its amplitude is greatest. The scattering of stars as each wave decays takes place over narrow ranges of angular momentum, causing abrupt changes to the impedance of the disk to subsequent traveling waves. Partial reflections of waves at these newly created features allows new standing-wave instabilities to appear that saturate and decay in their turn, scattering particles at new locations, creating a recurring cycle. The spiral activity causes the general level of random motion to rise, gradually decreasing the ability of the disk to support further activity unless the disk contains a dissipative gas component from which stars form on near-circular orbits. We also show that this interpretation is consistent with the behavior reported in other recent simulations with low-mass disks.

Motivated by recent observations of distinct parallel and perpendicular signatures in magnetic helicity measurements segregated by wave period and angle between the local magnetic field and the solar wind velocity, this paper undertakes a comparison of three intervals of Ulysses data with synthetic time series generated from a physically motivated turbulence model. From these comparisons, it is hypothesized that the observed signatures result from a perpendicular cascade of Alfvénic fluctuations and a local, non-turbulent population of ion-cyclotron or whistler waves generated by temperature anisotropy instabilities. By constraining the model's free parameters through comparison to in situ data, it is found that, on average, ∼95% of the power near dissipative scales is contained in a perpendicular Alfvénic cascade and that the parallel fluctuations are propagating nearly unidirectionally. The effects of aliasing on magnetic helicity measurements are considered and shown to be significant near the Nyquist frequency.

We study the nonlocal nonlinear coupling and generation of kinetic Alfvén waves (KAWs) and kinetic slow waves (KSWs) by magnetohydrodynamic Alfvén waves (MHD AWs) in conditions typical for the solar wind in the inner heliosphere. This cross-scale process provides an alternative to the turbulent energy cascade passing through many intermediate scales. The nonlinearities we study are proportional to the scalar products of wave vectors and hence are called "scalar" ones. Despite the strong Landau damping of kinetic waves, we found fast growing KAWs and KSWs at perpendicular wavelengths close to the ion gyroradius. Using the parametric decay formalism, we investigate two independent decay channels for the pump AW: forward decay (involving co-propagating product waves) and backward decay (involving counter-propagating product waves). The growth rate of the forward decay is typically 0.05 but can exceed 0.1 of the pump wave frequency. The resulting spectral transport is nonlocal and anisotropic, sharply increasing perpendicular wavenumbers but not parallel ones. AWs and KAWs propagating against the pump AW grow with about the same rate and contribute to the sunward wave flux in the solar wind. Our results suggest that the nonlocal decay of MHD AWs into KAWs and KSWs is a robust mechanism for the cross-scale spectral transport of the wave energy from MHD to dissipative kinetic scales in the solar wind and similar media.

We quantitatively assess, by means of comprehensive numerical simulations, the ability of broadband photometric surveys to recover the broad emission line region (BLR) size in quasars under various observing conditions and for a wide range of object properties. Focusing on the general characteristics of the Large Synoptic Survey Telescope (LSST), we find that the slope of the size–luminosity relation for the BLR in quasars can be determined with unprecedented accuracy, on the order of a few percent, over a broad luminosity range and out to z ∼ 3. In particular, major emission lines for which the BLR size can be reliably measured with LSST include Hα, Mg ii λ2799, C iii] λ1909, C iv λ1549, and Lyα, amounting to a total of ≳ 10⁵ time-delay measurements for all transitions. Combined with an estimate for the emission line velocity dispersion, upcoming photometric surveys will facilitate the estimation of black hole masses in active galactic nuclei over a broad range of luminosities and redshifts, allow for refined calibrations of BLR size–luminosity–redshift relations in different transitions, as well as lead to more reliable cross-calibration with other black hole mass estimation techniques.

The beaming pattern of radiation emitted by a relativistically moving source, such as jets in microquasars, active galactic nuclei, and gamma-ray bursts, is a key issue for understanding acceleration and radiation processes in these objects. In this paper, we introduce a formalism based on a solution of the photon transfer equation to study the beaming patterns for emission produced by electrons accelerated in the jet and the upscattering photons of low-energy radiation fields of external origin (the so-called external Compton scenario). The formalism allows us to treat non-stationary, non-homogeneous, and anisotropic distributions of electrons, but assuming homogeneous/isotropic and non-variable target photon fields. We demonstrate the non-negligible impact of the anisotropy in the electron distribution on angular and spectral characteristics of the EC radiation.

Black hole binaries exhibit a wide range of variability phenomena, from large-scale state changes to broadband noise and quasi-periodic oscillations, but the physical nature of much of this variability is poorly understood. We examine the variability properties of three GRMHD simulations of thin accretion disks around black holes of varying spin, producing light curves and power spectra as would be seen by observers. We find that the simulated power spectra show a broad feature at high frequency, which increases in amplitude with the inclination of the observer. We show that this high-frequency feature is a product of the Doppler effect and that its location is a function of the mass and spin of the black hole. This Doppler feature demonstrates that power spectral properties of the accretion disk can be tied to, and potentially used to determine, physical properties of the black hole.

We have measured the line-of-sight velocity distribution from integrated stellar light at two points in the outer halo of M87 (NGC 4486), the second-rank galaxy in the Virgo Cluster. The data were taken at R = 480'' (∼41.5 kpc) and R = 526'' (∼45.5 kpc) along the SE major axis. The second moment for a non-parametric estimate of the full velocity distribution is 420 ± 23 km s⁻¹ and 577 ± 35 km s⁻¹, respectively. There is intriguing evidence in the velocity profiles for two kinematically distinct stellar components at the position of our pointing. Under this assumption, we employ a two-Gaussian decomposition and find the primary Gaussian having rest velocities equal to M87 (consistent with zero rotation) and second moments of 383 ± 32 km s⁻¹ and 446 ± 43 km s⁻¹, respectively. The asymmetry seen in the velocity profiles suggests that the stellar halo of M87 is not in a relaxed state and confuses a clean dynamical interpretation. That said, either measurement (full or two component model) shows a rising velocity dispersion at large radii, consistent with previous integrated light measurements, yet significantly higher than globular cluster measurements at comparable radial positions. These integrated light measurements at large radii, and the stark contrast they make to the measurements of other kinematic tracers, highlight the rich kinematic complexity of environments like the center of the Virgo Cluster and the need for caution when interpreting kinematic measurements from various dynamical tracers.

We present RingFinder, a tool for finding galaxy-scale strong gravitational lenses in multi-band imaging data. By construction, the method is sensitive to configurations involving a massive foreground ETG and a faint, background, blue source. RingFinder detects the presence of blue residuals embedded in an otherwise smooth red light distribution by difference imaging in two bands. The method is automated for efficient application to current and future surveys, having originally been designed for the 150 deg² Canada–France–Hawaii Telescope Legacy Survey (CFHTLS). We describe each of the steps of RingFinder. We then carry out extensive simulations to assess completeness and purity. For sources with magnification μ > 4, RingFinder reaches 42% (25%) completeness and 29% (86%) purity before (after) visual inspection. The completeness of RingFinder is substantially improved in the particular range of Einstein radii 08 ⩽ R_Ein ⩽ 20 and lensed images brighter than g = 22.5, where it can be as high as ∼70%. RingFinder does not introduce any significant bias in the source or deflector population. We conclude by presenting the final catalog of RingFinder CFHTLS galaxy-scale strong lens candidates. Additional information obtained with Hubble Space Telescope and Keck adaptive optics high-resolution imaging, and with Keck and Very Large Telescope spectroscopy, is used to assess the validity of our classification and measure the redshift of the foreground and the background objects. From an initial sample of 640,000 ETGs, RingFinder returns 2500 candidates, which we further reduce by visual inspection to 330 candidates. We confirm 33 new gravitational lenses from the main sample of candidates, plus an additional 16 systems taken from earlier versions of RingFinder. First applications are presented in the Strong Lensing Legacy Survey galaxy-scale lens sample paper series.

Testing whether close-in massive exoplanets (hot Jupiters) can enhance the stellar activity in their host primary is crucial for the models of stellar and planetary evolution. Among systems with hot Jupiters, HD 189733 is one of the best studied because of its proximity, strong activity, and the presence of a transiting planet, which allows transmission spectroscopy and a measure of the planetary radius and its density. Here we report on the X-ray activity of the primary star, HD 189733 A, using a new XMM-Newton observation and a comparison with the previous X-ray observations. The spectrum in the quiescent intervals is described by two temperatures at 0.2 keV and 0.7 keV, while during the flares a third component at 0.9 keV is detected. With the analysis of the summed Reflection Grating Spectrometer spectra, we obtain estimates of the electron density in the range n_e = (1.6–13) × 10¹⁰ cm⁻³, and thus the corona of HD 189733 A appears denser than the solar one. For the third time, we observe a large flare that occurred just after the eclipse of the planet. Together with the flares observed in 2009 and 2011, the events are restricted to a small planetary phase range of ϕ = 0.55–0.65. Although we do not find conclusive evidence of a significant excess of flares after the secondary transits, we suggest that the planet might trigger such flares when it passes close to the locally high magnetic field of the underlying star at particular combinations of stellar rotational phases and orbital planetary phases. For the most recent flares, a wavelet analysis of the light curve suggests a loop of length of four stellar radii at the location of the bright flare, and a local magnetic field of the order of 40–100 G, in agreement with the global field measured in other studies. The loop size suggests an interaction of magnetic nature between planet and star, separated by only ∼8R_*. The X-ray variability of HD 189733 A is larger than the variability of field stars and young Pleiades of similar spectral type and X-ray luminosity. We also detect the stellar companion (HD 189733 B, ∼12'' from the primary star) in this XMM-Newton observation. Its very low X-ray luminosity (L_X = 3.4 × 10²⁶ erg s⁻¹) confirms the old age of this star and of the binary system. The high activity of the primary star is best explained by a transfer of angular momentum from the planet to the star.

Sakurai's object (V4334 Sgr), the final flash object discovered in the mid-1990s, underwent rapid cooling during the first decade of the 21st century becoming as faint as K ∼ 25. This stage of evolution has ceased. Between observations in 2010 September and 2013 April V4334 Sgr brightened >2 mag to K = 14.2 and the effective temperature increased to ∼590 K. AO images show a central source and two extended globules defining a 13° position angle. The globules span a spatial extent of ∼03 in 2013. This spatial extent is consistent with sizes derived from spectral energy distributions taken over the previous decade and a debris cloud expanding at 0.055 mas d⁻¹ since late 1998. Near-simultaneous 0.85–2.5 μm spectra reveal helium lines attributed to a wind-interaction shock. The He i 1.0830 μm emission has a spectral width of ∼1000 km s⁻¹ and a spatial extent of ∼14. The helium shell is fragmented, spatially asymmetric, and five times larger than the dust debris cloud. [C i] and [N i] forbidden lines are present in the 1 μm region spectrum. The forbidden line spectrum is similar to that of proto-planetary nebulae. The [C i] 9850 Å line is spatially extended. The expansion velocity and change of angular size limit the distance to 2.1–3.7 kpc.

We present an analysis of optical spectra from the globular cluster RZ2109 in NGC 4472, which hosts the first unambiguous globular cluster black hole. We use these spectra to determine the elemental composition of the emission line system associated with this source, and to constrain the age and metallicity of the host globular cluster. For the emission line system of RZ2109, our analysis indicates the [O iii] λ5007 equivalent width is 33.82 ± 0.39 Å and the Hβ equivalent width is 0.32 ± 0.32 Å, producing a formal [O iii] λ5007/Hβ emission line ratio of 106 for a 3200 km s⁻¹ measurement aperture covering the full velocity width of the [O iii] λ5007 line. Within a narrower 600 km s⁻¹ aperture covering the highest luminosity velocity structure in the line complex, we find [O iii] λ5007/Hβ = 62. The measured [O iii] λ5007/Hβ ratios are significantly higher than can be produced in radiative models of the emission line region with solar composition, and the confidence interval limits exclude all but models which have gas masses much larger than those for a single star. Therefore, we conclude that the region from which the [O iii] λ5007 emission originates is hydrogen-depleted relative to solar composition gas. This finding is consistent with emission from an accretion-powered outflow driven by a hydrogen-depleted donor star, such as a white dwarf, being accreted onto a black hole.

We probe the structure and composition of the atmospheres of five hot Jupiter exoplanets using the Hubble Space Telescope Wide Field Camera 3 (WFC3) instrument. We use the G141 grism (1.1–1.7 μm) to study TrES-2b, TrES-4b, and CoRoT-1b in transit; TrES-3b in secondary eclipse; and WASP-4b in both. This wavelength region includes a predicted absorption feature from water at 1.4 μm, which we expect to be nondegenerate with the other molecules that are likely to be abundant for hydrocarbon-poor (e.g., solar composition) hot Jupiter atmospheres. We divide our wavelength regions into 10 bins. For each bin we produce a spectrophotometric light curve spanning the time of transit or eclipse. We correct these light curves for instrumental systematics without reference to an instrument model. For our transmission spectra, our mean 1σ precision per bin corresponds to variations of 2.1, 2.8, and 3.0 atmospheric scale heights for TrES-2b, TrES-4b, and CoRoT-1b, respectively. We find featureless spectra for these three planets. We are unable to extract a robust transmission spectrum for WASP-4b. For our dayside emission spectra, our mean 1σ precision per bin corresponds to a planet-to-star flux ratio of 1.5 × 10⁻⁴ and 2.1 × 10⁻⁴ for WASP-4b and TrES-3b, respectively. We combine these estimates with previous broadband measurements and conclude that for both planets isothermal atmospheres are disfavored. We find no signs of features due to water. We confirm that WFC3 is suitable for studies of transiting exoplanets, but in staring mode multivisit campaigns are necessary to place strong constraints on water abundance.

We present the average rest-frame spectrum of high-redshift dusty, star-forming galaxies from 250 to 770 GHz. This spectrum was constructed by stacking Atacama Large Millimeter/submillimeter Array (ALMA) 3 mm spectra of 22 such sources discovered by the South Pole Telescope and spanning z = 2.0–5.7. In addition to multiple bright spectral features of ¹²CO, [C i], and H₂O, we also detect several faint transitions of ¹³CO, HCN, HNC, HCO⁺, and CN, and use the observed line strengths to characterize the typical properties of the interstellar medium of these high-redshift starburst galaxies. We find that the ¹³CO brightness in these objects is comparable to that of the only other z > 2 star-forming galaxy in which ¹³CO has been observed. We show that the emission from the high-critical density molecules HCN, HNC, HCO⁺, and CN is consistent with a warm, dense medium with T_kin ∼ 55 K and $n_\mathrm{H_2} \gtrsim 10^{5.5}$ cm⁻³. High molecular hydrogen densities are required to reproduce the observed line ratios, and we demonstrate that alternatives to purely collisional excitation are unlikely to be significant for the bulk of these systems. We quantify the average emission from several species with no individually detected transitions, and find emission from the hydride CH and the linear molecule CCH for the first time at high redshift, indicating that these molecules may be powerful probes of interstellar chemistry in high-redshift systems. These observations represent the first constraints on many molecular species with rest-frame transitions from 0.4 to 1.2 mm in star-forming systems at high redshift, and will be invaluable in making effective use of ALMA in full science operations.

Observations of the afterglows of long gamma-ray bursts (GRBs) allow the study of star-forming galaxies across most of cosmic history. Here we present observations of GRB 111008A, from which we can measure metallicity, chemical abundance patterns, dust-to-metals ratio (DTM), and extinction of the GRB host galaxy at z = 5.0. The host absorption system is a damped Lyα absorber with a very large neutral hydrogen column density of $\log N({{\rm H}\,{\scriptsize {I}}})/{\rm cm}^{-2}=22.30\pm 0.06$ and a metallicity of [S/H] = −1.70 ± 0.10. It is the highest-redshift GRB with such a precise metallicity measurement. The presence of fine-structure lines confirms the z = 5.0 system as the GRB host galaxy and makes this the highest redshift where Fe ii fine-structure lines have been detected. The afterglow is mildly reddened with A_V = 0.11 ± 0.04 mag, and the host galaxy has a DTM that is consistent with being equal to or lower than typical values in the Local Group.

We present a quantitative analysis of the low-resolution (∼4.5 Å) spectra of 12 late-B and early-A blue supergiants (BSGs) in the metal-poor dwarf galaxy NGC 3109. A modified method of analysis is presented which does not require use of the Balmer jump as an independent T_eff indicator, as used in previous studies. We determine stellar effective temperatures, gravities, metallicities, reddening, and luminosities, and combine our sample with the early-B-type BSGs analyzed by Evans et al. to derive the distance to NGC 3109 using the flux-weighted gravity–luminosity relation (FGLR). Using primarily Fe-group elements, we find an average metallicity of $[\bar{Z}]$ = −0.67 ± 0.13, and no evidence of a metallicity gradient in the galaxy. Our metallicities are higher than those found by Evans et al. based on the oxygen abundances of early-B supergiants ($[\bar{Z}]$ = −0.93 ± 0.07), suggesting a low α/Fe ratio for the galaxy. We adjust the position of NGC 3109 on the BSG-determined galaxy mass–metallicity relation accordingly and compare it to metallicity studies of H ii regions in star-forming galaxies. We derive an FGLR distance modulus of 25.55 ± 0.09 (1.27 Mpc) that compares well with Cepheid and tip of the red giant branch distances. The FGLR itself is consistent with those found in other galaxies, demonstrating the reliability of this method as a measure of extragalactic distances.

We describe a new paradigm for understanding both relativistic motions and particle acceleration in the M87 jet: a magnetically dominated relativistic flow that naturally produces four relativistic magnetohydrodynamic (MHD) shocks (forward/reverse fast and slow modes). We apply this model to a set of optical super- and subluminal motions discovered by Biretta and coworkers with the Hubble Space Telescope during 1994–1998. The model concept consists of ejection of a single relativistic Poynting jet, which possesses a coherent helical (poloidal + toroidal) magnetic component, at the remarkably flaring point HST-1. We are able to reproduce quantitatively proper motions of components seen in the optical observations of HST-1 with the same model we used previously to describe similar features in radio very long baseline interferometry observations in 2005–2006. This indicates that the quad relativistic MHD shock model can be applied generally to recurring pairs of super/subluminal knots ejected from the upstream edge of the HST-1 complex as observed from radio to optical wavelengths, with forward/reverse fast-mode MHD shocks then responsible for observed moving features. Moreover, we identify such intrinsic properties as the shock compression ratio, degree of magnetization, and magnetic obliquity and show that they are suitable to mediate diffusive shock acceleration of relativistic particles via the first-order Fermi process. We suggest that relativistic MHD shocks in Poynting-flux-dominated helical jets may play a role in explaining observed emission and proper motions in many active galactic nuclei.

We present results from near-infrared spectroscopy of 26 emission-line galaxies at z ∼ 2.2 and z ∼ 1.5 obtained with the Folded-port InfraRed Echellette (FIRE) spectrometer on the 6.5 m Magellan Baade telescope. The sample was selected from the WFC3 Infrared Spectroscopic Parallels survey, which uses the near-infrared grism of the Hubble Space Telescope Wide Field Camera 3 (WFC3) to detect emission-line galaxies over 0.3 ≲ z ≲ 2.3. Our FIRE follow-up spectroscopy (R ∼ 5000) over 1.0–2.5 μm permits detailed measurements of the physical properties of the z ∼ 2 emission-line galaxies. Dust-corrected star formation rates for the sample range from ∼5–100 M_☉ yr⁻¹ with a mean of 29 M_☉ yr⁻¹. We derive a median metallicity for the sample of 12 + log(O/H) = 8.34 or ∼0.45 Z_☉. The estimated stellar masses range from ∼10^8.5–10^9.5 M_☉, and a clear positive correlation between metallicity and stellar mass is observed. The average ionization parameter measured for the sample, log U ≈ −2.5, is significantly higher than what is found for most star-forming galaxies in the local universe, but similar to the values found for other star-forming galaxies at high redshift. We derive composite spectra from the FIRE sample, from which we measure typical nebular electron densities of ∼100–400 cm⁻³. Based on the location of the galaxies and composite spectra on diagnostic diagrams, we do not find evidence for significant active galactic nucleus activity in the sample. Most of the galaxies, as well as the composites, are offset diagram toward higher [O iii]/Hβ at a given [N ii]/Hα, in agreement with other observations of z ≳ 1 star-forming galaxies, but composite spectra derived from the sample do not show an appreciable offset from the local star-forming sequence on the [O iii]/Hβ versus [S ii]/Hα diagram. We infer a high nitrogen-to-oxygen abundance ratio from the composite spectrum, which may contribute to the offset of the high-redshift galaxies from the local star-forming sequence in the [O iii]/Hβ versus [N ii]/Hα diagram. We speculate that the elevated nitrogen abundance could result from substantial numbers of Wolf–Rayet stars in starbursting galaxies at z ∼ 2.

We present spectral energy distributions (SEDs) of 69 QSOs at z > 5, covering a rest frame wavelength range of 0.1 μm to ∼80 μm, and centered on new Spitzer and Herschel observations. The detection rate of the QSOs with Spitzer is very high (97% at λ_rest ≲ 4 μm), but drops toward the Herschel bands with 30% detected in PACS (rest frame mid-infrared) and 15% additionally in the SPIRE (rest frame far-infrared; FIR). We perform multi-component SED fits for Herschel-detected objects and confirm that to match the observed SEDs, a clumpy torus model needs to be complemented by a hot (∼1300 K) component and, in cases with prominent FIR emission, also by a cold (∼50 K) component. In the FIR-detected cases the luminosity of the cold component is of the order of 10¹³ L_☉ which is likely heated by star formation. From the SED fits we also determine that the active galactic nucleus (AGN) dust-to-accretion disk luminosity ratio declines with UV/optical luminosity. Emission from hot (∼1300 K) dust is common in our sample, showing that nuclear dust is ubiquitous in luminous QSOs out to redshift 6. However, about 15% of the objects appear under-luminous in the near infrared compared to their optical emission and seem to be deficient in (but not devoid of) hot dust. Within our full sample, the QSOs detected with Herschel are found at the high luminosity end in L_UV/opt and L_NIR and show low equivalent widths (EWs) in Hα and in Lyα. In the distribution of Hα EWs, as determined from the Spitzer photometry, the high-redshift QSOs show little difference to low-redshift AGN.

We present the first microlensing candidate for a free-floating exoplanet–exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ∼ 4 Jupiter masses hosting a sub-Earth mass moon. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M_Lπ_rel where M_L is the lens system mass and π_rel is the lens–source relative parallax. If the lens system is nearby (large π_rel), then M_L is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens–source relative proper motion, μ_rel = 19.6 ± 1.6 mas yr⁻¹, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, there is an alternate solution that has a lower proper motion and fits the data almost as well. This solution is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of $M_{\rm host} = 0.12^{+0.19}_{ -0.06}\,M_\odot$ and $m_{\rm comp} = 18^{+28}_{ -10}\,{M_\oplus }$, at a projected separation of $a_\perp = 0.84^{+0.25}_{ -0.14}$ AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

We develop a new method to measure source proper motions in microlensing events, which can partially overcome problems due to blending. It takes advantage of the fact that the source position is known precisely from the microlensing event itself. We apply this method to the event MOA-2011-BLG-262, which has a short timescale t_E = 3.8 day, a companion mass ratio q = 4.7 × 10⁻³, and a very high or high lens–source relative proper motion μ_rel = 20 mas yr⁻¹ or 12 mas yr⁻¹ (for two possible models). These three characteristics imply that the lens could be a brown dwarf or a massive planet with a roughly Earth-mass "moon." The probability of such an interpretation would be greatly increased if it could be shown that the high lens–source relative proper motion was primarily due to the lens rather than the source. Based on the long-term monitoring data of the Galactic bulge from the Optical Gravitational Lensing Experiment, we measure the source proper motion that is small, ${\boldsymbol {\mu }}_s = (-2.3, -0.9)\pm (2.8,2.6)\,{\rm mas\:yr^{-1}}$ in a (north, east) Galactic coordinate frame. These values are then important input into a Bayesian analysis of the event presented in a companion paper by Bennett et al.

We examine the exchange of angular momentum between the component spins and the orbit in semi-detached double white dwarf binaries undergoing mass transfer through direct impact of the transfer stream. We approximate the stream as a series of discrete massive particles ejected in the ballistic limit at the inner Lagrangian point of the donor toward the accretor. This work improves upon similar earlier studies in a number of ways. First, we self-consistently calculate the total angular momentum of the orbit at all times. This includes changes in the orbital angular momentum during the ballistic trajectory of the ejected mass, as well as changes during the ejection/accretion due to the radial component of the particle's velocity. Second, we calculate the particle's ballistic trajectory for each system, which allows us to determine the precise position and velocity of the particle upon accretion. We can then include specific information about the radius of the accretor as well as the angle of impact. Finally, we ensure that the total angular momentum is conserved, which requires the donor star spin to vary self-consistently. With these improvements, we calculate the angular momentum change of the orbit and each binary component across the entire parameter space of direct impact double white dwarf binary systems. We find a significant decrease in the amount of angular momentum removed from the orbit during mass transfer, as well as cases where this process increases the angular momentum of the orbit at the expense of the spin angular momentum of the donor. We conclude that, unlike earlier claims in the literature, mass transfer through direct impact need not destabilize the binary and that the quantity and sign of the orbital angular momentum transfer depends on the binary properties, particularly the masses of the double white dwarf binary component stars. This stabilization may significantly impact the population synthesis calculations of the expected numbers of events/systems for which double white dwarfs may be a progenitor, e.g., Type Ia supernovae, Type.Ia supernovae, and AM CVn.

A number of brown dwarfs are now known to be variable with observed amplitudes as large as 10%–30% at some wavelengths. While spatial inhomogeneities in cloud coverage and thickness are likely responsible for much of the observed variability, it is possible that some of the variations arise from atmospheric temperature fluctuations instead of, or in addition to, clouds. To better understand the role that thermal variability might play we present a case study of brown dwarf variability using a newly developed one-dimensional, time-stepping model of atmospheric thermal structure. We focus on the effects of thermal perturbations, intentionally simplifying the problem through omission of clouds and atmospheric circulation. Model results demonstrate that thermal perturbations occurring deep in the atmosphere (at pressures greater than 10 bar) of a model T-dwarf can be communicated to the upper atmosphere through radiative heating via the windows in near-infrared water opacity. The response time depends on where in the atmosphere a thermal perturbation is introduced. We show that, for certain periodic perturbations, the emission spectrum can have complex time- and wavelength-dependent behaviors, including phase shifts in times of maximum flux observed at different wavelengths. Since different wavelengths probe different levels in the atmosphere, these variations track a wavelength-dependent set of radiative exchanges happening between different atmospheric levels as a perturbation evolves in time. We conclude that thermal—as well as cloud—fluctuations must be considered as possible contributors to the observed brown dwarf variability.

SONYC—Substellar Objects in Nearby Young Clusters—is a survey program to investigate the frequency and properties of substellar objects in nearby star-forming regions. We present a new imaging and spectroscopic survey conducted in the young (∼1 Myr), nearby (∼200 pc) star-forming region Lupus 3. Deep optical and near-infrared images were obtained with MOSAIC-II and NEWFIRM at the CTIO 4 m telescope, covering ∼1.4 deg² on the sky. The i-band completeness limit of 20.3 mag is equivalent to 0.009–0.02 M_☉, for A_V ⩽ 5. Photometry and 11–12 yr baseline proper motions were used to select candidate low-mass members of Lupus 3. We performed a spectroscopic follow-up of 123 candidates, using VIMOS at the Very Large Telescope, and we identify 7 probable members, among which 4 have spectral type later than M6.0 and T_eff ⩽ 3000 K, i.e., are probably substellar in nature. Two of the new probable members of Lupus 3 appear underluminous for their spectral class and exhibit emission line spectrum with strong H_α or forbidden lines associated with active accretion. We derive a relation between the spectral type and effective temperature: T_eff = (4120 ± 175) − (172 ± 26) × SpT, where SpT refers to the M spectral subtype between 1 and 9. Combining our results with the previous works on Lupus 3, we show that the spectral type distribution is consistent with that in other star-forming regions, as well as the derived star-to-brown dwarf ratio of 2.0–3.3. We compile a census of all spectroscopically confirmed low-mass members with spectral type M0 or later.

In this work, we apply an updated version of the Neutron Monitor (NM) Based Anisotropic GLE Pure Power Law (NMBANGLE PPOLA) model, in order to derive the characteristics of the ground-level enhancement (GLE) on 2012 May 17 (GLE71), the spectral properties of the related solar energetic particle (SEP) event, the spatial distributions of the high-energy solar cosmic ray fluxes at the top of the atmosphere, and the time evolution of the location of the GLE source. Our modeling, based uniquely on the use of ground-level NM data, leads to the following main results. The SEP spectrum related to GLE71 was rather soft during the whole duration of the event, manifesting some weak acceleration episodes only during the initial phase (at ∼01:55–02:00 UT) and at ∼02:30–02:35 UT and ∼02:55–03:00 UT. The spectral index of the modeled SEP spectrum supports the coronal mass ejection–shock driven particle acceleration scenario, in agreement with past results based on the analysis of satellite measurements. During the initial phase of GLE71, the solar proton source at the top of the atmosphere was located above the northern hemisphere, implying that the asymptotic directions of viewing of the northern hemisphere NMs were more favorably located for registering the event than the southern ones. The spatial distribution of the solar proton fluxes at the top of the atmosphere during the main phase manifested a large variation along longitude and latitude. At the rigidity of 1 GV, the maximum primary solar proton flux resulted on the order of ∼3 × 10⁴ part. m⁻² s⁻¹ sr⁻¹ GV⁻¹.

We use a newly developed observing mode on the Hubble Space Telescope (HST) and Wide Field Camera 3, spatial scanning, to increase source sampling a thousand-fold and measure changes in source positions to a precision of 20–40 μas, more than an order of magnitude better than attainable in pointed observations. This observing mode can usefully measure the parallaxes of bright stars at distances of up to 5 kpc, a factor of 10 farther than achieved thus far with HST. The technique should also provide a unique crosscheck of future parallax measurements from Gaia. Long-period classical Cepheid variable stars in the Milky Way, nearly all of which reside beyond 1 kpc, are especially compelling targets for parallax measurements from scanning, as they may be used to anchor a determination of the Hubble constant to ∼1%. We illustrate the method by measuring to high precision the parallax of a classical Cepheid, SY Aurigae, at a distance of more than 2 kpc, using five epochs of spatial-scan data obtained at intervals of six months. Rapid spatial scans also enable photometric measurements of bright Milky Way Cepheids—which would otherwise saturate even in the shortest possible pointed observations—on the same flux scale as extragalactic Cepheids, which is a necessity for reducing a leading source of systematic error in the Hubble constant. We demonstrate this capability with photometric measurements of SY Aur on the same system used for Cepheids in Type Ia supernova host galaxies. While the technique and results presented here are preliminary, an ongoing program with HST is collecting such parallax measurements for another 18 Cepheids to produce a better anchor for the distance scale.

The Magellanic Bridge is the nearest low-metallicity, tidally stripped environment, offering a unique high-resolution view of physical conditions in merging and forming galaxies. In this paper, we present an analysis of candidate massive young stellar objects (YSOs), i.e., in situ, current massive star formation (MSF) in the Bridge using Spitzer mid-IR and complementary optical and near-IR photometry. While we definitely find YSOs in the Bridge, the most massive are ∼10 M_☉, ≪45 M_☉ found in the LMC. The intensity of MSF in the Bridge also appears to be decreasing, as the most massive YSOs are less massive than those formed in the past. To investigate environmental effects on MSF, we have compared properties of massive YSOs in the Bridge to those in the LMC. First, YSOs in the Bridge are apparently less embedded than in the LMC: 81% of Bridge YSOs show optical counterparts, compared to only 56% of LMC sources with the same range of mass, circumstellar dust mass, and line-of-sight extinction. Circumstellar envelopes are evidently more porous or clumpy in the Bridge's low-metallicity environment. Second, we have used whole samples of YSOs in the LMC and the Bridge to estimate the probability of finding YSOs at a given H i column density, N(H i). We found that the LMC has ∼3 × higher probability than the Bridge for N(H i) >12 × 10²⁰ cm⁻², but the trend reverses at lower N(H i). Investigating whether this lower efficiency relative to H i is due to less efficient molecular cloud formation or to less efficient cloud collapse, or to both, will require sensitive molecular gas observations.

We consider the evolution of supermassive black hole binaries at the center of spherical, axisymmetric, and triaxial galaxies, using direct N-body integrations as well as analytic estimates. We find that the rates of binary hardening exhibit a significant N-dependence in all the models, at least for N in the investigated range of 10⁵ ⩽ N ⩽ 10⁶. Binary hardening rates are also substantially lower than would be expected if the binary "loss cone" remained "full," as it would be if the orbits supplying stars to the binary were being efficiently replenished. The difference in binary hardening rates between the spherical and nonspherical models is less than a factor of two even in the simulations with the largest N. By studying the orbital populations of our models, we conclude that the rate of supply of stars to the binary via draining of centrophilic orbits is indeed expected to be much lower than the full-loss-cone rate, consistent with our simulations. We argue that the binary's evolution in the simulations is driven in roughly equal amounts by collisional and collisionless effects, even at the highest N-values currently accessible. While binary hardening rates would probably reach a limiting value for large N, our results suggest that we cannot approach that rate with currently available algorithms and computing hardware. The extrapolation of results from N-body simulations to real galaxies is therefore not straightforward, casting doubt on recent claims that triaxiality or axisymmetry alone are capable of solving the final-parsec problem in gas-free galaxies.

Massive young stellar clusters are strong sources of radiation and mechanical energy. Their powerful winds and radiation pressure sweep up interstellar gas into thin expanding shells that trap the ionizing radiation produced by the central clusters affecting the dynamics and the distribution of their ionized gas. Here we continue our comparison of the star cluster winds and radiation pressure effects on the dynamics of shells around young massive clusters. We calculate the impact that radiation pressure has on the distribution of matter and thermal pressure within such shells, as well as on the density-weighted ionization parameter U_w, and put our results on the diagnostic diagram, which allows one to discriminate between the wind-dominated and radiation-dominated regimes. We found that model-predicted values of the ionization parameter agree well with typical values found in local starburst galaxies. Radiation pressure may affect the inner structure and the dynamics of wind-driven shells, but only during the earliest stages of evolution (before ∼3 Myr) or if a major fraction of the star cluster mechanical luminosity is dissipated or radiated away within the star cluster volume and thus the star cluster mechanical energy output is significantly smaller than star cluster synthetic models predict. However, even in these cases radiation dominates over the wind dynamical pressure only if the exciting cluster is embedded into a high-density ambient medium.