We present an overview of the Large Program, “Early Planet Formation in Embedded Disks (eDisk),” conducted with the Atacama Large Millimeter/submillimeter Array (ALMA). The ubiquitous detections of substructures, particularly rings and gaps, in protoplanetary disks around T Tauri stars raise the possibility that at least some planet formation may have already started during the embedded stages of star formation. In order to address exactly how and when planet formation is initiated, the program focuses on searching for substructures in disks around 12 Class 0 and 7 Class I protostars in nearby (<200 pc) star-forming regions through 1.3 mm continuum observations at a resolution of ∼7 au (0.″04). The initial results show that the continuum emission, mostly arising from dust disks around the sample protostars, has relatively few distinctive substructures, such as rings and spirals, in marked contrast to Class II disks. The dramatic difference may suggest that substructures quickly develop in disks when the systems evolve from protostars to Class II sources, or alternatively that high optical depth of the continuum emission could obscure internal structures. Kinematic information obtained through CO isotopologue lines and other lines reveals the presence of Keplerian disks around protostars, providing us with crucial physical parameters, in particular, the dynamical mass of the central protostars. We describe the background of the eDisk program, the sample selection and their ALMA observations, and the data reduction, and we also highlight representative first-look results.

The Hall effect is recently shown to be efficient in magnetized dense molecular cores and could lead to a bimodal formation of rotationally supported discs (RSDs) in the first core phase. However, how such Hall dominated systems evolve in the protostellar accretion phase remains unclear. We carry out 2D axisymmetric simulations including Hall effect and ohmic dissipation, with realistic magnetic diffusivities computed from our equilibrium chemical network. We find that Hall effect only becomes efficient when the large population of very small grains (VSGs: ≲100 Å) is removed from the standard Mathis–Rumpl–Nordsieck size distribution. With such an enhanced Hall effect, however, the bimodality of disc formation does not continue into the main accretion phase. The outer part of the initial ∼40 au disc formed in the anti-aligned configuration ($\boldsymbol {\Omega \cdot B}\lt 0$) flattens into a thin rotationally supported Hall current sheet as Hall effect moves the poloidal magnetic field radially inward relative to matter, leaving only the inner ≲10–20 au RSD. In the aligned configuration ($\boldsymbol {\Omega \cdot B}\gt 0$), disc formation is suppressed initially but a counter-rotating disc forms subsequently due to efficient azimuthal Hall drift. The counter-rotating disc first grows to ∼30 au as Hall effect moves the magnetic field radially outward, but only the inner ≲10 au RSD is long lived like in the anti-aligned case. Besides removing VSGs, cosmic ray ionization rate should be below a few 10−16 s−1 for Hall effect to be efficient in disc formation. We conclude that Hall effect produces small ≲10–20 au discs regardless of the polarity of the magnetic field, and that radially outward diffusion of magnetic fields remains crucial for disc formation and growth.

Disks are essential to the formation of both stars and planets, but how they form in magnetized molecular cloud cores remains debated. This work focuses on how the disk formation is affected by turbulence and ambipolar diffusion (AD), both separately and in combination, with an emphasis on the protostellar mass accretion phase of star formation. We find that a relatively strong, sonic turbulence on the core scale strongly warps but does not completely disrupt the well-known magnetically-induced flattened pseudodisk that dominates the inner protostellar accretion flow in the laminar case, in agreement with previous work. The turbulence enables the formation of a relatively large disk at early times with or without ambipolar diffusion, but such a disk remains strongly magnetized and does not persist to the end of our simulation unless a relatively strong ambipolar diffusion is also present. The AD-enabled disks in laminar simulations tend to fragment gravitationally. The disk fragmentation is suppressed by initial turbulence. The ambipolar diffusion facilitates the disk formation and survival by reducing the field strength in the circumstellar region through magnetic flux redistribution and by making the field lines there less pinched azimuthally, especially at late times. We conclude that turbulence and ambipolar diffusion complement each other in promoting disk formation. The disks formed in our simulations inherit a rather strong magnetic field from its parental core, with a typical plasma-β of order a few tens or smaller, which is 2-3 orders of magnitude lower than the values commonly adopted in MHD simulations of protoplanetary disks. To resolve this potential tension, longer-term simulations of disk formation and evolution with increasingly more realistic physics are needed.

Understanding how material accretes onto the rotationally supported disk from the surrounding envelope of gas and dust in the youngest protostellar systems is important for describing how disks are formed. Magnetohydrodynamic simulations of magnetized, turbulent disk formation usually show spiral-like streams of material (accretion flows) connecting the envelope to the disk. However, accretion flows in these early stages of protostellar formation still remain poorly characterized, due to their low intensity, and possibly some extended structures are disregarded as being part of the outflow cavity. We use ALMA archival data of a young Class 0 protostar, Lupus 3-MMS, to uncover four extended accretion flow–like structures in C18O that follow the edges of the outflows. We make various types of position–velocity cuts to compare with the outflows and find the extended structures are not consistent with the outflow emission, but rather more consistent with a simple infall model. We then use a dendrogram algorithm to isolate five substructures in position–position–velocity space. Four out of the five substructures fit well (>95%) with our simple infall model, with specific angular momenta between 2.7–6.9 × 10−4 km s−1 pc and mass-infall rates of 0.5–1.1 × 10−6 M ⊙ yr−1. Better characterization of the physical structure in the supposed “outflow cavities” is important to disentangle the true outflow cavities and accretion flows.