Search

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., Accepted for publication in MNRAS, 20 pages

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., 11 pages, 9 figures, accepted for publication in MNRAS

Many critical decisions faced in early discovery require a thorough understanding of the dynamic behavior of pharmacological pathways following target engagement. From fundamental decisions on the optimal target to pursue and the ultimate drug product profile (combination of modality, potency, and compound properties) expected to elicit the desired clinical outcome to tactical program decisions such as what chemical series to pursue, what chemical properties require optimization, and what compounds to synthesize and progress, all demand detailed consideration of pharmacodynamics. Model-based target pharmacology assessment (mTPA) is a computational approach centered around large-scale virtual exploration of pharmacokinetic and pharmacodynamic models built early in discovery to guide these decisions. The present work summarizes several examples (use cases) from programs at GlaxoSmithKline that demonstrate the utility of mTPA throughout the drug discovery lifecycle.

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., Comment: Accepted for publication in MNRAS Letters. 5 pages, 5 figures

Planet-disc interactions build up local pressure maxima that may halt the radial drift of protoplanetary dust, and pile it up in rings and crescents. ALMA observations of the HD135344B disc revealed two rings in the thermal continuum stemming from ~mm-sized dust. At higher frequencies the inner ring is brighter relative to the outer ring, which is also shaped as a crescent rather than a full ring. In near-IR scattered light images, the disc is modulated by a 2-armed grand-design spiral originating inside the ALMA inner ring. Such structures may be induced by a massive companion evacuating the central cavity, and by a giant planet in the gap separating both rings, that channels the accretion of small dust and gas through its filamentary wakes while stopping the larger dust from crossing the gap. Here we present ALMA observations in the J=(2-1)CO isotopologue lines and in the adjacent continuum, with up to 12km baselines. Angular resolutions of 0.03" reveal the tentative detection of a filament connecting both rings, and which coincides with a local discontinuity in the pitch angle of the IR spiral, proposed previously as the location of the protoplanet driving this spiral. Line diagnostics suggest that turbulence, or superposed velocity components, is particularly strong in the spirals. The 12CO(2-1) 3-D rotation curve points at stellocentric accretion at radii within the inner dust ring, with a radial velocity of up to ~6%+-0.5% Keplerian, which corresponds to an excessively large accretion rate of ~2E-6M_sun/yr if all of the CO layer follows the 12CO(2-1) kinematics. This suggests that only the surface layers of the disc are undergoing accretion, and that the line broadening is due to superposed laminar flows., accepted to MNRAS

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., 16 pages, 7 figures, 1 table, 255 references. Invited review article accepted in EPJ+

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., 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

We present ALMA observations of the $^{12}$CO, $^{13}$CO, C$^{18}$O J=2-1 transitions and the 1.3\,mm continuum emission for the circumbinary disc around HD 142527, at an angular resolution of $\approx$\,0\farcs3. We observe multiple spiral structures in intensity, velocity and velocity dispersion for the $^{12}$CO and $^{13}$CO gas tracers. A newly detected $^{12}$CO spiral originates from the dust horseshoe, and is rotating at super-Keplerian velocity or vertically ascending, whilst the inter-spiral gas is rotating at sub-Keplerian velocities. This new spiral possibly connects to a previously identified spiral, thus spanning > 360$^\circ$. A spatial offset of ~30 au is observed between the $^{12}$CO and $^{13}$CO spirals, to which we hypothesize that the gas layers are propagating at different speeds (``surfing'') due to a non-zero vertical temperature gradient. Leveraging the varying optical depths between the CO isotopologues, we reconstruct temperature and column density maps of the outer disc. Gas surface density peaks at r\,$\approx$\,180\,au, coincident with the peak of continuum emission. Here the dust grains have a Stokes number of $\approx$\,1, confirming radial and azimuthal trapping in the horseshoe. We measure a cavity radius at half-maximum surface density of $\approx$\,100\,au, and a cavity eccentricity between 0.3 and 0.45.

We present 1.3 mm continuum ALMA long-baseline observations at 3-5 au resolution of 10 of the brightest discs from the Ophiuchus DIsc Survey Employing ALMA (ODISEA) project. We identify a total of 26 narrow rings and gaps distributed in 8 sources and 3 discs with small dust cavities (r $$20 au). We find that the 1.3 mm radial profiles of these objects are in good agreement with those produced by numerical simulations of dust evolution and planet-disc interactions, which predict the accumulation of mm-sized grains at the edges of planet-induced cavities. Our long-baseline observations resulted in the largest sample of discs observed at $\sim$3-5 au resolution in any given star-forming region (15 objects when combined with Ophiuchus objects in the DSHARP Large Program) and allow for a demographic study of the brightest $\sim5\%$ of the discs in Ophiuchus (i.e. the most likely formation sites of giant planets in the cloud). We use this unique sample to propose an evolutionary sequence and discuss a scenario in which the substructures observed in massive protoplanetary discs are mainly the result of planet formation and dust evolution. If this scenario is correct, the detailed study of disc substructures might provide a window to investigate a population of planets that remains mostly undetectable by other techniques., Comment: 21 pages, 10 figures. Appendix with 3 additional figures. Version 2, identical to previous one, but now accepted for publication (in MNRAS)

How does the cavity size in circumbinary discs depend on disc and binary properties? We investigate by simulating disc cavities carved by binary companions using smoothed particle hydrodynamics (SPH). We find that a cavity is quickly opened on the dynamical time, while the cavity size is set on the viscous time. In agreement with previous findings, we find long term cavity sizes of 2-5 times the binary semi-major axis, increasing with eccentricity and decreasing with disc aspect ratio. When considering binaries inclined with respect to the disc we find three regimes: i) discs that evolve towards a coplanar orbit have a large cavity, slightly smaller than that of an initially coplanar disc; ii) discs that evolve towards a polar orbit by breaking have a small cavity, equal in size to that of an initially polar disc; iii) discs that evolve towards a polar orbit via warping have an intermediate-sized cavity. We find typical gas depletions inside the cavity of $\gtrsim 2$ orders of magnitude in surface density., Comment: 12 pages, 13 figures, accepted to MNRAS

During the common envelope binary interaction, the expanding layers of the gaseous common envelope recombine and the resulting recombination energy has been suggested as a contributing factor to the ejection of the envelope. In this paper we perform a comparative study between simulations with and without the inclusion of recombination energy. We use two distinct setups, comprising 0.88-M$_{\odot}$ and 1.8-M$_{\odot}$ giants, that have been studied before and can serve as benchmarks. In so doing we conclude that (i) the final orbital separation is not affected by the choice of equation of state. In other words, simulations that unbind but a small fraction of the envelope result in similar final separations to those that, thanks to recombination energy, unbind a far larger fractions. (ii) The adoption of a tabulated equation of state results in a much greater fraction of unbound envelope and we demonstrate the cause of this to be the release of recombination energy. (iii) The fraction of hydrogen recombination energy that is allowed to do work should be about half of that which our adiabatic simulations use. (iv) However, for the heavier star simulation we conclude that it is helium and not hydrogen recombination energy that unbinds the gas and we determine that all helium recombination energy is thermalised in the envelope and does work. (v)~The outer regions of the expanding common envelope are likely to see the formation of dust. This dust would promote additional unbinding and shaping of the ejected envelope into axisymmetric morphologies., Comment: 18 pages, 11 figures. We welcome any comments or questions

Artificial intelligence and machine learning have demonstrated their potential role in predictive chemistry and synthetic planning of small molecules; there are at least a few reports of companies employing in silico synthetic planning into their overall approach to accessing target molecules. A data-driven synthesis planning program is one component being developed and evaluated by the Machine Learning for Pharmaceutical Discovery and Synthesis (MLPDS) consortium, comprising MIT and 13 chemical and pharmaceutical company members. Together, we wrote this perspective to share how we think predictive models can be integrated into medicinal chemistry synthesis workflows, how they are currently used within MLPDS member companies, and the outlook for this field.

To better predict the impacts of solar eruptions on Earth, understanding the low-corona evolution of CMEs is crucial because this influential early phase is highly dynamic. We therefore investigate the evolution of CME properties, such as the evolution of flux rope footpoints as well as the magnetic flux enclosed in the flux rope, during this stage of the eruption. To simulate the eruption we make use of the data-driven time-dependent magnetofrictional method (TMFM), which has been proven to accurately capture a flux rope's early evolution and lift-off. We then developed a semi-automatized method for identifying the flux rope and extracting these flux ropes from 3D data cubes and tracking their evolution in time. The extraction algorithm is based on the twist parameter Tw in a 2D plane close to the polarity inversion line as a proxy for the flux rope and its temporal evolution. It is then applied to TMFM simulations of the active region AR12473, which produced an eruption on 28th of December 2015 (see e.g., Price et al, 2020). This CME was also accompanied by an M1.9 flare, that peaked at about 12:45 UT. The extracted flux rope footpoints are then compared against observational data from SDO's AIA instrument in the 1600 Å wavelength. This comparison yields a very good match with coverage parameters (see Asvestari et al, 2019) in the range of 60-70 %. The magnetic flux is extracted from the footpoints that are rooted in one specific polarity region.

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., 19 pages, 17 figures, accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS, 15 pages, 13 figures, 1 tables, 2 appendices

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., Comment: Accepted to ApJL, 14 pages, 10 figures

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa)., 11 pages, 5 figures, accepted for publication in MNRAS

[Image: see text] We previously described the discovery of GSK5852 (1), a non-nucleoside polymerase (NS5B) inhibitor of hepatitis C virus (HCV), in which an N-benzyl boronic acid was essential for potent antiviral activity. Unfortunately, facile benzylic oxidation resulted in a short plasma half-life (5 h) in human volunteers, and a backup program was initiated to remove metabolic liabilities associated with 1. Herein, we describe second-generation NS5B inhibitors including GSK8175 (49), a sulfonamide-N-benzoxaborole analog with low in vivo clearance across preclinical species and broad-spectrum activity against HCV replicons. An X-ray structure of NS5B protein cocrystallized with 49 revealed unique protein-inhibitor interactions mediated by an extensive network of ordered water molecules and the first evidence of boronate complex formation within the binding pocket. In clinical studies, 49 displayed a 60–63 h half-life and a robust decrease in viral RNA levels in HCV-infected patients, thereby validating our hypothesis that reducing benzylic oxidation would improve human pharmacokinetics and lower efficacious doses relative to 1.

Nearly all young stars are initially surrounded by `protoplanetary' discs of gas and dust, and in the case of single stars at least 30\% of these discs go on to form planets. The process of protoplanetary disc formation can result in initial misalignments, where the disc orbital plane is different to the stellar equator in single star systems, or to the binary orbital plane in systems with two stars. A quirk of the dynamics means that initially misaligned `circumbinary' discs -- those that surround two stars -- are predicted to evolve to one of two possible stable configurations, one where the disc and binary orbital planes are coplanar, and one where they are perpendicular (a `polar' configuration). Prior work has found coplanar circumbinary discs, but no polar examples were known until now. Here we report the first discovery of a protoplanetary circumbinary disc in the polar configuration, supporting the predictions that such discs should exist. The disc shows some characteristics that are similar to discs around single stars, and that are attributed to dust growth. Thus, the first stages of planet formation appear able to proceed in polar circumbinary discs., Comment: Submitted version - do not read. Published in Nature Astronomy, free PDF at https://rdcu.be/bgNSO

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., 5 pages, 5 figures, accepted to MNRAS Letters

The concept of using precipitation inhibitors (PIs) to sustain supersaturation is well established for amorphous formulations but less in the case of lipid-based formulations (LBF). This study applied a systematic in silico-in vitro-in vivo approach to assess the merits of incorporating PIs in supersaturated LBFs (sLBF) using the model drug venetoclax. sLBFs containing hydroxypropyl methylcellulose (HPMC), hydroxypropyl methylcellulose acetate succinate (HPMCAS), polyvinylpyrrolidone (PVP), PVP-co-vinyl acetate (PVP/VA), Pluronic F108, and Eudragit EPO were assessed in silico calculating a drug-excipient mixing enthalpy, in vitro using a PI solvent shift test, and finally, bioavailability was assessed in vivo in landrace pigs. The estimation of pure interaction enthalpies of the drug and the excipient was deemed useful in determining the most promising PIs for venetoclax. The sLBF alone (i.e., no PI present) displayed a high initial drug concentration in the aqueous phase during in vitro screening. sLBF with Pluronic F108 displayed the highest venetoclax concentration in the aqueous phase and sLBF with Eudragit EPO the lowest. In vivo, the sLBF alone showed the highest bioavailability of 26.3 ± 14.2%. Interestingly, a trend toward a decreasing bioavailability was observed for sLBF containing PIs, with PVP/VA being significantly lower compared to sLBF alone. In conclusion, the ability of a sLBF to generate supersaturated concentrations of venetoclax in vitro was translated into increased absorption in vivo. While in silico and in vitro PI screening suggested benefits in terms of prolonged supersaturation, the addition of a PI did not increase in vivo bioavailability. The findings of this study are of particular relevance to pre-clinical drug development, where the high in vivo exposure of venetoclax was achieved using a sLBF approach, and despite the perceived risk of drug precipitation from a sLBF, including a PI may not be merited in all cases.

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., 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

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\sigma$ significance. The near-IR absolute magnitudes of the object point towards a planetary-mass companion. However, our analysis of the 1.0-3.8$\mu$m spectrum extracted for the companion suggests it is a young M5.5 dwarf, based on both the 1.13-$\mu$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., Comment: 24 pages, 14 figures, 4 tables, to be published in MNRAS

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., 22 pages, 27 figures, accepted for publication in MNRAS

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., Comment: 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

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$\alpha$ 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 $\theta\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., Comment: 19 pages, 14 figures, 1 table, accepted for publication in MNRAS

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., Comment: Accepted for publication in MNRAS, 17 pages excluding appendices, 15 figures

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., Comment: 26 pages, 15 figures; revised paper after reviewer's comments

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., Comment: 10 Pages, 6 figures. Accepted in ApJL