Table of contents

The amount of high-energy stellar radiation reaching the surface of protoplanetary disks is essential to determine their chemistry and physical evolution. Here, we use millimetric and centimetric radio data to constrain the extreme-ultraviolet (EUV) luminosity impinging on 14 disks around young (∼2–10 Myr) sun-like stars. For each object we identify the long-wavelength emission in excess to the dust thermal emission, attribute that to free–free disk emission, and thereby compute an upper limit to the EUV reaching the disk. We find upper limits lower than 10⁴² photons s⁻¹ for all sources without jets and lower than 5 × 10⁴⁰ photons s⁻¹ for the three older sources in our sample. These latter values are low for EUV-driven photoevaporation alone to clear out protoplanetary material in the timescale inferred by observations. In addition, our EUV upper limits are too low to reproduce the [Ne ii] 12.81 μm luminosities from three disks with slow [Ne ii]-detected winds. This indicates that the [Ne ii] line in these sources primarily traces a mostly neutral wind where Ne is ionized by 1 keV X-ray photons, implying higher photoevaporative mass loss rates than those predicted by EUV-driven models alone. In summary, our results suggest that high-energy stellar photons other than EUV may dominate the dispersal of protoplanetary disks around sun-like stars.

We present time series analyses of the full Kepler data set of Zw 229−15. This Kepler light curve—with a baseline greater than 3 yr, composed of virtually continuous, evenly sampled 30 minute measurements—is unprecedented in its quality and precision. We utilize two methods of power spectral analysis to investigate the optical variability and search for evidence of a bend frequency associated with a characteristic optical variability timescale. Each method yields similar results. The first interpolates across data gaps to use the standard Fourier periodogram. The second, using the CARMA-based time-domain modeling technique of Kelly et al., does not need evenly sampled data. Both methods find excess power at high frequencies that may be due to Kepler instrumental effects. More importantly, both also show strong bends (Δα ∼ 2) at timescales of ∼5 days, a feature similar to those seen in the X-ray power spectral densities of active galactic nuclei (AGNs) but never before in the optical. This observed ∼5 day timescale may be associated with one of several physical processes potentially responsible for the variability. A plausible association could be made with light-crossing dynamical or thermal timescales depending on the assumed value of the accretion disk size and on unobserved disk parameters such as α and H/R. This timescale is not consistent with the viscous timescale, which would be years in a ∼10⁷ M_☉ AGN such as Zw 229−15. However, there must be a second bend on long (≳ 1 yr) timescales and that feature could be associated with the viscous timescale.

Maximum likelihood estimation of the bulk flow from radial peculiar motions of galaxies generally assumes a constant velocity field inside the survey volume. This assumption is inconsistent with the definition of bulk flow as the average of the peculiar velocity field over the relevant volume. This follows from a straightforward mathematical relation between the bulk flow of a sphere and the velocity potential on its surface. This inconsistency also exists for ideal data with exact radial velocities and full spatial coverage. Based on the same relation, we propose a simple modification to correct for this inconsistency.

We present observations of compound flux rope formation, which occurred on 2014 January 1, via merging of two nearby filament channels, the associated dynamics, and its stability using multiwavelength data. We also discuss the dynamics of cool and hot plasma moving along the newly formed compound flux rope. The merging started after the interaction between the southern leg of the northward filament and the northern leg of the southward filament at ≈01:21 UT and continued until a compound flux rope formed at ≈01:33 UT. During the merging, the cool filament plasma heated up and started to move along both sides of the compound flux rope, i.e., toward the north (≈265 km s⁻¹) and south (≈118 km s⁻¹) from the point of merging. After traveling a distance of ≈150 Mm toward the north, the plasma cooled down and started to return back to the south (≈14 km s⁻¹) after ≈02:00 UT. The observations provide a clear example of compound flux rope formation via merging of two different flux ropes and the occurrence of a flare through tether cutting reconnection. However, the compound flux rope remained stable in the corona and had a confined eruption. The coronal magnetic field decay index measurements revealed that both the filaments and the compound flux rope axis lie within the stability domain (decay index <1.5), which may be the possible cause for their stability. The present study also deals with the relationship between the filament's chirality (sinistral) and the helicity (positive) of the surrounding flux rope.

We present the statistics of faint submillimeter/millimeter galaxies (SMGs) and serendipitous detections of a submillimeter/millimeter line emitter (SLE) with no multi-wavelength continuum counterpart revealed by the deep ALMA observations. We identify faint SMGs with flux densities of 0.1–1.0 mJy in the deep Band-6 and Band-7 maps of 10 independent fields that reduce cosmic variance effects. The differential number counts at 1.2 mm are found to increase with decreasing flux density down to 0.1 mJy. Our number counts indicate that the faint (0.1–1.0 mJy, or SFR_IR ∼ 30–300 M_☉ yr⁻¹) SMGs contribute nearly a half of the extragalactic background light (EBL), while the remaining half of the EBL is mostly contributed by very faint sources with flux densities of <0.1 mJy (SFR_IR ≲ 30 M_☉ yr⁻¹). We conduct counts-in-cells analysis with multifield ALMA data for the faint SMGs, and obtain a coarse estimate of galaxy bias, b_g < 4. The galaxy bias suggests that the dark halo masses of the faint SMGs are ≲ 7 × 10¹² M_☉, which is smaller than those of bright (>1 mJy) SMGs, but consistent with abundant high-z star-forming populations, such as sBzKs, LBGs, and LAEs. Finally, we report the serendipitous detection of SLE-1, which has no continuum counterparts in our 1.2 mm-band or multi-wavelength images, including ultra deep HST/WFC3 and Spitzer data. The SLE has a significant line at 249.9 GHz with a signal-to-noise ratio of 7.1. If the SLE is not a spurious source made by the unknown systematic noise of ALMA, the strong upper limits of our multi-wavelength data suggest that the SLE would be a faint galaxy at z ≳ 6.

A formula for the tidal dissipation rate in a spherical body is derived from first principles to correct some mathematical inaccuracies found in the literature. The development is combined with the Darwin–Kaula formalism for tides. Our intermediate results are compared with those by Zschau and Platzman. When restricted to the special case of an incompressible spherical planet spinning synchronously without libration, our final formula can be compared with the commonly used expression from Peale & Cassen. However, the two turn out to differ, as in our expression the contributions from all Fourier modes are positive-definite, which is not the case with the formula from Peale & Cassen. Examples of the application of our expression for the tidal damping rate are provided in the work by Makarov & Efroimsky (Paper II) published back to back with the current paper.

In Efroimsky & Makarov (Paper I), we derived from the first principles a formula for the tidal heating rate in a homogeneous sphere, compared it with the previously used formulae, and noted the differences. Now we present case studies: Mercury, Kepler-10 b, and a triaxial Io. A sharp frequency dependence of k₂/Q near spin–orbit resonances yields a sharp dependence of k₂/Q (and, therefore, of tidal heating) upon the spin rate. Thereby physical libration plays a major role in tidal heating of synchronously rotating planets. The magnitude of libration in the spin rate being defined by the planet's triaxiality, the latter becomes a factor determining the dissipation rate. Other parameters equal, a strongly triaxial synchronized body generates more heat than a similar body of a more symmetrical shape. After an initially triaxial object melts and loses its triaxiality, dissipation becomes less intensive; the body can solidify, with the tidal bulge becoming a new figure with triaxiality lower than the original. We derive approximate expressions for the dissipation rate in a Maxwell planet with the Maxwell time longer than the inverse tidal frequency. The expressions derived pertain to the 1:1 and 3:2 resonances and a nonresonant case; so they are applicable to most close-in super-Earths detected. In these planets, the heating outside synchronism is weakly dependent on the eccentricity and obliquity, provided both these parameters's values are moderate. According to our calculation, Kepler-10 b could hardly survive the intensive tidal heating without being synchronized, circularized, and reshaped through a complete or partial melt-down.

Pre-outburst and post-outburst light curves and post-outburst eclipse timings are analyzed to measure any period (P) change related to nova CI Aql's outburst of early 2000 and a mean post-outburst dP/dt, which then lead to estimates of the accreting component's rate of mass (M) change and its overall outburst-related change of mass over roughly a decade of observations. We apply a recently developed procedure for unified analysis of three timing-related data types (light curves, radial velocities, and eclipse timings), although with only light curves and timings in this case. Fits to the data are reasonably good without need for a disk in the light-curve model, although the disk certainly exists and has an important role in our post-outburst mass flow computations. Initial experiments showed that, although there seems to be an accretion hot spot, it has essentially no effect on derived outburst-related ΔP or on post-outburst dP/dt. Use of atomic time (HJED) in place of HJD also has essentially nil effect on ΔP and dP/dt. We find ΔP consistently negative in various types of solutions, although at best only marginally significant statistically in any one experiment. Pre-outburst HJD₀ and P results are given, as are post-outburst HJD₀, P, and dP/dt, with light curves and eclipse times as joint input, and also with only eclipse time input. Post-outburst dP/dt is negative at about 2.4σ. Explicit formulae for mass transfer rates and epoch-to-epoch mass change are developed and applied. A known offset in the magnitude zero point for 1991–1994 is corrected.

Spectropolarimetric temporal series from Fe i λ6301.5 Å and Ca ii infrared triplet lines are obtained by applying the Stokes synthesis code NICOLE to a numerical simulation of wave propagation in a sunspot umbra from MANCHA code. The analysis of the phase difference between Doppler velocity and intensity core oscillations of the Fe i λ6301.5 Å line reveals that variations in the intensity are produced by opacity fluctuations rather than intrinsic temperature oscillations, except for frequencies between 5 and 6.5 mHz. On the other hand, the photospheric magnetic field retrieved from the weak field approximation provides the intrinsic magnetic field oscillations associated to wave propagation. Our results suggest that this is due to the low magnetic field gradient of our sunspot model. The Stokes parameters of the chromospheric Ca ii infrared triplet lines show striking variations as shock waves travel through the formation height of the lines, including emission self-reversals in the line core and highly abnormal Stokes V profiles. Magnetic field oscillations inferred from the Ca ii infrared lines using the weak field approximation appear to be related with the magnetic field strength variation between the photosphere and the chromosphere.

Magnetic energy released in the corona by solar flares reaches the chromosphere where it drives characteristic upflows and downflows known as evaporation and condensation. These flows are studied here for the case where energy is transported to the chromosphere by thermal conduction. An analytic model is used to develop relations by which the density and velocity of each flow can be predicted from coronal parameters including the flare's energy flux F. These relations are explored and refined using a series of numerical investigations in which the transition region (TR) is represented by a simplified density jump. The maximum evaporation velocity, for example, is well approximated by v_e ≃ 0.38(F/ρ_{co, 0})^1/3, where ρ_{co, 0} is the mass density of the pre-flare corona. This and the other relations are found to fit simulations using more realistic models of the TR both performed in this work, and taken from a variety of previously published investigations. These relations offer a novel and efficient means of simulating coronal reconnection without neglecting entirely the effects of evaporation.

The Galactic cosmic ray spectrum exhibits subtle variations over the 22 yr solar magnetic cycle in addition to the more dramatic variations over the 11 yr sunspot cycle. Neutron monitors are large ground-based detectors that provide accurate measurements of variations in the cosmic ray flux at the top of the atmosphere above the detector. At any given location the magnetic field of the Earth excludes particles below a well-defined rigidity (momentum per unit charge) known as the cutoff rigidity, which can be accurately calculated using detailed models of the geomagnetic field. By carrying a neutron monitor to different locations, e.g., on a ship, the Earth itself serves as a magnet spectrometer. By repeating such latitude surveys with identical equipment, a sensitive measurement of changes in the spectrum can be made. In this work, we analyze data from the 1994 through 2007 series of latitude surveys conducted by the Bartol Research Institute, the University of Tasmania, and the Australian Antarctic Division. We confirm the curious "crossover" in spectra measured near solar minima during epochs of opposite solar magnetic polarity, and show that it is directly related to a sudden change in the spectral behavior of solar modulation at the time of the polarity reversal, as revealed from contemporaneous variations in the survey data and a fixed station. We suggest that the spectral change and crossover result from the interaction of effects due to gradient/curvature drifts with a systematic change in the interplanetary diffusion coefficient caused by turbulent magnetic helicity.

In contrast to photometric transits, whose peak signal occurs at mid-transit due to occultation of the brightest region of the disk, polarimetric transits provide a signal upon ingress and egress due to occultation of the polarized stellar limb. Limb polarization, the bright corollary to limb darkening, arises from the 90° scattering angle and low optical depth experienced by photons at the limb. In addition to the ratio R_p/R_*, the amplitude of a polarimetric transit is expected to be controlled by the strength and width of the stellar limb polarization profile, which depend on the scattering-to-total opacity ratio at the stellar limb. We present a short list of the systems providing the highest expected signal-to-noise ratio for detection of this effect, and we draw particular attention to HD 80606b. This planet is spin/orbit misaligned, has a three-hour ingress, and has a bright parent star, which make it an attractive target. We report on test observations of an HD 80606b ingress with the POLISH2 polarimeter at the Lick Observatory Shane 3 m telescope. We conclude that unmodeled telescope systematic effects prevented polarimetric detection of this event. We outline a roadmap for further refinements of exoplanet polarimetry, whose eventual success will require a further factor of ten reduction in systematic noise.

Using the IRAM 30 m telescope, we have surveyed a 1 × 08 part of the Orion molecular cloud in the ¹²CO and ¹³CO (2–1) lines with a maximal spatial resolution of ∼11'' and spectral resolution of ∼0.4 km s⁻¹. The cloud appears filamentary, clumpy, and with a complex kinematical structure. We derive an estimated mass of the cloud of 7700 M_☉ (half of which is found in regions with visual extinctions A_V below ∼10) and a dynamical age for the nebula of the order of 0.2 Myr. The energy balance suggests that magnetic fields play an important role in supporting the cloud, at large and small scales. According to our analysis, the turbulent kinetic energy in the molecular gas due to outflows is comparable to turbulent kinetic energy resulting from the interaction of the cloud with the H ii region. This latter feedback appears negative, i.e., the triggering of star formation by the H ii region is inefficient in Orion. The reduced data as well as additional products such as the column density map are made available online (http://userpages.irap.omp.eu/~oberne/Olivier_Berne/Data).

We report ground-based, low-frequency (<100 MHz) radio imaging, spectral, and polarimeter observations of the type II radio burst associated with the solar coronal mass ejection (CME) that occurred on 2013 May 2. The spectral observations indicate that the burst has fundamental (F) and harmonic (H) emission components with split-band and herringbone structures. The imaging observations at 80 MHz indicate that the H component of the burst was located close to leading edge of the CME at a radial distance of r ≈ 2 R_☉ in the solar atmosphere. The polarimeter observations of the type II burst, also at 80 MHz, indicate that the peak degree of circular polarization (dcp) corresponding to the emission generated in the corona ahead of and behind the associated MHD shock front are ≈0.05 ± 0.02 and ≈0.1 ± 0.01, respectively. We calculated the magnetic field B in the above two coronal regions by adopting the empirical relationship between the dcp and B for the harmonic plasma emission and the values are ≈(0.7–1.4) ± 0.2 G and ≈(1.4–2.8) ± 0.1 G, respectively.

Using data from the STEREO and SOHO spacecraft, we show that temporal organization of energy release events in the quiet solar corona is close to random, in contrast to the clustered behavior of flaring times in solar active regions. The locations of the quiet-Sun events follow the meso- and supergranulation pattern of the underling photosphere. Together with earlier reports of the scale-free event size statistics, our findings suggest that quiet solar regions responsible for bulk coronal heating operate in a driven self-organized critical state, possibly involving long-range Alfvénic interactions.

Hydromagnetic turbulence affects the evolution of large-scale magnetic fields through mean-field effects like turbulent diffusion and the α effect. For stronger fields, these effects are usually suppressed or quenched, and additional anisotropies are introduced. Using different variants of the test-field method, we determine the quenching of the turbulent transport coefficients for the forced Roberts flow, isotropically forced non-helical turbulence, and rotating thermal convection. We see significant quenching only when the mean magnetic field is larger than the equipartition value of the turbulence. Expressing the magnetic field in terms of the equipartition value of the quenched flows, we obtain for the quenching exponents of the turbulent magnetic diffusivity about 1.3, 1.1, and 1.3 for Roberts flow, forced turbulence, and convection, respectively. However, when the magnetic field is expressed in terms of the equipartition value of the unquenched flows, these quenching exponents become about 4, 1.5, and 2.3, respectively. For the α effect, the exponent is about 1.3 for the Roberts flow and 2 for convection in the first case, but 4 and 3, respectively, in the second. In convection, the quenching of turbulent pumping follows the same power law as turbulent diffusion, while for the coefficient describing the ${\boldsymbol{\Omega}} \times \boldsymbol {{J}}$ effect nearly the same quenching exponent is obtained as for α. For forced turbulence, turbulent diffusion proportional to the second derivative along the mean magnetic field is quenched much less, especially for larger values of the magnetic Reynolds number. However, we find that in corresponding axisymmetric mean-field dynamos with dominant toroidal field the quenched diffusion coefficients are the same for the poloidal and toroidal field constituents.

Photospheric electric fields, estimated from sequences of vector magnetic field and Doppler measurements, can be used to estimate the flux of magnetic energy (the Poynting flux) into the corona and as time-dependent boundary conditions for dynamic models of the coronal magnetic field. We have modified and extended an existing method to estimate photospheric electric fields that combines a poloidal–toroidal decomposition (PTD) of the evolving magnetic field vector with Doppler and horizontal plasma velocities. Our current, more comprehensive method, which we dub the "PTD-Doppler-FLCT Ideal" (PDFI) technique, can now incorporate Doppler velocities from non-normal viewing angles. It uses the FISHPACK software package to solve several two-dimensional Poisson equations, a faster and more robust approach than our previous implementations. Here, we describe systematic, quantitative tests of the accuracy and robustness of the PDFI technique using synthetic data from anelastic MHD (ANMHD) simulations, which have been used in similar tests in the past. We find that the PDFI method has less than 1% error in the total Poynting flux and a 10% error in the helicity flux rate at a normal viewing angle (θ = 0) and less than 25% and 10% errors, respectively, at large viewing angles (θ < 60°). We compare our results with other inversion methods at zero viewing angle and find that our method's estimates of the fluxes of magnetic energy and helicity are comparable to or more accurate than other methods. We also discuss the limitations of the PDFI method and its uncertainties.

In the last few years, it has been found that transverse waves are present at all times in coronal loops or spicules. Their energy has been estimated with an expression derived for bulk Alfvén waves in homogeneous media, with correspondingly uniform wave energy density and flux. The kink mode, however, is localized in space with the energy density and flux dependent on the position in the cross-sectional plane. The more relevant quantities for the kink mode are the integrals of the energy density and flux over the cross-sectional plane. The present paper provides an approximation to the energy propagated by kink modes in an ensemble of flux tubes by means of combining the analysis of single flux tube kink oscillations with a filling factor for the tube cross-sectional area. This finally allows one to compare the expressions for energy flux of Alfvén waves with an ensemble of kink waves. We find that the correction factor for the energy in kink waves, compared to the bulk Alfvén waves, is between f and 2f, where f is the density filling factor of the ensemble of flux tubes.

We present the results of the longest yet undertaken search for apparently extragalactic radio bursts at the Bleien Radio Observatory covering 21,000 hr (898 days). The data were searched for events of less than 50 ms FWHM duration showing a ν⁻² drift in the spectrogram characteristic of the delay of radio waves in plasma. We have found five cases suggesting dispersion measures between 350 and 400 cm⁻³ pc while searching in the range of 75–2000 cm⁻³ pc. Four of the five events occurred between 10:27 and 11:24 a.m. local civil time. The only exception occurred at night with the full Moon in the beam. It was an event that poorly fits plasma dispersion, but had the characteristics of a solar Type III burst. However, we were not able to confirm that it was a lunar reflection. All events were observed with a log-periodic dipole within 6800 hr, but none with a more directional horn antenna observing the rest of the time. These properties suggest a terrestrial origin of the "peryton" type reported before. However, the cause of these events remains ambiguous.

The redshift-dependent fraction of color-selected galaxies revealing Lyman alpha (Lyα) emission, x_Lyα  has become the most valuable constraint on the evolving neutrality of the early intergalactic medium. However, in addition to resonant scattering by neutral gas, the visibility of Lyα is also dependent on the intrinsic properties of the host galaxy, including its stellar population, dust content, and the nature of outflowing gas. Taking advantage of significant progress we have made in determining the line-emitting properties of z ≃ 4–6 galaxies, we propose an improved method, based on using the measured slopes of the rest-frame ultraviolet continua of galaxies, to interpret the growing body of near-infrared spectra of z > 7 galaxies in order to take into account these host galaxy dependencies. In a first application of our new method, we demonstrate its potential via a new spectroscopic survey of 7 < z < 8 galaxies undertaken with the Keck MOSFIRE spectrograph. Together with earlier published data, our data provide improved estimates of the evolving visibility of Lyα, particularly at redshift z ≃ 8. As a by-product, we also present a promising new line-emitting galaxy candidate, detected at 4.0σ at redshift z = 7.62. We discuss the improving constraints on the evolving neutral fraction over 6 < z < 8 and the implications for cosmic reionization.

We report photometric results from a subset of a Herschel–PACS program to observe cool dust in disks around low-mass stars as a complement to our earlier program to measure far-infrared emission from brown dwarfs. In this latest study we observed five low-mass objects in the nearby R Corona Australis region and detected at least three at 70 μm. Using a Monte Carlo radiative transfer code we have investigated the disk masses and geometry based on detailed spectral energy distribution (SED) modeling, and we compare these new results to those from our earlier larger sample of brown dwarfs. In particular, our SED analysis for these five objects shows again that disk geometries of brown dwarfs or low-mass stars are generally similar to their higher mass counterparts like T Tauri disks, but the range of disk mass extends to well below the value found in T Tauri stars.

We present new determinations of the iron abundance for 220 stars belonging to the stellar system Terzan 5 in the Galactic bulge. The spectra have been acquired with FLAMES at the Very Large Telescope of the European Southern Observatory and DEIMOS at the Keck II Telescope. This is by far the largest spectroscopic sample of stars ever observed in this stellar system. From this data set, a subsample of targets with spectra unaffected by TiO bands was extracted and statistically decontaminated from field stars. Once combined with 34 additional stars previously published by our group, a total sample of 135 member stars covering the entire radial extent of the system has been used to determine the metallicity distribution function of Terzan 5. The iron distribution clearly shows three peaks: a super-solar component at [Fe/H] ≃ 0.25 dex, accounting for ∼29% of the sample, a dominant sub-solar population at [Fe/H] ≃ −0.30 dex, corresponding to ∼62% of the total, and a minor (6%) metal-poor component at [Fe/H] ≃ −0.8 dex. Such a broad, multi-modal metallicity distribution demonstrates that Terzan 5 is not a genuine globular cluster but the remnant of a much more complex stellar system.

The achievable level of precision on photospheric abundances of stars is a major limiting factor on investigations of exoplanet host star characteristics, the chemical histories of star clusters, and the evolution of the Milky Way and other galaxies. While model-induced errors can be minimized through the differential analysis of spectrally similar stars, the maximum achievable precision of this technique has been debated. As a test, we derive differential abundances of 19 elements from high-quality asteroid-reflected solar spectra taken using a variety of instruments and conditions. We treat the solar spectra as being from unknown stars and use the resulting differential abundances, which are expected to be zero, as a diagnostic of the error in our measurements. Our results indicate that the relative resolution of the target and reference spectra is a major consideration, with use of different instruments to obtain the two spectra leading to errors up to 0.04 dex. Use of the same instrument at different epochs for the two spectra has a much smaller effect (∼0.007 dex). The asteroid used to obtain the solar standard also has a negligible effect (∼0.006 dex). Assuming that systematic errors from the stellar model atmospheres have been minimized, as in the case of solar twins, we confirm that differential chemical abundances can be obtained at sub-0.01 dex precision with due care in the observations, data reduction, and abundance analysis.

We present high resolution transmission spectra of giant planet atmospheres from a coupled three-dimensional (3D) atmospheric dynamics and transmission spectrum model that includes Doppler shifts which arise from winds and planetary motion. We model Jovian planets covering more than two orders of magnitude in incident flux, corresponding to planets with 0.9–55 day orbital periods around solar-type stars. The results of our 3D dynamical models reveal certain aspects of high resolution transmission spectra that are not present in simple one-dimensional (1D) models. We find that the hottest planets experience strong substellar to anti-stellar (SSAS) winds, resulting in transmission spectra with net blueshifts of up to 3 km s⁻¹, whereas less irradiated planets show almost no net Doppler shifts. We find only minor differences between transmission spectra for atmospheres with temperature inversions and those without. Compared to 1D models, peak line strengths are significantly reduced for the hottest atmospheres owing to Doppler broadening from a combination of rotation (which is faster for close-in planets under the assumption of tidal locking) and atmospheric winds. Finally, high resolution transmission spectra may be useful in studying the atmospheres of exoplanets with optically thick clouds since line cores for very strong transitions should remain optically thick to very high altitude. High resolution transmission spectra are an excellent observational test for the validity of 3D atmospheric dynamics models, because they provide a direct probe of wind structures and heat circulation. Ground-based exoplanet spectroscopy is currently on the verge of being able to verify some of our modeling predictions, most notably the dependence of SSAS winds on insolation. We caution that interpretation of high resolution transmission spectra based on 1D atmospheric models may be inadequate, as 3D atmospheric motions can produce a noticeable effect on the absorption signatures.

In most theories of planet formation, the snow-line represents a boundary between the emergence of the interior rocky planets and the exterior ice giants. The wide separation of the snow-line makes the discovery of transiting worlds challenging, yet transits would allow for detailed subsequent characterization. We present the discovery of Kepler-421b, a Uranus-sized exoplanet transiting a G9/K0 dwarf once every 704.2 days in a near-circular orbit. Using public Kepler photometry, we demonstrate that the two observed transits can be uniquely attributed to the 704.2 day period. Detailed light curve analysis with BLENDER validates the planetary nature of Kepler-421b to >4σ confidence. Kepler-421b receives the same insolation as a body at ∼2 AU in the solar system, as well as a Uranian albedo, which would have an effective temperature of ∼180 K. Using a time-dependent model for the protoplanetary disk, we estimate that Kepler-421b's present semi-major axis was beyond the snow-line after ∼3 Myr, indicating that Kepler-421b may have formed at its observed location.

We present [C i] and [C ii] observations of a linear edge region in the Taurus molecular cloud, and model this region as a cylindrically symmetric photon-dominated region (PDR) exposed to a low-intensity UV radiation field. The sharp, long profile of the linear edge makes it an ideal case to test PDR models and determine cloud parameters. We compare observations of the [C i], ³P₁ → ³P₀ (492 GHz), [C i] ³P₂ → ³P₁ (809 GHz), and [C ii] ²P_3/2 → ²P_1/2 (1900 GHz) transitions, as well as the lowest rotational transitions of ¹²CO and ¹³CO, with line intensities produced by the RATRAN radiative transfer code from the results of the Meudon PDR code. We constrain the density structure of the cloud by fitting a cylindrical density function to visual extinction data. We study the effects of variation of the FUV field, ¹²C/¹³C isotopic abundance ratio, sulfur depletion, cosmic ray ionization rate, and inclination of the filament relative to the sky-plane on the chemical network of the PDR model and resulting line emission. We also consider the role of suprathermal chemistry and density inhomogeneities. We find good agreement between the model and observations, and that the integrated line intensities can be explained by a PDR model with an external FUV field of 0.05 G₀, a low ratio of ¹²C to ¹³C ∼43, a highly depleted sulfur abundance (by a factor of at least 50), a cosmic ray ionization rate (3–6) × 10⁻¹⁷ s⁻¹, and without significant effects from inclination, clumping or suprathermal chemistry.

We apply the anisotropic resonant scattering model developed to explain the presence of non-Zeeman circular polarization signals recently detected in the ¹²CO (J = 2 → 1) and (J = 1 → 0) transitions in molecular clouds to Stokes V spectra of SiO v = 1 and v = 2, (J = 1 → 0) masers commonly observed in evolved stars. It is found that the observed antisymmetric "S"- and symmetric "∪"- or "∩"-shaped spectral profiles naturally arise when the maser radiation scatters off populations of foreground molecules located outside the velocity range covered by the background maser radiation. Using typical values for the relevant physical parameters, it is estimated that magnetic field strengths on the order of a few times 15 mG are sufficient to explain the observational results found in the literature.

ALMA HCO⁺ observations of the infrared dark cloud G0.253+0.016 located in the central molecular zone of the Galaxy are presented. The 89 GHz emission is area-filling, optically thick, and sub-thermally excited. Two types of filaments are seen in absorption against the HCO⁺ emission. Broad-line absorption filaments (BLAs) have widths of less than a few arcseconds (0.07–0.14 pc), lengths of 30–50 arcsec (1.2–1.8 pc), and absorption profiles extending over a velocity range larger than 20 km s⁻¹. The BLAs are nearly parallel to the nearby G0.18 non-thermal filaments and may trace HCO⁺ molecules gyrating about highly ordered magnetic fields located in front of G0.253+0.016 or edge-on sheets formed behind supersonic shocks propagating orthogonal to our line of sight in the foreground. Narrow-line absorption filaments (NLAs) have line widths less than 20 km s⁻¹. Some NLAs are also seen in absorption in other species with high optical depth, such as HCN, and occasionally in emission where the background is faint. The NLAs, which also trace low-density, sub-thermally excited HCO⁺ molecules, are mostly seen on the blueshifted side of the emission from G0.253+0.016. If associated with the surface of G0.253+0.016, the kinematics of the NLAs indicate that the cloud surface is expanding. The decompression of entrained, milli-Gauss magnetic fields may be responsible for the re-expansion of the surface layers of G0.253+0.016 as it recedes from the Galactic center following a close encounter with Sgr A.

We report on the determination of astrometric, spin, and orbital parameters for PSR J1518+0204C, a "black widow" binary millisecond pulsar in the globular cluster (GC) M5. The accurate position and orbital parameters obtained from radio timing allowed us to search for the optical companion. By using WFC3/Hubble Space Telescope images, we identified a very faint variable star (m_F390W ≳ 24.8, m_F606W ≳ 24.3, m_F814W ≳ 23.1) located at only 025 from the pulsar's timing position. Due to its strong variability, this star is visible only in a subsample of images. However, the light curve obtained folding the available data with the orbital parameters of the pulsar shows a maximum at the pulsar inferior conjunction and a possible minimum at the pulsar superior conjunction. Furthermore, the shape of the optical modulation indicates a heating process possibly due to the pulsar wind. This is the first identification of an optical companion to a black widow pulsar in the dense stellar environment of a GC.

We perform a census of ionized gas outflows using a sample of ∼23,000 type 2 active galactic nuclei (AGNs) out to z ∼ 0.1. By measuring the velocity offset of narrow emission lines, i.e., [O iii] λ5007 and Hα, with respect to the systemic velocity measured from the stellar absorption lines, we find that 47% of AGNs display an [O iii] line-of-sight velocity offset ⩾ 20 km s⁻¹. The fraction of the [O iii] velocity offset in type 2 AGNs is comparable to that in type 1 AGNs after considering the projection effect. AGNs with a large [O iii] velocity offset preferentially have a high Eddington ratio, implying that the detected velocity offsets are related to black hole activity. The distribution of the host galaxy inclination is clearly different between the AGNs with blueshifted [O iii] and the AGNs with redshifted [O iii], supporting the combined model of the biconical outflow and dust obscuration. In addition, for ∼3% of AGNs, [O iii] and Hα show comparable large velocity offsets, indicating a more complex gas kinematics than decelerating outflows in a stratified narrow-line region.

We explore the distribution of cool (∼10⁴ K) gas around galaxies and its dependence on galaxy properties. By cross-correlating about 50,000 Mg ii absorbers with millions of sources from the SDSS (optical), WISE (IR), and GALEX (UV) surveys we effectively extract about 2000 galaxy–absorber pairs at z ∼ 0.5 and probe relations between absorption strength and galaxy type, impact parameter and azimuthal angle. We find that cool gas traced by Mg ii absorbers exists around both star-forming and passive galaxies with a similar incidence rate on scales greater than 100 kpc but each galaxy type exhibits a different behavior on smaller scales: Mg ii equivalent width does not correlate with the presence of passive galaxies whereas stronger Mg ii absorbers tend to be found in the vicinity of star-forming galaxies. This effect is preferentially seen along the minor axis of these galaxies, suggesting that some of the gas is associated with outflowing material. In contrast, the distribution of cool gas around passive galaxies is consistent with being isotropic on the same scales. We quantify the average excess Mg ii equivalent width $\langle \delta W_{0}^{\rm Mg\,\scriptsize{II}}\rangle$ as a function of galaxy properties and find $\langle \delta W_{0}^{\rm Mg\,\scriptsize{II}}\rangle \propto {\rm SFR}^{1.2}$, sSFR^0.5, and $M_\ast ^{0.4}$ for star-forming galaxies. This work demonstrates that the dichotomy between star-forming and passive galaxies is reflected in the circumgalactic medium traced by low-ionized gas. We also measure the covering fraction of Mg ii absorption and find it to be about 2–10 times higher for star-forming galaxies than passive ones within 50 kpc. We estimate the amount of neutral gas in the halo of 〈log M_*/M_☉〉 ∼ 10.8 galaxies to be a few × 10⁹ M_☉ for both types of galaxies. Finally, we find that correlations between absorbers and sources detected in the UV and IR lead to physical trends consistent with those measured in the optical.

A significant fraction of Kepler systems are closely packed, largely coplanar, and circular. We study the stability of a six-planet system, Kepler-11, to gain insights on the dynamics and formation history of such systems. Using a technique called "frequency maps" as fast indicators of long-term stability, we explore the stability of the Kepler-11 system by analyzing the neighborhood space around its orbital parameters. Frequency maps provide a visual representation of chaos and stability, and their dependence on orbital parameters. We find that the current system is stable, but lies within a few percent of several dynamically dangerous two-body mean-motion resonances. Planet eccentricities are restricted below a small value, ∼0.04, for long-term stability, but planet masses can be more than twice their reported values (thus allowing for the possibility of mass loss by past photoevaporation). Based on our frequency maps, we speculate on the origin of instability in closely packed systems. We then proceed to investigate how the system could have been assembled. The stability constraints on Kepler-11 (mainly eccentricity constraints) suggest that if the system were assembled in situ, a dissipation mechanism must have been at work to neutralize the eccentricity excitation. On the other hand, if migration was responsible for assembling the planets, there has to be little differential migration among the planets to avoid them either getting trapped into mean motion resonances, or crashing into each other.

We study the Lyα profiles of 36 spectroscopically detected Lyα-emitters (LAEs) at z ∼ 2–3, using Keck MOSFIRE to measure systemic redshifts and velocity dispersions from rest-frame optical nebular emission lines. The sample has a median optical magnitude ${\cal R}=26.0$, and ranges from ${\cal R}\simeq 23$ to ${\cal R}\gt27$, corresponding to rest-frame UV absolute magnitudes M_UV ≃ −22 to M_UV > −18.2. Dynamical masses range from M_dyn < 1.3 × 10⁸ M_☉ to M_dyn = 6.8 × 10⁹ M_☉, with a median value of M_dyn = 6.3 × 10⁸ M_☉. Thirty of the 36 Lyα emission lines are redshifted with respect to the systemic velocity with at least 1σ significance, and the velocity offset with respect to systemic Δv_Lyα is correlated with the ${\cal R}$-band magnitude, M_UV, and the velocity dispersion measured from nebular emission lines with >3σ significance: brighter galaxies with larger velocity dispersions tend to have larger values of Δv_Lyα. We also make use of a comparison sample of 122 UV-color-selected ${\cal R}<25.5$ galaxies at z ∼ 2, all with Lyα emission and systemic redshifts measured from nebular emission lines. Using the combined LAE and comparison samples for a total of 158 individual galaxies, we find that Δv_Lyα is anti-correlated with the Lyα equivalent width with 7σ significance. Our results are consistent with a scenario in which the Lyα profile is determined primarily by the properties of the gas near the systemic redshift; in such a scenario, the opacity to Lyα photons in lower mass galaxies may be reduced if large gaseous disks have not yet developed and if the gas is ionized by the harder spectrum of young, low metallicity stars.

We investigate the enrichment in elements produced by the slow neutron-capture process (s-process) in the globular clusters M4 (NGC 6121) and M22 (NGC 6656). Stars in M4 have homogeneous abundances of Fe and neutron-capture elements, but the entire cluster is enhanced in s-process elements (Sr, Y, Ba, Pb) relative to other clusters with a similar metallicity. In M22, two stellar groups exhibit different abundances of Fe and s-process elements. By subtracting the mean abundances of s-poor from s-rich stars, we derive s-process residuals or empirical s-process distributions for M4 and M22. We find that the s-process distribution in M22 is more weighted toward the heavy s-peak (Ba, La, Ce) and Pb than M4, which has been enriched mostly with light s-peak elements (Sr, Y, Zr). We construct simple chemical evolution models using yields from massive star models that include rotation, which dramatically increases s-process production at low metallicity. We show that our massive star models with rotation rates of up to 50% of the critical (break-up) velocity and changes to the preferred ¹⁷O(α, γ)²¹Ne rate produce insufficient heavy s-elements and Pb to match the empirical distributions. For models that incorporate asymptotic giant branch yields, we find that intermediate-mass yields (with a ²²Ne neutron source) alone do not reproduce the light-to-heavy s-element ratios for M4 and M22, and that a small contribution from models with a ¹³C pocket is required. With our assumption that ¹³C pockets form for initial masses below a transition range between 3.0 and 3.5 M_☉, we match the light-to-heavy s-element ratio in the s-process residual of M22 and predict a minimum enrichment timescale of between 240 and 360 Myr. Our predicted value is consistent with the 300 Myr upper limit age difference between the two groups derived from isochrone fitting.

This paper presents an analysis of the dependence of the opacity for high-energy γ-rays to γγ absorption by low-energy photons on the polarization of the γ-ray and target photons. This process has so far only been considered using the polarization-averaged γγ absorption cross section. It is demonstrated that in the case of polarized γ-ray emission, subject to source-intrinsic γγ absorption by polarized target photons, this may lead to a slight overestimation of the γγ opacity by up to ∼10% in the case of a perfectly ordered magnetic field. Thus, for realistic astrophysical scenarios with partially ordered magnetic fields, the use of the polarization-averaged γγ cross section is justified for practical purposes, such as estimates of minimum Doppler factors inferred for γ-ray bursts and blazars, based on γγ transparency arguments; this paper quantifies the small error incurred by the unpolarized-radiation approximation. Furthermore, it is shown that polarization-dependent γγ absorption of initially polarized γ-rays can lead to a slight increase in the polarization beyond the spectral break caused by γγ absorption. This amount is distinctly different from the change in polarization expected if the same spectral break were produced by a break in the underlying electron distribution. This may serve as a diagnostic of whether γγ absorption is relevant in sources such as γ-ray bursts and blazars where the γ-ray emission may be intrinsically highly polarized.

The high polarization observed in the prompt phase of some gamma-ray bursts invites extensive study of the emission mechanism. In this paper, we investigate the polarization properties of the synchrotron self-Compton (SSC) process from a highly relativistic jet. A magnetic-dominated, baryon-loaded jet ejected from the central engine travels with a large Lorentz factor. Shells with slightly different velocities collide with each other and produce shocks. The shocks accelerate electrons to a power-law distribution and, at the same time, magnify the magnetic field. Electrons move in the magnetic field and produce synchrotron photons. Synchrotron photons suffer from the Compton scattering (CS) process and then are detected by an observer located slightly off-axis. We analytically derive the formulae of photon polarization in the SSC process in two magnetic configurations: a magnetic field in the shock plane and perpendicular to the shock plane. We show that photons induced by the SSC process can be highly polarized, with the maximum polarization Π ∼ 24% in the energy band [0.5, 5] MeV. The polarization depends on the viewing angles, peaking in the plane perpendicular to the magnetic field. In the energy band [0.05, 0.5] MeV, in which most γ-ray polarimeters are active, the polarization is about twice that in the Thomson limit, reaching Π ∼ 20%. This implies that the Klein–Nishina effect, which is often neglected in the literature, should be carefully considered.

The effect of line overlap in the Lyman and Werner bands, often ignored in galactic studies of the atomic-to-molecular transition, greatly enhances molecular hydrogen self-shielding in low metallicity environments and dominates over dust shielding for metallicities below about 10% solar. We implement that effect in cosmological hydrodynamics simulations with an empirical model, calibrated against the observational data, and provide fitting formulae for the molecular hydrogen fraction as a function of gas density on various spatial scales and in environments with varied dust abundance and interstellar radiation field. We find that line overlap, while important for detailed radiative transfer in the Lyman and Werner bands, has only a minor effect on star formation on galactic scales, which, to a much larger degree, is regulated by stellar feedback.

KA1858+4850 is a narrow-line Seyfert 1 galaxy at redshift 0.078 and is among the brightest active galaxies monitored by the Kepler mission. We have carried out a reverberation mapping campaign designed to measure the broad-line region size and estimate the mass of the black hole in this galaxy. We obtained 74 epochs of spectroscopic data using the Kast Spectrograph at the Lick 3 m telescope from 2012 February to November, and obtained complementary V-band images from five other ground-based telescopes. We measured the Hβ light curve lag with respect to the V-band continuum light curve using both cross-correlation techniques (CCF) and continuum light curve variability modeling with the JAVELIN method and found rest-frame lags of $\tau _\mathrm{CCF} = 13.53^{+2.03}_{-2.32}$ days and τ_JAVELIN$= 13.15^{+1.08}_{-1.00}$ days. The Hβ rms line profile has a width of σ_line = 770  ±  49 km s⁻¹. Combining these two results and assuming a virial scale factor of f = 5.13, we obtained a virial estimate of $M_{\mathrm{BH}} = 8.06^{+1.59}_{-1.72} \times 10^6 {M}_{\odot }$ for the mass of the central black hole and an Eddington ratio of L/L_Edd ≈ 0.2. We also obtained consistent but slightly shorter emission-line lags with respect to the Kepler light curve. Thanks to the Kepler mission, the light curve of KA1858+4850 has among the highest cadences and signal-to-noise ratios ever measured for an active galactic nucleus; thus, our black hole mass measurement will serve as a reference point for relations between black hole mass and continuum variability characteristics in active galactic nuclei.

We present theoretical C iv λλ1548,1550 absorption line profiles for QSOs calculated assuming the accretion disk wind (ADW) scenario. The results suggest that the multiple absorption troughs seen in many QSOs may be due to the discontinuities in the ion balance of the wind (caused by X-rays), rather than discontinuities in the density/velocity structure. The profiles are calculated from a 2.5-dimensional time-dependent hydrodynamic simulation of a line-driven disk wind for a typical QSO black hole mass, a typical QSO luminosity, and for a standard Shakura–Sunyaev disk. We include the effects of ionizing X-rays originating from within the inner disk radius by assuming that the wind is shielded from the X-rays from a certain viewing angle up to 90° ("edge on"). In the shielded region, we assume constant ionization equilibrium, and thus constant line-force parameters. In the non-shielded region, we assume that both the line-force and the C iv populations are nonexistent. The model can account for P-Cygni absorption troughs (produced at edge on viewing angles), multiple absorption troughs (produced at viewing angles close to the angle that separates the shielded region and the non-shielded region), and for detached absorption troughs (produced at an angle in between the first two absorption line types); that is, the model can account for the general types of broad absorption lines seen in QSOs as a viewing angle effect. The steady nature of ADWs, in turn, may account for the steady nature of the absorption structure observed in multiple-trough broad absorption line QSOs. The model parameters are M_bh = 10⁹ M_☉ and L_disk = 10⁴⁷ erg s⁻¹.

We report the tentative detection in space of the nitrosylium ion, NO⁺. The observations were performed toward the cold dense core Barnard 1-b. The identification of the NO⁺J = 2–1 line is supported by new laboratory measurements of NO⁺ rotational lines up to the J = 8–7 transition (953207.189 MHz), which leads to an improved set of molecular constants: B₀ = 59597.1379(62) MHz, D₀ = 169.428(65) kHz, and eQq₀(N) = −6.72(15) MHz. The profile of the feature assigned to NO⁺ exhibits two velocity components at 6.5 and 7.5 km s⁻¹, with column densities of 1.5 × 10¹² and 6.5 × 10¹¹ cm⁻², respectively. New observations of NO and HNO, also reported here, allow us to estimate the following abundance ratios: X(NO)/X(NO⁺) ≃ 511, and X(HNO)/X(NO⁺) ≃ 1. This latter value provides important constraints on the formation and destruction processes of HNO. The chemistry of NO⁺ and other related nitrogen-bearing species is investigated by the means of a time-dependent gas phase model which includes an updated chemical network according to recent experimental studies. The predicted abundance for NO⁺ and NO is found to be consistent with the observations. However, that of HNO relative to NO is too high. No satisfactory chemical paths have been found to explain the observed low abundance of HNO. HSCN and HNCS are also reported here with an abundance ratio of ≃ 1. Finally, we have searched for NNO, NO₂, HNNO⁺, and NNOH⁺, but only upper limits have been obtained for their column density, except for the latter for which we report a tentative 3σ detection.

We have used high-cadence radial velocity measurements from the Hobby–Eberly Telescope with published velocities from the Lick 3 m Shane Telescope, combined with astrometric data from the Hubble Space Telescope (HST) Fine Guidance Sensors to refine the orbital parameters of the HD 128311 system, and determine an inclination of 5595 ± 1455 and true mass of 3.789 $^{+0.924}_{ -0.432}$M_JUP for HD 128311 c. The combined radial velocity data also reveal a short period signal which could indicate a third planet in the system with an Msin i of 0.133 ± 0.005 M_JUP or stellar phenomena. Photometry from the T12 0.8 m automatic photometric telescope at the Fairborn Observatory and HST are used to determine a photometric period close to, but not within the errors of the radial velocity signal. We performed a cross-correlation bisector analysis of the radial velocity data to look for correlations with the photometric period and found none. Dynamical integrations of the proposed system show long-term stability with the new orbital parameters of over 10 million years. Our new orbital elements do not support the claims of HD 128311 b and c being in mean motion resonance.

We present the gravitational microlensing discovery of a 4 M_Uranus planet that orbits a 0.7 M_☉ star at ≈18 AU. This is the first known analog of Uranus. Similar planets, i.e., cold ice giants, are inaccessible to either radial velocity or transit methods because of the long orbital periods, while low reflected light prevents direct imaging. We discuss how similar planets may contaminate the sample of the very short microlensing events that are interpreted as free-floating planets with an estimated rate of 1.8 per main-sequence star. Moreover, the host star has a nearby stellar (or brown dwarf) companion. The projected separation of the planet is only about three times smaller than that of the companion star, suggesting significant dynamical interactions.

We present a general framework for incorporating astrophysical information into Bayesian parameter estimation techniques used by gravitational wave data analysis to facilitate multi-messenger astronomy. Since the progenitors of transient gravitational wave events, such as compact binary coalescences, are likely to be associated with a host galaxy, improvements to the source sky location estimates through the use of host galaxy information are explored. To demonstrate how host galaxy properties can be included, we simulate a population of compact binary coalescences and show that for ∼8.5% of simulations within 200 Mpc, the top 10 most likely galaxies account for a ∼50% of the total probability of hosting a gravitational wave source. The true gravitational wave source host galaxy is in the top 10 galaxy candidates ∼10% of the time. Furthermore, we show that by including host galaxy information, a better estimate of the inclination angle of a compact binary gravitational wave source can be obtained. We also demonstrate the flexibility of our method by incorporating the use of either the B or K band into our analysis.

We present griz_P1 light curves of 146 spectroscopically confirmed Type Ia supernovae (SNe Ia; 0.03 < z < 0.65) discovered during the first 1.5 yr of the Pan-STARRS1 Medium Deep Survey. The Pan-STARRS1 natural photometric system is determined by a combination of on-site measurements of the instrument response function and observations of spectrophotometric standard stars. We find that the systematic uncertainties in the photometric system are currently 1.2% without accounting for the uncertainty in the Hubble Space Telescope Calspec definition of the AB system. A Hubble diagram is constructed with a subset of 113 out of 146 SNe Ia that pass our light curve quality cuts. The cosmological fit to 310 SNe Ia (113 PS1 SNe Ia + 222 light curves from 197 low-z SNe Ia), using only supernovae (SNe) and assuming a constant dark energy equation of state and flatness, yields $w=-1.120^{+0.360}_{-0.206}\hbox{(Stat)} ^{+0.269}_{-0.291}\hbox{(Sys)}$. When combined with BAO+CMB(Planck)+H₀, the analysis yields $\Omega _{\rm M}=0.280^{+0.013}_{-0.012}$ and $w=-1.166^{+0.072}_{-0.069}$ including all identified systematics. The value of w is inconsistent with the cosmological constant value of −1 at the 2.3σ level. Tension endures after removing either the baryon acoustic oscillation (BAO) or the H₀ constraint, though it is strongest when including the H₀ constraint. If we include WMAP9 cosmic microwave background (CMB) constraints instead of those from Planck, we find $w=-1.124^{+0.083}_{-0.065}$, which diminishes the discord to <2σ. We cannot conclude whether the tension with flat ΛCDM is a feature of dark energy, new physics, or a combination of chance and systematic errors. The full Pan-STARRS1 SN sample with ∼three times as many SNe should provide more conclusive results.

We probe the systematic uncertainties from the 113 Type Ia supernovae (SN Ia) in the Pan-STARRS1 (PS1) sample along with 197 SN Ia from a combination of low-redshift surveys. The companion paper by Rest et al. describes the photometric measurements and cosmological inferences from the PS1 sample. The largest systematic uncertainty stems from the photometric calibration of the PS1 and low-z samples. We increase the sample of observed Calspec standards from 7 to 10 used to define the PS1 calibration system. The PS1 and SDSS-II calibration systems are compared and discrepancies up to ∼0.02 mag are recovered. We find uncertainties in the proper way to treat intrinsic colors and reddening produce differences in the recovered value of w up to 3%. We estimate masses of host galaxies of PS1 supernovae and detect an insignificant difference in distance residuals of the full sample of 0.037 ± 0.031 mag for host galaxies with high and low masses. Assuming flatness and including systematic uncertainties in our analysis of only SNe measurements, we find w =$-1.120^{+0.360}_{-0.206}{\rm (Stat)} ^{+0.269}_{-0.291}{\rm (Sys)}$. With additional constraints from Baryon acoustic oscillation, cosmic microwave background (CMB) (Planck) and H₀ measurements, we find $w=-1.166^{+0.072}_{-0.069}$ and $\Omega _m=0.280^{+0.013}_{-0.012}$ (statistical and systematic errors added in quadrature). The significance of the inconsistency with w = −1 depends on whether we use Planck or Wilkinson Microwave Anisotropy Probe measurements of the CMB: $w_{{\rm BAO+H0+SN+{WMAP}}}=-1.124^{+0.083}_{-0.065}$.

A comprehensive spectral analysis of both the solar background magnetic field (SBMF) in cycles 21–23 and the sunspot magnetic field in cycle 23 reported in our recent paper showed the presence of two principal components (PCs) of SBMF having opposite polarity, e.g., originating in the northern and southern hemispheres, respectively. Over a duration of one solar cycle, both waves are found to travel with an increasing phase shift toward the northern hemisphere in odd cycles 21 and 23 and to the southern hemisphere in even cycle 22. These waves were linked to solar dynamo waves assumed to form in different layers of the solar interior. In this paper, for the first time, the PCs of SBMF in cycles 21–23 are analyzed with the symbolic regression technique using Hamiltonian principles, allowing us to uncover the underlying mathematical laws governing these complex waves in the SBMF presented by PCs and to extrapolate these PCs to cycles 24–26. The PCs predicted for cycle 24 very closely fit (with an accuracy better than 98%) the PCs derived from the SBMF observations in this cycle. This approach also predicts a strong reduction of the SBMF in cycles 25 and 26 and, thus, a reduction of the resulting solar activity. This decrease is accompanied by an increasing phase shift between the two predicted PCs (magnetic waves) in cycle 25 leading to their full separation into the opposite hemispheres in cycle 26. The variations of the modulus summary of the two PCs in SBMF reveals a remarkable resemblance to the average number of sunspots in cycles 21–24 and to predictions of reduced sunspot numbers compared to cycle 24: 80% in cycle 25 and 40% in cycle 26.

Hybrid simulations are performed to investigate the dynamics of both solar wind protons and interplanetary pickup ions (PUIs) around the corotating interaction region (CIR). The one-dimensional system is applied in order to focus on processes in the direction of CIR propagation. The CIR is bounded by forward and reverse shocks, which are responsible for particle acceleration. The effective acceleration of solar wind protons takes place when the reverse shock (fast wind side) favors a quasi-parallel regime. The diffusive process accounts for this acceleration, and particles can gain energy in a suprathermal range (on the order of 10 keV). In contrast, the PUI acceleration around the shock differs from the conventional model in which the motional electric field along the shock surface accelerates particles. Owing to their large gyroradius, PUIs can gyrate between the upstream and downstream, several proton inertial lengths away from the shock. This "cross-shock" gyration results in a net velocity increase in the field-aligned component, indicating that the magnetic mirror force is responsible for acceleration. The PUIs that remain in the vicinity of the shock for a long duration (tens of gyroperiods) gain much energy and are reflected back toward the upstream. These reflected energetic PUIs move back and forth along the magnetic field between a pair of CIRs that are magnetically connected. The PUIs are repeatedly accelerated in each reflection, leading to a maximum energy gain close to 100 keV. This mechanism can be evaluated in terms of "preacceleration" for the generation of anomalous cosmic rays.

We investigate intensity variations and energy deposition in five coronal loops in active region cores. These were selected for their strong variability in the AIA/SDO 94 Å intensity channel. We isolate the hot Fe xviii and Fe xxi components of the 94 Å and 131 Å by modeling and subtracting the "warm" contributions to the emission. HMI/SDO data allow us to focus on "inter-moss" regions in the loops. The detailed evolution of the inter-moss intensity time series reveals loops that are impulsively heated in a mode compatible with a nanoflare storm, with a spike in the hot 131 Å signals leading and the other five EUV emission channels following in progressive cooling order. A sharp increase in electron temperature tends to follow closely after the hot 131 Å signal confirming the impulsive nature of the process. A cooler process of growing emission measure follows more slowly. The Fourier power spectra of the hot 131 Å signals, when averaged over the five loops, present three scaling regimes with break frequencies near 0.1 min⁻¹ and 0.7 min⁻¹. The low frequency regime corresponds to 1/f noise; the intermediate indicates a persistent scaling process and the high frequencies show white noise. Very similar results are found for the energy dissipation in a 2D "hybrid" shell model of loop magneto-turbulence, based on reduced magnetohydrodynamics, that is compatible with nanoflare statistics. We suggest that such turbulent dissipation is the energy source for our loops.

Solar eruptions are usually associated with a variety of phenomena occurring in the low corona before, during, and after the onset of eruption. Though easily visible in coronagraph observations, so-called stealth coronal mass ejections (CMEs) do not obviously exhibit any of these low-coronal signatures. The presence or absence of distinct low-coronal signatures can be linked to different theoretical models to establish the mechanisms by which the eruption is initiated and driven. In this study, 40 CMEs without low-coronal signatures occurring in 2012 are identified. Their observational and kinematic properties are analyzed and compared to those of regular CMEs. Solar eruptions without clear on-disk or low-coronal signatures can lead to unexpected space weather impacts, since many early warning signs for significant space weather activity are not present in these events. A better understanding of their initiation mechanism(s) will considerably improve the ability to predict such space weather events.

We present the mass distribution at the core of SDSS J1531+3414, a strong-lensing cluster at z = 0.335. We find that the mass distribution is well described by two cluster-scale halos with a contribution from cluster-member galaxies. New Hubble Space Telescope observations of SDSS J1531+3414 reveal a signature of ongoing star formation associated with the two central galaxies at the core of the cluster, in the form of a chain of star forming regions at the center of the cluster. Using the lens model presented here, we place upper limits on the contribution of a possible lensed image to the flux at the central region, and rule out that this emission is coming from a background source.

Young stars accrete mass from circumstellar disks and, in many cases, the accretion coincides with a phase of massive outflows, which can be highly collimated. Those jets emit predominantly in the optical and IR wavelength range. However, in several cases, X-ray and UV observations reveal a weak but highly energetic component in those jets. X-rays are observed both from stationary regions close to the star and from knots in the jet several hundred AU from the star. In this article, we show semianalytically that a fast stellar wind that is recollimated by the pressure from a slower, more massive disk wind can have the right properties to power stationary X-ray emission. The size of the shocked regions is compatible with observational constraints. Our calculations support a wind–wind interaction scenario for the high-energy emission near the base of young stellar object jets. For the specific case of DG Tau, a stellar wind with a mass-loss rate of 5 × 10⁻¹⁰ M_☉ yr⁻¹ and a wind speed of 800 km s⁻¹ reproduces the observed X-ray spectrum. We conclude that a stellar wind recollimation shock is a viable scenario to power stationary X-ray emission close to the jet launching point.

We have obtained a sample of ∼22,000 red giant branch (RGB) stars based on stellar parameters, provided by the ninth data release of the Sloan Digital Sky Survey, and the CH(G)/MgH indices, measured from the included spectra. The Galactic rest-frame velocity of V_gsr versus longitude for the sample shows the existence of several groups of stars from globular clusters and known streams. Excluding these substructures, a sample of ∼16,000 RGB stars from the general field is used to investigate the properties of the thick disk, the inner halo, and the outer halo of our Galaxy. The metallicity and rotational velocity distributions are investigated for stars at 0 kpc < |Z| < 10 kpc. It is found that the canonical thick disk dominates at 0 kpc < |Z| < 2 kpc and its contribution becomes negligible at |Z| > 3 kpc. The MWTD is present and overlaps with the inner halo at 1 kpc < |Z| < 3 kpc. The inner halo starts at 2 kpc < |Z| < 3 kpc and becomes the dominated population for 4 kpc < |Z| < 10 kpc. For halo stars with |Z| > 5 kpc, bimodal metallicity distributions are found for 20 kpc < |Z| < 25 kpc and 35 kpc < RR < 45 kpc, which suggests a dual halo, the inner and the outer halo, as reported in Carollo et al. at low |Z| values. The peak of metallicity for the inner halo is at [Fe/H] ∼ −1.6 and appears to be at [Fe/H] ∼ −2.3 for the outer halo. The transition point from the inner to the outer halo is located at |Z| ∼ 20 kpc and RR ∼ 35 kpc.

We study particle trapping at the edge of a gap opened by a planet in a protoplanetary disk. In particular, we explore the effects of turbulence driven by the magnetorotational instability on particle trapping, using global three-dimensional magnetohydrodynamic (MHD) simulations including Lagrangian dust particles. We study disks either in the ideal MHD limit or dominated by ambipolar diffusion (AD) which plays an essential role at the outer regions of a protoplanetary disk. With ideal MHD, strong turbulence (the equivalent viscosity parameter α ∼ 10⁻²) in disks prevents vortex formation at the edge of the gap opened by a 9 M_J planet, and most particles (except the particles that drift fastest) pile up at the outer gap edge almost axisymmetrically. When AD is considered, turbulence is significantly suppressed (α ≲ 10⁻³), and a large vortex forms at the edge of the planet induced gap, which survives ∼1000 orbits. The vortex can efficiently trap dust particles that span 3 orders of magnitude in size within 100 planetary orbits. We have also carried out two-dimensional hydrodynamical (HD) simulations using viscosity as an approximation to MHD turbulence. These HD simulations can reproduce vortex generation at the gap edge as seen in MHD simulations. Finally, we use our simulation results to generate synthetic images for ALMA dust continuum observations on Oph IRS 48 and HD 142527, which show good agreement with existing observations. Predictions for future ALMA cycle 2 observations have been made. We conclude that the asymmetry in ALMA observations can be explained by dust trapping vortices and the existence of vortices could be the evidence that the outer protoplanetary disks are dominated by AD with α < 10⁻³ at the disk midplane.

We have obtained deep images of the highly isolated (d = 1 Mpc) Aquarius dwarf irregular galaxy (DDO 210) with the Hubble Space Telescope Advanced Camera for Surveys. The resulting color–magnitude diagram (CMD) reaches more than a magnitude below the oldest main-sequence turnoff, allowing us to derive the star formation history (SFH) over the entire lifetime of the galaxy with a timing precision of ≈10% of the lookback time. Using a maximum likelihood fit to the CMD we find that only ≈10% of all star formation in Aquarius took place more than 10 Gyr ago (lookback time equivalent to redshift z ≈ 2). The star formation rate increased dramatically ≈6–8 Gyr ago (z ≈ 0.7–1.1) and then declined until the present time. The only known galaxy with a more extreme confirmed delay in star formation is Leo A, a galaxy of similar M_H i/M_⋆, dynamical mass, mean metallicity, and degree of isolation. The delayed stellar mass growth in these galaxies does not track the mean dark matter accretion rate from CDM simulations. The similarities between Leo A and Aquarius suggest that if gas is not removed from dwarf galaxies by interactions or feedback, it can linger for several gigayears without cooling in sufficient quantity to form stars efficiently. We discuss possible causes for the delay in star formation including suppression by reionization and late-time mergers. We find reasonable agreement between our measured SFHs and select cosmological simulations of isolated dwarfs. Because star formation and merger processes are both stochastic in nature, delayed star formation in various degrees is predicted to be a characteristic (but not a universal) feature of isolated small galaxies.

The spatial morphology and dynamical status of a young, still-forming stellar cluster provide valuable clues to the conditions during the star formation event and the processes that regulated it. We analyze the Orion Nebula Cluster (ONC), utilizing the latest censuses of its stellar content and membership estimates over a large wavelength range. We determine the center of mass of the ONC and study the radial dependence of angular substructure. The core appears rounder and smoother than the outskirts, which is consistent with a higher degree of dynamical processing. At larger distances, the departure from circular symmetry is mostly driven by the elongation of the system, with very little additional substructure, indicating a somewhat evolved spatial morphology or an expanding halo. We determine the mass density profile of the cluster, which is well fitted by a power law that is slightly steeper than a singular isothermal sphere. Together with the interstellar medium density, which is estimated from average stellar extinction, the mass content of the ONC is insufficient by a factor ∼1.8 to reproduce the observed velocity dispersion from virialized motions, in agreement with previous assessments that the ONC is moderately supervirial. This may indicate recent gas dispersal. Based on the latest estimates for the age spread in the system and our density profiles, we find that at the half-mass radius, 90% of the stellar population formed within ∼5–8 free-fall times (t_ff). This implies a star formation efficiency per t_ff of _ff ∼ 0.04–0.07 (i.e., relatively slow and inefficient star formation rates during star cluster formation).

The spectrum of methanol (CH₃OH) has been characterized between 214.6 and 265.4 GHz for astrophysically significant temperatures. Four hundred and eighty-six spectra with absolute intensity calibration recorded between 240 and 389 K provided a means for the calculation of the complete experimental spectrum (CES) of methanol as a function of temperature. The CES includes contributions from v_t = 3 and other higher states that are difficult to model quantum mechanically (QM). It also includes the spectrum of the ¹³C isotopologue in terrestrial abundance. In general the QM models provide frequencies that are within 1 MHz of their experimental values, but there are several outliers that differ by tens of MHz. As in our recent work on methanol in the 560–654 GHz region, significant intensity differences between our experimental intensities and cataloged values were found. In this work these differences are explored in the context of several QM analyses. The experimental results presented here are analyzed to provide a frequency point-by-point catalog that is well suited for the simulation of crowded and overlapped spectra. Additionally, a catalog in the usual line frequency, line strength, and lower state energy format is provided.

Weak spectral features in BL Lacertae objects (BL Lacs) often provide a unique opportunity to probe the inner region of this rare type of active galactic nucleus. We present a Hubble Space Telescope/Cosmic Origins Spectrograph observation of the BL Lac H 2356–309. A weak Lyα emission line was detected. This is the fourth detection of a weak Lyα emission feature in the ultraviolet (UV) band in the so-called high-energy peaked BL Lacs, after Stocke et al. Assuming the line-emitting gas is located in the broad line region (BLR) and the ionizing source is the off-axis jet emission, we constrain the Lorentz factor (Γ) of the relativistic jet to be ⩾8.1 with a maximum viewing angle of 36. The derived Γ is somewhat larger than previous measurements of Γ ≈ 3–5, implying a covering factor of ∼3% of the line-emitting gas. Alternatively, the BLR clouds could be optically thin, in which case we constrain the BLR warm gas to be ∼10⁻⁵ M_☉. We also detected two H i and one O vi absorption lines that are within |Δv| < 150 km s⁻¹ of the BL Lac object. The O vi and one of the H i absorbers likely coexist due to their nearly identical velocities. We discuss several ionization models and find a photoionization model where the ionizing photon source is the BL Lac object that can fit the observed ion column densities with reasonable physical parameters. This absorber can either be located in the interstellar medium of the host galaxy or in the BLR.

We present optical and near-infrared (NIR) imaging data of the radio-loud, narrow-line Seyfert 1 galaxy 1H 0323+342, which shows intense and variable gamma-ray activity discovered by the Fermi satellite with the Large Area Telescope. Near-infrared and optical images are used to investigate the structural properties of the host galaxy of 1H 0323+342; this together with optical spectroscopy allows us to examine its black hole mass. Based on two-dimensional (2D) multiwavelength surface-brightness modeling, we find that statistically, the best model fit is a combination of a nuclear component and a Sérsic profile (n ∼ 2.8). However, the presence of a disk component (with a small bulge n ∼ 1.2) also remains a possibility and cannot be ruled out with the present data. Although at first glance a spiral-arm-like structure is revealed in our images, a 2D Fourier analysis of the imagery suggests that this structure corresponds to an asymmetric ring, likely associated with a recent violent dynamical interaction. We discuss our results in the context of relativistic jet production and galaxy evolution.

We introduce the equations of magneto-quantum-radiative hydrodynamics. By rewriting them in a dimensionless form, we obtain a set of parameters that describe scale-dependent ratios of characteristic hydrodynamic quantities. We discuss how these dimensionless parameters relate to the scaling between astrophysical observations and laboratory experiments.

Using the known detection limits for high-resolution imaging observations and the statistical properties of true binary and line-of-sight companions, we estimate the binary fraction of Kepler exoplanet host stars. Our speckle imaging programs at the WIYN 3.5 m and Gemini North 8.1 m telescopes have observed over 600 Kepler objects of interest and detected 49 stellar companions within ∼1 arcsec. Assuming binary stars follow a log-normal period distribution for an effective temperature range of 3000–10,000 K, then the model predicts that the vast majority of detected sub-arcsecond companions are long period (P > 50 yr), gravitationally bound companions. In comparing the model predictions to the number of real detections in both observational programs, we conclude that the overall binary fraction of host stars is similar to the 40%–50% rate observed for field stars.

We improve on our previous treatments of the long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e., PdV work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead zone. The GI-induced spiral waves propagate well inside of the gravitationally unstable region before they trigger outbursts at R ≲ 1 AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local α treatments for GI. In our standard model where zero dead-zone residual viscosity (α_rd) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as 0.01 in terms of the generic α parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite α_rd triggers thermally driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when α_rd = 0. The inclusion of non-zero residual viscosity in the dead zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a "quiescent" accretion phase with $\dot{M}_{\rm acc}\sim 10^{-8}\hbox{--}10^{-7} \,M_{\odot }\,{\rm yr}^{-1}$ over 60% of the time and spend less than 15% of the infall phase in accretion outbursts. While our models indicate that episodic mass accretion during protostellar evolution can qualitatively help explain the low accretion luminosities seen in most low-mass protostars, detailed tests of the mechanism will require model calculations for a range of protostellar masses with some constraint on the initial core angular momentum, which affects the length of time spent in a quasi-steady disk accretion phase.

We study how active galactic nucleus (AGN) activity changes across environments from galaxy pairs to clusters using 143,843 galaxies with z < 0.2 from the Sloan Digital Sky Survey. Using a refined technique, we apply a continuous measure of AGN activity, characteristic of the ionization state of the narrow-line emitting gas. Changes in key emission-line ratios ([N ii] λ6548/Hα, [O iii] λ5007/Hβ) between different samples allow us to disentangle different environmental effects while removing contamination. We confirm that galaxy interactions enhance AGN activity. However, conditions in the central regions of clusters are inhospitable for AGN activity even if galaxies are in pairs. These results can be explained through models of gas dynamics in which pair interactions stimulate the transfer of gas to the nucleus and clusters suppress gas availability for accretion onto the central black hole.

Strong singly ionized magnesium (Mg ii) absorption lines in quasar spectra typically serve as a proxy for intervening galaxies along the line of sight. Previous studies have found a correlation between the number of these Mg ii absorbers and the Faraday rotation measure (RM) at ≈5 GHz. We cross-match a sample of 35,752 optically identified non-intrinsic Mg ii absorption systems with 25,649 polarized background radio sources for which we have measurements of both the spectral index and RM at 1.4 GHz. We use the spectral index to split the resulting sample of 599 sources into flat-spectrum and steep-spectrum subsamples. We find that our flat-spectrum sample shows significant (∼3.5σ) evidence for a correlation between Mg ii absorption and RM at 1.4 GHz, while our steep-spectrum sample shows no such correlation. We argue that such an effect cannot be explained by either luminosity or other observational effects, by evolution in another confounding variable, by wavelength-dependent polarization structure in an active galactic nucleus, by the Galactic foreground, by cosmological expansion, or by partial coverage models. We conclude that our data are most consistent with intervenors directly contributing to the Faraday rotation along the line of sight, and that the intervening systems must therefore have coherent magnetic fields of substantial strength ($\bar{B}=1.8\pm 0.4$ μG). Nevertheless, the weak nature of the correlation will require future high-resolution and broadband radio observations in order to place it on a much firmer statistical footing.

No true extrasolar Earth analog is known. Hundreds of planets have been found around Sun-like stars that are either Earth-sized but on shorter periods, or else on year-long orbits but somewhat larger. Under strong assumptions, exoplanet catalogs have been used to make an extrapolated estimate of the rate at which Sun-like stars host Earth analogs. These studies are complicated by the fact that every catalog is censored by non-trivial selection effects and detection efficiencies, and every property (period, radius, etc.) is measured noisily. Here we present a general hierarchical probabilistic framework for making justified inferences about the population of exoplanets, taking into account survey completeness and, for the first time, observational uncertainties. We are able to make fewer assumptions about the distribution than previous studies; we only require that the occurrence rate density be a smooth function of period and radius (employing a Gaussian process). By applying our method to synthetic catalogs, we demonstrate that it produces more accurate estimates of the whole population than standard procedures based on weighting by inverse detection efficiency. We apply the method to an existing catalog of small planet candidates around G dwarf stars. We confirm a previous result that the radius distribution changes slope near Earth's radius. We find that the rate density of Earth analogs is about 0.02 (per star per natural logarithmic bin in period and radius) with large uncertainty. This number is much smaller than previous estimates made with the same data but stronger assumptions.

We apply hydrodynamic evaporation models to different synthetic planet populations that were obtained from a planet formation code based on the core-accretion paradigm. We investigated the evolution of the planet populations using several evaporation models, which are distinguished by the driving force of the escape flow (X-ray or EUV), the heating efficiency in energy-limited evaporation regimes, or both. Although the mass distribution of the planet populations is barely affected by evaporation, the radius distribution clearly shows a break at approximately 2 R_⊕. We find that evaporation can lead to a bimodal distribution of planetary sizes and to an "evaporation valley" running diagonally downward in the orbital distance—planetary radius plane, separating bare cores from low-mass planets that have kept some primordial H/He. Furthermore, this bimodal distribution is related to the initial characteristics of the planetary populations because low-mass planetary cores can only accrete small primordial H/He envelopes and their envelope masses are proportional to their core masses. We also find that the population-wide effect of evaporation is not sensitive to the heating efficiency of energy-limited description. However, in two extreme cases, namely without evaporation or with a 100% heating efficiency in an evaporation model, the final size distributions show significant differences; these two scenarios can be ruled out from the size distribution of Kepler candidates.

We report the discovery of star formation activity in perhaps the most distant molecular cloud in the extreme outer galaxy. We performed deep near-infrared imaging with the Subaru 8.2 m telescope, and found two young embedded clusters at two CO peaks of "Digel Cloud 1" at the kinematic distance of D = 16 kpc (Galactocentric radius R_G = 22 kpc). We identified 18 and 45 cluster members in the two peaks, and the estimated stellar densities are ∼5 and ∼3 pc⁻², respectively. The observed K-band luminosity function suggests that the age of the clusters is less than 1 Myr and also that the distance to the clusters is consistent with the kinematic distance. On the sky, Cloud 1 is located very close to the H i peak of high-velocity cloud Complex H, and there are some H i intermediate velocity structures between the Complex H and the Galactic disk, which could indicate an interaction between them. We suggest the possibility that Complex H impacting on the Galactic disk has triggered star formation in Cloud 1 as well as the formation of the Cloud 1 molecular cloud.

The Earth's obliquity is stabilized by the Moon, which facilitates a rapid precession of the Earth's spin axis, detuning the system away from resonance with orbital modulation. It is, however, likely that the architecture of the solar system underwent a dynamical instability-driven transformation, where the primordial configuration was more compact. Hence, the characteristic frequencies associated with orbital perturbations were likely faster in the past, potentially allowing for secular resonant encounters. In this work, we examine if, at any point in the Earth's evolutionary history, the obliquity varied significantly. Our calculations suggest that even though the orbital perturbations were different, the system nevertheless avoided resonant encounters throughout its evolution. This indicates that the Earth obtained its current obliquity during the formation of the Moon.

Using multiwavelength imaging observations, in EUV, white light and radio, and radio spectral data over a large frequency range, we analyzed the triggering and development of a complex eruptive event. This one includes two components, an eruptive jet and a coronal mass ejection (CME), which interact during more than 30 minutes, and can be considered as physically linked. This was an unusual event. The jet is generated above a typical complex magnetic configuration that has been investigated in many former studies related to the build-up of eruptive jets; this configuration includes fan-field lines originating from a corona null point above a parasitic polarity, which is embedded in one polarity region of a large active region. The initiation and development of the CME, observed first in EUV, does not show usual signatures. In this case, the eruptive jet is the main actor of this event. The CME appears first as a simple loop system that becomes destabilized by magnetic reconnection between the outer part of the jet and the ambient medium. The progression of the CME is closely associated with the occurrence of two successive type II bursts from a distinct origin. An important part of this study is the first radio type II burst for which the joint spectral and imaging observations were allowed: (1) to follow, step by step, the evolution of the spectrum and of the trajectory of the radio burst, in relationship with the CME evolution and (2) to obtain, without introducing an electronic density model, the B field and the Alfvén speed.

The ratio of the rms electron density fluctuations to the background density in the solar wind (density modulation index, _N ≡ ΔN/N) is of vital importance for understanding several problems in heliospheric physics related to solar wind turbulence. In this paper, we have investigated the behavior of _N in the inner heliosphere from 0.26 to 0.82 AU. The density fluctuations ΔN have been deduced using extensive ground-based observations of interplanetary scintillation at 327 MHz, which probe spatial scales of a few hundred kilometers. The background densities (N) have been derived using near-Earth observations from the Advanced Composition Explorer. Our analysis reveals that 0.001 ≲ _N ≲ 0.02 and does not vary appreciably with heliocentric distance. We also find that _N declines by 8% from 1998 to 2008. We discuss the impact of these findings on problems ranging from our understanding of Forbush decreases to the behavior of the solar wind dynamic pressure over the recent peculiar solar minimum at the end of cycle 23.

We report here our latest search for molecular outflows from young brown dwarfs and very low-mass stars in nearby star-forming regions. We have observed three sources in Taurus with the Submillimeter Array and the Combined Array for Research in Millimeter-wave Astronomy at 230 GHz frequency to search for CO J = 2 → 1 outflows. We obtain a tentative detection of a redshifted and extended gas lobe at about 10 arcsec from the source GM Tau, a young brown dwarf in Taurus with an estimated mass of 73 M_J, which is right below the hydrogen-burning limit. No blueshifted emission around the brown dwarf position is detected. The redshifted gas lobe that is elongated in the northeast direction suggests a possible bipolar outflow from the source with a position angle of about 36°. Assuming that the redshifted emission is outflow emission from GM Tau, we then estimate a molecular outflow mass in the range from 1.9 × 10⁻⁶ M_☉ to 2.9 × 10⁻⁵ M_☉ and an outflow mass-loss rate from 2.7 × 10⁻⁹ M_☉ yr⁻¹ to 4.1 × 10⁻⁸ M_☉ yr⁻¹. These values are comparable to those we have observed in the young brown dwarf ISO-Oph 102 of 60 M_J in ρ Ophiuchi and the very low-mass star MHO 5 of 90 M_J in Taurus. Our results suggest that the outflow process in very low-mass objects is episodic with a duration of a few thousand years and the outflow rate of active episodes does not significantly change for different stages of the formation process of very low-mass objects. This may provide us with important implications that clarify the formation process of brown dwarfs.

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165–2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

Compact binary millisecond pulsars with main-sequence donors, often referred to as "redbacks," constitute the long-sought link between low-mass X-ray binaries and millisecond radio pulsars and offer a unique probe of the interaction between pulsar winds and accretion flows. We present a systematic study of eight nearby redbacks, using more than 100 observations obtained with Swift's X-ray Telescope. We distinguish between three main states: pulsar, disk, and outburst states. We find X-ray mode switching in the disk state of PSR J1023+0038 and XSS J12270–4859, similar to what was found in the other redback that showed evidence for accretion: rapid, recurrent changes in X-ray luminosity (0.5–10 keV, L_X), between (6–9) × 10³² erg s⁻¹ (disk-passive state) and (3–5) × 10³³ erg s⁻¹ (disk-active state). This strongly suggests that mode switching—which has not been observed in quiescent low-mass X-ray binaries—is universal among redback millisecond pulsars in the disk state. We briefly explore the implications for accretion disk truncation and find that the inferred magnetospheric radius in the disk state of PSR J1023+0038 and XSS J12270–4859 lies outside the light cylinder. Finally, we note that all three redbacks that have developed accretion disks have relatively high L_X in the pulsar state (>10³² erg s⁻¹).

We present new deep, high-resolution radio images of the diffuse minihalo in the cool core of the galaxy cluster RX J1720.1+2638. The images have been obtained with the Giant Metrewave Radio Telescope at 317, 617, and 1280 MHz and with the Very Large Array at 1.5, 4.9, and 8.4 GHz, with angular resolutions ranging from 1'' to 10''. This represents the best radio spectral and imaging data set for any minihalo. Most of the radio flux of the minihalo arises from a bright central component with a maximum radius of ∼80 kpc. A fainter tail of emission extends out from the central component to form a spiral-shaped structure with a length of ∼230 kpc, seen at frequencies 1.5 GHz and below. We find indication of a possible steepening of the total radio spectrum of the minihalo at high frequencies. Furthermore, a spectral index image shows that the spectrum of the diffuse emission steepens with increasing distance along the tail. A striking spatial correlation is observed between the minihalo emission and two cold fronts visible in the Chandra X-ray image of this cool core. These cold fronts confine the minihalo, as also seen in numerical simulations of minihalo formation by sloshing-induced turbulence. All these observations favor the hypothesis that the radio-emitting electrons in cluster cool cores are produced by turbulent re-acceleration.

We investigated the relation between compact jet emission and X-ray variability properties of all black hole transients with multiwavelength coverage during their outburst decays. We studied the evolution of all power spectral components (including low-frequency quasi-periodic oscillations; QPOs), and related this evolution to changes in jet properties tracked by radio and infrared observations. We grouped sources according to their tracks in the radio/X-ray luminosity relation and show that the standards show stronger broadband X-ray variability than outliers at a given X-ray luminosity when the compact jet turns on. This trend is consistent with the internal shock model and can be important for the understanding of the presence of tracks in the radio/X-ray luminosity relation. We also observed that the total and the QPO rms amplitudes increase together during the earlier part of the outburst decay, but after the compact jet turns, either the QPO disappears or its rms amplitude decreases significantly while the total rms amplitudes remain high. We discuss these results with a scenario including a variable corona and a non-variable disk with a mechanism for the QPO separate from the mechanism that creates broad components. Finally, we evaluated the timing predictions of the magnetically dominated accretion flow model that can explain the presence of tracks in the radio/X-ray luminosity relation.

The sensitivity of blind gamma-ray pulsar searches in multiple years worth of photon data, as from the Fermi LAT, is primarily limited by the finite computational resources available. Addressing this "needle in a haystack" problem, here we present methods for optimizing blind searches to achieve the highest sensitivity at fixed computing cost. For both coherent and semicoherent methods, we consider their statistical properties and study their search sensitivity under computational constraints. The results validate a multistage strategy, where the first stage scans the entire parameter space using an efficient semicoherent method and promising candidates are then refined through a fully coherent analysis. We also find that for the first stage of a blind search incoherent harmonic summing of powers is not worthwhile at fixed computing cost for typical gamma-ray pulsars. Further enhancing sensitivity, we present efficiency-improved interpolation techniques for the semicoherent search stage. Via realistic simulations we demonstrate that overall these optimizations can significantly lower the minimum detectable pulsed fraction by almost 50% at the same computational expense.

Our knowledge of the diffuse emission that encompasses active regions is very limited. In this paper we investigate two off-limb active regions, namely, AR 10939 and AR 10961, to probe the underlying heating mechanisms. For this purpose, we have used spectral observations from Hinode/EIS and employed the emission measure (EM) technique to obtain the thermal structure of these diffuse regions. Our results show that the characteristic EM distributions of the diffuse emission regions peak at log  T = 6.25 and the coolward slopes are in the range 1.4–3.3. This suggests that both low- as well as high-frequency nanoflare heating events are at work. Our results provide additional constraints on the properties of these diffuse emission regions and their contribution to the background/foreground when active region cores are observed on-disk.

We present a numerical study in resistive magnetohydrodynamics (MHD) where the initial equilibrium configuration contains adjacent, oppositely directed, parallel current channels. Since oppositely directed current channels repel, the equilibrium is liable to an ideal magnetohydrodynamic tilt instability. This tilt evolution, previously studied in planar settings, involves two magnetic islands or flux ropes, which on Alfvénic timescales undergo a combined rotation and separation. This in turn leads to the creation of (near) singular current layers, posing severe challenges to numerical approaches. Using our open-source grid-adaptive MPI-AMRVAC software, we revisit the planar evolution case in compressible MHD, as well as its extension to two-and-a-half-dimensional (2.5D) and full three-dimensional (3D) scenarios. As long as the third dimension can be ignored, pure tilt evolutions result that are hardly affected by out of plane magnetic field components. In all 2.5D runs, our simulations do show secondary tearing type disruptions throughout the near singular current sheets in the far nonlinear saturation regime. In full 3D runs, both current channels can be liable to additional ideal kink deformations. We discuss the effects of having both tilt and kink instabilities acting simultaneously in the violent, reconnection-dominated evolution. In 3D, both the tilt and the kink instabilities can be stabilized by tension forces. As a concrete space plasma application, we argue that interacting tilt-kink instabilities in repelling current channels provide a novel route to initiate solar coronal mass ejections, distinctly different from the currently favored pure kink or torus instability routes.

The standard, or fast, solutions of m-CAK line-driven wind theory cannot account for slowly outflowing disks like the ones that surround Be stars. It has been previously shown that there exists another family of solutions—the Ω-slow solutions—that is characterized by much slower terminal velocities and higher mass-loss rates. We have solved the one-dimensional m-CAK hydrodynamical equation of rotating radiation-driven winds for this latter solution, starting from standard values of the line force parameters (α, k, and δ), and then systematically varying the values of α and k. Terminal velocities and mass-loss rates that are in good agreement with those found in Be stars are obtained from the solutions with lower α and higher k values. Furthermore, the equatorial densities of such solutions are comparable to those that are typically assumed in ad hoc models. For very high values of k, we find that the wind solutions exhibit a new kind of behavior.

We consider the lifetimes of spots on the Sun and other stars from the standpoint of magnetic diffusion. While normal magnetic diffusivity predicts lifetimes of sunspots that are too large by at least two orders of magnitude, turbulent magnetic diffusivity accounts for both the functional form of the solar empirical spot–lifetime relation and for the observed sunspot lifetimes, provided that the relevant diffusion length is the supergranule size. Applying this relation to other stars, the value of turbulent diffusivity depends almost entirely on supergranule size, with very weak dependence on other variables such as magnetic field strength and density. Overall, the best observational data for other stars is consistent with the extension of the solar relation, provided that stellar supergranule sizes for some stars are significantly larger than they are on the Sun.

In this paper we derive some first order differential equations which model the classical and the relativistic thin layer approximations. The circumstellar medium is assumed to follow a density profile of the Plummer type, the Lane–Emden (n = 5) type, or a power law. The first order differential equations are solved analytically, numerically, by a series expansion, or by recursion. The initial conditions are chosen in order to model the temporal evolution of SN 1993J over 10 yr and a smaller chi-squared is obtained for the Plummer case with η = 6. The stellar mass ejected by the SN progenitor prior to the explosion, expressed in solar mass, is identified with the total mass associated with the selected density profile and varies from 0.217 to 0.402 when the central number density is 10⁷ particles per cubic centimeter. The FWHM of the three density profiles, which can be identified with the size of the pre-SN 1993J envelope, varies from 0.0071 pc to 0.0092 pc.

Accurate mass-loss rate estimates are crucial keys in the study of wind properties of massive stars and for testing different evolutionary scenarios. From a theoretical point of view, this implies solving a complex set of differential equations in which the radiation field and the hydrodynamics are strongly coupled. The use of an analytical expression to represent the radiation force and the solution of the equation of motion has many advantages over numerical integrations. Therefore, in this work, we present an analytical expression as a solution of the equation of motion for radiation-driven winds in terms of the force multiplier parameters. This analytical expression is obtained by employing the line acceleration expression given by Villata and the methodology proposed by Müller & Vink. On the other hand, we find useful relationships to determine the parameters for the line acceleration given by Müller & Vink in terms of the force multiplier parameters.

A new set of theoretical Hα emission line profiles of Be stars has been computed using the code Beray, which solves the transfer equation along a series of rays passing through the star+disk system, representing an improved treatment over earlier work done by the authors. The new profiles were compared with the previous work, and general trends (such as line profile shapes and correlations between line equivalent widths as a function of initial density ρ₀ and power law index n) were recovered. Additionally, Beray was employed to model the spectra of eight well-known Be shell stars. Some degeneracy was found in the choice of model parameters, highlighting the need to employ alternate observables to constrain the models. However, the inclination angle of the model seemed relatively insensitive to the choices of other parameters, and we show that, with our models, only a very small range of inclination angles can adequately reproduce the observations. Five of our eight targets were found to have inclination angles of 70° or higher, and two more were found to have inclination angles of 67° and 65°. The observation of one target—4 Aquilae—could only be reproduced by models created at an inclination angle of approximately 45°.

We report deep long-slit emission-line spectra, the line flux ratios, and Doppler profile shapes of various bright optical lines. The low-ionization lines (primarily [N i], [O i], [S ii], and [N ii]) originate in shocked knots, as reported by many previous observers. Dust-scattered lines of higher ionization are seen throughout the lobes but do not peak in the knots. Our analysis of these line profiles and the readily discernible stellar continuum shows that (1) the central star is an active symbiotic (whose spectrum resembles the central stars of highly bipolar and young planetary nebulae such as M2–9 and Hen2–437) whose compact companion shows a WC8-type spectrum, (2) extended nebular lines of [O iii] and He i originate in the heavily obscured nuclear H ii region, and (3) the Balmer lines observed throughout the lobes are dominated by reflected Hα emission from the symbiotic star. Comparing our line ratios with those observed historically shows that (1) the [O iii]/Hβ and He i/Hβ ratios have been steadily rising by large amounts throughout the nebula, (2) the Hα/Hβ ratio is steadily decreasing while Hγ/Hβ remains nearly constant, and (3) the low-ionization line ratios formed in the shocked knots have been in decline in different ways at various locations. We show that the first two of these results might be expected if the symbiotic central star has been active and if its bright Hα line has faded significantly in the past 20 years.

We show the importance of γ and positron transport for the formation of late-time spectra in Type Ia supernovae (SNe Ia). The goal is to study the imprint of magnetic fields (B) on late-time IR line profiles, particularly the [Fe ii] feature at 1.644 μm, which becomes prominent two to three months after the explosion. As a benchmark, we use the explosion of a Chandrasekhar mass (M_Ch) white dwarf (WD) and, specifically, a delayed detonation model that can reproduce the light curves and spectra for a Branch-normal SN Ia. We assume WDs with initial magnetic surface fields between 1 and 10⁹ G. We discuss large-scale dipole and small-scale magnetic fields. We show that positron transport effects must be taken into account for the interpretation of emission features starting at about one to two years after maximum light, depending on the size of B. The [Fe ii] line profile and its evolution with time can be understood in terms of the overall energy input by radioactive decay and the transition from a γ-ray to a positron-dominated regime. We find that the [Fe ii] line at 1.644 μm can be used to analyze the overall chemical and density structure of the exploding WD up to day 200 without considering B. At later times, positron transport and magnetic field effects become important. After about day 300, the line profile allows one to probe the size of the B-field. The profile becomes sensitive to the morphology of B at about day 500. In the presence of a large-scale dipole field, a broad line is produced in M_Ch mass explosions that may appear flat-topped or rounded depending on the inclination at which the SN is observed. Small or no directional dependence of the spectra is found for small-scale B. We note that narrow-line profiles require central ⁵⁶Ni as shown in our previous studies. Persistent broad-line, flat-topped profiles require high-density burning, which is the signature of a WD close to M_Ch. Good time coverage is required to separate the effects of optical depth, the size and morphology of B, and the aspect angle of the observer. The spectra require a resolution of about 500 km s⁻¹ and a signal-to-noise ratio of about 20%. Two other strong near-IR spectral features at about 1.5 and 1.8 μm are used to demonstrate the importance of line blending, which may invalidate a kinematic interpretation of emission lines. Flat-topped line profiles between 300 and 400 days have been observed and reported in the literature. They lend support for M_Ch mass explosions in at least some cases and require magnetic fields equal to or in excess of 10⁶ G. We briefly discuss the effects of the size and morphology of B on light curves, as well as limitations. We argue that line profiles are a more direct measurement of B than light curves because they measure both the distribution of ⁵⁶Ni and the redistribution of the energy input by positrons rather than the total energy input. Finally, we discuss possible mechanisms for the formation of high B-fields and the limitations of our analysis.

The Kepler mission has released ∼4229 transiting planet candidates. There are approximately 222 candidate systems with three planets. Among them, the period ratios of planet pairs near 1.5 and 2.0 reveal that two peaks exist for which the proportions of the candidate systems are ∼7.0% and 18.0%, respectively. In this work, we study the formation of mean motion resonance (MMR) systems, particularly for the planetary configurations near 3:2 and 2:1 MMRs, and we concentrate on the interplay between the resonant configuration and the combination of stellar accretion rate, stellar magnetic field, speed of migration, and additional planets. We perform more than 1000 runs by assuming a system with a solar-like star and three surrounding planets. From the statistical results, we find that under the formation scenario, the proportions near 1.5 and 2.0 can reach 14.5% and 26.0%, respectively. In addition, $\dot{M}=0.1\times 10^{-8} \,M_\odot \ {\rm yr^{-1}}$ is propitious toward the formation of 3:2 resonance, whereas $\dot{M}=2\times 10^{-8} \,M_\odot \ {\rm yr^{-1}}$ contributes to the formation of 2:1 resonance. The speed-reduction factor of type I migration f₁ ⩾ 0.3 facilitates 3:2 MMRs, whereas f₁ ⩾ 0.1 facilitates 2:1 MMRs. If additional planets are present in orbits within the innermost or beyond the outermost planet in a three-planet system, 3:2:1 MMRs can be formed, but the original systems trapped in 4:2:1 MMRs are not affected by the supposed planets. In summary, we conclude that this formation scenario will provide a likely explanation for Kepler candidates involved in 2:1 and 3:2 MMRs.

Magnetic interactions between close-in planets and their host star can play an important role in the secular orbital evolution of the planets, as well as the rotational evolution of their host. As long as the planet orbits inside the Alfvén surface of the stellar wind, the magnetic interaction between the star and the planet can modify the wind properties and also lead to direct angular momentum transfers between the two. We model these star–planet interactions using compressible magnetohydrodynamic (MHD) simulations, and quantify the angular momentum transfers between the star, the planet, and the stellar wind. We study the cases of magnetized and non-magnetized planets and vary the orbital radius inside the Alfvén surface of the stellar wind. Based on a grid of numerical simulations, we propose general scaling laws for the modification of the stellar wind torque, for the torque between the star and the planet, and for the planet migration associated with the star–planet magnetic interactions. We show that when the coronal magnetic field is large enough and the star is rotating sufficiently slowly, the effect of the magnetic star–planet interaction is comparable to tidal effects and can lead to a rapid orbital decay.

We present the time-resolved spectral analysis of the XMM-Newton data of NGC 1365  collected during one XMM-Newton observation, which caught this "changing-look" active galactic nucleus in a high flux state characterized also by a low column density (N_H ∼ 10²² cm⁻²) of the X-ray absorber. During this observation, the low-energy photoelectric cut-off is at about ∼1 keV and the primary continuum can be investigated with the XMM-Newton-RGS data, which show strong spectral variability that can be explained as a variable low N_H that decreased from N_H ∼ 10²³ cm⁻² to 10²² cm⁻² in a 100 ks timescale. The spectral analysis of the last segment of the observation revealed the presence of several absorption features that can be associated with an ionized (log ξ ∼ 2 erg cm s⁻¹) outflowing wind (v_out ∼ 2000 km s⁻¹). We detected for the first time a possible P-Cygni profile of the Mg xii Lyα line associated with this mildly ionized absorber indicative of a wide angle outflowing wind. We suggest that this wind is a low ionization zone of the highly ionized wind present in NGC 1365, which is responsible for the iron K absorption lines and is located within the variable X-ray absorber. At the end of the observation, we detected a strong absorption line at E ∼ 0.76 keV most likely associated with a lower ionization zone of the absorber (log ξ ∼ 0.2 erg cm s⁻¹, N_H ∼ 10²² cm⁻²), which suggests that the variable absorber in NGC 1365 could be a low ionization zone of the disk wind.

The Fermi source 2FGL J0523.3−2530 has recently been identified as a candidate millisecond pulsar binary with an orbital period of 16.5 hr. We have carried out detailed studies of the source's emission properties by analyzing data taken with the Fermi Large Area Telescope in the 0.2–300 GeV energy range. Long-term, yearly variability from the source has been found, with a factor of four flux variations in 1–300 GeV. From spectral analysis, we find an extra spectral component at 2–3 GeV that causes the source brightening. While no orbital modulations have been found from the Fermi data over the whole period of 2008–2014, orbital modulation in the source's >2 GeV emission is detected during the last 1.5 yr of the Fermi observation. Our results support the millisecond pulsar binary nature of 2FGL J0523.3−2530. Multi-wavelength observations of the source are warranted in order to find any correlated flux variations and thus help determine the origin of the long-term variability, which currently is not understood.

NGC 6872, hereafter the Condor, is a large spiral galaxy that is interacting with its closest companion, the S0 galaxy IC 4970. The extent of the Condor provides an opportunity for detailed investigation of the impact of the interaction on the current star formation rate and its history across the galaxy, on the age and spatial distribution of its stellar population, and on the mechanism that drives the star formation activity. To address these issues we analyzed the far-ultraviolet (FUV) to near-infrared (near-IR) spectral energy distribution of seventeen 10 kpc diameter regions across the galaxy, and derived their star formation history, current star formation rate, and stellar population and mass. We find that most of the star formation takes place in the extended arms, with very little star formation in the central 5 kpc of the galaxy, in contrast to what was predicted from previous numerical simulations. There is a trend of increasing star formation activity with distance from the nucleus of the galaxy, and no evidence for a recent increase in the current star formation rate due to the interaction. The nucleus itself shows no significant current star formation activity. The extent of the Condor also provides an opportunity to test the applicability of a single standard prescription for conversion of the FUV + IR (22 μm) intensities to a star formation rate for all regions. We find that the conversion factor differs from region to region, arising from regional differences in the stellar populations.

We investigate the use of spiral arm pitch angles as a probe of disk galaxy mass profiles. We confirm our previous result that spiral arm pitch angles (P) are well correlated with the rate of shear (S) in disk galaxy rotation curves. We use this correlation to argue that imaging data alone can provide a powerful probe of galactic mass distributions out to large look-back times. We then use a sample of 13 galaxies, with Spitzer 3.6 μm imaging data and observed Hα rotation curves, to demonstrate how an inferred shear rate coupled with a bulge–disk decomposition model and a Tully–Fisher-derived velocity normalization can be used to place constraints on a galaxy's baryon fraction and dark matter halo profile. Finally, we show that there appears to be a trend (albeit a weak correlation) between spiral arm pitch angle and halo concentration. We discuss implications for the suggested link between supermassive black hole (SMBH) mass and dark halo concentration, using pitch angle as a proxy for SMBH mass.

The very high-energy (VHE) gamma-ray spectral index of high-energy peaked blazars correlates strongly with its corresponding redshift, whereas no such correlation is observed in the X-ray or GeV bands. We attribute this correlation to photon–photon absorption of TeV photons with the extragalactic background light (EBL), and utilizing this we compute the allowed flux range for the EBL, which is independent of previous estimates. The observed VHE spectrum of the sources in our sample can be well approximated by a power law, and if the de-absorbed spectrum is also assumed to be a power law, then we show that the spectral shape of EBL will be n() ∼ klog (/_p). We estimate the range of values for the parameters defining the EBL spectrum, k and _p, such that the correlation of the intrinsic VHE spectrum with redshift is nullified. The estimated EBL depends only on the observed correlation and the assumption of a power-law source spectrum. Specifically, it does not depend on the spectral modeling or radiative mechanism of the sources or on any theoretical shape of the EBL spectrum obtained through cosmological calculations. The estimated EBL spectrum is consistent with the upper and lower limits imposed by different observations. Moreover, it also agrees closely with the theoretical estimates obtained through cosmological evolution models.

We apply the mathematical formalism of vector spherical harmonics decomposition to convective stellar velocity fields from multidimensional hydrodynamics simulations and show that the resulting power spectra furnish a robust and stable statistical description of stellar convective turbulence. Analysis of the power spectra helps identify key physical parameters of the convective process such as the dominant scale of the turbulent motions that influence the structure of massive evolved pre-supernova stars. We introduce the numerical method that can be used to calculate vector spherical harmonics power spectra from two-dimensional (2D) and three-dimensional (3D) convective shell simulation data. Using this method we study the properties of oxygen shell burning and convection for a 15 M_☉ star simulated by the hydrodynamics code FLASH in 2D and 3D. We discuss the importance of realistic initial conditions to achieving successful core-collapse supernova explosions in multidimensional simulations. We show that the calculated power spectra can be used to generate realizations of the velocity fields of presupernova convective shells. We find that the slope of the solenoidal mode power spectrum remains mostly constant throughout the evolution of convection in the oxygen shell in both 2D and 3D simulations. We also find that the characteristic radial scales of the convective elements are smaller in 3D than in 2D, while the angular scales are larger in 3D.

We report the discovery of 24 Lyman-break candidates at 7 ≲ z ≲ 10.5, in the Hubble Frontier Fields (HFF) imaging data of A2744 (z = 0.308), plus Spitzer/IRAC data and archival ACS data. The sample includes a triple image system with a photometric redshift of z ≃ 7.4. This high redshift is geometrically confirmed by our lens model corresponding to deflection angles that are 12% larger than the lower-redshift systems used to calibrate the lens model at z = 2.019. The majority of our high-redshift candidates are not expected to be multiply lensed given their locations in the image plane and the brightness of foreground galaxies, but are magnified by factors of ∼1.3–15, so that we are seeing further down the luminosity function than comparable deep-field imaging. It is apparent that the redshift distribution of these sources does not smoothly extend over the full redshift range accessible at z < 12, but appears to break above z = 9. Nine candidates are clustered within a small region of 20'' across, representing a potentially unprecedented concentration. Given the poor statistics, however, we must await similar constraints from the additional HFF clusters to properly examine this trend. The physical properties of our candidates are examined using the range of lens models developed for the HFF program by various groups including our own, for a better estimate of underlying systematics. Our spectral-energy-distribution fits for the brightest objects suggest stellar masses of ≃ 10⁹M_☉, star formation rates of ≃ 4 M_☉ yr⁻¹, and a typical formation redshift of z ≲ 19.

We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in six-dimensional phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams that evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5%–20%, depending on the choice of analytic parameterization. Collectively, the mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1- and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.

I present a new framework for modeling the dynamics of tidal streams. The framework consists of simple models for the initial action-angle distribution of tidal debris, which can be straightforwardly evolved forward in time. Taking advantage of the essentially one-dimensional nature of tidal streams, the transformation to position–velocity coordinates can be linearized and interpolated near a small number of points along the stream, thus allowing for efficient computations of a stream's properties in observable quantities. I illustrate how to calculate the stream's average location (its "track") in different coordinate systems, how to quickly estimate the dispersion around its track, and how to draw mock stream data. As a generative model, this framework allows one to compute the full probability distribution function and marginalize over or condition it on certain phase-space dimensions as well as convolve it with observational uncertainties. This will be instrumental in proper data analysis of stream data. In addition to providing a computationally efficient practical tool for modeling the dynamics of tidal streams, the action-angle nature of the framework helps elucidate how the observed width of the stream relates to the velocity dispersion or mass of the progenitor, and how the progenitors of "orphan" streams could be located. The practical usefulness of the proposed framework crucially depends on the ability to calculate action-angle variables for any orbit in any gravitational potential. A novel method for calculating actions, frequencies, and angles in any static potential using a single orbit integration is described in the Appendix.

The use of a high precision pulsar timing array is a promising approach to detecting gravitational waves in the very low frequency regime (10⁻⁶–10⁻⁹ Hz) that is complementary to ground-based efforts (e.g., LIGO, Virgo) at high frequencies (∼10–10³ Hz) and space-based ones (e.g., LISA) at low frequencies (10⁻⁴–10⁻¹ Hz). One of the target sources for pulsar timing arrays is individual supermassive black hole binaries which are expected to form in galactic mergers. In this paper, a likelihood-based method for detection and parameter estimation is presented for a monochromatic continuous gravitational wave signal emitted by such a source. The so-called pulsar terms in the signal that arise due to the breakdown of the long-wavelength approximation are explicitly taken into account in this method. In addition, the method accounts for equality and inequality constraints involved in the semi-analytical maximization of the likelihood over a subset of the parameters. The remaining parameters are maximized over numerically using Particle Swarm Optimization. Thus, the method presented here solves the monochromatic continuous wave detection and parameter estimation problem without invoking some of the approximations that have been used in earlier studies.

In this paper we report on a two-year study to estimate the Ne/O abundance ratio in the gas phase of the local interstellar cloud (LIC). Based on the first two years of observations with the Interstellar Boundary Explorer, we determined the fluxes of interstellar neutral (ISN) O and Ne atoms at the Earth's orbit in spring 2009 and 2010. A temporal variation of the Ne/O abundance ratio at the Earth's orbit could be expected due to solar cycle-related effects such as changes of ionization. However, this study shows that there is no significant change in the Ne/O ratio at the Earths orbit from 2009 to 2010. We used time-dependent survival probabilities of the ISNs to calculate the Ne/O abundance ratio at the termination shock. Then we estimated the Ne/O abundance ratio in the gas phase of the LIC with the use of filtration factors and the ionization fractions. From our analysis, the Ne/O abundance ratio in the LIC is 0.33 ± 0.07, which is in agreement with the abundance ratio inferred from pickup-ion measurements.

We present high-resolution H i 21 cm observations with the Karl G. Jansky Very Large Array for three H i rich galaxies in absorption against radio quasars. Our sample contains six sightlines with impact parameters from 2.6 to 32.4 kpc. We detected a narrow H i absorber of FWHM 1.1 km s⁻¹ at 444.5 km s⁻¹ toward SDSS J122106.854+454852.16 probing the dwarf galaxy UCG 7408 at an impact parameter of 2.8 kpc. The absorption feature was barely resolved and its width corresponds to a maximum kinetic temperature, T_k ≈ 26 K. We estimate a limiting peak optical depth of 1.37 and a column density of 6 × 10¹⁹ cm⁻². The physical extent of the absorber is 0.04 kpc² and covers ∼25%–30% of the background source. A comparison between the emission and absorption strengths suggests the cold-to-total H i column density in the absorber is ∼30%. Folding in the covering fraction, the cold-to-total H i mass is ∼10%. This suggest that condensation of warm H i (T_s ∼ 1000 K) to cold phase (T_s < 100 K) is suppressed in UGC 7408. The unusually low temperature of the H i absorber also indicates inefficiency in condensation of atomic gas into molecular gas. The suppression in condensation is likely to be the result of low metal content in this galaxy. The same process might explain the low efficiency of star formation in dwarf galaxies despite their huge gas reservoirs. We report the non-detection of H i in absorption in five other sightlines. This indicates that either the cold gas distribution is highly patchy or the gas is much warmer (T_s > 1000 K) toward these sightlines.

In the Galactic halo, metal-bearing Galactic halo material mixes into high velocity clouds (HVCs) as they hydrodynamically interact. This interaction begins long before the clouds completely dissipate and long before they slow to the velocity of the Galactic material. In order to make quantitative estimates of the mixing efficiency and resulting metal enrichment of HVCs, we made detailed two- and three-dimensional simulations of cloud–interstellar medium interactions. Our simulations track the hydrodynamics and time-dependent ionization levels. They assume that the cloud originally has a warm temperature and extremely low metallicity while the surrounding medium has a high temperature, low density, and substantial metallicity, but our simulations can be generalized to other choices of initial metallicities. In our simulations, mixing between cloud and halo gas noticeably raises the metallicity of the high velocity material. We present plots of the mixing efficiency and metal enrichment as a function of time.

The origin of the knee in the cosmic ray spectrum remains to be an unsolved fundamental problem. There are various kinds of models that predict different break positions and the compositions of the knee. In this work, we suggest the use of diffuse γ-rays and neutrinos as probes to test these models. Based on several typical types of composition models, the diffuse γ-ray and neutrino spectra are calculated and show distinctive cutoff behaviors at energies from tens of TeV to multi-PeV. The expected flux will be observable by the newly upgraded Tibet-ASγ+MD (muon detector) experiment as well as more sensitive future projects, such as LHAASO and HiSCORE. By comparing the neutrino spectrum with the recent observations by the IceCube experiment, we find that the diffuse neutrinos from interactions between the cosmic rays and the interstellar medium may not be responsible to the majority of the IceCube events. Future measurements of the neutrinos may be able to identify the Galactic diffuse component and shed further light on the problem of the knee of cosmic rays.