Circumplanetary disks (CPDs) may be essential to the formation of planets, regulating their spin and accretion evolution. We perform a series of 3D hydrodynamics simulations in both the isothermal and adiabatic limits to systematically measure the rotation rates, sizes, and masses of CPDs as functions of $q_{\rm thermal}$, the ratio of the planet mass to the disk thermal mass. Our $q_{\rm thermal}$ ranges from 0.1 to 4; for our various disk temperatures, this corresponds to planet masses between 1 Earth mass and 4 Jupiter masses. Within this parameter space, we find that isothermal CPDs are disky and bound within $\sim$10\% of the planet's Bondi radius $r_{\rm B}$, with the innermost $\sim0.05\,r_{\rm B}$ in full rotational support. Adiabatic CPDs are spherical (therefore not actually "disks"), bound within $\sim0.2\,r_{\rm B}$, and mainly pressure-supported with rotation rates scaling linearly with $q_{\rm thermal}$; extrapolation suggests full rotational support of adiabatic envelopes at $\sim10\,q_{\rm thermal}$. Fast rotation and 3D supersonic flow render isothermal CPDs significantly different in structure from --- and orders of magnitude less massive than --- their 1D isothermal hydrostatic counterparts. Inside a minimum-mass solar nebula, even a maximally cooled, isothermal CPD around a 10 Earth-mass core may have less than 1 Earth mass, suggesting that gas giant formation may hinge on angular momentum transport processes in CPDs. Our CPD sizes and masses appear consistent with the regular satellites orbiting solar system giants.

Numerical simulations of dust-gas dynamics are one of the fundamental tools in astrophysical research, such as the study of star and planet formation. It is common to find tightly coupled dust and gas in astrophysical systems, which demands that any practical integration method be able to take time steps $\Delta t$ much longer than the stopping time $t_{\rm s}$ due to drag. A number of methods have been developed to ensure stability in this stiff ($\Delta t\gg t_{\rm s}$) regime, but there remains large room for improvement in terms of accuracy. In this paper, we describe an easy-to-implement method, the "staggered semi-analytic method" (SSA), and conduct numerical tests to compare it to other implicit and semi-analytic methods, including the $2^{\rm nd}$ order implicit method and the Verlet method. SSA makes use of a staggered step to better approximate the terminal velocity in the stiff regime. In applications to protoplanetary disks, this not only leads to orders-of-magnitude higher accuracy than the other methods, but also provides greater stability, making it possible to take time steps 100 times larger in some situations. SSA is also $2^{\rm nd}$ order accurate and symplectic when $\Delta t \ll t_{\rm s}$. More generally, the robustness of SSA makes it applicable to linear dust-gas drag in virtually any context.

The wide, deep cavities of transition disks are often believed to have been hollowed out by nascent planetary systems. PDS 70, a ${\sim}5$ Myr old transition disk system in which a multi-Jupiter-mass planet candidate at 22 au coexists with a ${\sim}30$ au gas and ${\sim}60$ au dust-continuum gap, provides a valuable case study for this hypothesis. Using the PEnGUIn hydrodynamics code, we simulate the orbital evolution and accretion of PDS 70b in its natal disk over the lifetime of the system. When the accreting planet reaches about 2.5 Jupiter masses, it spontaneously grows in eccentricity and consumes material from a wide swathe of the PDS 70 disk; radiative transfer post-processing with DALI shows that this accurately reproduces the observed gap profile. Our results demonstrate that super-Jupiter planets can single-handedly carve out transition disk cavities, and indicate that the high eccentricities measured for such giants may be a natural consequence of disk-planet interaction.

Protoplanets can produce structures in protoplanetary disks via gravitational disk-planet interactions. Once detected, such structures serve as signposts of planet formation. Here we investigate the kinematic signatures in disks produced by multi-Jupiter mass ($M_{\rm J}$) planets using 3D hydrodynamics and radiative transfer simulations. Such a planet opens a deep gap, and drives transonic vertical motions inside. Such motions include both a bulk motion of the entire half-disk column, and turbulence on scales comparable to and smaller than the scale height. They significantly broaden molecular lines from the gap, producing double-peaked line profiles at certain locations, and a kinematic velocity dispersion comparable to thermal after azimuthal averaging. The same planet does not drive fast vertical motions outside the gap, except at the inner spiral arms and the disk surface. Searching for line broadening induced by multi-$M_{\rm J}$ planets inside gaps requires an angular resolution comparable to the gap width, an assessment of the gap gas temperature to within a factor of 2, and a high sensitivity needed to detect line emission from the gap.

We present $H$- and $K_{\rm s}$-bands observations of the LkH$\alpha$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investigated the azimuthal profiles at the ring and the outer-disk regions as well as radial profiles in the directions of the spiral-like structures and semi-major axis. Azimuthal analysis shows a large variety in wavelength and implies that the disk has non-axisymmetric dust distributions. The radial profiles in the major-axis direction (PA=$271^\circ$) suggest that the outer region ($r\geq0\farcs25$) may be influenced by shadows of the inner region of the disk. The spiral-like directions (PA=10$^\circ$ and 230$^\circ$) show different radial profiles, which suggests that the surfaces of the spiral-like structures are highly flared and/or have different dust properties. Finally, a color-map of the disk shows a lack of an outer eastern region in the $H$-band disk, which may hint the presence of an inner object that casts a directional shadow onto the disk.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

The majority of gas giants (planets of masses $\gtrsim10^2 M_\oplus$) are found to reside at distances beyond $\sim1$ au from their host stars. Within 1 au, the planetary population is dominated by super-Earths of $2-20 M_\oplus$. We show that this dichotomy between inner super-Earths and outer gas giants can be naturally explained should they form in nearly inviscid disks. In laminar disks, a planet can more easily repel disk gas away from its orbit. The feedback torque from the pile-up of gas inside the planet's orbit slows down and eventually halts migration. A pressure bump outside the planet's orbit traps pebbles and solids, starving the core. Gas giants are born cold and stay cold: more massive cores are preferentially formed at larger distances, and they barely migrate under disk feedback. We demonstrate this using 2D hydrodynamical simulations of disk-planet interaction lasting up to $10^5$ years: we track planet migration and pebble accretion until both come to an end by disk feedback. Whether cores undergo runaway gas accretion to become gas giants or not is determined by computing 1D gas accretion models. Our simulations show that in an inviscid minimum mass solar nebula, gas giants do not form inside $\sim$0.5 au, nor can they migrate there while the disk is present. We also explore the dependence on disk mass, and find that gas giants form further out in less massive disks.

HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD 100453AB between 2015-2017 utilizing extreme adaptive optics systems on the Very Large Telescope and Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a=1.06"+/-0.09", e=0.17+/-0.07, and i=32.5+/- 6.5 degrees. We utilized publicly available ALMA CO data to constrain the inclination of the disk to i~28 degrees, which is relatively co-planar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamical and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-dirven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material.

We explain the fast-moving, ripple-like features in the edge-on debris disk orbiting the young M dwarf AU Mic. The bright features are clouds of sub-micron dust repelled by the host star's wind. The clouds are produced by avalanches: radial outflows of dust that gain exponentially more mass as they shatter background disk particles in collisional chain reactions. The avalanches are triggered from a region a few AU across -- the "avalanche zone" -- located on AU Mic's primary "birth" ring, at a true distance of $\sim$35 AU from the star but at a projected distance more than a factor of 10 smaller: the avalanche zone sits directly along the line of sight to the star, on the side of the ring nearest Earth, launching clouds that disk rotation sends wholly to the southeast, as observed. The avalanche zone marks where the primary ring intersects a secondary ring of debris left by the catastrophic disruption of a progenitor up to Varuna in size, less than tens of thousands of years ago. Only where the rings intersect are particle collisions sufficiently violent to spawn the sub-micron dust needed to seed the avalanches. We show that this picture works quantitatively, reproducing the masses, sizes, and velocities of the observed escaping clouds. The Lorentz force exerted by the wind's magnetic field, whose polarity reverses periodically according to the stellar magnetic cycle, promises to explain the observed vertical undulations. The timescale between avalanches, about 10 yr, might be set by time variability of the wind mass-loss rate or, more speculatively, by some self-regulating limit cycle.

Recently, high angular resolution imaging instruments such as SPHERE and GPI have discovered many spiral-arm-like features in near-infrared scattered light images of protoplanetary disks. Theory and simulations have suggested that these arms are most likely excited by planets forming in the disks; however, a quantitative relation between the arm-to-disk brightness contrast and planet mass is still missing. Using 3D hydrodynamics and radiative transfer simulations, we examine the morphology and contrast of planet-induced arms in disks. We find a power-law relation for the face-on arm contrast ($\delta_{\rm max}$) as a function of planet mass ($M_{\rm p}$) and disk aspect ratio ($h/r$): $\delta_{\rm max}\approx((M_{\rm p}/M_{\rm J})/(h/r)^{1.38})^{0.22}$. With current observational capability, at a 30 AU separation, the minimum planet mass for driving detectable arms in a disk around a 1 Myr 1$M_\odot$ star at 140 pc at low inclinations is around Saturn mass. For planets more massive than Neptune masses, they typically drive multiple arms. Therefore in observed disks with spirals, it is unlikely that each spiral arm originates from a different planet. We also find only massive perturbers with at least multi-Jupiter masses are capable of driving bright arms with $\delta_{\rm max}\gtrsim2$ as found in SAO 206462, MWC 758, and LkH$\alpha$~330, and these arms do not follow linear wave propagation theory. Additionally, we find the morphology and contrast of the primary and secondary arms are largely unaffected by a modest level of viscosity with $\alpha\lesssim0.01$. Finally, the contrast of the arms in the SAO 206462 disk suggests that the perturber SAO 206462 b at $\sim100$ AU is about $5-10M_{\rm J}$ in mass.

High contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their masses, from the observed gap properties, such as their depths and widths. Combing 2D and 3D hydrodynamics simulations with 3D radiative transfer simulations, we investigate the morphology of planet-opened gaps in near-infrared scattered light images. Quantitatively, we obtain correlations that directly link intrinsic gap depths and widths in the gas surface density to observed depths and widths in images of disks at modest inclinations under finite angular resolution. Subsequently, the properties of the surface density gaps enable us to derive the disk scale height at the location of the gap $h$, and to constrain the quantity $M_{\rm p}^2/\alpha$, where $M_{\rm p}$ is the mass of the gap-opening planet and $\alpha$ characterizes the viscosity in the gap. As examples, we examine the gaps recently imaged by VLT/SPHERE, Gemini/GPI, and Subaru/HiCIAO in HD 97048, TW Hya, HD 169142, LkCa 15, and RX J1615.3-3255. Scale heights of the disks and possible masses of the gap-opening planets are derived assuming each gap is opened by a single planet. Assuming $\alpha=10^{-3}$, the derived planet mass in all cases are roughly between 0.1-1 $M_{\rm J}$.

We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at H{\alpha}. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

Two longstanding problems in planet formation include (1) understanding how planets survive migration, and (2) articulating the process by which protoplanetary disks disperse---and in particular how they accrete onto their central stars. We can go a long way toward solving both problems if the disk gas surrounding planets has no intrinsic diffusivity ("viscosity"). In inviscid, laminar disks, a planet readily repels gas away from its orbit. On short timescales, zero viscosity gas accumulates inside a planet's orbit to slow Type I migration by orders of magnitude. On longer timescales, multiple super-Earths (distributed between, say, $\sim$0.1--10 AU) can torque inviscid gas out of interplanetary space, either inward to feed their stars, or outward to be blown away in a wind. We explore this picture with 2D hydrodynamics simulations of Earths and super-Earths embedded in inviscid disks, confirming their slow/stalled migration even under gas-rich conditions, and showing that disk transport rates range up to $\sim$$10^{-7} M_\odot~{\rm yr^{-1}}$ and scale as $\dot{M} \propto \Sigma M_{\rm p}^{3/2}$, where $\Sigma$ is the disk surface density and $M_{\rm p}$ is the planet mass. Gas initially sandwiched between two planets is torqued past both into the inner and outer disks. In sum, sufficiently compact systems of super-Earths can clear their natal disk gas, in a dispersal history that may be complicated and non-steady, but which conceivably leads over Myr timescales to large gas depletions similar to those characterizing transition disks.

Planet migration is inherently a three-dimensional (3D) problem, because Earth-size planetary cores are deeply embedded in protoplanetary disks. Simulations of these 3D disks remain challenging due to the steep requirement in resolution. Using two different hydrodynamics code, FARGO3D and PEnGUIn, we simulate disk-planet interaction for a 1 to 5 Earth-mass planet embedded in an isentropic disk. We measure the torque on the planet and ensure that the measurements are converged both in resolution and between the two codes. We find that the torque is independent of the smoothing length of the planet's potential ($r_{\rm s}$), and that it has a weak dependence on the adiabatic index of the gaseous disk ($\gamma$). The torque values correspond to an inward migration rate qualitatively similar to previous linear calculations. We perform additional simulations with explicit radiative transfer using FARGOCA, and again find agreement between 3D simulations and existing torque formulae. We also present the flow pattern around the planets, and show that active flow is present within the planet's Hill sphere, and meridional vortices are shed downstream. The vertical flow speed near the planet is faster for a smaller $r_{\rm s}$ or $\gamma$, up to supersonic speeds for the smallest $r_{\rm s}$ and $\gamma$ in our study.

Giant planets can clear deep gaps when embedded in 2D (razor-thin) viscous circumstellar disks. We show by direct simulation that giant planets are just as capable of carving out gaps in 3D. Surface density maps are similar between 2D and 3D, even in detail. In particular, the scaling $\Sigma_{\rm gap} \propto q^{-2}$ of gap surface density with planet mass, derived from a global "zero-dimensional" balance of Lindblad and viscous torques, applies equally well to results obtained at higher dimensions. Our 3D simulations reveal extensive, near-sonic, meridional flows both inside and outside the gaps; these large-scale circulations might bear on disk compositional gradients, in dust or other chemical species. At high planet mass, gap edges are mildly Rayleigh unstable and intermittently shed streams of material into the gap - less so in 3D than in 2D.

Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a disk by a companion varies with inclination and position angle in NIR scattered light. We compare the cases of a $3M_{\rm J}$ and a $0.1M_{\odot}$ companion, and make predictions suitable for testing with Gemini/GPI, VLT/NACO/SPHERE, and Subaru/HiCIAO/SCExAO. We find that the two trailing arms produced by an external perturber can have a variety of morphologies in inclined systems - they may appear as one trailing arm; two trailing arms on the same side of the disk; or two arms winding in opposite directions. The disk ring outside a planetary gap may also mimic spiral arms when viewed at high inclinations. We suggest potential explanations for the features observed in HH 30, HD 141569 A, AK Sco, HD 100546, and AB Aur. We emphasize that inclined views of companion-induced features cannot be converted into face-on views using simple and commonly practiced image deprojections.

Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.

Recent observations of protoplanetary disk have reported spiral structures that are potential signatures of embedded planets, and modeling efforts have shown that a single planet can excite multiple spiral arms, in contrast to conventional disk-planet interaction theory. Using two and three-dimensional hydrodynamics simulations to perform a systematic parameter survey, we confirm the existence of multiple spiral arms in disks with a single planet, and discover a scaling relation between the azimuthal separation of the primary and secondary arm, $\phi_{\rm sep}$, and the planet-to-star mass ratio $q$: $\phi_{\rm sep} = 102^{\circ} (q/0.001)^{0.2}$ for companions between Neptune mass and 16 Jupiter masses around a 1 solar mass star, and $\phi_{\rm sep} = 180^{\circ}$ for brown dwarf mass companions. This relation is independent of the disk's temperature, and can be used to infer a planet's mass to within an accuracy of about 30% given only the morphology of a face-on disk. Combining hydrodynamics and Monte-Carlo radiative transfer calculations, we verify that our numerical measurements of $\phi_{\rm sep}$ are accurate representations of what would be measured in near-infrared scattered light images, such as those expected to be taken by Gemini/GPI, VLT/SPHERE, or Subaru/SCExAO in the future. Finally, we are able to infer, using our scaling relation, that the planet responsible for the spiral structure in SAO 206462 has a mass of about 6 Jupiter masses.

3D modifications to the well-studied 2D flow topology around an embedded planet have the potential to resolve long-standing problems in planet formation theory. We present a detailed analysis of the 3D isothermal flow field around a 5 Earth-mass planet on a fixed circular orbit, simulated using our multi-GPU hydrodynamics code PEnGUIn. We find that, overall, the horseshoe region has a columnar structure extending vertically much beyond the Hill sphere of the planet. This columnar structure is only broken for some of the widest horseshoe streamlines, along which high altitude fluid descends rapidly into the planet's Bondi sphere, performs one horseshoe turn, and exits the Bondi sphere radially in the midplane. A portion of this flow exits the horseshoe region altogether, which we refer to as the "transient" horseshoe flow. The flow continues as it rolls up into a pair of up-down symmetric horizontal vortex lines shed into the wake of the planet. This flow, unique to 3D, affects both planet accretion and migration. It prevents the planet from sustaining a hydrostatic atmosphere due to its intrusion into the Bondi sphere, and leads to a significant corotation torque on the planet, unanticipated by 2D analysis. In the reported simulation, starting with a $\Sigma\sim r^{-3/2}$ radial surface density profile, this torque is positive and partially cancels with the negative differential Lindblad torque, resulting in a factor of 3 slower planet migration rate. Finally, we report 3D effects can be suppressed by a sufficiently large disk viscosity, leading to results similar to 2D.

An instability can potentially operate in highly irradiated disks where the disk sharply transitions from being radially transparent to opaque (the 'transition region'). Such conditions may exist at the inner edges of transitional disks around T Tauri stars and accretion disks around AGNs. We derive the criterion for this instability, which we term the 'irradiation instability', or IRI. We also present the linear growth rate as a function of beta, the ratio between radiation force and gravity, and c_s, the sound speed of the disk, obtained using two methods: a semi-analytic analysis of the linearized equations and a numerical simulation using the GPU-accelerated hydrodynamical code PEnGUIn. In particular, we find that IRI occurs at beta~0.1 if the transition region extends as wide as ~0.05r, and at higher beta values if it is wider. This threshold value applies to c_s ranging from 3% of the Keplerian orbital speed to 5%, and becomes higher if c_s is lower. Furthermore, in the nonlinear evolution of the instability, disks with a large beta and small c_s exhibit 'clumping', extreme local surface density enhancements that can reach over ten times the initial disk surface density.

Gap clearing by giant planets has been proposed to explain the optically thin cavities observed in many protoplanetary disks. How much material remains in the gap determines not only how detectable young planets are in their birth environments, but also how strong corotation torques are, which impacts how planets can survive fast orbital migration. We determine numerically how the average surface density inside the gap, sigma_gap, depends on planet-to-star mass ratio q, Shakura-Sunyaev viscosity parameter alpha, and disk height-to-radius aspect ratio h/r. Our results are derived from our new GPU-accelerated Lagrangian hydrodynamical code PEnGUIn, and are verified by independent simulations with ZEUS90. For Jupiter-like planets, we find sigma_gap \propto q^-2.2 alpha^1.4 (h/r)^6.6, and for near brown dwarf masses, sigma_gap \propto q^-1 alpha^1.3 (h/r)^6.1. Surface density contrasts inside and outside gaps can be as large as 10^4, even when the planet does not accrete. We derive a simple analytic scaling, sigma_gap \propto q^-2 alpha^1 (h/r)^5, that compares reasonably well to empirical results, especially at low Neptune-like masses, and use discrepancies to highlight areas for progress.

In this paper we refine the calculation of primordial helium recombination, accounting for several additional effects that were neglected or treated more approximately in previous studies. These include consideration of (i) time-dependent radiative transfer interaction between the 2^1 P_1 - 1^1 S_0 and 2^3 P_1 - 1^1 S_0 resonances; (ii) time-dependent radiative transfer for the partially overlapping n^1 P_1 - 1^1 S_0, n^1 D_2 - 1^1 S_0 and n^3 P_1 - 1^1 S_0 series with 3 <= n <= 10; (iii) electron scattering within a kernel approach. We also briefly discuss the effect of electron scattering and HI quadrupole lines on the recombination of hydrogen. Although the physics of all the considered processes is interesting and subtle, for the standard cosmology, with Y_p \sim 0.24, the overall correction to the ionization history during helium recombination with respect to the previous implementation of CosmoRec remains smaller than |DeltaNe/Ne| \sim 0.05%. The dominant improvement is caused by consistent inclusion of resonance scattering for the 2^1 P_1 - 1^1 S_0 resonance. For cosmologies with a large helium fraction, Y_p \sim 0.4, the difference reaches |D N_e/N_e|\sim 0.22% at z\sim1800, however, the overall correction to the CMB power spectra is small, exceeding |D C_l/C_l| \sim 0.05 % only at l>\sim 3000. In comparison to the current version of Recfast the difference reaches |D C_l/C_l|\sim0.4% at l\sim 3000 for Y_p \sim 0.4, and also for the standard value Y_p \sim 0.24 we find differences |DC_l/C_l|>\sim 0.1% at l\sim 2500. The new processes are now included by the cosmological recombination code CosmoRec and can be activated as needed for most settings without affecting its runtime significantly.