Source author record

W. Lyra

W. Lyra appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

16works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

16 published item(s)

preprint2022arXiv

A Predicted Dearth of Majority Hypervolatile Ices in Oort Cloud Comets

We present new, ice species-specific New Horizons/Alice upper gas coma production limits from the 01 Jan 2019 MU69/Arrokoth flyby of Gladstone et al. (2021) and use them to make predictions about the rarity of majority hypervolatile (CO, N$_2$, CH$_4$) ices in KBOs and Oort Cloud comets. These predictions have a number of important implications for the study of the Oort Cloud, including: determination of hypervolatile rich comets as the first objects emplaced into the Oort Cloud; measurement of CO/N$_2$/CH$_4$ abundance ratios in the proto-planetary disk from hypervolatile rich comets; and population statistical constraints on early (< 20 Myr) planetary aggregation driven versus later (> 50 Myr) planetary migration driven emplacement of objects into the Oort Cloud. They imply that the phenomenon of ultra-distant active comets like C/2017K2 (Jewitt et al. 2017, Hui et al. 2018) should be rare, and thus not a general characteristic of all comets. They also suggest that interstellar object 2I/Borisov did not originate in a planetary system that was inordinately CO rich (Bodewits et al. 2020), but rather could have been ejected onto an interstellar trajectory very early in its natal system's history.

preprint2022arXiv

Hydro-, Magnetohydro-, and Dust-Gas Dynamics of Protoplanetary Disks

The building of planetary systems is controlled by the gas and dust dynamics of protoplanetary disks. While the gas is simultaneously accreted onto the central star and dissipated away by winds, dust grains aggregate and collapse to form planetesimals and eventually planets. This dust and gas dynamics involves instabilities, turbulence and complex non-linear interactions which ultimately control the observational appearance and the secular evolution of these disks. This chapter is dedicated to the most recent developments in our understanding of the dynamics of gaseous and dusty disks, covering hydrodynamic and magnetohydrodynamic turbulence, gas-dust instabilities, dust clumping and disk winds. We show how these physical processes have been tested from observations and highlight standing questions that should be addressed in the future.

preprint2022arXiv

The thermodynamics of stellar multiplicity: dynamical evolution of binary star populations in dense stellar environments

We recently derived, using the density-of-states approximation, analytic distribution functions for the outcomes of direct single-binary scatterings (Stone & Leigh 2019). Using these outcome distribution functions, we present in this paper a self-consistent statistical mechanics-based analytic model obtained using the Fokker-Planck limit of the Boltzmann equation. Our model quantifies the dominant gravitational physics, combining both strong and weak single-binary interactions, that drives the time evolution of binary orbital parameter distributions in dense stellar environments. We focus in particular the distributions of binary orbital energies and eccentricities. We find a novel steady state distribution of binary eccentricities, featuring strong depletions of both the highest and the lowest eccentricity binaries. In energy space, we compare the predictions of our analytic model to the results of numerical N-body simulations, and find that the agreement is good for the initial conditions considered here. This work is a first step toward the development of a fully self-consistent semi-analytic model for dynamically evolving binary star populations in dense stellar environments due to direct few-body interactions.

preprint2020arXiv

On the Origin and Thermal Stability of Arrokoths and Plutos Ices

We discuss in a thermodynamic, geologically empirical way the long-term nature of the stable majority ices that could be present in Kuiper Belt Object 2014 MU69 after its 4.6 Gyr residence in the EKB as a cold classical object. Considering the stability versus sublimation into vacuum for the suite of ices commonly found on comets, Centaurs, and KBOs at the average ~40K sunlit surface temperature of MU69 over Myr to Gyr, we find only 3 common ices that are truly refractory: HCN, CH3OH, and H2O (in order of increasing stability). NH3 and H2CO ices are marginally stable and may be removed by any positive temperature excursions in the EKB, as produced every 1e8 - 1e9 yrs by nearby supernovae and passing O/B stars. To date the NH team has reported the presence of abundant CH3OH and evidence for H2O on MU69s surface (Lisse et al. 2017, Grundy et al. 2020). NH3 has been searched for, but not found. We predict that future absorption feature detections will be due to an HCN or poly-H2CO based species. Consideration of the conditions present in the EKB region during the formation era of MU69 lead us to infer that it formed "in the dark", in an optically thick mid-plane, unable to see the nascent, variable, highly luminous Young Stellar Object-TTauri Sun, and that KBOs contain HCN and CH3OH ice phases in addition to the H2O ice phases found in their Short Period comet descendants. Finally, when we apply our ice thermal stability analysis to bodies/populations related to MU69, we find that methanol ice may be ubiquitous in the outer solar system; that if Pluto is not a fully differentiated body, then it must have gained its hypervolatile ices from proto-planetary disk sources in the first few Myr of the solar systems existence; and that hypervolatile rich, highly primordial comet C/2016 R2 was placed onto an Oort Cloud orbit on a similar timescale.

preprint2020arXiv

The Pencil Code, a modular MPI code for partial differential equations and particles: multipurpose and multiuser-maintained

The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of add-ons ranging from astrophysical magnetohydrodynamics (MHD) to meteorological cloud microphysics and engineering applications in combustion. Nevertheless, the framework is general and can also be applied to situations not related to hydrodynamics or even PDEs, for example when just the message passing interface or input/output strategies of the code are to be used. The code can also evolve Lagrangian (inertial and noninertial) particles, their coagulation and condensation, as well as their interaction with the fluid.

preprint2015arXiv

Compact dust concentration in the MWC 758 protoplanetary disk

The formation of planetesimals requires that primordial dust grains grow from micron- to km-sized bodies. Dust traps caused by gas pressure maxima have been proposed as regions where grains can concentrate and grow fast enough to form planetesimals, before radially migrating onto the star. We report new VLA Ka & Ku observations of the protoplanetary disk around the Herbig Ae/Be star MWC 758. The Ka image shows a compact emission region in the outer disk indicating a strong concentration of big dust grains. Tracing smaller grains, archival ALMA data in band 7 continuum shows extended disk emission with an intensity maximum to the north-west of the central star, which matches the VLA clump position. The compactness of the Ka emission is expected in the context of dust trapping, as big grains are trapped more easily than smaller grains in gas pressure maxima. We develop a non-axisymmetric parametric model inspired by a steady state vortex solution with parameters adequately selected to reproduce the observations, including the spectral energy distribution. Finally, we compare the radio continuum with SPHERE scattered light data. The ALMA continuum spatially coincides with a spiral-like feature seen in scattered light, while the VLA clump is offset from the scattered light maximum. Moreover, the ALMA map shows a decrement that matches a region devoid of scattered polarised emission. Continuum observations at a different wavelength are necessary to conclude if the VLA-ALMA difference is an opacity or a real dust segregation.

preprint2014arXiv

Intermediate mass black holes in AGN disks II. Model predictions & observational constraints

If intermediate mass black holes (IMBHs) grow efficiently in gas disks around supermassive black holes, their host active galactic nucleus (AGN) disks should exhibit myriad observational signatures. Gap-opening IMBHs in AGN disks can exhibit spectral features and variability analagous to gapped protoplanetary disks. A gap-opening IMBH in the innermost disk imprints ripples and oscillations on the broad Fe K$α$ line which may be detectable with future X-ray missions. A non-gap-opening IMBH will accrete and produce a soft X-ray excess relative to continuum emission. An IMBH on a retrograde orbit in an AGN disk will not open a gap and will generate soft X-rays from a bow-shock 'headwind'. Accreting IMBH in a large cavity can generate ULX-like X-ray luminosities and LINER-like optical line ratios from local ionized gas. We propose that many LINERs house a weakly accreting MBH binary in a large central disk cavity and will be luminous sources of gravitational waves (GW). IMBHs in galactic nuclei may also be detected via intermittent observational signatures including: UV/X-ray flares due to tidal disruption events, asymmetric X-ray intensity distributions as revealed by AGN transits, quasi-periodic oscillations and underluminous Type Ia supernovae. GW emitted during IMBH inspiral and collisions may be detected with eLISA and LIGO, particularly from LINERs. We summarize observational signatures and compare to current data where possible or suggest future observations.

preprint2014arXiv

Rossby wave instability does not require sharp resistivity gradients

Rossby wave instability (RWI) at dead zone boundaries may play an important role in planet formation. Viscous hydrodynamics results suggest RWI is excited only when the viscosity changes over a radial distance less than two density scale heights. However in the disks around Solar-mass T Tauri stars, it is not viscosity but magnetic forces that provide the accretion stress beyond about 10 AU, where surface densities are low enough so stellar X-rays and interstellar cosmic rays can penetrate. Here we aim to explore the conditions for RWI in the smooth transition with increasing distance, from resistive and magnetically-dead to conducting and magnetically-active. We perform 3D unstratified MHD simulations with the Pencil Code, using static resistivity profiles. As a result, we find that in MHD, contrary to viscous models, the RWI is triggered even with a gradual change in resistivity extending from 10 to 40 AU (i.e., spanning 15 scale heights for aspect ratio 0.1). This is because magneto-rotational turbulence sets in abruptly when the resistivity reaches a threshold level. At higher resistivities the longest unstable wavelength is quenched, resulting in a sharp decline of the Maxwell stress towards the star. The sharp gradient in the magnetic forces leads to a localized density bump, that is in turn Rossby wave unstable. We conclude that even weak gradients in the resistivity can lead to sharp transitions in the Maxwell stress. The upshot is that the RWI is more easily activated in the outer disk than previously thought. Rossby vortices at the outer dead zone boundary thus could underlie the dust asymmetries seen in the outer reaches of transition disks.

preprint2013arXiv

A Parameter Study for Baroclinic Vortex Amplification

Recent studies have shown that baroclinic vortex amplification is strongly dependent on certain factors, namely, the global entropy gradient, the efficiency of thermal diffusion and/or relaxation as well as numerical resolution. We conduct a comprehensive study of a broad range and combination of various entropy gradients, thermal diffusion and thermal relaxation time-scales via local shearing sheet simulations covering the parameter space relevant for protoplanetary disks. We measure the Reynolds stresses as a function of our control parameters and see that there is angular momentum transport even for entropy gradients as low as $β=-{d\ln s}/{d\ln r}={1}/{2}$, which corresponds to values observed in protoplanetary accretion disks. The amplification-rate of the perturbations, $Γ$, appears to be proportional to $β^2$ and thus proportional to the square of the \BV ($Γ\propto β^2 \propto N^2$). The saturation level of Reynolds stresses on the other hand seems to be proportional to $β^{1/2}$. This highlights the importance of baroclinic effects even for the low entropy gradients expected in protoplanetary disks.

preprint2013arXiv

Formation of sharp eccentric rings in debris disks with gas but without planets

"Debris disks" around young stars (analogues of the Kuiper Belt in our Solar System) show a variety of non-trivial structures attributed to planetary perturbations and used to constrain the properties of the planets. However, these analyses have largely ignored the fact that some debris disks are found to contain small quantities of gas, a component that all such disks should contain at some level. Several debris disks have been measured with a dust-to-gas ratio around unity at which the effect of hydrodynamics on the structure of the disk cannot be ignored. Here we report linear and nonlinear modelling that shows that dust-gas interactions can produce some of the key patterns attributed to planets. We find a robust clumping instability that organizes the dust into narrow, eccentric rings, similar to the Fomalhaut debris disk. The conclusion that such disks might contain planets is not necessarily required to explain these systems.

preprint2012arXiv

Intermediate mass black holes in AGN disks: I. Production & Growth

Here we propose a mechanism for efficiently growing intermediate mass black holes (IMBH) in disks around supermassive black holes. Stellar mass objects can efficiently agglomerate when facilitated by the gas disk. Stars, compact objects and binaries can migrate, accrete and merge within disks around supermassive black holes. While dynamical heating by cusp stars excites the velocity dispersion of nuclear cluster objects (NCOs) in the disk, gas in the disk damps NCO orbits. If gas damping dominates, NCOs remain in the disk with circularized orbits and large collision cross-sections. IMBH seeds can grow extremely rapidly by collisions with disk NCOs at low relative velocities, allowing for super-Eddington growth rates. Once an IMBH seed has cleared out its feeding zone of disk NCOs, growth of IMBH seeds can become dominated by gas accretion from the AGN disk. However, the IMBH can migrate in the disk and expand its feeding zone, permitting a super-Eddington accretion rate to continue. Growth of IMBH seeds via NCO collisions is enhanced by a pile-up of migrators. We highlight the remarkable parallel between the growth of IMBH in AGN disks with models of giant planet growth in protoplanetary disks. If an IMBH becomes massive enough it can open a gap in the AGN disk. IMBH migration in AGN disks may stall, allowing them to survive the end of the AGN phase and remain in galactic nuclei. Our proposed mechanisms should be more efficient at growing IMBH in AGN disks than the standard model of IMBH growth in stellar clusters. Dynamical heating of disk NCOs by cusp stars is transferred to the gas in a AGN disk helping to maintain the outer disk against gravitational instability. Model predictions, observational constraints and implications are discussed in a companion paper (Paper II).

preprint2011arXiv

On rapid migration and accretion within disks around supermassive black holes

Galactic nuclei should contain a cluster of stars and compact objects in the vicinity of the central supermassive black hole due to stellar evolution, minor mergers and gravitational dynamical friction. By analogy with protoplanetary migration, nuclear cluster objects (NCOs) can migrate in the accretion disks that power active galactic nuclei by exchanging angular momentum with disk gas. Here we show that an individual NCO undergoing runaway outward migration comparable to Type III protoplanetary migration can generate an accretion rate corresponding to Seyfert AGN or quasar luminosities. Multiple migrating NCOs in an AGN disk can dominate traditional viscous disk accretion and at large disk radii, ensemble NCO migration and accretion could provide sufficient heating to prevent the gravitational instability from consuming disk gas in star formation. The magnitude and energy of the X-ray soft excess observed at ~0.1-1keV in Seyfert AGN could be explained by a small population of ~10^{2}-10^{3} accreting stellar mass black holes or a few ULXs. NCO migration and accretion in AGN disks are therefore extremely important mechanisms to add to realistic models of AGN disks.

preprint2011arXiv

The baroclinic instability in the context of layered accretion. Self-sustained vortices and their magnetic stability in local compressible unstratified models of protoplanetary disks

Turbulence and angular momentum transport in accretion disks remains a topic of debate. With the realization that dead zones are robust features of protoplanetary disks, the search for hydrodynamical sources of turbulence continues. A possible source is the baroclinic instability (BI), which has been shown to exist in unmagnetized non-barotropic disks. We present shearing box simulations of baroclinicly unstable, magnetized, 3D disks, in order to assess the interplay between the BI and other instabilities, namely the magneto-rotational instability (MRI) and the magneto-elliptical instability. We find that the vortices generated and sustained by the baroclinic instability in the purely hydrodynamical regime do not survive when magnetic fields are included. The MRI by far supersedes the BI in growth rate and strength at saturation. The resulting turbulence is virtually identical to an MRI-only scenario. We measured the intrinsic vorticity profile of the vortex, finding little radial variation in the vortex core. Nevertheless, the core is disrupted by an MHD instability, which we identify with the magneto-elliptic instability. This instability has nearly the same range of unstable wavelengths as the MRI, but has higher growth rates. In fact, we identify the MRI as a limiting case of the magneto-elliptic instability, when the vortex aspect ratio tends to infinity (pure shear flow). We conclude that vortex excitation and self-sustenance by the baroclinic instability in protoplanetary disks is viable only in low ionization, i.e., the dead zone. Our results are thus in accordance with the layered accretion paradigm. A baroclinicly unstable dead zone should be characterized by the presence of large-scale vortices whose cores are elliptically unstable, yet sustained by the baroclinic feedback. As magnetic fields destroy the vortices and the MRI outweighs the BI, the active layers are unmodified.

preprint2010arXiv

Orbital migration of low-mass planets in evolutionary radiative models: Avoiding catastrophic infall

Outward migration of low-mass planets has recently been shown to be a possibility in non-barotropic disks. We examine the consequences of this result in evolutionary models of protoplanetary disks. Planet migration occurs towards equilibrium radii with zero torque. These radii themselves migrate inwards because of viscous accretion and photoevaporation. We show that as the surface density and temperature fall, the planet orbital migration and disk depletion timescales eventually become comparable, with the precise timing depending on the mass of the planet. When this occurs, the planet decouples from the equilibrium radius. At this time, however, the gas surface density is already too low to drive substantial further migration. A higher mass planet, of 10 Earth masses, can open a gap during the late evolution of the disk, and stops migrating. Low mass planets, with 1 or 0.1 Earth masses, released beyond 1 AU in our models, avoid migrating into the star. Our results provide support for the reduced migration rates adopted in recent planet population synthesis models.

preprint2009arXiv

Fine structure of the chromospheric activity in Solar-type stars - The Halpha Line

A calibration of H-alpha as both a chromospheric diagnostic and an age indicator is presented, complementing the works previously done on this subject (Herbig 1985, Pasquini & Pallavicini 1991. The chromospheric diagnostic was built with a statistically significant sample, covering nine years of observations, and including 175 solar neighborhood stars. Regarding the age indicator, the presence of stars for which very accurate ages are determined, such as those belonging to clusters and kinematic groups, lends confidence to our analysis. We also investigate the possibility that stars of the same age might have gone through different tracks of chromospheric decay, identifying - within the same age range - effects of metallicity and mass. These parameters, however, as well as age, seem to be significant only for dwarf stars, losing their meaning when we analyze stars in the subgiant branch. This result suggests that, in these evolved stars, the emission mechanism cannot be magnetohydrodynamical in nature, in agreement with recent models (Fawzy et al. 2002c, and references therein). The Sun is found to be a typical star in its H-alpha chromospheric flux, for its age, mass and metallicity. As a byproduct of this work, we developed an automatic method to determine temperatures from the wings of H-alpha, which means the suppression of the error inherent to the visual procedure used in the literature.

preprint2006arXiv

On the difference between nuclear and contraction ages

Context:Ages derived from the low mass stars still contracting onto the main sequence often differ from ages derived from the high mass ones that have already evolved away from it. Aims:We investigate the general claim of disagreement between these two independent age determinations by presenting UBVRI photometry of the young galactic open clusters NGC 2232, NGC 2516, NGC 2547 and NGC 4755, spanning the age range ~10-150 Myr Methods: We derive reddenings, distances and nuclear ages by fitting ZAMS and isochrones to color-magnitudes and color-color diagrams. To derive contraction ages, we use four different pre-main sequence models, with an empirically calibrated color-temperature relation to match the Pleiades cluster sequence. Results: When using exclusively the V vs. V-I color-magnitude diagram and empirically calibrated isochrones, there is consistency between nuclear and contraction ages for the studied clusters. Although the contraction ages seem systematically underestimated, in none of cases they deviate by more than one standard deviation from the nuclear ages.