Showing 1–50 of 167 results for author: Price, D J

Dust formation during the interaction of binary stars by common envelope

Authors: Luis C. Bermúdez-Bustamante, Orsola De Marco, Lionel Siess, Daniel J. Price, Miguel González-Bolívar, Mike Y. M. Lau, Chunliang Mu, Ryosuke Hirai, Taïssa Danilovich, Mansi Kasliwal

Abstract: We performed numerical simulations of the common envelope (CE) interaction between two intermediate-mass asymptotic giant branch (AGB) stars and their low-mass companions. For the first time, formation and growth of dust in the envelope is calculated explicitly. We find that the first dust grains appear as early as $\sim$1-3 yrs after the onset of the CE, and are smaller than grains formed later.… ▽ More We performed numerical simulations of the common envelope (CE) interaction between two intermediate-mass asymptotic giant branch (AGB) stars and their low-mass companions. For the first time, formation and growth of dust in the envelope is calculated explicitly. We find that the first dust grains appear as early as $\sim$1-3 yrs after the onset of the CE, and are smaller than grains formed later. As the simulations progress, a high-opacity dusty shell forms, resulting in the CE photosphere being up to an order of magnitude larger than it would be without the inclusion of dust. At the end of the simulations, the total dust yield is $0.0082~M_{\odot}$ ($0.022~M_{\odot}$) for a CE with a $1.7~M_{\odot}$ ($3.7~M_{\odot}$) AGB star. Dust formation does not substantially lead to more mass unbinding or substantially alter the orbital evolution. △ Less

Submitted 10 July, 2024; originally announced July 2024.

Comments: 6 pages, 6 figures, contributed talk to the Proceedings IAU Symposium No. 384

Compaction during fragmentation and bouncing produces realistic dust grain porosities in protoplanetary discs

Authors: Stéphane Michoulier, Jean-François Gonzalez, Daniel J. Price

Abstract: Context: In protoplanetary discs, micron-sized dust grows to form millimetre- to centimetre-sized pebbles but encounters several barriers during its evolution. Collisional fragmentation and radial drift impede further dust growth to planetesimal size. Fluffy grains have been hypothesised to solve these problems. While porosity leads to faster grain growth, the implied porosity values obtained from… ▽ More Context: In protoplanetary discs, micron-sized dust grows to form millimetre- to centimetre-sized pebbles but encounters several barriers during its evolution. Collisional fragmentation and radial drift impede further dust growth to planetesimal size. Fluffy grains have been hypothesised to solve these problems. While porosity leads to faster grain growth, the implied porosity values obtained from previous simulations were larger than suggested by observations. Aims: In this paper, we study the influence of porosity on dust evolution taking into account growth, bouncing, fragmentation, compaction, rotational disruption and snow lines, in order to understand their impact on dust evolution. Methods: We develop a module for porosity evolution for the 3D Smoothed Particle Hydrodynamics (SPH) code Phantom that accounts for dust growth and fragmentation. This mono-disperse model is integrated into both a 1D code and the 3D code to capture the overall evolution of dust and gas. Results: We show that porosity helps dust growth and leads to the formation of larger solids than when considering compact grains, as predicted by previous work. Our simulations taking into account compaction during fragmentation show that large millimetre grains are still formed, but are 10 to 100 times more compact. Thus, mm sizes with typical filling factors of ~0.1 match the values measured on comets or via polarimetric observations of protoplanetary discs. △ Less

Submitted 21 June, 2024; originally announced June 2024.

Comments: 24 pages, 21 figures, accepted for publication in A&A

V892 Tau: A tidally perturbed circumbinary disc in a triple stellar system

Authors: Antoine Alaguero, Nicolás Cuello, François Ménard, Simone Ceppi, Álvaro Ribas, Rebecca Nealon, Miguel Vioque, Andrés Izquierdo, James Miley, Enrique Macías, Daniel J. Price

Abstract: V892 Tau is a young binary star surrounded by a circumbinary disc which show hints of interaction with the low-mass nearby star V892 Tau NE. The goal of this paper is to constrain the orbit of V892 Tau NE and to determine the resulting circumbinary disc dynamics. We present new ALMA observations of the V892 Tau circumbinary disc at a twice higher angular and spectral resolution. We model the data… ▽ More V892 Tau is a young binary star surrounded by a circumbinary disc which show hints of interaction with the low-mass nearby star V892 Tau NE. The goal of this paper is to constrain the orbit of V892 Tau NE and to determine the resulting circumbinary disc dynamics. We present new ALMA observations of the V892 Tau circumbinary disc at a twice higher angular and spectral resolution. We model the data with V892 Tau as a triple system and perform a grid of hydrodynamical simulations testing several orbits of the companion. The simulation outputs are then post-processed to build synthetic maps that we compare to the observations. The 12CO emission of the disc shows clear non-Keplerian features such as spiral arms. When comparing the data with our synthetic observations, we interpret these features as ongoing interactions with the companion. Our simulations indicate that an eccentricity of 0.5 of the companion is needed to reproduce the observed disc extent and that a mutual inclination of approximately 60° with the inner binary reproduces the measured disc tilt. In order to explain most of the features of the circumbinary disc, we propose that V892 Tau NE follows a misaligned eccentric orbit, with an eccentricity between 0.2 and 0.5 and a mutual inclination between 30° and 60°. Such a misaligned companion suggests the disc is oscillating and precessing with time, stabilising in an intermediate plane with a non-zero mutual inclination with the inner binary. Given that orbital configuration, we show that the stability of future planets is compromised in the second half of the disc once the gas has dissipated. △ Less

Submitted 21 May, 2024; originally announced May 2024.

Comments: Accepted for publication in A&A, 29 pages, 15 figures

Decretion disc size in Be/X-ray binaries depends upon the disc aspect ratio

Authors: Rebecca G. Martin, Stephen H. Lubow, Philip J. Armitage, Daniel J. Price

Abstract: With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation… ▽ More With three-dimensional hydrodynamical simulations we show that the size of the decretion disc and the structure of the accretion flow onto the neutron star in a Be/X-ray binary strongly depends upon the disc aspect ratio, $H/R$. We simulate a Be star disc that is coplanar to the orbit of a circularly or moderately eccentric neutron star companion, thereby maximising the effects of tidal truncation. For low disc aspect ratio, $H/R\lesssim 0.1$, the disc is efficiently tidally truncated by the neutron star. Most material that escapes the Roche lobe of the Be star is accreted by the neutron star through tidal streams. For larger disc aspect ratio, the outflow rate through the Be star disc is higher, tidal truncation becomes inefficient, the disc fills the Roche lobe and extends to the orbit of the companion. Some material escapes the binary as a gas stream that begins near the L2 point. While the accretion rate onto the neutron star is higher, the fraction of the outflow that is accreted by the neutron star is smaller. Low density Be star discs are expected to be approximately isothermal, such that $H/R$ increases with radius. Tidal truncation is therefore weaker for larger separation binaries, and lower mass primaries. △ Less

Submitted 27 April, 2024; originally announced April 2024.

Comments: Accepted for publication in MNRAS

Eddington envelopes: The fate of stars on parabolic orbits tidally disrupted by supermassive black holes

Authors: Daniel J. Price, David Liptai, Ilya Mandel, Joanna Shepherd, Giuseppe Lodato, Yuri Levin

Abstract: Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the retur… ▽ More Stars falling too close to massive black holes in the centres of galaxies can be torn apart by the strong tidal forces. Simulating the subsequent feeding of the black hole with disrupted material has proved challenging because of the range of timescales involved. Here we report a set of simulations that capture the relativistic disruption of the star, followed by one year of evolution of the returning debris stream. These reveal the formation of an expanding asymmetric bubble of material extending to hundreds of astronomical units -- an outflowing Eddington envelope with an optically thick inner region. Such envelopes have been hypothesised as the reprocessing layer needed to explain optical/UV emission in tidal disruption events, but never produced self-consistently in a simulation. Our model broadly matches the observed light curves with low temperatures, faint luminosities, and line widths of 10,000--20,000 km/s. △ Less

Submitted 9 July, 2024; v1 submitted 14 April, 2024; originally announced April 2024.

Comments: 12 pages + appendices, 15 figures, accepted to ApJL. Data+movies on zenodo at https://doi.org/10.5281/zenodo.11438154 . Movies on YouTube... watch this one first: https://www.youtube.com/watch?v=KlffYCtYxFU then: https://www.youtube.com/watch?v=bPurbeQNgvI&t=19s and finally watch the expanding "bubble" with rho and T at photosphere: https://www.youtube.com/watch?v=D4dyZeUbfvM

Partial tidal disruption events: The elixir of life

Authors: Megha Sharma, Daniel J. Price, Alexander Heger

Abstract: In our Galactic Center, about 10,000 to 100,000 stars are estimated to have survived tidal disruption events, resulting in partially disrupted remnants. These events occur when a supermassive black hole (SMBH) tidally interacts with a star, but not enough to completely disrupt the star. We use the 1D stellar evolution code Kepler and the 3D smoothed particle hydrodynamics code Phantom to model the… ▽ More In our Galactic Center, about 10,000 to 100,000 stars are estimated to have survived tidal disruption events, resulting in partially disrupted remnants. These events occur when a supermassive black hole (SMBH) tidally interacts with a star, but not enough to completely disrupt the star. We use the 1D stellar evolution code Kepler and the 3D smoothed particle hydrodynamics code Phantom to model the tidal disruption of 1, 3, and 10 solar mass stars at zero-age (ZAMS), middle-age (MAMS), and terminal-age main-sequence (TAMS). We map the disruption remnants into Kepler in order to understand their post-distribution evolution. We find distinct characteristics in the remnants, including increased radius, rapid core rotation, and differential rotation in the envelope. The remnants undergo composition mixing that affects their stellar evolution. Whereas the remnants formed by disruption of ZAMS models evolve similarly to unperturbed models of the same mass, for MAMS and TAMS stars, the remnants have higher luminosity and effective temperature. Potential observational signatures include peculiarities in nitrogen and carbon abundances, higher luminosity, rapid rotation, faster evolution, and unique tracks in the Hertzsprung-Russell diagram △ Less

Submitted 6 March, 2024; originally announced March 2024.

Comments: 22 pages, 27 figures, submitted to MNRAS, comments welcome

The orbit of HD 142527 B is too compact to explain many of the disc features

Authors: M. Nowak, S. Rowther, S. Lacour, F. Meru, R. Nealon, D. J. Price

Abstract: HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13… ▽ More HD 142527 A is a young and massive Herbig Ae/Be star surrounded by a highly structured disc. The disc shows numerous morphological structures, such as spiral arms, a horseshoe region of dust emission, a set of shadows cast by an inner disc on the outer disc, and a large cavity extending from $\simeq{}$30 au to $\simeq{}$130 au. HD 142527 A also has a lower mass companion, HD 142527 B (M = 0.13 $\pm$ 0.03 $M_\odot{}$), which is thought to be responsible for most of the structures observed in the surrounding disc. We gathered VLTI/GRAVITY observations of HD 142527, either from our own programmes or from the ESO archive. We used this inhomogeneous set of data to extract a total of seven high-precision measurements of the relative astrometry between HD 142527 A and B, spread from mid-2017 to early 2021. Combined with what is available in the literature, we now have 9 yr of astrometric monitoring on HD 142527. We used orbit fitting tools to determine the orbital parameters of HD 142527 B, and used them as inputs for a 3D hydrodynamical model of the disc to determine whether or not the binary is able to create the structures observed in the disc. Our VLTI/GRAVITY astrometry gives excellent constraints on the orbit of HD 142527 B. We show that the secondary is following an orbit of semi-major axis a = 10.80 $\pm$ 0.22 au, with moderate eccentricity (e = 0.47 $\pm$ 0.01). With such a compact orbit, we show that HD 142527 B can only generate a gap and spiral arms of $\sim$30 au in the disc, which is much smaller than what is revealed by observations. Even from a theoretical standpoint, the observed cavity size of $\sim$100 au far exceeds even the most generous predictions for a companion like HD 142527 B on such a compact orbit. Thus, we conclude that the low-mass companion cannot be solely responsible for the observed morphology of the disc surrounding the system. △ Less

Submitted 5 February, 2024; originally announced February 2024.

Comments: 11 pages, 6 figures. Accepted for publication in A&A

Probing initial distributions of orbital eccentricity and disc misalignment via polar discs

Authors: Simone Ceppi, Nicolás Cuello, Giuseppe Lodato, Cristiano Longarini, Daniel J. Price, Daniel Elsender, Matthew R. Bate

Abstract: In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the di… ▽ More In a population of multiple protostellar systems with discs, the sub-population of circumbinary discs whose orbital plane is highly misaligned with respect to the binary's orbital plane constrains the initial distribution of orbital parameters of the whole population. We show that by measuring the polar disc fraction and the average orbital eccentricity in the polar discs, one can constrain the distributions of initial eccentricity and mutual inclination in multiple stellar systems at birth. △ Less

Submitted 8 January, 2024; originally announced January 2024.

Comments: accepted for publication in A&A

Dust formation in common envelope binary interactions -- II: 3D simulations with self-consistent dust formation

Authors: Luis C. Bermúdez-Bustamante, Orsola De Marco, Lionel Siess, Daniel J. Price, Miguel González-Bolívar, Mike Y. M. Lau, Chunliang Mu, Ryosuke Hirai, Taïssa Danilovich, Mansi M. Kasliwal

Abstract: We performed numerical simulations of the common envelope (CE) interaction between two thermally-pulsing asymptotic giant branch (AGB) stars of 1.7 $M_\odot$ and 3.7 $M_\odot$, and their 0.6 $M_\odot$ compact companion. We use tabulated equations of state to take into account recombination energy. For the first time, formation and growth of dust in the envelope is calculated explicitly, using a ca… ▽ More We performed numerical simulations of the common envelope (CE) interaction between two thermally-pulsing asymptotic giant branch (AGB) stars of 1.7 $M_\odot$ and 3.7 $M_\odot$, and their 0.6 $M_\odot$ compact companion. We use tabulated equations of state to take into account recombination energy. For the first time, formation and growth of dust in the envelope is calculated explicitly, using a carbon dust nucleation network with a gas phase C/O number ratio of 2.5. By the end of the simulations, the total dust yield are $\sim8.2\times10^{-3}~M_\odot$ and $\sim2.2\times10^{-2}~M_\odot$ for the CE with a 1.7 $M_\odot$ and a 3.7 $M_\odot$ AGB star, respectively, close to the theoretical limit. Dust formation does not substantially lead to more mass unbinding or substantially alter the orbital evolution. The first dust grains appear as early as $\sim1-3$ yrs after the onset of the CE rapidly forming an optically thick shell at $\sim10-20$ au, growing in thickness and radius to values of $\sim400-500$ au by $\sim40$ yrs. These large objects have approximate temperatures of 400 K. While dust yields are commensurate with those of single AGB stars of comparable mass, the dust in CE ejections forms over decades as opposed to tens of thousands of years. It is likely that these rapidly evolving IR objects correspond to the post-optically-luminous tail of the lightcurve of some luminous red novae. The simulated characteristics of dusty CEs also lend further support to the idea that extreme carbon stars and the so called ``water fountains" may be objects observed in the immediate aftermath of a CE event. △ Less

Submitted 7 January, 2024; originally announced January 2024.

Comments: 17 pages, 17 figures, 1 table

Optical Appearance of Eccentric Tidal Disruption Events

Authors: Fangyi, Hu, Daniel J. Price, Ilya Mandel

Abstract: Stars approaching supermassive black holes can be tidally disrupted. Despite being expected to emit X-rays, TDEs have been largely observed in optical bands, which is poorly understood. In this Letter, we simulate the tidal disruption of a $1~M_\odot$ main sequence star on an eccentric ($e=0.95$) orbit with a periapsis distance one or five times smaller than the tidal radius ($β= 1$ or $5$) using… ▽ More Stars approaching supermassive black holes can be tidally disrupted. Despite being expected to emit X-rays, TDEs have been largely observed in optical bands, which is poorly understood. In this Letter, we simulate the tidal disruption of a $1~M_\odot$ main sequence star on an eccentric ($e=0.95$) orbit with a periapsis distance one or five times smaller than the tidal radius ($β= 1$ or $5$) using general relativistic smoothed particle hydrodynamics. We follow the simulation for up to a year post-disruption. We show that accretion disks in eccentric TDEs are masked by unbound material outflowing at $\sim10,000~$km/s. Assuming electron scattering opacity, this material would be visible as a $\sim100~$au photosphere at $\sim10^4~$K, in line with observations of candidate TDEs. △ Less

Submitted 5 December, 2023; originally announced December 2023.

Comments: 11 pages, 5 figures, submitted to ApJL

Journal ref: ApJL 963 L27 (2024)

The Automatic Identification and Tracking of Coronal Flux Ropes -- Part II: New Mathematical Morphology-based Flux Rope Extraction Method and Deflection Analysis

Authors: Andreas Wagner, Slava Bourgeois, Emilia K. J. Kilpua, Ranadeep Sarkar, Daniel J. Price, Anshu Kumari, Jens Pomoell, Stefaan Poedts, Teresa Barata, Robertus Erdélyi, Orlando Oliveira, Ricardo Gafeira

Abstract: We present a magnetic flux rope (FR) extraction tool for solar coronal magnetic field modelling data, which builds upon the methodology from Wagner et al. (2023). We apply the scheme to magnetic field simulations of active regions AR12473 and AR11176. We compare the method to its predecessor and study the 3D movement of the newly extracted FRs up to heights of 200 and 300 Mm, respectively. The ext… ▽ More We present a magnetic flux rope (FR) extraction tool for solar coronal magnetic field modelling data, which builds upon the methodology from Wagner et al. (2023). We apply the scheme to magnetic field simulations of active regions AR12473 and AR11176. We compare the method to its predecessor and study the 3D movement of the newly extracted FRs up to heights of 200 and 300 Mm, respectively. The extraction method is based on the twist parameter and a variety of mathematical morphology algorithms, including the opening transform and the morphological gradient. We highlight the differences between the methods by investigating the circularity of the FRs in the plane we extract from. The simulations for the active regions are carried out with a time-dependent data-driven magnetofrictional model (TMFM; Pomoell et al. (2019)). We investigate the FR trajectories by tracking their apex throughout the full simulation time span. We demonstrate that this upgraded methodology provides the user with more tools and less a-priori assumptions about the FR shape that, in turn, leads to a more accurate set of field lines. The propagation analysis yields that the erupting FR from AR12473 showcases stronger dynamics than the AR11176 FR and a significant deflection during its ascent through the domain. The AR11176 FR appears more stable, though there still is a notable deflection. This confirms that at these low coronal heights, FRs do undergo significant changes in the direction of their propagation even for less dynamic cases. The modelling results are also verified with observations, with AR12473 being indeed dynamic and eruptive, while AR11176 only features an eruption outside of our simulation time window. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Inhomogeneous Cosmology using General Relativistic Smoothed Particle Hydrodynamics coupled to Numerical Relativity

Authors: Spencer J. Magnall, Daniel J. Price, Paul D. Lasky, Hayley J. Macpherson

Abstract: We perform three-dimensional simulations of homogeneous and inhomogeneous cosmologies via the coupling of a numerical relativity code for spacetime evolution and smoothed particle hydrodynamics (SPH) code. Evolution of a flat dust and radiation dominated Friedmann-Lemaître-Roberston-Walker (FLRW) spacetime shows an agreement of exact solutions with residuals on the order $10^{-6}$ and $10^{-3}$ re… ▽ More We perform three-dimensional simulations of homogeneous and inhomogeneous cosmologies via the coupling of a numerical relativity code for spacetime evolution and smoothed particle hydrodynamics (SPH) code. Evolution of a flat dust and radiation dominated Friedmann-Lemaître-Roberston-Walker (FLRW) spacetime shows an agreement of exact solutions with residuals on the order $10^{-6}$ and $10^{-3}$ respectively, even at low grid resolutions. We demonstrate evolution of linear perturbations of density, velocity and metric quantities to the FLRW with residuals of only $10^{-2}$ compared to exact solutions. Finally, we demonstrate the evolution of non-linear perturbations of the metric past shell-crossing, such that dark matter halo formation is possible. We show that numerical relativistic smoothed particle hydrodynamics is a viable method for understanding non-linear effects in cosmology. △ Less

Submitted 27 July, 2023; originally announced July 2023.

Comments: 15 pages, 12 figures, submitted to PRD

Tidal truncation of circumplanetary disks fails above a critical disk aspect ratio

Authors: Rebecca G. Martin, Philip J. Armitage, Stephen H. Lubow, Daniel J. Price

Abstract: We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodyn… ▽ More We use numerical simulations of circumplanetary disks to determine the boundary between disks that are radially truncated by the tidal potential, and those where gas escapes the Hill sphere. We consider a model problem, in which a coplanar circumplanetary disk is resupplied with gas at an injection radius smaller than the Hill radius. We evolve the disk using the PHANTOM Smoothed Particle Hydrodynamics code until a steady-state is reached. We find that the most significant dependence of the truncation boundary is on the disk aspect ratio $H/R$. Circumplanetary disks are efficiently truncated for $H/R \lesssim 0.2$. For $H/R \simeq 0.3$, up to about half of the injected mass, depending on the injection radius, flows outwards through the decretion disk and escapes. As expected from analytic arguments, the conditions ($H/R$ and Shakura-Sunyaev $α$) required for tidal truncation are independent of planet mass. A simulation with larger $α=0.1$ shows stronger outflow than one with $α=0.01$, but the dependence on transport efficiency is less important than variations of $H/R$. Our results suggest two distinct classes of circumplanetary disks: tidally truncated thin disks with dust-poor outer regions, and thicker actively decreting disks with enhanced dust-to-gas ratios. Applying our results to the PDS 70c system, we predict a largely truncated circumplanetary disk, but it is possible that enough mass escapes to support an outward flow of dust that could explain the observed disk size. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Comments: Accepted for publication in ApJ

Dust Formation in Common Envelope Binary Interaction -- I: 3D Simulations Using the Bowen Approximation

Authors: Miguel González-Bolívar, Luis C. Bermúdez-Bustamante, Orsola De Marco, Lionel Siess, Daniel J. Price, Mansi Kasliwal

Abstract: We carried out 3D smoothed particle hydrodynamics simulations of the common envelope binary interaction using the approximation of Bowen to calculate the dust opacity in order to investigate the resulting dust-driven accelerations. We have simulated two types of binary star: a 1.7 and a 3.7 $M_{\odot}$ thermally-pulsating, asymptotic giant branch stars with a 0.6 $M_{\odot}$ companion. We carried… ▽ More We carried out 3D smoothed particle hydrodynamics simulations of the common envelope binary interaction using the approximation of Bowen to calculate the dust opacity in order to investigate the resulting dust-driven accelerations. We have simulated two types of binary star: a 1.7 and a 3.7 $M_{\odot}$ thermally-pulsating, asymptotic giant branch stars with a 0.6 $M_{\odot}$ companion. We carried out simulations using both an ideal gas and a tabulated equations of state, with the latter considering the recombination energy of the envelope. We found that the dust-driven wind leads to a relatively small increase in the unbound gas, with the effect being smaller for the tabulated equation of state simulations and for the more massive primary. Dust acceleration does contribute to envelope expansion with only a slightly elongated morphology, if we believe the results from the tabulated equation of state as more reliable. The Bowen opacities in the outer envelopes of the two models, at late times, are large enough that the photosphere of the post-inspiral object is about ten times larger compared to the same without accounting for the dust opacities. As such, the prediction of the appearance of the transient would change substantially if dust is included. △ Less

Submitted 13 April, 2024; v1 submitted 28 June, 2023; originally announced June 2023.

Comments: 14 Pages, 13 Figures. Accepted for publication in MNRAS

An Integral-based Approach for the Vector Potential in Smoothed Particle Magnetohydrodynamics

Authors: Terrence S. Tricco, Daniel J. Price

Abstract: A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge choice. This new discretisation is tested using the Orszag-Tang MHD vortex in a 3D configuration. The corresponding conservative equations of motion are derive… ▽ More A new implementation for the time evolution of the magnetic vector potential is obtained for smoothed particle magnetohydrodynamics by considering the induction equation in integral form. Galilean invariance is achieved through proper gauge choice. This new discretisation is tested using the Orszag-Tang MHD vortex in a 3D configuration. The corresponding conservative equations of motion are derived, but are not found to solve the MHD equations in the continuum limit. Tests are performed using a hybrid approach instead, whereby the equations of motion based on the magnetic field instead of vector potential are used. Test results experience the same numerical instability as with the Price (2010) formulation. We conclude that this new formulation is non-viable. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: 7 pages, 5 figures, Proceedings of the "17th SPHERIC International Workshop", Rhodes, Greece, Jun 27-29, 2023

The Automatic Identification and Tracking of Coronal Flux Ropes -- Part I: Footpoints and Fluxes

Authors: Andreas Wagner, Emilia K. J. Kilpua, Ranadeep Sarkar, Daniel J. Price, Anshu Kumari, Farhad Daei, Jens Pomoell, Stefaan Poedts

Abstract: Investigating the early-stage evolution of an erupting flux rope from the Sun is important to understand the mechanisms of how it looses its stability and its space weather impacts. Our aim is to develop an efficient scheme for tracking the early dynamics of erupting solar flux ropes and use the algorithm to analyse its early-stage properties. The algorithm is tested on a data-driven simulation of… ▽ More Investigating the early-stage evolution of an erupting flux rope from the Sun is important to understand the mechanisms of how it looses its stability and its space weather impacts. Our aim is to develop an efficient scheme for tracking the early dynamics of erupting solar flux ropes and use the algorithm to analyse its early-stage properties. The algorithm is tested on a data-driven simulation of an eruption that took place in active region AR12473. We investigate the modelled flux rope's footpoint movement and magnetic flux evolution and compare with observational data from the Solar Dynamics Observatory's Atmospheric Imaging Assembly in the 211 $\unicode{x212B}$ and 1600 $\unicode{x212B}$ channels. To carry out our analysis, we use the time-dependent data-driven magnetofrictional model (TMFM). We also perform another modelling run, where we stop the driving of the TMFM midway through the flux rope's rise through the simulation domain and evolve it instead with a zero-beta magnetohydrodynamic (MHD) approach. The developed algorithm successfully extracts a flux rope and its ascend through the simulation domain. We find that the movement of the modelled flux rope footpoints showcases similar trends in both TMFM and relaxation MHD run: they recede from their respective central location as the eruption progresses and the positive polarity footpoint region exhibits a more dynamic behaviour. The ultraviolet brightenings and extreme ultraviolet dimmings agree well with the models in terms of their dynamics. According to our modelling results, the toroidal magnetic flux in the flux rope first rises and then decreases. In our observational analysis, we capture the descending phase of toroidal flux. In conclusion, the extraction algorithm enables us to effectively study the flux rope's early dynamics and derive some of its key properties such as footpoint movement and toroidal magnetic flux. △ Less

Submitted 26 June, 2023; originally announced June 2023.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 677, A81 (2023)

Observational Signatures of Circumbinary Discs I: Kinematics

Authors: Josh Calcino, Daniel J. Price, Christophe Pinte, Himanshi Garg, Brodie J. Norfolk, Valentin Christiaens, Hui Li, Richard Teague

Abstract: We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by i) a central cavity, ii) spiral arms both in and outside of their central cavities, iii) non-localised perturbations… ▽ More We present five morphological and kinematic criteria to aid in asserting the binary nature of a protoplanetary disc, based on 3D hydrodynamical simulations of circumbinary discs post-processed with Monte Carlo radiative transfer. We find that circumbinary discs may be identified by i) a central cavity, ii) spiral arms both in and outside of their central cavities, iii) non-localised perturbations in their iso-velocity curves, iv) asymmetry between the lines of maximum speed of the blue and red-shifted wings and v) asymmetry between the area of the blue and red-shifted wings. We provide quantitative metrics for the last two criteria that can be used, in conjunction with the morphological criteria, to signal whether a protoplanetary disc is likely to be a circumbinary disc. △ Less

Submitted 13 June, 2023; originally announced June 2023.

Comments: 26 pages, 26 figures, accepted to MNRAS

Effects of optimisation parameters on data-driven magnetofrictional modelling of active regions

Authors: A. Kumari, D. J. Price, F. Daei, J. Pomoell, E. K. J. Kilpua

Abstract: Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing th… ▽ More Data-driven time-dependent magnetofrictional modelling (TMFM) of active region magnetic fields has been proven to be a useful tool to study the corona. The input to the model is the photospheric electric field that is inverted from a time series of the photospheric magnetic field. Constraining the complete electric field, that is, including the non-inductive component, is critical for capturing the eruption dynamics. We present a detailed study of the effects of optimisation of the non-inductive electric field on the TMFM of AR12473. We aim to study the effects of varying the non-inductive electric field on the data-driven coronal simulations, for two alternative parametrisations. By varying parameters controlling the strength of the non-inductive electric field, we wish to explore the changes in flux rope formation and their early evolution and other parameters, for instance, axial flux and magnetic field magnitude.The non-inductive electric field component in the photosphere is critical for energising and introducing twist to the coronal magnetic field, thereby allowing unstable configurations to be formed. We estimated this component using an approach based on optimising the injection of magnetic energy. However, the flux rope formation, evolution and eruption time varies depending on the values of the optimisation parameters. The flux rope is formed and has overall similar evolution and properties with a large range of non-inductive electric fields needed to determine the non-inductive electric field component that is critical for energising and introducing twist to the coronal magnetic field. This study shows that irrespective of non-inductive electric field values, flux ropes are formed and erupted, which indicates that data-driven TMFM can be used to estimate flux rope properties early in their evolution without employing a lengthy optimisation process. △ Less

Submitted 25 May, 2023; originally announced May 2023.

Comments: 13 pages, 9 figures, 1 table. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 675, A80 (2023)

The role of the drag force in the gravitational stability of dusty planet-forming disc -- II. Numerical simulations

Authors: Cristiano Longarini, Philip J. Armitage, Giuseppe Lodato, Daniel J. Price, Simone Ceppi

Abstract: Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed p… ▽ More Young protostellar discs are likely to be both self-gravitating, and to support grain growth to sizes where the particles decoupled from the gas. This combination could lead to short-wavelength fragmentation of the solid component in otherwise non-fragmenting gas discs, forming Earth-mass solid cores during the Class 0/I stages of Young Stellar Object evolution. We use three-dimensional smoothed particle hydrodynamics simulations of two-fluid discs, in the regime where the Stokes number of the particles St>1, to study how the formation of solid clumps depends on the disc-to-star mass ratio, the strength of gravitational instability, and the Stokes number. Gravitational instability of the simulated discs is sustained by local cooling. We find that the ability of the spiral structures to concentrate solids increases with the cooling time, and decreases with the Stokes number, while the relative dynamical temperature between gas and dust of the particles decreases with the cooling time and the disc-to-star mass ratio, and increases with the Stokes number. Dust collapse occurs in a subset of high disc mass simulations, yielding clumps whose mass is close to linear theory estimates, namely 1-10 Earth masses. Our results suggest that if planet formation occurs via this mechanism, the best conditions correspond to near the end of the self-gravitating phase, when the cooling time is long and the Stokes number close to unity. △ Less

Submitted 5 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS, 20 pages

Confirmation and Keplerian motion of the gap-carving protoplanet HD 169142 b

Authors: Iain Hammond, Valentin Christiaens, Daniel J. Price, Claudia Toci, Christophe Pinte, Sandrine Juillard, Himanshi Garg

Abstract: We present the re-detection of a compact source in the face-on protoplanetary disc surrounding HD 169142, using VLT/SPHERE data in YJH bands. The source is found at a separation of 0.''319 ($\sim$37 au) from the star. Three lines of evidence argue in favour of the signal tracing a protoplanet: (i) it is found in the annular gap separating the two bright rings of the disc, as predicted by theory; (… ▽ More We present the re-detection of a compact source in the face-on protoplanetary disc surrounding HD 169142, using VLT/SPHERE data in YJH bands. The source is found at a separation of 0.''319 ($\sim$37 au) from the star. Three lines of evidence argue in favour of the signal tracing a protoplanet: (i) it is found in the annular gap separating the two bright rings of the disc, as predicted by theory; (ii) it is moving at the expected Keplerian velocity for an object at $\sim$37 au in the 2015, 2017 and 2019 datasets; (iii) we also detect a spiral-shaped signal whose morphology is consistent with the expected outer spiral wake triggered by a planet in the gap, based on dedicated hydrodynamical simulations of the system. The YJH colours we extracted for the object are consistent with tracing scattered starlight, suggesting that the protoplanet is enshrouded in a significant amount of dust, as expected for a circumplanetary disc or envelope surrounding a gap-clearing Jovian-mass protoplanet. △ Less

Submitted 23 February, 2023; v1 submitted 22 February, 2023; originally announced February 2023.

Comments: Accepted for publication in MNRAS Letters. 5 pages, 5 figures

Kinematic and thermal signatures of the directly imaged protoplanet candidate around Elias 2-24

Authors: Christophe Pinte, Iain Hammond, Daniel J. Price, Valentin Christiaens, Sean M. Andrews, Gaël Chauvin, Laura M. Pérez, Sebastián Jorquera, Himanshi Garg, Brodie J. Norfolk, Josh Calcino, Mickaël Bonnefoy

Abstract: We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require loc… ▽ More We report kinematic and thermal signatures associated with the directly imaged protoplanet candidate in the Elias 2-24 disc. Using the DSHARP ALMA observations of the $^{12}$CO J=2-1 line, we show that the disc kinematics are perturbed, with a detached CO emission spot at the location of the planet candidate and traces of spiral wakes, and also that the observed CO emission intensities require local heating. While the foreground extinction hides the velocity channels associated with the planet, preventing a planet mass estimate, the level of gas heating implied by the CO emission indicates the presence of a warm, embedded giant planet. Comparison with models show this could either be a $\gtrsim 5$M$_\mathrm{Jup}$, or a lower mass ( $\gtrsim 2$M$_\mathrm{Jup}$) but accreting proto-planet. △ Less

Submitted 20 January, 2023; originally announced January 2023.

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters

Dynamics of dust grains in turbulent molecular clouds. Conditions for decoupling and limits of different numerical implementations

Authors: Benoît Commerçon, Ugo Lebreuilly, Daniel J. Price, Francesco Lovascio, Guillaume Laibe, Patrick Hennebelle

Abstract: Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of… ▽ More Dust grain dynamics in molecular clouds is regulated by its interplay with supersonic turbulent gas motions. The conditions under which dust grains decouple from the dynamics of gas remain poorly constrained. We first aim to investigate the critical dust grain size for dynamical decoupling, using both analytical predictions and numerical experiments. Second, we aim to set the range of validity of two fundamentally different numerical implementations for the evolution of dust and gas mixtures in turbulent molecular clouds. We carried out a suite of numerical experiments using two different schemes. First, we used a monofluid formalism in the terminal velocity approximation (TVA) on a Eulerian grid. Second, we used a two-fluid scheme, in which the dust dynamics is handled with Lagrangian super-particles, and the gas dynamics on a Eulerian grid. The monofluid results are in good agreement with the theoretical critical size for decoupling. We report dust dynamics decoupling for Stokes number St>0.1, that is, dust grains of $s>4~μ$m in size. We find that the TVA is well suited for grain sizes of 10 $μ$m in molecular clouds, in particular in the densest regions. However, the maximum dust enrichment measured in the low-density material where St>1 is questionable. In the Lagrangian dust experiments, we show that the results are affected by the numerics for all dust grain sizes. At St<<1, the dust dynamics is largely affected by artificial trapping in the high-density regions, leading to spurious variations of the dust concentration. At St>1, the maximum dust enrichment is regulated by the grid resolution used for the gas dynamics. The results of previous similar numerical work should therefore be revisited with respect to the limitations we highlight in this study. Dust enrichment of submicron dust grains is unlikely to occur in the densest parts of molecular clouds. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 22 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 671, A128 (2023)

Hot Jupiter engulfment by a red giant in 3D hydrodynamics

Authors: Mike Y. M. Lau, Matteo Cantiello, Adam S. Jermyn, Morgan MacLeod, Ilya Mandel, Daniel J. Price

Abstract: Transit and radial-velocity surveys over the past two decades have uncovered a significant population of short-period exoplanets. Among them are hot Jupiters, which are gas giant planets with orbital periods of a few days and found in 0.1-1% of Sun-like stars. Hot Jupiters are expected to be engulfed during their host star's radial expansion on the red giant branch. Planetary engulfment has been s… ▽ More Transit and radial-velocity surveys over the past two decades have uncovered a significant population of short-period exoplanets. Among them are hot Jupiters, which are gas giant planets with orbital periods of a few days and found in 0.1-1% of Sun-like stars. Hot Jupiters are expected to be engulfed during their host star's radial expansion on the red giant branch. Planetary engulfment has been studied extensively as it may account for observed rapidly rotating and chemically enriched giant stars. We perform 3D hydrodynamical simulations of hot Jupiter engulfment by a 1 solar mass, 4 solar radii early red giant. Our "global" simulations simultaneously resolve the stellar envelope and planetary structure, modelling the hot Jupiter as a polytropic gas sphere. We find that approximately 90% of the hot Jupiter's mass is ablated in the convective part of the giant envelope, which would enhance the surface lithium abundance by 0.1 dex. The hot Jupiter is disrupted by a combination of ram pressure and tidal forces near the base of the convective envelope, with the deepest material penetrating to the radiative zone. The star experiences modest spin-up (~1 km/s), although engulfing a more massive companion could produce a rapidly rotating giant. Drag heating near the surface could exceed the unperturbed stellar luminosity and power an optical transient. For the amount of unbound ejecta recorded in the simulation, H-recombination could also power a transient that is around ten times the pre-engulfment luminosity, for several days. △ Less

Submitted 27 October, 2022; originally announced October 2022.

Comments: 14 pages, 11 figures

Solar Energetic Particle Time Series Analysis with Python

Authors: Christian Palmroos, Jan Gieseler, Nina Dresing, Diana E. Morosan, Eleanna Asvestari, Aleksi Yli-Laurila, Daniel J. Price, Saku Valkila, Rami Vainio

Abstract: Solar Energetic Particles (SEPs) are charged particles accelerated within the solar atmosphere or the interplanetary space by explosive phenomena such as solar flares or Coronal Mass Ejections (CMEs). Once injected into the interplanetary space, they can propagate towards Earth, causing space weather related phenomena. For their analysis, interplanetary in-situ measurements of charged particles ar… ▽ More Solar Energetic Particles (SEPs) are charged particles accelerated within the solar atmosphere or the interplanetary space by explosive phenomena such as solar flares or Coronal Mass Ejections (CMEs). Once injected into the interplanetary space, they can propagate towards Earth, causing space weather related phenomena. For their analysis, interplanetary in-situ measurements of charged particles are key. The recently expanded spacecraft fleet in the heliosphere not only provides much-needed additional vantage points, but also increases the variety of missions and instruments for which data loading and processing tools are needed. This manuscript introduces a series of Python functions that will enable the scientific community to download, load, and visualize charged particle measurements of the current space missions that are especially relevant to particle research as time series or dynamic spectra. In addition, further analytical functionality is provided that allows the determination of SEP onset times as well as their inferred injection times. The full workflow, which is intended to be run within Jupyter Notebooks and can also be approachable for Python laymen, will be presented with scientific examples. All functions are written in Python, with the source code publicly available at GitHub under a permissive license. Where appropriate, available Python libraries are used, and their application is described. △ Less

Submitted 14 December, 2022; v1 submitted 19 October, 2022; originally announced October 2022.

Journal ref: Front. Astron. Space Sci. 9:1073578 (2022)

Sustained FU Orionis-type outbursts from colliding discs in stellar flybys

Authors: Elisabeth M. A. Borchert, Daniel J. Price, Christophe Pinte, Nicolás Cuello

Abstract: We perform 3D hydrodynamics simulations of disc-disc stellar flybys with on-the-fly Monte Carlo radiative transfer. We show that pre-existing circumstellar discs around both stars result in fast rising ($\sim$yrs) outbursts lasting 2-5 times longer than for a star-disc flyby. The perturber always goes into outburst ($\dot{M}>10^{-5}~{\rm M_{\odot}~ yr^{-1}}$). Whereas we find that the primary goes… ▽ More We perform 3D hydrodynamics simulations of disc-disc stellar flybys with on-the-fly Monte Carlo radiative transfer. We show that pre-existing circumstellar discs around both stars result in fast rising ($\sim$yrs) outbursts lasting 2-5 times longer than for a star-disc flyby. The perturber always goes into outburst ($\dot{M}>10^{-5}~{\rm M_{\odot}~ yr^{-1}}$). Whereas we find that the primary goes into a decades long outburst only when the flyby is retrograde to the circumprimary disc rotation. High accretion rates during the outburst are triggered by angular momentum cancellation in misaligned material generated by the encounter. A large fraction of accreted material is alien. △ Less

Submitted 3 October, 2022; originally announced October 2022.

Comments: 11 pages, 9 figures, accepted for publication in MNRAS

Solar-MACH: An open-source tool to analyze solar magnetic connection configurations

Authors: Jan Gieseler, Nina Dresing, Christian Palmroos, Johan L. Freiherr von Forstner, Daniel J. Price, Rami Vainio, Athanasios Kouloumvakos, Laura Rodríguez-García, Domenico Trotta, Vincent Génot, Arnaud Masson, Markus Roth, Astrid Veronig

Abstract: The Solar MAgnetic Connection HAUS tool (Solar-MACH) is an open-source tool completely written in Python that derives and visualizes the spatial configuration and solar magnetic connection of different observers (i.e., spacecraft or planets) in the heliosphere at different times. For doing this, the magnetic connection in the interplanetary space is obtained by the classic Parker Heliospheric Magn… ▽ More The Solar MAgnetic Connection HAUS tool (Solar-MACH) is an open-source tool completely written in Python that derives and visualizes the spatial configuration and solar magnetic connection of different observers (i.e., spacecraft or planets) in the heliosphere at different times. For doing this, the magnetic connection in the interplanetary space is obtained by the classic Parker Heliospheric Magnetic Field (HMF). In close vicinity of the Sun, a Potential Field Source Surface (PFSS) model can be applied to connect the HMF to the solar photosphere. Solar-MACH is especially aimed at providing publication-ready figures for the analyses of Solar Energetic Particle events (SEPs) or solar transients such as Coronal Mass Ejections (CMEs). It is provided as an installable Python package (listed on PyPI and conda-forge), but also as a web tool at solar-mach.github.io that completely runs in any web browser and requires neither Python knowledge nor installation. The development of Solar-MACH is open to everyone and takes place on GitHub, where the source code is publicly available under the BSD 3-Clause License. Established Python libraries like sunpy and pfsspy are utilized to obtain functionalities when possible. In this article, the Python code of Solar-MACH is explained, and its functionality is demonstrated using real science examples. In addition, we introduce the overarching SERPENTINE project, the umbrella under which the recent development took place. △ Less

Submitted 12 January, 2023; v1 submitted 3 October, 2022; originally announced October 2022.

Comments: 14 pages, 5 figures

Journal ref: Front. Astron. Space Sci. 9:1058810 (2023)

3D simulations of AGB stellar winds -- I. Steady winds and dust formation

Authors: L. Siess, W. Homan, S. Toupin, D. J. Price

Abstract: Aims. We present the implementation of the treatment of particle ejection and dust nucleation in the smoothed particle hydrodynamics (SPH) code phantom. These developments represent the first step toward a more complete modeling of dust-driven winds emanating from AGB stars. Methods. The AGB outflow is modeled by injecting the SPH particles from a spherical inner boundary. This boundary is a serie… ▽ More Aims. We present the implementation of the treatment of particle ejection and dust nucleation in the smoothed particle hydrodynamics (SPH) code phantom. These developments represent the first step toward a more complete modeling of dust-driven winds emanating from AGB stars. Methods. The AGB outflow is modeled by injecting the SPH particles from a spherical inner boundary. This boundary is a series of concentric shells, with the AGB star at its center, and the particles are positioned on these shells on the vertices of an isocahedron geodesic surface. The outermost shell is ejected with a predefined radial velocity, and subsequent lower shells replenish the ejected ones, all rotated randomly to improve the isotropy of the outflow. The physical properties of the particles on these shells are set by solving the 1D analytic steady wind equations. The formation of dust is calculated starting from a compact chemical network for carbon-rich material, which creates the building blocks of the solid-state particles. Subsequently, the theory of the moments is used to obtain dust growth rates, without requiring knowledge on the grain size distribution. Results. We tested our implementation against a series of 1D reference solutions. We demonstrate that our method is able to reproduce Parker-type wind solutions. For the trans-sonic solution, small oscillations are present in the vicinity of the sonic point, but these do not impact the trans-sonic passage or terminal wind velocity. Supersonic solutions always compare nicely with 1D analytic profiles. We also tested our implementation of dust using two formalisms: an analytic prescription for the opacity devised by Bowen and the full treatment of carbon-dust formation. Both simulations reproduce the 1D analytic solution displaying the expected additional acceleration when the gas temperature falls below the condensation temperature. △ Less

Submitted 29 August, 2022; originally announced August 2022.

Comments: 16 pages, 12 figures

Journal ref: A&A 667, A75 (2022)

The Origin of the Doppler-flip in HD 100546: a large scale spiral arm generated by an inner binary companion

Authors: Brodie J. Norfolk, Christophe Pinte, Josh Calcino, Iain Hammond, Nienke van der Marel, Daniel J. Price, Sarah T. Maddison, Valentin Christiaens, Jean-Francois Gonzalez, Dori Blakely, Giovanni Rosotti, Christian Ginski

Abstract: Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner sp… ▽ More Companions at sub-arcsecond separation from young stars are difficult to image. However their presence can be inferred from the perturbations they create in the dust and gas of protoplanetary disks. Here we present a new interpretation of SPHERE polarised observations that reveal the previously detected inner spiral in the disk of HD 100546. The spiral coincides with a newly detected 12CO inner spiral and the previously reported CO emission Doppler-flip, which has been interpreted as the signature of an embedded protoplanet. Comparisons with hydrodynamical models indicate that this Doppler-flip is instead the kinematic counterpart of the spiral, which is likely generated by an inner companion inside the disk cavity. △ Less

Submitted 4 August, 2022; originally announced August 2022.

Comments: Accepted to ApJL, 14 pages, 10 figures

Close encounters: How stellar flybys shape planet-forming discs

Authors: Nicolás Cuello, François Ménard, Daniel J. Price

Abstract: We review the role of stellar flybys and encounters in shaping planet-forming discs around young stars, based on the published literature on this topic in the last 30 years. Since most stars $\leq~2$ Myr old harbour protoplanetary discs, tidal perturbations affect planet formation. First, we examine the probability of experiencing flybys or encounters: More than 50\% of stars with planet-forming d… ▽ More We review the role of stellar flybys and encounters in shaping planet-forming discs around young stars, based on the published literature on this topic in the last 30 years. Since most stars $\leq~2$ Myr old harbour protoplanetary discs, tidal perturbations affect planet formation. First, we examine the probability of experiencing flybys or encounters: More than 50\% of stars with planet-forming discs in a typical star forming environment should experience a close stellar encounter or flyby within 1000 au. Second, we detail the dynamical effects of flybys on planet-forming discs. Prograde, parabolic, disc-penetrating flybys are the most destructive. Grazing and penetrating flybys in particular lead to the capture of disc material by the secondary to form a highly misaligned circumsecondary disc with respect to the disc around the primary. One or both discs may undergo extreme accretion and outburst events, similar to the ones observed in FU Orionis-type stars. Warps and broken discs are distinct signatures of retrograde flybys. Third, we review some recently observed stellar systems with discs where a stellar flyby or an encounter is suspected -- including UX Tau, RW Aur, AS 205, Z CMa, and FU Ori. Finally, we discuss the implications of stellar flybys for planet formation and exoplanet demographics, including possible imprints of a flyby in the Solar System in the orbits of trans-Neptunian objects and the Sun's obliquity. △ Less

Submitted 22 December, 2022; v1 submitted 20 July, 2022; originally announced July 2022.

Comments: 16 pages, 7 figures, 1 table, 255 references. Invited review article accepted in EPJ+

External or internal companion exciting the spiral arms in CQ Tau?

Authors: Iain Hammond, Valentin Christiaens, Daniel J. Price, Maria Giulia Ubeira-Gabellini, Jennifer Baird, Josh Calcino, Myriam Benisty, Giuseppe Lodato, Leonardo Testi, Christophe Pinte, Claudia Toci, Davide Fedele

Abstract: We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M… ▽ More We present new high-contrast images in near-infrared wavelengths ($λ$ = 1.04, 1.24, 1.62, 2.18 and 3.78$μ$m) of the young variable star CQ Tau, aiming to constrain the presence of companions in the protoplanetary disc. We reached a Ks-band contrast of 14 magnitudes with SPHERE/IRDIS at separations greater than 0."4 from the star. Our mass sensitivity curve rules out giant planets above 4 M$_{\rm Jup}$ immediately outside the spiral arms at $\sim$60 au and above 2-3 M$_{\rm Jup}$ beyond 100 au to 5$σ$ confidence assuming hot-start models. We do, however, detect four spiral arms, a double-arc and evidence for shadows in scattered light cast by a misaligned inner disc. Our observations may be explained by an unseen close-in companion on an inclined and eccentric orbit. Such a hypothesis would also account for the disc CO cavity and disturbed kinematics. △ Less

Submitted 31 July, 2022; v1 submitted 18 July, 2022; originally announced July 2022.

Comments: 10 figures, 14 pages. Accepted for publication in MNRAS on 25 July 2022

Report number: MN-22-0627-MJ.R1

Kinematic evidence for an embedded planet in the IM Lupi disc

Authors: Harrison J. Verrios, Daniel J. Price, Christophe Pinte, Thomas Hilder, Josh Calcino

Abstract: We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disc can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2-3 times th… ▽ More We test the hypothesis that an embedded giant planet in the IM Lupi protostellar disc can produce velocity kinks seen in CO line observations as well as the spiral arms seen in scattered light and continuum emission. We inject planets into 3D hydrodynamics simulations of IM Lupi, generating synthetic observations using Monte Carlo radiative transfer. We find that an embedded planet of 2-3 times the mass of Jupiter can reproduce non-Keplerian velocity perturbations, or `kinks', in the 12CO J=2-1 channel maps. Such a planet can also explain the spiral arms seen in 1.25mm dust continuum emission and 1.6 micron scattered light images. We show that the wake of the planet can be traced in the observed peak velocity map, which appears to closely follow the morphology expected from our simulations and from analytic models of planet-disc interaction. △ Less

Submitted 6 July, 2022; originally announced July 2022.

Comments: 10 pages, 6 figures, accepted to ApJL. Radiative transfer models and a movie showing change with azimuth available from https://doi.org/10.26180/20145620.v3

Common envelopes in massive stars II: The distinct roles of hydrogen and helium recombination

Authors: Mike Y. M. Lau, Ryosuke Hirai, Daniel J. Price, Ilya Mandel

Abstract: The role of recombination during a common-envelope event has been long debated. Many studies have argued that much of hydrogen recombination energy, which is radiated in relatively cool and optically-thin layers, might not thermalise in the envelope. On the other hand, helium recombination contains 30% of the total recombination energy, and occurs much deeper in the stellar envelope. We investigat… ▽ More The role of recombination during a common-envelope event has been long debated. Many studies have argued that much of hydrogen recombination energy, which is radiated in relatively cool and optically-thin layers, might not thermalise in the envelope. On the other hand, helium recombination contains 30% of the total recombination energy, and occurs much deeper in the stellar envelope. We investigate the distinct roles played by hydrogen and helium recombination in a common-envelope interaction experienced by a 12 solar mass red supergiant donor. We perform adiabatic, 3D hydrodynamical simulations that (i) include hydrogen, helium, and molecular hydrogen recombination, (ii) include hydrogen and helium recombination, (iii) include only helium recombination, and (iv) do not include recombination energy. By comparing these simulations, we find that the addition of helium recombination energy alone ejects 30% more envelope mass, and leads to a 16% larger post-plunge-in separation. Under the adiabatic assumption, adding hydrogen recombination energy increases the amount of ejected mass by a further 40%, possibly unbinding the entire envelope, but does not affect the post-plunge separation. Most of the ejecta becomes unbound at relatively high (>70%) degrees of hydrogen ionisation, where the hydrogen recombination energy is likely to expand the envelope instead of being radiated away. △ Less

Submitted 31 August, 2022; v1 submitted 13 June, 2022; originally announced June 2022.

Comments: 10 pages and 6 figures. Accepted for publication in MNRAS

Common envelope binary interaction simulations between a thermally-pulsating AGB star and a low mass companion

Authors: Miguel Gonzalez-Bolivar, Orsola De Marco, Mike Y. M. Lau, Ryosuke Hirai, Daniel J. Price

Abstract: At least one in five of all planetary nebulae are the product of a common envelope (CE) interaction, where the companion in-spirals into the envelope of an asymptotic giant branch (AGB) star ejecting the nebula and leaving behind a compact binary. In this work we carry out 3D smoothed particle hydrodynamics simulations of the CE interaction between a 1.7 $M_{\odot}$ AGB star and a 0.6 $M_{\odot}$… ▽ More At least one in five of all planetary nebulae are the product of a common envelope (CE) interaction, where the companion in-spirals into the envelope of an asymptotic giant branch (AGB) star ejecting the nebula and leaving behind a compact binary. In this work we carry out 3D smoothed particle hydrodynamics simulations of the CE interaction between a 1.7 $M_{\odot}$ AGB star and a 0.6 $M_{\odot}$ companion. We model the AGB structure using a 1D stellar model taken at the seventh thermal pulse. The interaction takes place when the giant is on the expanding phase of the seventh thermal pulse and has a radius of 250 $R_{\odot}$. The post-CE orbital separations varies between 20 and 31 $R_{\odot}$, with the inclusion of recombination energy resulting in wider separations. Based on the observed short in-spiral time-scales, we suggest that thermal pulses can trigger CEs, extending the ability of AGB stars to capture companions into CEs, that would lead to the prediction of a larger population of post-AGB, post-CE binaries. Simulations that include a tabulated equation of state unbind a great deal more gas, likely unbinding the entire envelope on short time-scales. The shape of the CE after the in-spiral is more spherical for AGB than red giant branch stars, and even more so if recombination energy is included. We expect that the planetary nebula formed from this CE will have different features from those predicted by Zou et al. 2020. △ Less

Submitted 24 October, 2022; v1 submitted 19 May, 2022; originally announced May 2022.

Comments: 17 pages excluding appendices, 17 figures, Accepted on MNRAS

Journal ref: Monthly Notices of the Royal Astronomical Society, Volume 517, Issue 3, December 2022, Pages 3181-3199

Accretion rates in hierarchical triple systems with discs

Authors: Simone Ceppi, Nicolás Cuello, Giuseppe Lodato, Cathie Clarke, Claudia Toci, Daniel J. Price

Abstract: Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical… ▽ More Young multiple systems accrete most of their final mass in the first few Myr of their lifetime, during the protostellar and protoplanetary phases. Previous studies showed that in binary systems the majority of the accreted mass falls onto the lighter star, thus evolving to mass equalisation. However, young stellar systems often comprise more than two stars, which are expected to be in hierarchical configurations. Despite its astrophysical relevance, differential accretion in hierarchical systems remains to be understood. In this work, we investigate whether the accretion trends expected in binaries are valid for higher order multiples. We performed a set of 3D Smoothed Particle Hydrodynamics simulations of binaries and of hierarchical triples (HTs) embedded in an accretion disc, with the code Phantom. We identify for the first time accretion trends in HTs and their deviations compared to binaries. These deviations, due to the interaction of the small binary with the infalling material from the circum-triple disc, can be described with a semi-analytical prescription. Generally, the smaller binary of a HT accretes more mass than a single star of the same mass as the smaller binary. We found that in a HT, if the small binary is heavier than the third body, the standard differential accretion scenario (whereby the secondary accretes more of the mass) is hampered. Reciprocally, if the small binary is lighter than the third body, the standard differential accretion scenario is enhanced. The peculiar differential accretion mechanism we find in HTs is expected to affect their mass ratio distribution. △ Less

Submitted 18 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS, 15 pages, 13 figures, 1 tables, 2 appendices

On the origin of magnetic fields in stars II: The effect of numerical resolution

Authors: James Wurster, Matthew R. Bate, Daniel J. Price, Ian A. Bonnell

Abstract: Are the kG-strength magnetic fields observed in young stars a fossil field left over from their formation or are they generated by a dynamo? Our previous numerical study concluded that magnetic fields must originate by a dynamo process. Here, we continue that investigation by performing even higher numerical resolution calculations of the gravitational collapse of a 1~M$_\odot$ rotating, magnetise… ▽ More Are the kG-strength magnetic fields observed in young stars a fossil field left over from their formation or are they generated by a dynamo? Our previous numerical study concluded that magnetic fields must originate by a dynamo process. Here, we continue that investigation by performing even higher numerical resolution calculations of the gravitational collapse of a 1~M$_\odot$ rotating, magnetised molecular cloud core through the first and second collapse phases until stellar densities are reached. Each model includes Ohmic resistivity, ambipolar diffusion, and the Hall effect. We test six numerical resolutions, using between $10^5$ and $3\times10^7$ particles to model the cloud. At all but the lowest resolutions, magnetic walls form in the outer parts of the first hydrostatic core, with the maximum magnetic field strength located within the wall rather than at the centre of the core. At high resolution, this magnetic wall is disrupted by the Hall effect, producing a magnetic field with a spiral-shaped distribution of intensity. As the second collapse occurs, this field is dragged inward and grows in strength, with the maximum field strength increasing with resolution. As the second core forms, the maximum field strength exceeds 1~kG in our highest resolution simulations, and the stellar core field strength exceeds this threshold at the highest resolution. Our resolution study suggests that kG-strength magnetic fields may be implanted in low-mass stars during their formation, and may persist over long timescales given that the diffusion timescale for the magnetic field exceeds the age of the Universe. △ Less

Submitted 18 January, 2022; originally announced January 2022.

Comments: 19 pages, 17 figures, accepted for publication in MNRAS

Multi-spacecraft observations of the structure of the sheath of an interplanetary coronal mass ejection and related energetic ion enhancement

Authors: E. K. J. Kilpua, S. W. Good, N. Dresing, R. Vainio, E. E. Davies, R. J. Forsyth, J. Gieseler, B. Lavraud, E. Asvestari, D. E. Morosan, J. Pomoell, D. J. Price, D. Heyner, T. S. Horbury, V. Angelini, H. O'Brien, V. Evans, J. Rodriguez-Pacheco, R. Gómez Herrero, G. C. Ho, R. Wimmer-Schweingruber

Abstract: Sheaths ahead of coronal mass ejections (CMEs) are large heliospheric structures that form with CME expansion and propagation. Turbulent and compressed sheaths contribute to the acceleration of particles in the corona and in interplanetary space, but the relation of their internal structures to particle energization is still relatively little studied. In particular, the role of sheaths in accelera… ▽ More Sheaths ahead of coronal mass ejections (CMEs) are large heliospheric structures that form with CME expansion and propagation. Turbulent and compressed sheaths contribute to the acceleration of particles in the corona and in interplanetary space, but the relation of their internal structures to particle energization is still relatively little studied. In particular, the role of sheaths in accelerating particles when the shock Mach number is low is a significant open problem. This work seeks to provide new insights on the internal structure of CME sheaths with regard to energetic particle enhancements. A good opportunity to achieve this aim was provided by observations of a sheath made by radially aligned spacecraft at 0.8 and $\sim$ 1 AU (Solar Orbiter, Wind, ACE and BepiColombo) on 19-21 April 2020. The sheath was preceded by a weak shock. Energetic ion enhancements occurred at different locations within the sheath structure at Solar Orbiter and L1. Magnetic fluctuation amplitudes at inertial-range scales increased in the sheath relative to the upstream wind. However, when normalised to the local mean field, fluctuation amplitudes did not increase significantly; magnetic compressibility of fluctuation also did not increase. Various substructures were embedded within the sheath at the different spacecraft, including multiple heliospheric current sheet (HCS) crossings and a small-scale flux rope. At L1, the ion flux enhancement was associated with the HCS crossings, while at Solar Orbiter, the enhancement occurred within the rope. Substructures that are swept from the upstream solar wind and compressed in the sheath can act as particularly effective acceleration sites. A possible acceleration mechanism is betatron acceleration associated with the small-scale flux rope and the warped HCS in the sheath. △ Less

Submitted 17 December, 2021; originally announced December 2021.

Comments: 14 pages, 12 figures; published in Astronomy & Astrophysics, Solar Orbiter First Results (Cruise Phase) special issue

The protoplanetary disc around HD 169142: circumstellar or circumbinary?

Authors: P. P. Poblete, N. Cuello, S. Pérez, S. Marino, J. Calcino, E. Macías, Á. Ribas, A. Zurlo, J. Cuadra, M. Montesinos, S. Zúñiga-Fernández, A. Bayo, C. Pinte, F. Ménard, D. J. Price

Abstract: Stellar binaries represent a substantial fraction of stellar systems, especially among young stellar objects. Accordingly, binaries play an important role in setting the architecture of a large number of protoplanetary discs. Binaries in coplanar and polar orientations with respect to the circumbinary disc are stable configurations and could induce non-axisymmetric structures in the dust and gas d… ▽ More Stellar binaries represent a substantial fraction of stellar systems, especially among young stellar objects. Accordingly, binaries play an important role in setting the architecture of a large number of protoplanetary discs. Binaries in coplanar and polar orientations with respect to the circumbinary disc are stable configurations and could induce non-axisymmetric structures in the dust and gas distributions. In this work, we suggest that the structures shown in the central region of the protoplanetary disc HD 169142 are produced by the presence of an inner stellar binary and a circumbinary (P-type) planet. We find that a companion with a mass-ratio of 0.1, semi-major axis of 9.9 au, eccentricity of 0.2, and inclination of 90°, together with a 2 Jupiter Mass coplanar planet on a circular orbit at 45 au reproduce the structures at the innermost ring observed at 1.3 mm and the shape of spiral features in scattered light observations. The model predicts changes in the disc's dust structure, and star's astrometric parameters, which would allow testing its veracity by monitoring this system over the next 20 years. △ Less

Submitted 14 December, 2021; v1 submitted 26 November, 2021; originally announced November 2021.

Comments: 11 pages, 7 figures, accepted for publication in MNRAS

On the rise times in FU Orionis events

Authors: Elisabeth M. A. Borchert, Daniel J. Price, Christophe Pinte, Nicolás Cuello

Abstract: We examine whether stellar flyby simulations can initiate FU Orionis outbursts using 3D hydrodynamics simulations coupled to live Monte Carlo radiative transfer. We find that a flyby where the secondary penetrates the circumprimary disc triggers a 1-2 year rise in the mass accretion rate to $10^{-4}~{\rm M_{\odot}~ yr^{-1}}$ that remains high ($\gtrsim 10^{-5}~{\rm M_{\odot}~yr^{-1}}$) for more th… ▽ More We examine whether stellar flyby simulations can initiate FU Orionis outbursts using 3D hydrodynamics simulations coupled to live Monte Carlo radiative transfer. We find that a flyby where the secondary penetrates the circumprimary disc triggers a 1-2 year rise in the mass accretion rate to $10^{-4}~{\rm M_{\odot}~ yr^{-1}}$ that remains high ($\gtrsim 10^{-5}~{\rm M_{\odot}~yr^{-1}}$) for more than a hundred years, similar to the outburst observed in FU Ori. Importantly, we find that the less massive star becomes the dominant accretor, as observed. △ Less

Submitted 13 December, 2021; v1 submitted 24 November, 2021; originally announced November 2021.

Comments: 5 pages, 4 figures, accepted for publication in MNRAS, movies available at https://emborchert.github.io/videos, press release available at https://www.monash.edu/science/news/current/astronomers-find-clue-to-solar-system-formation-through-well-known-star

Mapping the Planetary Wake in HD 163296 with Kinematics

Authors: Josh Calcino, Thomas Hilder, Daniel J. Price, Christophe Pinte, Francesco Bollati, Giuseppe Lodato, Brodie J. Norfolk

Abstract: We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD~163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk major axis confirm that the wake induces predominantly radial motion, as predicted by models of planet-disk interaction. These resu… ▽ More We map the planetary wake associated with the embedded protoplanet creating the CO kink in the disk of HD~163296. We show that the wake can be traced by a series of correlated perturbations in the peak velocity map. The sign change of the perturbations across the disk major axis confirm that the wake induces predominantly radial motion, as predicted by models of planet-disk interaction. These results provide the first direct confirmation of planet wakes generated by Lindblad resonances. Mapping the wake provides a constraint on the disk aspect ratio, which is required to measure the mass of the planet. △ Less

Submitted 6 April, 2022; v1 submitted 14 November, 2021; originally announced November 2021.

Comments: 10 Pages, 6 figures. Accepted in ApJL

Gravitational waves from tidal disruption events: an open and comprehensive catalogue

Authors: Martina Toscani, Giuseppe Lodato, Daniel J. Price, David Liptai

Abstract: We present an online, open and comprehensive template library of gravitational waveforms produced during the tidal disruptions of stars by massive black holes, spanning a broad space of parameters. We build this library thanks to a new feature that we implement in the general relativistic version of PHANTOM, a smoothed particle hydrodynamics code for three dimensional simulations in general relati… ▽ More We present an online, open and comprehensive template library of gravitational waveforms produced during the tidal disruptions of stars by massive black holes, spanning a broad space of parameters. We build this library thanks to a new feature that we implement in the general relativistic version of PHANTOM, a smoothed particle hydrodynamics code for three dimensional simulations in general relativity. We first perform a series of numerical tests to show that the gravitational wave (GW) signal obtained is in excellent agreement with the one expected from theory. This benchmark is done for well studied scenarios (such as binary stellar systems). We then apply our code to calculate the GW signals from tidal disruption events (TDEs), finding that our results are consistent with the theoretical estimates obtained in previous studies for selected parameters. We illustrate interesting results from the catalogue, where we stress how the gravitational signal is affected by variations of some parameters (like black hole spin, stellar orbital eccentricity and inclination). The full catalogue is available online. It is intended to be a living catalogue. △ Less

Submitted 9 November, 2021; originally announced November 2021.

Comments: Accepted for publication on MNRAS

Common envelopes in massive stars: Towards the role of radiation pressure and recombination energy in ejecting red supergiant envelopes

Authors: Mike Y. M. Lau, Ryosuke Hirai, Miguel González-Bolívar, Daniel J. Price, Orsola De Marco, Ilya Mandel

Abstract: We perform 3D hydrodynamical simulations of a common-envelope event involving a 12 solar mass red supergiant donor. Massive stars are expected to be qualitatively different from low-mass stars as their envelopes have significant support from radiation pressure, which increases both the final separation and amount of mass ejected through the common-envelope interaction. We perform adiabatic simulat… ▽ More We perform 3D hydrodynamical simulations of a common-envelope event involving a 12 solar mass red supergiant donor. Massive stars are expected to be qualitatively different from low-mass stars as their envelopes have significant support from radiation pressure, which increases both the final separation and amount of mass ejected through the common-envelope interaction. We perform adiabatic simulations that include radiation energy through the equation of state, which results in ejecting 60 per cent more mass (up to two thirds of the total envelope mass becoming unbound, or more) and yield a 10 per cent larger final separation compared to simulations that assume an ideal gas. When also including recombination energy, we find that at least three quarters of the envelope, and possibly the entire envelope, may be unbound. The final separation further increases by almost 20 per cent. The additional amount of ejected material is mainly due to energy injected from helium recombination. Hydrogen recombination plays a comparatively small role, as it mainly occurs in gas that has already become unbound. We conclude that the internal energy of the envelope can be a significant energy source for ejecting the common envelope, but ultimately radiation transport and convection need to be included. △ Less

Submitted 12 January, 2022; v1 submitted 1 November, 2021; originally announced November 2021.

Comments: Accepted for publication in MNRAS, 17 pages excluding appendices, 15 figures

Dust Traffic Jams in Inclined Circumbinary Protoplanetary Discs I. Morphology and Formation Theory

Authors: Hossam Aly, Jean-François Gonzalez, Rebecca Nealon, Cristiano Longarini, Giuseppe Lodato, Daniel J. Price

Abstract: Gas and dust in inclined orbits around binaries experience precession induced by the binary gravitational torque. The difference in precession between gas and dust alters the radial drift of weakly coupled dust and leads to density enhancements where the radial drift is minimised. We explore this phenomenon using 3D hydrodynamical simulations to investigate the prominence of these `dust traffic ja… ▽ More Gas and dust in inclined orbits around binaries experience precession induced by the binary gravitational torque. The difference in precession between gas and dust alters the radial drift of weakly coupled dust and leads to density enhancements where the radial drift is minimised. We explore this phenomenon using 3D hydrodynamical simulations to investigate the prominence of these `dust traffic jams' and the evolution of the resulting dust sub-structures at different disc inclinations and binary eccentricities. We then derive evolution equations for the angular momentum of warped dust discs and implement them in a 1D code and present calculations to further explain these traffic jams. We find that dust traffic jams in inclined circumbinary discs provide significant dust density enhancements that are long lived and can have important consequences for planetesimal formation. △ Less

Submitted 27 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS, 15 pages, 13 figures

Dust growth, fragmentation and self-induced dust traps in PHANTOM

Authors: Arnaud Vericel, Jean-François Gonzalez, Daniel J. Price, Guillaume Laibe, Christophe Pinte

Abstract: We present the implementation of a dust growth and fragmentation module in the public Smoothed Particle Hydrodynamics (SPH) code PHANTOM. This module is made available for public use with this paper. The coagulation model considers locally monodisperse dust size distributions around single values that are carried by the SPH particles. Along with the presentation of the model, implementation and te… ▽ More We present the implementation of a dust growth and fragmentation module in the public Smoothed Particle Hydrodynamics (SPH) code PHANTOM. This module is made available for public use with this paper. The coagulation model considers locally monodisperse dust size distributions around single values that are carried by the SPH particles. Along with the presentation of the model, implementation and tests, we showcase growth and fragmentation in a few typical circumstellar disc simulations and revisit previous results. The module is also interfaced with the radiative transfer code MCFOST, which facilitates the comparison between simulations and ALMA observations by generating synthetic maps. Circumstellar disc simulations with growth and fragmentation reproduce the `self-induced dust trap' mechanism first proposed by Gonzalez et al., which supports its existence. Synthetic images of discs featuring this mechanism suggest it would be detectable by ALMA as a bright axisymmetric ring at several tens of au from the star. With this paper, our aim is to provide a public tool to be able to study and explore dust growth in a variety of applications related to planet formation. △ Less

Submitted 2 August, 2021; originally announced August 2021.

Comments: 22 pages, 27 figures, accepted for publication in MNRAS

Circumbinary and circumstellar discs around the eccentric binary IRAS 04158+2805 -- a testbed for binary-disc interaction

Authors: Enrico Ragusa, Daniele Fasano, Claudia Toci, Gaspard Duchêne, Nicolás Cuello, Marion Villenave, Gerrit van der Plas, Giuseppe Lodato, François Ménard, Daniel J. Price, Christophe Pinte, Karl Stapelfeldt, Schuyler Wolff

Abstract: IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$α$ imaging and Keck ast… ▽ More IRAS~04158+2805 has long been thought to be a very low mass T-Tauri star (VLMS) surrounded by a nearly edge-on, extremely large disc. Recent observations revealed that this source hosts a binary surrounded by an extended circumbinary disc with a central dust cavity. In this paper, we combine ALMA multi-wavelength observations of continuum and $^{12}$CO line emission, with H$α$ imaging and Keck astrometric measures of the binary to develop a coherent dynamical model of this system. The system features an azimuthal asymmetry detected at the western edge of the cavity in Band~7 observations and a wiggling outflow. Dust emission in ALMA Band 4 from the proximity of the individual stars suggests the presence of marginally resolved circumstellar discs. We estimate the binary orbital parameters from the measured arc of the orbit from Keck and ALMA astrometry. We further constrain these estimates using considerations from binary-disc interaction theory. We finally perform three SPH gas + dust simulations based on the theoretical constraints; we post-process the hydrodynamic output using radiative transfer Monte Carlo methods and directly compare the models with observations. Our results suggest that a highly eccentric $e\sim 0.5\textrm{--}0.7$ equal mass binary, with a semi-major axis of $\sim 55$ au, and small/moderate orbital plane vs. circumbinary disc inclination $θ\lesssim 30^\circ$ provides a good match with observations. A dust mass of $\sim 1.5\times 10^{-4} {\rm M_\odot}$ best reproduces the flux in Band 7 continuum observations. Synthetic CO line emission maps qualitatively capture both the emission from the central region and the non-Keplerian nature of the gas motion in the binary proximity. △ Less

Submitted 28 July, 2021; originally announced July 2021.

Comments: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

Formation of eccentric gas discs from sublimating or partially disrupted asteroids orbiting white dwarfs

Authors: David Trevascus, Daniel J. Price, Rebecca Nealon, David Liptai, Christopher J. Manser, Dimitri Veras

Abstract: Of the 21 known gaseous debris discs around white dwarfs, a large fraction of them display observational features that are well described by an eccentric distribution of gas. In the absence of embedded objects or additional forces, these discs should not remain eccentric for long timescales, and should instead circularise due to viscous spreading. The metal pollution and infrared excess we observe… ▽ More Of the 21 known gaseous debris discs around white dwarfs, a large fraction of them display observational features that are well described by an eccentric distribution of gas. In the absence of embedded objects or additional forces, these discs should not remain eccentric for long timescales, and should instead circularise due to viscous spreading. The metal pollution and infrared excess we observe from these stars is consistent with the presence of tidally disrupted sub-stellar bodies. We demonstrate, using smoothed particle hydrodynamics, that a sublimating or partially disrupting planet on an eccentric orbit around a white dwarf will form and maintain a gas disc with an eccentricity within 0.1 of, and lower than, that of the orbiting body. We also demonstrate that the eccentric gas disc observed around the white dwarf SDSS J1228+1040 can be explained by the same hypothesis. △ Less

Submitted 3 May, 2021; originally announced May 2021.

Comments: 5 pages, 5 figures, accepted to MNRAS Letters

Uncovering Erosion Effects on Magnetic Flux Rope Twist

Authors: Sanchita Pal, Emilia Kilpua, Simon Good, Jens Pomoell, Daniel J. Price

Abstract: Magnetic clouds (MCs) are transient structures containing large-scale magnetic flux ropes from solar eruptions. The twist of magnetic field lines around the rope axis reveals information about flux rope formation processes and geoeffectivity. During propagation, MC flux ropes may erode via reconnection with the ambient solar wind. Any erosion reduces the magnetic flux and helicity of the ropes, an… ▽ More Magnetic clouds (MCs) are transient structures containing large-scale magnetic flux ropes from solar eruptions. The twist of magnetic field lines around the rope axis reveals information about flux rope formation processes and geoeffectivity. During propagation, MC flux ropes may erode via reconnection with the ambient solar wind. Any erosion reduces the magnetic flux and helicity of the ropes, and changes their cross-sectional twist profiles. This study relates twist profiles in MC flux ropes observed at 1 AU to the amount of erosion undergone by the MCs in interplanetary space. The twist profiles of two well-identified MC flux ropes associated with the clear appearance of post eruption arcades in the solar corona are analysed. To infer the amount of erosion, the magnetic flux content of the ropes in the solar atmosphere is estimated, and compared to estimates at 1 AU. The first MC shows a monotonically decreasing twist from the axis to periphery, while the second displays high twist at the axis, rising twist near the edges, and lower twist in between. The first MC displays a larger reduction in magnetic flux between the Sun and 1 AU, suggesting more erosion, than that seen in the second MC. In the second cloud, rising twist at the rope edges may have been due to an envelope of overlying coronal field lines with relatively high twist, formed by reconnection beneath the erupting flux rope in the low corona. This high-twist envelope remained almost intact from the Sun to 1 AU due to the low erosion levels. In contrast, the high-twist envelope of the first cloud may have been entirely peeled away via erosion by the time it reaches 1 AU. △ Less

Submitted 8 April, 2021; originally announced April 2021.

Comments: 10 pages, 5 figures, 1 table

Journal ref: A&A 650, A176 (2021)

HD 143006: circumbinary planet or misaligned disc?

Authors: G. Ballabio, R. Nealon, R. D. Alexander, N. Cuello, C. Pinte, D. J. Price

Abstract: Misalignments within protoplanetary discs are now commonly observed, and features such as shadows in scattered light images indicate departure from a co-planar geometry. VLT/SPHERE observations of the disc around HD 143006 show a large-scale asymmetry, and two narrow dark lanes which are indicative of shadowing. ALMA observations also reveal the presence of rings and gaps in the disc, along with a… ▽ More Misalignments within protoplanetary discs are now commonly observed, and features such as shadows in scattered light images indicate departure from a co-planar geometry. VLT/SPHERE observations of the disc around HD 143006 show a large-scale asymmetry, and two narrow dark lanes which are indicative of shadowing. ALMA observations also reveal the presence of rings and gaps in the disc, along with a bright arc at large radii. We present new hydrodynamic simulations of HD 143006, and show that a configuration with both a strongly inclined binary and an outer planetary companion is the most plausible to explain the observed morphological features. We compute synthetic observations from our simulations, and successfully reproduce both the narrow shadows and the brightness asymmetry seen in IR scattered light. Additionally, we reproduce the large dust observed in the mm continuum, due to a 10 Jupiter mass planet detected in the CO kinematics. Our simulations also show the formation of a circumplanetary disc, which is misaligned with respect to the outer disc. The narrow shadows cast by the inner disc and the planet-induced "kink" in the disc kinematics are both expected to move on a time-scale of $\sim$ 5-10 years, presenting a potentially observable test of our model. If confirmed, HD 143006 would be the first known example of a circumbinary planet on a strongly misaligned orbit. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: Accepted for publication in MNRAS, 10 pages, 8 figures. Movies of the simulation available at the following links: https://www.youtube.com/watch?v=F8GAw-Fzpyg&list=PLgaPAkHEP_RqxHJUvj9DIcsjEpWqbTPwQ&index=3 and https://www.youtube.com/watch?v=HCZQdS3t4l8&list=PLgaPAkHEP_RqxHJUvj9DIcsjEpWqbTPwQ&index=1

Electromagnetic Signatures from the Tidal Tail of a Black Hole - Neutron Star Merger

Authors: Siva Darbha, Daniel Kasen, Francois Foucart, Daniel J. Price

Abstract: Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave (GW) observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial sta… ▽ More Black hole - neutron star (BH-NS) mergers are a major target for ground-based gravitational wave (GW) observatories. A merger can also produce an electromagnetic counterpart (a kilonova) if it ejects neutron-rich matter that assembles into heavy elements through r-process nucleosynthesis. We study the kilonova signatures of the unbound dynamical ejecta of a BH-NS merger. We take as our initial state the results from a numerical relativity simulation, and then use a general relativistic hydrodynamics code to study the evolution of the ejecta with parameterized r-process heating models. The unbound dynamical ejecta is initially a flattened, directed tidal tail largely confined to a plane. Heating from the r-process inflates the ejecta into a more spherical shape and smooths its small-scale structure, though the ejecta retains its bulk directed motion. We calculate the electromagnetic signatures using a 3D radiative transfer code and a parameterized opacity model for lanthanide-rich matter. The light curve varies with viewing angle due to two effects: asphericity results in brighter emission for orientations with larger projected areas, while Doppler boosting results in brighter emission for viewing angles more aligned with the direction of bulk motion. For typical r-process heating rates, the peak bolometric luminosity varies by a factor of $\sim 3$ with orientation while the peak in the optical bands varies by $\sim 3$ magnitudes. The spectrum is blue-shifted at viewing angles along the bulk motion, which increases the $V$-band peak magnitude to $\sim -14$ despite the lanthanide-rich composition. △ Less

Submitted 24 August, 2021; v1 submitted 4 March, 2021; originally announced March 2021.

Comments: 26 pages, 15 figures; revised paper after reviewer's comments

The theory of kinks -- I. A semi-analytic model of velocity perturbations due to planet-disc interaction

Authors: Francesco Bollati, Giuseppe Lodato, Daniel J. Price, Christophe Pinte

Abstract: A new technique to detect protoplanets is by observing the kinematics of the surrounding gas. Gravitational perturbations from a planet produce peculiar `kinks' in channel maps of different gas species. In this paper, we show that such kinks can be reproduced using semi-analytic models for the velocity perturbation induced by a planet. In doing so we i) confirm that the observed kinks are consiste… ▽ More A new technique to detect protoplanets is by observing the kinematics of the surrounding gas. Gravitational perturbations from a planet produce peculiar `kinks' in channel maps of different gas species. In this paper, we show that such kinks can be reproduced using semi-analytic models for the velocity perturbation induced by a planet. In doing so we i) confirm that the observed kinks are consistent with the planet-induced wake; ii) show how to quantify the planet mass from the kink amplitude; in particular, we show that the kink amplitude scales with the square root of the planet mass for channels far from the planet velocity, steepening to linear as the channels approach the planet; iii) show how to extend the theory to include the effect of damping, which may be needed in order to have localized kinks. △ Less

Submitted 19 April, 2021; v1 submitted 26 February, 2021; originally announced March 2021.

Comments: Submission to MNRAS Accepted

A faint companion around CrA-9: protoplanet or obscured binary?

Authors: V. Christiaens, M. -G. Ubeira-Gabellini, H. Cánovas, P. Delorme, B. Pairet, O. Absil, S. Casassus, J. H. Girard, A. Zurlo, Y. Aoyama, G-D. Marleau, L. Spina, N. van der Marel, L. Cieza, G. Lodato, S. Pérez, C. Pinte, D. J. Price, M. Reggiani

Abstract: Understanding how giant planets form requires observational input from directly imaged protoplanets. We used VLT/NACO and VLT/SPHERE to search for companions in the transition disc of 2MASS J19005804-3645048 (hereafter CrA-9), an accreting M0.75 dwarf with an estimated age of 1-2 Myr. We found a faint point source at $\sim$0.7'' separation from CrA-9 ($\sim$108 au projected separation). Our 3-epoc… ▽ More Understanding how giant planets form requires observational input from directly imaged protoplanets. We used VLT/NACO and VLT/SPHERE to search for companions in the transition disc of 2MASS J19005804-3645048 (hereafter CrA-9), an accreting M0.75 dwarf with an estimated age of 1-2 Myr. We found a faint point source at $\sim$0.7'' separation from CrA-9 ($\sim$108 au projected separation). Our 3-epoch astrometry rejects a fixed background star with a $5σ$ significance. The near-IR absolute magnitudes of the object point towards a planetary-mass companion. However, our analysis of the 1.0-3.8$μ$m spectrum extracted for the companion suggests it is a young M5.5 dwarf, based on both the 1.13-$μ$m Na index and comparison with templates of the Montreal Spectral Library. The observed spectrum is best reproduced with high effective temperature ($3057^{+119}_{-36}$K) BT-DUSTY and BT-SETTL models, but the corresponding photometric radius required to match the measured flux is only $0.60^{+0.01}_{-0.04}$ Jovian radius. We discuss possible explanations to reconcile our measurements, including an M-dwarf companion obscured by an edge-on circum-secondary disc or the shock-heated part of the photosphere of an accreting protoplanet. Follow-up observations covering a larger wavelength range and/or at finer spectral resolution are required to discriminate these two scenarios. △ Less

Submitted 20 February, 2021; originally announced February 2021.

Comments: 24 pages, 14 figures, 4 tables, to be published in MNRAS