SED modeling of the Class 0 protostar L1527 IRS

We model the spectral energy distribution (SED) of the Class 0 protostar L1527 IRS using a dust continuum radiative transfer code RADMC-3D to study the initial condition of gravitational collapse. To constrain the envelope structure, we use the data obtained by Herschel /PACS, which covers the far-infrared regime ($55-190{\mu}m$) where the SED of L1527 IRS peaks. According to our modeling, a more flattened density profile fits the far-infrared SED of L1527 IRS better than the density profile of a rotating and infalling envelope. Thus, we employ the density structure of a Bonnor-Ebert sphere, which consists of the inner flat-topped and the outer power-law regions and is often used for describing the density structure of the youngest sources in the low mass star formation process. A Bonnor-Ebert sphere fits very well the observed SED at ${\lambda}$ > $10{\mu}m$, suggesting that L1527 IRS might collapse from an unstable Bonnor-Ebert sphere rather than a singular isothermal sphere.