Search

Relativistically broadened and redshifted 6.4 – 6.9 keV iron K lines are observed from many accretion powered objects, including X-ray binaries and active galactic nuclei(AGN). Existence of gas close to the central engine implies large radiation intensities and correspondingly large gas densities if the gas is to remain partially ionized. Simple estimates indicate that high gas densities are needed to allow survival of iron against ionization. These are high enough that rates for many atomic processes are affected by mechanisms related to interactions with nearby ions and electrons. Radiation intensities are high enough that stimulated processes can be important. Most models currently in use for interpreting relativistic lines use atomic rate coefficients designed for use at low densities and neglect stimulated processes. In our work so far we have presented atomic structure calculations with the goal of providing physically appropriate models at densities consistent with line-emitting gas near compact objects. In this paper we apply these rates to photoionization calculations, and produce ionization balance curves and X-ray emissivities and opacities which are appropriate for high densities and high radiation intensities. The final step in our program will be presented in a subsequent paper: Model atmosphere calculations which incorporate these rates into synthetic spectra.

This is the second in a series of papers presenting the results of fully general relativistic simulations of stellar tidal disruptions in which the stars' initial states are realistic main-sequence models. In the first paper, we gave an overview of this program and discussed the principal observational implications of our work. Here we describe our calculational method, which includes a new method for calculating fully relativistic stellar self-gravity, and provide details about the outcomes of full disruptions, focusing on the stellar mass dependence of the outcomes for a black hole of mass 106 M(⊙). We consider eight different stellar masses, from 0.15 M(⊙) to 10 M(⊙). We find that, relative to the traditional order-of-magnitude estimate r t, the physical tidal radius of low-mass stars (M(*) ≲ 0.7 M(⊙)) is larger by tens of percent, while for high-mass stars (M(*) ≲ 1 M(⊙)) it is smaller by a factor of 2–2.5. The traditional estimate of the range of energies found in the debris is ≈1.4× too large for low-mass stars, but is a factor of ~2 too small for high-mass stars; in addition, the energy distribution for high-mass stars has significant wings. For all stars undergoing tidal encounters, we find that mass loss continues for many stellar vibration times because the black hole's tidal gravity competes with the instantaneous stellar gravity at the star's surface until the star has reached a distance from the black hole ~O(10)r(t).

This paper introduces a series of papers presenting a quantitative theory for the tidal disruption of main-sequence stars by supermassive black holes. Using fully general relativistic hydrodynamics simulations and MESA-model initial conditions, we explore the pericenter-dependence of tidal disruption properties for eight stellar masses(MM0.1510)and six black hole masses( MM105 105BH7). We present here the results most relevant to observations. The effects of internal stellar structure and relativity decouple for both the disruption cross section and the characteristic energy width of the debris. Moreover, the full disruption cross section is almost independent of Må for Må/Me3. Independent of Må, relativistic effects increase the critical pericenter distance for full disruption events by up to a factor of∼3 relative to the Newtonian prediction. The probability of a direct capture is also independent of Må; at MBH/Me;5×106 this probability is equal to the probability of a complete disruption. The breadth of the debris energy distribution ΔE can differ from the standard estimate by factors of0.35−2, depending on Må and MBH, implying a corresponding change(∝(ΔE)−3/2)in the characteristic mass-return timescale. We provide analytic forms, suitable for use in both event rate estimates and parameter inference, to describe all these trends. For partial disruptions, we find a nearly universal relation between the star’s angular momentum and the fraction of Må remaining. Within the “empty loss-cone” regime, partial disruptions must precede full disruptions. These partial disruptions can drastically affect the rate and appearance of subsequent total disruptions.

To adequately constrain the frequency of energy deposition in active region cores in the solar corona, systematic comparisons between detailed models and observational data are needed. In this paper, we describe a pipeline for forward modeling active region emission using magnetic field extrapolations and field-aligned hydrodynamic models. We use this pipeline to predict time-dependent emission from active region NOAA 1158 for low-, intermediate-, and high-frequency nanoflares. In each pixel of our predicted multi-wavelength, time-dependent images, we compute two commonly used diagnostics: the emission measure slope and the time lag. We find that signatures of the heating frequency persist in both of these diagnostics. In particular, our results show that the distribution of emission measure slopes narrows and the mean decreases with decreasing heating frequency and that the range of emission measure slopes is consistent with past observational and modeling work. Furthermore, we find that the time lag becomes increasingly spatially coherent with decreasing heating frequency while the distribution of time lags across the whole active region becomes more broad with increasing heating frequency. In a follow-up paper, we train a random forest classifier on these predicted diagnostics and use this model to classify real observations of NOAA 1158 in terms of the underlying heating frequency.

The 3D characterization of magnetic flux ropes observed in the heliosphere has been a challenging task for decades. This is mainly due to the limitations on inferring the 3D global topology and physical properties from the 1D time series from any spacecraft. To advance our understanding of magnetic flux ropes whose configuration departs from the typical stiff geometries, here we present an analytical solution for a 3D flux rope model with an arbitrary cross section and a toroidal global shape. This constitutes the next level of complexity following the elliptic-cylindrical (EC) geometry. The mathematical framework was established by Nieves-Chinchilla et al. with the EC flux rope model, which describes a magnetic topology with an elliptical cross section as a first approach to changes in the cross section. In the distorted-toroidal flux rope model, the cross section is described by a general function. The model is completely described by a nonorthogonal geometry and the Maxwell equations can be consistently solved to obtain the magnetic field and relevant physical quantities. As a proof of concept, this model is generalized in terms of the radial dependence of current density components. The last part of this paper is dedicated to a specific function, F(φ) = δ (1 - λ cosφ), to illustrate possibilities of the model. This model paves the way toward the investigation of complex distortions of magnetic structures in the solar wind. Future investigations will explore these distortions in depth by analyzing specific events; studying implications for physical quantities, such as magnetic fluxes, helicity, or energy; and evaluating the force balance with the ambient solar wind that allows such distortions.

In this paper, we explore the power of the cosmic microwave background (CMB) polarization (E-mode) data to corroborate four potential anomalies in CMB temperature data: the lack of large angular-scale correlations, the alignment of the quadrupole and octupole (Q–O), the point-parity asymmetry, and the hemispherical power asymmetry. We use CMB simulations with noise representative of three experiments—the Planck satellite, the Cosmology Large Angular Scale Surveyor (CLASS), and the LiteBIRD satellite—to test how current and future data constrain the anomalies. We find the correlation coefficients ρ between temperature and E-mode estimators to be less than 0.1, except for the point-parity asymmetry (ρ = 0.17 for cosmic-variance-limited simulations), confirming that E-modes provide a check on the anomalies that is largely independent of temperature data. Compared to Planck component-separated CMB data (smica), the putative LiteBIRD survey would reduce errors on E-mode anomaly estimators by factors of ∼3 for hemispherical power asymmetry and point-parity asymmetry, and by ∼26 for lack of large-scale correlation. The improvement in Q–O alignment is not obvious due to large cosmic variance, but we found the ability to pin down the estimator value will be improved by a factor ≳100. Improvements with CLASS are intermediate to these.

We present and analyze the optical and X-ray catalogs of moderate-redshift cluster candidates from the ROSA TOptical X-Ray Survey, or ROXS. The survey covers the sky area contained in the fields of view of 23 deep archival ROSA T PSPC pointings, 4.8 square degrees. The cross-correlated cluster catalogs were con- structed by comparing two independent catalogs extracted from the optical and X-ray bandpasses, using a matched-filter technique for the optical data and a wavelet technique for the X-ray data. We cross-identified cluster candidates in each catalog. As reported in Paper 1, the matched-filter technique found optical counter- parts for at least 60% (26 out of 43) of the X-ray cluster candidates; the estimated redshifts from the matched filter algorithm agree with at least 7 of 1 1 spectroscopic confirmations (Az 5 0.10). The matched filter technique. with an imaging sensitivity of ml N 23, identified approximately 3 times the number of candidates (155 candidates, 142 with a detection confidence >3 u) found in the X-ray survey of nearly the same area. There are 57 X-ray candidates, 43 of which are unobscured by scattered light or bright stars in the optical images. Twenty-six of these have fairly secure optical counterparts. We find that the matched filter algorithm, when applied to images with galaxy flux sensitivities of mI N 23, is fairly well-matched to discovering z 5 1 clusters detected by wavelets in ROSAT PSPC exposures of 8000-60,000 s. The difference in the spurious fractions between the optical and X-ray (30%) and IO%, respectively) cannot account for the difference in source number. In Paper I, we compared the optical and X-ray cluster luminosity functions and we found that the luminosity functions are consistent if the relationship between X-ray and optical luminosities is steep (Lx o( L&f). Here, in Paper 11, we present the cluster catalogs and a numerical simulation of the ROXS. We also present color-magnitude plots for several of the cluster candidates, and examine the prominence of the red sequence in each. We find that the X-ray clusters in our survey do not all have a prominent red sequence. We conclude that while the red sequence may be a distinct feature in the color-magnitude plots for virialized massive clusters, it may be less distinct in lower mass clusters of galaxies at even moderate redshifts. Multiple, complementary methods of selecting and defining clusters may be essential, particularly at high redshift where all methods start to run into completeness limits, incomplete understanding of physical evolution, and projection effects.

The Galactic low-mass X-ray binary AT2019wey (ATLAS19bcxp, SRGAJ043520.9+552226, SRGEJ043523.3+552234, ZTF19acwrvzk) was discovered as a new optical transient in Dec 2019, and independently as an X-ray transient in Mar 2020. In this paper, we present comprehensive NICER, NuSTAR, Chandra, Swift, and MAXI observations of AT2019wey from ~1 year prior to the discovery to the end of September 2020. AT2019wey appeared as a ~ 1 mCrab source and stayed at this flux density for several months, displaying a hard X-ray spectrum that can be modeled as a power-law with photon index Γ ~ 1.8. In June 2020 it started to brighten and reached ~20 mCrab in ~2 months. The inclination of this system can be constrained to I ≲ 30° by modelling the reflection spectrum. Starting from late-August (~ 59082 MJD), AT2019wey entered into the hard-intermediate state (HIMS), and underwent a few week-long timescale outbursts, where the brightening in soft X-rays is correlated with the enhancement of a thermal component. Low-frequency quasi-periodic oscillation (QPO) was observed in the HIMS. We detect no pulsation and in timing analysis of the NICER and NuSTAR data. The X-ray states and power spectra of AT2019wey are discussed against the landscape of low-mass X-ray binaries.

We present timing analysis for a Rossi X-ray Timing Explorer (RXTE) observation of Cygnus X-1 in its hard/low state. This was the first RXTE observation of Cyg X-1 taken after it transited back to this state from its soft/high state. RXTE's large effective area, superior timing capabilities, and ability to obtain long, uninterrupted observations have allowed us to obtain measurements of the power spectral density (PSD), coherence function, and Fourier time lags to a decade lower in frequency and half a decade higher in frequency than typically was achieved with previous instruments. Notable aspects of our observations include a weak 0.005 Hz feature in the PSD coincident with a coherence recovery; a 'hardening' of the high-frequency PSD with increasing energy; a broad frequency range measurement of the coherence function, revealing rollovers from unity coherence at both low and high frequency; and an accurate determination of the Fourier time lags over two and a half decades in frequency. As has been noted in previous similar observations, the time delay is approximately proportional to f(exp -0.7), and at a fixed Fourier frequency the time delay of the hard X-rays compared to the softest energy channel tends to increase logarithmically with energy. Curiously, the 0.01-0.2 Hz coherence between the highest and lowest energy bands is actually slightly greater than the coherence between the second highest and lowest energy bands. We carefully describe all of the analysis techniques used in this paper, and we make comparisons of the data to general theoretical expectations. In a companion paper, we make specific comparisons to a Compton corona model that we have successfully used to describe the energy spectral data from this observation.

The luminosity structure of spiral galaxies is studied using the technique of principal component analysis. It is found that approximately 94% of the variation in the luminosity distribution of galaxies can be accounted for by just two principal components. The principal luminosity components may contain valuable information about star formation history or whatever luminosity-regulating process occurs in galaxies. Practically, these principal components provide a new approach for the investigation of the luminosity structures of galaxies and their dependence on other properties. They also serve as an excellent objective classification system for galaxies. We introduce in this paper such a classification scheme and explore its various properties. The new system shows a number of very impressive characteristics. Most important, it can well segregate virtually all the important galactic properties we tested and does so much better than the conventional morphological classification systems. Of particular interest is that some distance-dependent parameters can also be determined to a surprisingly good accuracy; for example, absolute magnitude may be determined to an accuracy of approximately 0.6 mag (yet further improvement is believed to be highly possible). Second, the system is objective, and the classification procedure can be automated to a large degree; also the new system can apply to much smaller and fainter images than do eye-based clasification systems. These properties make the new system suitable for practical application, especially on very large (and deeper) digital image catalogs. Third, the classification is expressed in dimensionless numbers, yet the simple notation bears significant and easily understandable meaning, making it easy and convenient to use. Finally, the new system has another extremely useful feature: it provides a very powerful and convenient platform not only for classification, but also for easily recording, examining, and studying the variations and correlations of galaxy properties-all these may be carried out graphically by using the C-vectors and the C-diagrams introduced in the paper. We wil also give an example to demonstrate the use of the classification system for the study of the internal extinction problem in spiral galaxies.

Two-dimensional, time-dependent drift models have done exceptionally well in explaining major modulation features, especially during the A less than 0 magnetic polarity cycle of the heliospheric magnetic field when positively charge particles are drifting in along the heliospheric neutral sheet (HNS). These models were found to do well when the heliospheric 'tilt angles' alpha less than approximately 30 deg (le Roux & Potgeiter). However, they seem to do less well when alpha greater than approximately 30 deg during A less than 0 cycles seem to fail when this happens in A greater than 0 cycles. Progress was made in understanding these phases of the modulation cycle when merged interaction regions (MIRs) were incorporated in time-dependent drift models (Potgieter et al.). It was also explicitly shown that in obtaining large step decreases in cosmic rays, the MIRs had to be global, i.e., having a latitudinal extent of more than approximately 60 deg. Other classes of MIRs, such as local MIRs and co-rotating MIRs were found to be of secondary importance for establishing long-term modulation. In a previous paper we studied the effects of two consecutive, identical global MIRs, together with a changing wavy HNS, on long-term modulation (Potegieter & le Roux). This approach gave a very natural and convincing explanation for the observed step decreases in cosmic-ray modulation. Emphasis was placed in the declining and recovery phases of the 11 yr modulation cycle. In this paper, four consecutive, nonidentical global MIRs, in combination with a varying wavy HNS, were included in our time-dependent drift model in order to do simulations closer to what was observed between 1977 and 1987. By doing this we could model, for the first time, complete 11 and 22 yr cycles in the heliospheric modulation of galactic cosmic-rays, including the solar magnetic polarity reversals.

This paper presents the gravitational-wave measurement of the Hubble constant (H(0)) using the detections from the first and second observing runs of the Advanced LIGO and Virgo detector network. The presence of the transient electromagnetic counterpart of the binary neutron star GW170817 led to the first standard-siren measurement of H(0). Here we additionally use binary black hole detections in conjunction with galaxy catalogs and report a joint measurement. Our updated measurement is H(0) = 69(+16,-8) km/s Mpc (68.3% of the highest density posterior interval with a flat-in-log prior) which is an improvement by a factor of 1.04 (about 4%) over the GW170817-only value of 6(+17,-8) km/s Mpc. A significant additional contribution currently comes from GW170814, a loud and well-localized detection from a part of the sky thoroughly covered by the Dark Energy Survey. With numerous detections anticipated over the upcoming years, an exhaustive understanding of other systematic effects are also going to become increasingly important. These results establish the path to cosmology using gravitational-wave observations with and without transient electromagnetic counterparts.

The Nuclear Spectroscopic Telescope Array (NuSTAR) mission is the first focusing X-ray telescope in the hard X-ray (3–79 keV) band. Among the phenomena that can be studied in this energy band, some require high time resolution and stability: rotation-powered and accreting millisecond pulsars, fast variability from black holes and neutron stars, X-ray bursts, and more. Moreover, a good alignment of the timestamps of X-ray photons to UTC is key for multi-instrument studies of fast astrophysical processes. In this paper, we describe the timing calibration of the NuSTAR mission. In particular, we present a method to correct the temperature-dependent frequency response of the on-board temperature-compensated crystal oscillator. Together with measurements of the spacecraft clock offsets obtained during downlinks passes, this allows a precise characterization of the behavior of the oscillator. The calibrated NuSTAR event timestamps for a typical observation are shown to be accurate to a precision of ∼65 μs.

Solar coronal plumes long seemed to possess a simple geometry supporting spatially coherent, stable outflow without significant fine structure. Recent high-resolution observations have challenged this picture by revealing numerous transient, small-scale, collimated out-flows (“jetlets”) at the base of plumes. The dynamic filamentary structure of solar plumes above these outflows, and its relationship with the overall plume structure, have remained largely unexplored. We analyzed the statistics of continuously observed fine structure in-side a single representative bright plume within a mid-latitude coronal hole during 2016 July2-3. By applying advanced edge-enhancement and spatiotemporal analysis techniques to ex-tended series of high-resolution images from the Solar Dynamics Observatory’s Atmospheric Imaging Assembly, we determined that the plume was composed of numerous time-evolving filamentary substructures, referred to as “plumelets” in this paper, that accounted for most of the plume emission. The number of simultaneously identifiable plumelets was positively correlated with plume brightness, peaked in the fully formed plume, and remained saturated thereafter. The plumelets had transverse widths of 10 Mm and intermittently supported upwardly propagating periodic disturbances with phase speeds of 190-260 km s−1and longitudinal wavelengths of 55-65 Mm. The characteristic frequency (3.5 mHz) is commensurate with that of solar p-modes. Oscillations in neighboring plumelets are uncorrelated, indicating that the waves could be driven by p-mode flows at spatial scales smaller than the plumelet separation. Multiple independent sources of outflow within a single coronal plume should impart significant fine structure to the solar wind that may be detectable by Parker Solar Probe and Solar Orbiter.

High-cadence transient surveys are able to capture supernovae closer to their first light than ever before. Applying analytical models to such early emission, we can constrain the progenitor stars' properties. In this paper, we present observations of SN 2018fif (ZTF 18abokyfk). The supernova was discovered close to first light and monitored by the Zwicky Transient Facility (ZTF) and the Neil Gehrels Swift Observatory. Early spectroscopic observations suggest that the progenitor of SN 2018fif was surrounded by relatively small amounts of circumstellar material compared to all previous cases. This particularity, coupled with the high-cadence multiple-band coverage, makes it a good candidate to investigate using shock-cooling models. We employ the SOPRANOS code, an implementation of the model by Sapir & Waxman and its extension to early times by Morag et al. Compared with previous implementations, SOPRANOS has the advantage of including a careful account of the limited temporal validity domain of the shock-cooling model as well as allowing usage of the entirety of the early UV data. We find that the progenitor of SN 2018fif was a large red supergiant with a radius of R = 744.0 (sup +183.0) (sub -128.0) R⨀ and an ejected mass of M(sub ej) = 9.3 (sup +0.4)(sub -5.8) M⨀. Our model also gives information on the explosion epoch, the progenitor's inner structure, the shock velocity, and the extinction. The distribution of radii is double-peaked, with smaller radii corresponding to lower values of the extinction, earlier recombination times, and a better match to the early UV data. If these correlations persist in future objects, denser spectroscopic monitoring constraining the time of recombination, as well as accurate UV observations (e.g., with ULTRASAT), will help break the extinction/radius degeneracy and independently determine both.

When galaxies merge, the supermassive black holes in their centers may form binaries and emit low-frequency gravitational radiation in the process. In this paper, we consider the galaxy 3C 66B, which was used as the target of the first multimessenger search for gravitational waves. Due to the observed periodicities present in the photometric and astrometric data of the source, it has been theorized to contain a supermassive black hole binary. Its apparent 1.05-year orbital period would place the gravitational-wave emission directly in the pulsar timing band. Since the first pulsar timing array study of 3C 66B, revised models of the source have been published, and timing array sensitivities and techniques have improved dramatically. With these advances, we further constrain the chirp mass of the potential supermassive black hole binary in 3C 66B to less than (1.65±0.02)×109Me using data from the NANOGrav 11-year data set. This upper limit provides a factor of 1.6 improvement over previous limits and a factor of 4.3 over the first search done. Nevertheless, the most recent orbital model for the source is still consistent with our limit from pulsar timing array data. In addition, we are able to quantify the improvement made by the inclusion of source properties gleaned from electromagnetic data over “blind” pulsar timing array searches. With these methods, it is apparent that it is not necessary to obtain exact a priori knowledge of the period of a binary to gain meaningful astrophysical inferences.

We present model calculations of the transport processes of solar energetic particles in the corona and interplanetary medium for two events detected by Parker Solar Probe near its second perihelion on 2019 April 2 and April 4. In the 2019 April 2 event, the <100 keV proton differential intensity measured by the ISOIS EPI-Lo instrument increased by more than a factor of 10 above the pre-event intensity, whereas the ~1 MeV proton differential intensity detected by the EPI-Hi instrument did not show any intensity enhancement. In the 2019 April 4 event, the ~1 MeV proton intensity showed an increase of more than a factor of 100 above the pre-event intensity, but the <100 keV proton intensity enhancement was rather low and gradual. The observational fact that only the 2019 April 4 event displayed a high-energy intensity enhancement indicates that the associated acceleration process in the second event was more energetic than the first event. However, the gradual and low enhancement of the low-energy proton intensity in the second event is at odds with the scenario. In this paper, we investigate the injection and transport processes of protons in the corona and interplanetary magnetic fields with numerical model simulations. Our model calculations suggest that the gradual and low rise of the low-energy protons in the April 4 event was probably due to the different diffusion conditions in these two events and the fact that the spacecraft was moving away from the source region during these two events.

The Cosmology Large Angular Scale Surveyor (CLASS) is a telescope array that observes the cosmic microwave background (CMB) over 75% of the sky from the Atacama Desert, Chile, at frequency bands centered near 40, 90, 150, and 220 GHz. CLASS measures the large angular scale (1°~< θ >=90°) CMB polarization to constrain the tensor-to-scalar ratio at the r ∼ 0.01 level and the optical depth to last scattering to the sample variance limit. This paper presents the optical characterization of the 40 GHz telescope during its first observation era, from 2016 September to 2018 February. High signal-to-noise observations of the Moon establish the pointing and beam calibration. The telescope boresight pointing variation is <0°. 023 (<1.6% of the beam’s full width at half maximum (FWHM)). We estimate beam parameters per detector and in aggregate, as in the CMB survey maps. The aggregate beam has an FWHM of 1°.579 ± 0°.001 and a solid angle of 838 ± 6 μsr, consistent with physical optics simulations. The corresponding beam window function has a sub-percent error per multipole at ℓ < 200. An extended 90° beam map reveals no significant far sidelobes. The observed Moon polarization shows that the instrument polarization angles are consistent with the optical model and that the temperature-to-polarization leakage fraction is <10(exp −4) (95% C.L.). We find that the Moon-based results are consistent with measurements of M42, RCW 38, and Tau A from CLASS’s CMB survey data. In particular, Tau A measurements establish degree- level precision for instrument polarization angles.

An ensemble of inspiraling supermassive black hole binaries should produce a stochastic background of very low frequency gravitational waves. This stochastic background is predicted to be a power law, with a gravitational-wave strain spectral index of −2/3, and it should be detectable by a network of precisely timed millisecond pulsars, widely distributed on the sky. This paper reports a new “time slicing” analysis of the 11 yr data release from the North American Nanohertz Observatory for Gravitational Waves (NANOGrav) using 34 millisecond pulsars. Methods to flag potential “false-positive” signatures are developed, including techniques to identify responsible pulsars. Mitigation strategies are then presented. We demonstrate how an incorrect noise model can lead to spurious signals, and we show how independently modeling noise across 30 Fourier components, spanning NANOGrav’s frequency range, effectively diagnoses and absorbs the excess power in gravitational-wave searches. This results in a nominal, and expected, progression of our gravitational-wave statistics. Additionally, we show that the first interstellar medium event in PSR J1713+0747 pollutes the common red-noise process with low spectral index noise, and we use a tailored noise model to remove these effects.

Early-time observations of Type Ia supernovae (SNe Ia) are essential to constrain the properties of their progenitors. In this paper, we present high-quality light curves of 127 SNe Ia discovered by the Zwicky Transient Facility (ZTF) in2018. We describe our method to perform forced point-spread function photometry, which can be applied to other types of extragalactic transients. With a planned cadence of six observations per night (three g +three r), all of the 127 SNe Ia are detected in both g and r bands more than 10 days (in the rest frame) prior to the epoch of g-band maximum light. The redshifts of these objects range from z = 0.0181 to 0.165; the median redshift is 0.074. Among the 127 SNe, 50 are detected at least 14 days prior to maximum light (in the rest frame), with a subset of nine objects being detected more than 17 days before g-band peak. This is the largest sample of young SNe Ia collected to date; it can be used to study the shape and color evolution of the rising light curves in unprecedented detail. We discuss six peculiar events in this sample: one 02cx-like event ZTF18abclfee (SN 2018crl), one Ia-CSM SN ZTF18aaykjei (SN 2018cxk), and four objects with possible super-Chandrasekhar mass progenitors: ZTF18abhpgje (SN 2018eul), ZTF18abdpvnd (SN 2018dvf), ZTF18aawpcel (SN 2018cir), and ZTF18abddmrf (SN 2018dsx).

A comprehensive understanding of the structure of Doppler motions in the transition region including the center-to-limb variation and its relationship with the magnetic field structure is vital for the understanding of mass and energy transfer in the solar atmosphere. In this paper, we have performed such a study in an active region using the Si IV 1394 Å emission line recorded by the Interface Region Imaging Spectrograph and the line-of-sight photospheric magnetic field obtained by the Helioseismic and Magnetic Imager onboard the Solar Dynamics Observatory. The active region has two opposite polarity strong field regions separated by a weak field corridor, which widened as the active region evolved. On average, the strong field regions (corridor) show(s) redshifts of 5–10 (3–9) km s(exp -1) (depending on the date of observation). There is, however, a narrow lane in the middle of the corridor with near-zero Doppler shifts at all disk positions, suggesting that any flows there are very slow. The Doppler velocity distributions in the corridor seem to have two components—a low velocity component centered near 0 km s(exp -1) and a high-velocity component centered near 10 km s(exp -1). The high-velocity component is similar to the velocity distributions in the strong field regions, which have just one component. Both exhibit a small center-to-limb variation and seem to come from the same population of flows. To explain these results, we suggest that the emission from the lower transition region comes primarily from warm type II spicules, and we introduce the idea of a “chromospheric wall”—associated with classical cold spicules—to account for a diminished center-to-limb variation.

Warm absorbers and related phenomena are some of the observable manifestations of outflows or winds from active galactic nuclei (AGNs). Warm absorbers are common in low-luminosity AGNs. They have been extensively studied observationally and are well described by simple phenomenological models. However, major open questions remain. What is the driving mechanism? What is the density and geometrical distribution? How much associated fully ionized gas is there? What is the relation to the quasi-relativistic "ultrafast outflows"? In this paper we present synthetic spectra for the observable properties of warm absorber flows and associated quantities. We use ab initio dynamical models, i.e., solutions of the equations of motion for gas in finite difference form. The models employ various plausible assumptions for the origin of the warm absorber gas and the physical mechanisms affecting its motion. The synthetic spectra are presented as an observational test of these models. In this way we explore various scenarios for warm absorber dynamics. We show that observed spectra place certain requirements on the geometrical distribution of the warm absorber gas, and that not all dynamical scenarios are equally successful at producing spectra similar to what is observed.

The zodiacal dust complex, a population of dust and small particles that pervades the solar system, provides important insight into the formation and dynamics of planets, comets, asteroids, and other bodies. We present a new set of data obtained from direct measurements of momentum transfer to a spacecraft from individual particle impacts. This technique is made possible by the extreme precision of the instruments flown on the LISA Pathfinder spacecraft, a technology demonstrator for a future space-based gravitational wave observatory. Pathfinder employed a technique known as drag-free control that achieved rejection of external disturbances, including particle impacts, using a micropropulsion system. Using a simple model of the impacts and knowledge of the control system, we show that it is possible to detect impacts and measure properties such as the transferred momentum, direction of travel, and location of impact on the spacecraft. In this paper, we present the results of a systematic search for impacts during 4348 hr of Pathfinder data. We report a total of 54 candidates with transferred momenta ranging from 0.2 to 230 μNs. We furthermore make a comparison of these candidates with models of micrometeoroid populations in the inner solar system, including those resulting from Jupiter-family comets (JFCs), Oort Cloud comets, Halley-type comets, and asteroids. We find that our measured population is consistent with a population dominated by JFCs, with some evidence for a smaller contribution from Halley-type comets, in agreement with consensus models of the zodiacal dust complex in the momentum range sampled by LISA Pathfinder.

Ultraluminous X-ray (ULX) pulsars are a new class of object powered by apparent super-critical accretion onto magnetized neutron stars. Three sources in this class have been identified so far; M82X-2, NGC 5907 ULX-1, and NGC 7793 P13 have been found to have two properties in common; ∼1 s spin periods, and for NGC 5907 ULX-1 and NGC 7793 P13 periodic X-ray flux modulations on timescales of ∼60–80 days. M82 X-2 resides in a crowded field that includes the ULX M82 X-1 separated from X-2 by 5″, as well as other bright point sources. A 60 day modulation has been observed from the region, but the origin has been difficult to identify; both M82 X-1 and X-2 have been suggested as the source. In this paper we present the analysis of a systematic monitoring campaign by Chandra, the only X-ray telescope capable of resolving the crowded field. From a simple Lomb–Scargle periodogram analysis and a more sophisticated Gaussian Process analysis we find that only X-2 exhibits a periodic signal around 60 days, supporting previous claims that it is the origin. We also construct a phase-averaged flux profile of the modulations from higher-cadence Swift/XRT data and find that the flux variations in the Chandra data are fully consistent with the flux profile. Since the orbit of the neutron star and its companion is known to be 2.5 days, the ∼60 day period must be super-orbital in origin. The flux of the modulations varies by a factor of ∼100 from the minimum to the maximum, with no evidence for spectral variations, making the origin difficult to explain.

We study electron–positron pair production in polar caps of energetic pulsars to determine the maximum multiplicity of pair plasma a pulsar can produce under the most favorable conditions. This paper complements and updates our study of pair cascades presented in Timokhin & Harding (2015) with a more accurate treatment of the effects of ultrastrong B approximately greater than 3 x 10 (exp 12) G magnetic fields and emission processes of primary and secondary particles. We include pairs produced by curvature and synchrotron radiation photons as well as resonant Compton-scattered photons. We develop a semianalytical model of electron–positron cascades that can efficiently simulate pair cascades with an arbitrary number of microphysical processes and use it to explore cascade properties for a wide range of pulsar parameters. We argue that the maximum cascade multiplicity cannot exceed approximately a few x 10 (exp 5) and that the multiplicity has a rather weak dependence on pulsar period. The highest multiplicity is achieved in pulsars with magnetic field 4 x 10 (exp 12) is approximately greater than B is approximately greater than 10 (exp 13) G and hot surfaces, with T is approximately greater than 10 (exp 6)K. We also derive analytical expressions for several physical quantities relevant for electromagnetic cascade in pulsars, which may be useful in future works on pulsar cascades, including the upper limit on cascade multiplicity and various approximations for the parameter χ, the exponential factor in the expression for photon attenuation in strong magnetic fields.

As part of our work on nitrogen-rich ices, the IR spectra and band strengths used in a recent paper to identify and quantify radiation-induced changes in an N2+H2O ice near 15 K are examined, along with reports of (i) a chemical tracer for N2+H2O ices, (ii) a new IR feature of solid N2, and (iii) a striking 15N isotopic enrichment. Problems are found for each IR band strength used and for each of the three claims made, to the extent that none are supported by the results presented to date. In contrast, new work presented here, combined with several older investigations, strongly supports the formation of di- and triatomic nitrogen oxides in irradiated N2-rich ices. Observations and trends in the chemistry of N2-rich icy solids are described, and conclusions are drawn. A considerable amount of material from previous chemical studies of N2-rich systems, spanning more than a century, is brought together for the rst time and used to examine the chemistry of N2-rich ices in extraterrestrial environments. Needs are identied and suggestions made for future studies of N2-rich interstellar and planetary ice analogs.

This paper reports on the re-analysis of solar flares in which the hard X-rays (HXRs) come predominantly from the corona rather than from the more usual chromospheric footpoints. All of the 26 previously analyzed event time intervals, over 13 flares, are re-examined for consistency with a flare model in which electrons are accelerated near the top of a magnetic loop which has a sufficiently high density to stop most of the electrons by Coulomb collisions before they can reach the footpoints. Of particular importance in the previous analysis was the finding that the length of the coronal HXR source increased with energy in the 20–30 keV range. However, after allowing for the possibility that footpoint emission at the higher energies affects the inferred length of the coronal HXR source, and using analysis techniques that suppress the possible influence of such footpoint emission, we conclude that there is no longer evidence that the length of the HXR coronal sources increase with increasing energy. In fact, for the six flares and 12 time intervals that satisfied our selection criteria, the loop lengths decreased on average by 1.0 ± 0.2 arcsec between 20 and 30 keV, with a standard deviation of 3.5 arcsec. We find strong evidence that the peak of the coronal HXR source increases in altitude with increasing energy. For the thermal component of the emission, this is consistent with the standard CHSKP flare model in which magnetic reconnection in a coronal current sheet results in new hot loops being formed at progressively higher altitudes. The explanation for the nonthermal emission is not so clear.

A recent laboratory experiment of ideal magnetohydrodynamic instabilities revealed four distinct eruption regimes readily distinguished by the torus instability (TI) and helical kink instability (KI) parameters. To establish its observational counterpart, we collected 38 solar flares (stronger than GOES-class M5 in general) that took place within 45° of disk center during 2011–2017, 26 of which are associated with a halo or partial halo coronal mass ejection (CME; i.e., ejective events), while the others are CME-less (i.e., confined events). This is a complete sample of solar events satisfying our selection criteria detailed in the paper. For each event, we calculate a decay index n of the potential strapping field above the magnetic flux rope (MFR) in and around the flaring magnetic polarity inversion line (a TI parameter) and the unsigned twist number T(sub w) of the nonlinear force-free field lines forming the same MFR (a KI parameter). We then construct an n–T(sub w) diagram to investigate how the eruptiveness depends on these parameters. We find that (1) T(sub w) appears to play little role in discriminating between confined and ejective events; (2) the events with n ≳ 0.8 are all ejective, and all confined events have n ≲ 0.8. However, n ≳ 0.8 is not a necessary condition for eruption because some events with n ≲ 0.8 also erupted. In addition, we investigate the MFR’s geometrical parameters, apex height, and distance between footpoints, as a possible factors for the eruptiveness. We briefly discuss the difference of the present result for solar eruptions with that of the laboratory result in terms of the role played by magnetic reconnection.

High-time-resolution in situ wave observations show that Langmuir waves associated with solar type III radio bursts often occur as coherent localized one-dimensional magnetic-field-aligned wave packets with short durations of a few milliseconds and peak intensities well above the strong turbulence thresholds. In this paper, we report observations of a wave packet obtained by the time domain sampler of the STEREO WAVES experiment, which is unique in the sense that it is the most intense wave packet ever detected in association with a solar type III radio burst, with a peak intensity Et ∼ 107 mVm(sup-1). We show that this wave packet provides evidence for (1) oscillating two-stream instability (OTSI), (2) a collapsing soliton formed as a result of OTSI, (3) the formation of a soliton-caviton pair, and (4) excitation of second and third harmonic electromagnetic waves. We also show that the peak intensity and spatial width satisfy the threshold condition for this wave packet to be the collapsing Langmuir wave packet formed as a result of nucleation processes even when δn(sub b) > δn(sub p), where δn(sub b) and δn(sub p) are the levels of background and ponderomotive-force-induced density fluctuations, respectively. Thus, these observations provide unambiguous evidence for the spatial collapse of Langmuir waves in the source region of a type III radio burst, and the observed spectral evidence for OTSI and the ponderomotive-force-induced density cavity strongly suggest that OTSI is mostly likely responsible for the collapse of the observed wave packet.

With 7 years of Interstellar Boundary Explorer (IBEX) measurements of energetic neutral atoms (ENAs), IBEX has shown a clear correlation between dynamic changes in the solar wind and the heliosphere’s response in the formation of ENAs. In this paper, we investigate temporal variations in the latitudinal-dependent ENA spectrum from IBEX and their relationship to the solar wind speed observed at 1 au. We find that the variation in latitude of the transition in ENA spectral indices between low (≲1.8) and high (≳1.8) values, as well as the distribution of ENA spectral indices at high and low latitudes, correlates well with the evolution of the fast and slow solar wind latitudinal structure observed near 1 au. This correlation includes a delay due to the time it takes the solar wind to propagate to the termination shock and into the inner heliosheath, and for ENAs to be generated via charge-exchange and travel back toward 1 au. Moreover, we observe a temporal asymmetry in the steepening of the ENA spectrum in the northern and southern hemispheres, consistent with asymmetries observed in the solar wind and polar coronal holes. While this asymmetry is observed near the upwind direction of the heliosphere, it is not yet observed in the tail direction, suggesting a longer line-of-sight integration distance or different processing of the solar wind plasma downstream of the termination shock.

Condensible substances are nearly ubiquitous in planetary atmospheres. For the most familiar case—water vapor in Earth’s present climate—the condensible gas is dilute, in the sense that its concentration is everywhere small relative to the noncondensible background gases. A wide variety of important planetary climate problems involve nondilute condensible substances. These include planets near or undergoing a water vapor runaway and planets near the outer edge of the conventional habitable zone, for which CO2 is the condensible. Standard representations of convection in climate models rely on several approximations appropriate only to the dilute limit, while nondilute convection differs in fundamental ways from dilute convection. In this paper, a simple parameterization of convection valid in the nondilute as well as dilute limits is derived and used to discuss the basic character of nondilute convection. The energy conservation properties of the scheme are discussed in detail and are verified in radiative-convective simulations. As a further illustration of the behavior of the scheme, results for a runaway greenhouse atmosphere for both steady instellation and seasonally varying instellation corresponding to a highly eccentric orbit are presented. The latter case illustrates that the high thermal inertia associated with latent heat in nondilute atmospheres can damp out the effects of even extreme seasonal forcing.

This is the first of a series of papers dealing with transfer of radiation in dusty nebulae. The formulae for the general problem are presented. This paper discusses nebulae containing dust and hydrogen. Approximate solutions, accurate to within 50 percent, are obtained for spherical nebulae containing purely absorbing dust. Less accurate solutions are presented for the general case of scattering and absorbing dust. The results are presented in a form convenient for calculating the required number of ionizing photons, the fractional ionization of hydrogen, and the fractions of ionizing photons absorbed by dust and gas in terms of size, hydrogen density, dust albedo, and optical depth. Explicit relations between these quantities are obtained for uniform nebulae. Simple integrals are presented for calculating the above quantities for spherically symmetric nebulae with nonuniform gas and dust distribution.

A set of restrictions upon the magnetic field and visible angular size of a canonical nonthermal source (one whose emission is incoherent electron synchrotron) has been obtained in another paper. The degree of electron anisotropy and magnetic field disordering were left as free parameters. In the present paper these restrictions are applied to the variable compact radio sources 3C 84, 3C 120, 3C 273, 3C 279, 3C 454.3, CTA 102, 3C 446, PKS 2134 + 004, VRO 42.22.01, and OJ 287. The theoretically derived angular sizes for these sources are in reasonable agreement with published VLBI sizes, provided that the characteristic pitch angles are not extremely small. The magnetic field strengths determined from the degree of circular polarization are generally compatible with field strengths derived from the self-absorption and synchrotron self-Compton conditions, again provided that the pitch angles are not too small. It is argued that the self-absorption optical depth cannot be much less than unity near the low-frequency turnovers; otherwise, serious discrepancies with observation result.

In a previous paper, the PPN (parametrized post-Newtonian) formalism was used to analyze relativistic influences on stellar stability in nearly all metric theories of gravity. That analysis omitted all preferred-frame terms. In this paper, possible preferred-frame effects on stellar stability are examined and no new instabilities are found. Although terms linear in the preferred-frame velocity w (time-odd terms, analogous to viscosity and energy generation) change the shapes of the normal modes, their symmetry properties prevent them from changing the characteristic frequencies. Thus, no new vibrational or secular instabilities can occur. Terms quadratic in w do not change either the shapes of the normal modes or the characteristic frequencies for radial pulsations (except for the effects due to the renormalization of the gravitation constant which does not affect stability). Thus, they have no influence on radial stability. Terms quadratic in w do change both the normal modes and the characteristic frequencies of nonradial pulsations; but in the limit of a neutral mode these changes vanish.

It is shown that the observed data about the apsidal motion of gaseous rings around some Be stars, like Beta(1) Mon and Pi Aqr, cannot be explained solely by the influence of the equatorial bulge. This fact has prompted the development of a theory of apsidal motion and slow drift of orbital elements that is based on the mixing of slowly escaping gases from the central star with gases in the emission ring around it. Thus, the basic physical idea in this paper follows closely what has been advanced in the first paper of this series but the mathematical treatment has been broadened in order to deal with the actual orbits, instead of osculating orbits, of particles in the ring. In accordance with this theory, it is believed that in cases where no mixing of gases takes place, the apsidal motion is controlled solely by the equatorial bulge. But in cases where mixing takes place, the apsidal motion depends upon the manner of gas mixing as well.

Dicke has recently contested our statement in an earlier paper that faculae could account for a large part, if not all, of the solar oblateness signal measured by Dicke and Goldenberg in 1966. Using the facular oblateness signal published in our earlier paper and some hitherto unpublished data from his 1966 observations, he concludes that faculae account for only a small part (11%) of the observed excess oblateness. His analysis considers data only from a restricted 48-day sample and is based on the assumption that only the observed oblateness signal is subject to error. Our analysis considers data from all 64 days on which observations were made, and is based on the assumption that both the observed oblateness signal and the facular signal are subject to error. We find that faculae account for at least one-third to one-half of the observed excess oblateness, depending on whether 48 days or 64 days are used in the analysis.

This paper argues that the bridges and tails seen in some multiple galaxies are just tidal relics of close encounters. These consequences of the brief but violent tidal forces are here studied in a deliberately simple-minded fashion. Each encounter is considered to involve only two galaxies and to be roughly parabolic; each galaxy is idealized as just a disk of noninteracting test particles which initially orbit a central mass point. As shown here, the two-sided distortions provoked by gravity alone in such circumstances can indeed evolve kinematically into some remarkably narrow and elongated features. Besides extensive pictorial survey of tidal damage, this paper offers reconstructions of the orbits and outer shapes of four specific interacting pairs: Arp 295, M51 + NGC 5195, NGC 4676, and NGC 4038/9.

We present results from the first three years of observations of moving group (MG) targets in the Strategic Exploration of Exoplanets and Disks with Subaru (SEEDS) high-contrast imaging survey of exoplanets and disks using the Subaru telescope. We achieve typical contrasts of (is) approximately10(exp 5) at 1" and (is) approximately 10(exp 6) beyond 2" around 63 proposed members of nearby kinematic MGs. We review each of the kinematic associations to which our targets belong, concluding that five, beta Pictoris ((is) approximately 20 Myr), AB Doradus ((is) approximately 100 Myr), Columba ((is) approximately 30 Myr), Tucana-Horogium ((is) approximately 30 Myr), and TW Hydrae ((is) approximately 10 Myr), are sufficiently well-defined to constrain the ages of individual targets. Somewhat less than half of our targets are high-probability members of one of these MGs. For all of our targets, we combine proposed MG membership with other age indicators where available, including Ca ii HK emission, X-ray activity, and rotation period, to produce a posterior probability distribution of age. SEEDS observations discovered a substellar companion to one of our targets, kappa And, a late B star. We do not detect any other substellar companions, but do find seven new close binary systems, of which one still needs to be confirmed. A detailed analysis of the statistics of this sample, and of the companion mass constraints given our age probability distributions and exoplanet cooling models, will be presented in a forthcoming paper.

The Arches cluster is a young, densely packed massive star cluster in our Galaxy that shows a high level of star formation activity. The nature of the extended non-thermal X-ray emission around the cluster remains unclear. The observed bright Fe K(alpha) line emission at 6.4 keV from material that is neutral or in a low ionization state can be produced either by X-ray photoionization or by cosmic-ray particle bombardment or both. In this paper, we report on the first detection of the extended emission around the Arches cluster above 10 keV with the NuSTAR mission, and present results on its morphology and spectrum. The spatial distribution of the hard X-ray emission is found to be consistent with the broad region around the cluster where the 6.4 keV line is observed. The interpretation of the hard X-ray emission within the context of the X-ray reflection model puts a strong constraint on the luminosity of the possible illuminating hard X-ray source. The properties of the observed emission are also in broad agreement with the low-energy cosmic-ray proton excitation scenario. Key words: cosmic rays - Galaxy: center - ISM: general - X-rays: individual (Arches cluster)

One possible diagnostic of planet formation, orbital migration, and tidal evolution is the angle between a planet's orbital axis and the spin axis of its parent star. In general, Psi cannot be measured, but for transiting planets one can measure the angle lambda between the sky projections of the two axes via the Rossiter-McLaughlin effect. In this paper we showed how to combine measurements of lambda in different systems to derive statistical constraints on Psi, using a Bayesian analysis. The results provided evidence for two distinct modes of planet migration.

The Sloan Lens ACS (SLACS) Survey is an efficient Hubble Space Telescope (HST) Snapshot imaging survey for new galaxy-scale strong gravitational lenses. The targeted lens candidates are selected spectroscopically from the Sloan Digital Sky Survey (SDSS) database of galaxy spectra for having multiple nebular emission lines at a redshift significantly higher than that of the SDSS target galaxy. The SLACS survey is optimized to detect bright early-type lens galaxies with faint lensed sources in order to increase the sample of known gravitational lenses suitable for detailed lensing, photometric, and dynamical modeling. In this paper, the first in a series on the current results of our HST Cycle 13 imaging survey, we present a catalog of 19 newly discovered gravitational lenses, along with nine other observed candidate systems that are either possible lenses, nonlenses, or nondetections. The survey efficiency is thus >=68%. We also present Gemini 8 m and Magellan 6.5 m integral-field spectroscopic data for nine of the SLACS targets, which further support the lensing interpretation. A new method for the effective subtraction of foreground galaxy images to reveal faint background features is presented. We show that the SLACS lens galaxies have colors and ellipticities typical of the spectroscopic parent sample from which they are drawn (SDSS luminous red galaxies and quiescent MAIN sample galaxies), but are somewhat brighter and more centrally concentrated. Several explanations for the latter bias are suggested. The SLACS survey provides the first statistically significant and homogeneously selected sample of bright early-type lens galaxies, furnishing a powerful probe of the structure of early-type galaxies within the half-light radius. The high confirmation rate of lenses in the SLACS survey suggests consideration of spectroscopic lens discovery as an explicit science goal of future spectroscopic galaxy surveys.

A challenge to standard leptonic synchrotron self-Compton (SSC) models is the so-called orphan TeV flares, i.e., enhanced very high energy (VHE) gamma-ray emission without any contemporaneous X-ray flaring activity, that have recently been observed in TeV blazars (e.g., 1ES 1959+650). In order to explain the orphan TeV flare of 1 ES 1959+650 observed in 2002 June, the so-called hadronic synchrotron mirror model has been developed. Here relativistic protons are proposed to exist in the jet and interact with reflected electron synchrotron radiation of the precursor SSC flare. If the reflector is located in the cloud region, time shifts of several days are possible between the precursor and the orphan flare. The external photons, blue-shifted in the comoving jet frame, are able to excite the Delta(1232) resonance when interacting with protons of Lorentz factors gamma(sub p) approx. 10(exp 3)-10(exp 4). The decay products of this resonance include charged pions, which, on decay, give rise to neutrino production during the orphan flare. In this paper we calculate the expected neutrino emission for the 2002 June 4 orphan TeV flare of 1ES 1959+650. We compare our results with the recent observations of AMANDA-II of a neutrino event in spatial and temporal coincidence with the orphan flare of this blazar. We find that the expected neutrino signal from the hadronic synchrotron mirror model is insufficient to explain the claimed neutrino signal from the direction of 1ES 1959+650.

Shock acceleration is a ubiquitous phenomenon in astrophysical plasmas. Plasma waves and their associated instabilities (e.g., Buneman, Weibel, and other two-stream instabilities) created in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a three-dimensional relativistic electromagnetic particle (REMP) code, we have investigated particle acceleration associated with a relativistic electron-positron jet front propagating into an ambient electron-positron plasma with and without initial magnetic fields. We find small differences in the results for no ambient and modest ambient magnetic fields. New simulations show that the Weibel instability created in the collisionless shock front accelerates jet and ambient particles both perpendicular and parallel to the jet propagation direction. Furthermore, the nonlinear fluctuation amplitudes of densities, currents, and electric and magnetic fields in the electron-positron shock are larger than those found in the electron-ion shock studied in a previous paper at a comparable simulation time. This comes from the fact that both electrons and positrons contribute to generation of the Weibel instability. In addition, we have performed simulations with different electron skin depths. We find that growth times scale inversely with the plasma frequency, and the sizes of structures created by tine Weibel instability scale proportionally to the electron skin depth. This is the expected result and indicates that the simulations have sufficient grid resolution. While some Fermi acceleration may occur at the jet front, the majority of electron and positron acceleration takes place behind the jet front and cannot be characterized as Fermi acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying nonuniform, small-scale magnetic fields, which contribute to the electron s (positron s) transverse deflection behind the jet head. This small- scale magnetic field structure is appropriate to the generation of "jitter" radiation from deflected electrons (positrons) as opposed to synchrotron radiation. The jitter radiation has different properties than synchrotron radiation calculated assuming a uniform magnetic field. The jitter radiation resulting from small-scale magnetic field structures may be important for understanding the complex time structure and spectral evolution observed in gamma-ray bursts or other astrophysical sources containing relativistic jets and relativistic collisionless shocks.

We presented infrared photometry obtained with the IRAC camera on the Spitzer Space Telescope of a sample of 82 pre-main-sequence stars and brown dwarfs in the Taurus starforming region. We find a clear separation in some IRAC color-color diagrams between objects with and without disks. A few "transition" objects are noted, which correspond to systems in which the inner disk has been evacuated of small dust. Separating pure disk systems from objects with remnant protostellar envelopes is more difficult at IRAC wavelengths, especially for objects with infall at low rates and large angular momenta. Our results generally confirm the IRAC color classification scheme used in previous papers by Allen et al. and Megeath et al. to distinguish between protostars, T Tauri stars with disks, and young stars without (inner) disks. The observed IRAC colors are in good agreement with recent improved disk models, and in general accord with models for protostellar envelopes derived from analyzing a larger wavelength region. We also comment on a few Taurus objects of special interest. Our results should be useful for interpreting IRAC results in other, less well studied star-forming regions.

Our program consisted of an observation of a single source, 4U 0513-40, which we had previously identified as a candidate ultracompact binary (a system with an orbital period below 1 hour). Several other known or suspected ultracompact binaries have shown unusual elemental abundance ratios in their X-ray spectra. In this program, however, our observation found no unusual abundance ratios in the spectrum of 4U 0513-40. This result was included, along with results from a separate Chandra program, in a paper submitted for publication in the Astrophysical Journal.

This paper reports the ISO results on hard X-ray selected AGN which are less biased against red/obscured objects than other selection wavebands. We find that, as predicted, the IR continuum of these sources extends to redder sources than in optically/radio selected sample. This indicates that the latter samples miss a portion of the population which is fainter in the optical but can be easily picked up in the hard X-ray. The range of IR SEDs is roughly consistent with reddening of the IR continuum up to column densities of around 10(exp 23)/cu cm. Modeling of the full SED using dusty disk models demonstrated that varying the viewing angle can explain the observed SEDs, though rather large disks are required to fit the cooler, long wavelength emission. From the fits we can obtain estimates of the mass and inclination of the system.

In this paper we present optical long-slit spectroscopy and imaging of the protoplanetary nebula CRL 618. The optical lobes of CRL 618 consist of shock-excited gas, which emits many recombination and forbidden lines, and dust, which scatters light from the innermost regions.

We present results of numerical simulations carried out with a 2D radiation hydrodynamics code in order to study the impact of massive stars on their surrounding interstellar medium. This first paper deals with the evolution of the circumstellar gas around an isolated 60 M star.