No Arabic abstract
We model Spitzer Space Telescope observations of the Taurus Class 0 protostar L1527 IRS (IRAS 04368+2557) to provide constraints on its protostellar envelope structure. The nearly edge-on inclination of L1527 IRS, coupled with the highly spatially-resolved near to mid-infrared images of this object and the detailed IRS spectrum, enable us to constrain the outflow cavity geometry quite well, reducing uncertainties in the other derived parameters. The mid-infrared scattered light image shows a bright central source within a dark lane; the aspect ratio of this dark lane is such that it appears highly unlikely to be a disk shadow. In modeling this dark lane, we conclude that L1527 IRS is probably not described by a standard TSC envelope with simple bipolar cavities. We find it necessary to model the dark lane and central source as a modified inner envelope structure. This structure may be due either to a complex wind-envelope interaction or induced by the central binary. To fit the overall SED, we require the central source to have a large near to mid-infrared excess, suggesting substantial disk accretion. Our model reproduces the overall morphology and surface brightness distribution of L1527 IRS fairly well, given the limitations of using axisymmetric models to fit the non-axisymmetric real object, and the derived envelope infall rates are in reasonable agreement with some other investigations. IRAC observations of L1527 IRS taken 12 months apart show variability in total flux and variability in the opposing bipolar cavities, suggesting asymmetric variations in accretion. We also provide model images at high resolution for comparison to future observations with current ground-based instrumentation and future space-based telescopes.
Context. The bipolar outflows associated with low-mass protostars create cavities in the infalling envelope. These cavities are illuminated by the central protostar and inner disk, creating a bipolar scattered light nebula at near-infrared and mid-infrared wavelengths. The variability of the scattered light nebula in both total intensity and intensity as a function of position in the scattered light nebula can provide important insights into the structure of the inner disk that cannot be spatially resolved. Aims. We aim to determine the likelihood that a warped inner disk is the origin of the surface brightness variability in the bipolar scattered light nebula associated with L1527 IRS. Methods. We present results from near-IR imaging conducted over the course of seven years, with periods of monthly cadence monitoring. We used Monte Carlo radiative transfer models to interpret the observations. Results. We find a time varying, asymmetrical brightness in the scattered light nebulae within the outflow cavities of the protostar. Starting in 2007, the surface brightnesses of the eastern and western outflow cavities were roughly symmetric. Then, in 2009, the surface brightnesses of the cavities were found to be asymmetric, with a substantial increase in surface brightness and a larger increase in the eastern outflow cavity. More regular monitoring was conducted from 2011 to 2014, revealing a rotating pattern of surface brightness variability in addition to a slow change of the eastern and western outflow cavities toward symmetry, but still not as symmetric as observed in 2007. We find that an inner disk warp is a feasible mechanism to produce the rotating pattern of surface brightness variability.
Sub-millimeter spectral line and continuum emission from the protoplanetary disks and envelopes of protostars are powerful probes of their structure, chemistry, and dynamics. Here we present a benchmark study of our modeling code, RadChemT, that for the first time uses a chemical model to reproduce ALMA C$^{18}$O (2-1) and CARMA $^{12}$CO (1-0) and N$_{2}$H$^{+}$ (1-0) observations of L1527, that allow us to distinguish the disk, the infalling envelope and outflow of this Class 0/I protostar. RadChemT combines dynamics, radiative transfer, gas chemistry and gas-grain reactions to generate models which can be directly compared with observations for individual protostars. Rather than individually fit abundances to a large number of free parameters, we aim to best match the spectral line maps by (i) adopting a physical model based on density structure and luminosity derived primarily from previous work that fit SED and 2D imaging data, updating it to include a narrow jet detected in CARMA and ALMA data near ($leq 75$au) the protostar, and then (ii) computing the resulting astrochemical abundances for 292 chemical species. Our model reproduces the C$^{18}$O and N$_{2}$H$^{+}$ line strengths within a factor of 3.0; this is encouraging considering the pronounced abundance variation (factor $> 10^3$) between the outflow shell and CO snowline region near the midplane. Further, our modeling confirms suggestions regarding the anti-correlation between N$_{2}$H$^{+}$ and the CO snowline between 400 au to 2,000 au from the central star. Our modeling tools represent a new and powerful capability with which to exploit the richness of spectral line imaging provided by modern submillimeter interferometers.
We present the first resolved near infrared imagery of the transition disk Oph IRS 48 (WLY 2-48), which was recently observed with ALMA to have a strongly asymmetric sub-millimeter flux distribution. H-band polarized intensity images show a $sim$60AU radius scattered light cavity with two pronounced arcs of emission, one from Northeast to Southeast and one smaller, fainter and more distant arc in the Northwest. K-band scattered light imagery reveals a similar morphology, but with a clear third arc along the Southwestern rim of the disk cavity. This arc meets the Northwestern arc at nearly a right angle, revealing the presence of a spiral arm or local surface brightness deficit in the disk, and explaining the East-West brightness asymmetry in the H-band data. We also present 0.8-5.4$mu$m IRTF SpeX spectra of this object, which allow us to constrain the spectral class to A0$pm$1 and measure a low mass accretion rate of 10$^{-8.5}$M$_{odot}$/yr, both consistent with previous estimates. We investigate a variety of reddening laws in order to fit the mutliwavelength SED of Oph IRS 48 and find a best fit consistent with a younger, higher luminosity star than previous estimates.
Subarcsecond-resolution images of the rotational line emissions of CS and c-C$_3$H$_2$ obtained toward the low-mass protostar IRAS 04368$+$2557 in L1527 with the Atacama Large Millimeter/submillimeter Array are investigated to constrain the orientation of the outflow/envelope system. The distribution of CS consists of an envelope component extending from north to south and a faint butterfly-shaped outflow component. The kinematic structure of the envelope is well reproduced by a simple ballistic model of an infalling rotating envelope. Although the envelope has a nearly edge-on configuration, the inclination angle of the rotation axis from the plane of the sky is found to be 5$^circ$, where we find that the western side of the envelope faces the observer. This configuration is opposite to the direction of the large-scale ($sim$ 10$^4$ AU) outflow suggested previously from the $^{12}$CO ($J$=3$-$2) observation, and to the morphology of infrared reflection near the protostar ($sim$ 200 AU). The latter discrepancy could originate from high extinction by the outflow cavity of the western side, these discrepancies or may indicate that the outflow axis is not parallel to the rotation axis of the envelope. Position-velocity diagrams show the accelerated outflow cavity wall, and its kinematic structure in the 2000 AU scale is explained by a standard parabolic model with the inclination angle derived from the analysis of the envelope. The different orientation of the outflow between the small and large scale implies a possibility of precession of the outflow axis. The shape and the velocity of the outflow in the vicinity of the protostar are compared with those of other protostars.
The low-mass protostellar region L1527 is unusual because it contains observable abundances of unsaturated carbon-chain molecules including CnH radicals, H2Cn carbenes, cyanopolyynes, and the negative ions C4H- and C6H-, all of which are more associated with cold cores than with protostellar regions. Sakai et al. suggested that these molecules are formed in L1527 from the chemical precursor methane, which evaporates from the grains during the heat-up of the region. With the gas-phase osu.03.2008 network extended to include negative ions of the families Cn-, and CnH-, as well as the newly detected C3N-, we modeled the chemistry that occurs following methane evaporation at T~ 25-30 K. We are able to reproduce most of the observed molecular abundances in L1527 at a time of ~5000 yr. At later times, the overall abundance of anions become greater than that of electrons, which has an impact on many organic species and ions. The anion-to-neutral ratio in our calculation is in good agreement with observation for C6H- but exceeds the observed ratio by more than three orders of magnitude for C4H-. In order to explain this difference, further investigation is needed on the rate coefficients for electron attachment and other reactions regarding anions.