Source author record

Adam Frank

Adam Frank 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

63works
14topics
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

63 published item(s)

preprint2024arXiv

How negative feedback and the ambient environment limit the influence of recombination in common envelope evolution

We perform 3D hydrodynamical simulations to study recombination and ionization during the common envelope (CE) phase of binary evolution, and develop techniques to track the ionic transitions in time and space. We simulate the interaction of a $2\,M_\odot$ red giant branch primary and a $1\,M_\odot$ companion modeled as a particle. We compare a run employing a tabulated equation of state (EOS) that accounts for ionization and recombination, with a run employing an ideal gas EOS. During the first half of the simulations, $\sim15$ per cent more mass is unbound in the tabulated EOS run due to the release of recombination energy, but by simulation end the difference has become negligible. We explain this as being a consequence of (i) the tabulated EOS run experiences a shallower inspiral and hence smaller orbital energy release at late times because recombination energy release expands the envelope and reduces drag, and (ii) collision and mixing between expanding envelope gas, ejecta and circumstellar ambient gas assists in unbinding the envelope, but does so less efficiently in the tabulated EOS run where some of the energy transferred to bound envelope gas is used for ionization. The rate of mass unbinding is approximately constant in the last half of the simulations and the orbital separation steadily decreases at late times. A simple linear extrapolation predicts a CE phase duration of $\sim2\,\mathrm{yr}$, after which the envelope would be unbound.

preprint2023arXiv

The messy death of a multiple star system and the resulting planetary nebula as observed by JWST

Planetary nebulae (PNe), the ejected envelopes of red giant stars, provide us with a history of the last, mass-losing phases of 90 percent of stars initially more massive than the Sun. Here, we analyse James Webb Space Telescope (JWST) Early Release Observation (ERO) images of the PN NGC3132. A structured, extended H2 halo surrounding an ionised central bubble is imprinted with spiral structures, likely shaped by a low-mass companion orbiting the central star at 40-60 AU. The images also reveal a mid-IR excess at the central star interpreted as a dusty disk, indicative of an interaction with another, closer companion. Including the previously known, A-type visual companion, the progenitor of the NGC3132 PN must have been at least a stellar quartet. The JWST images allow us to generate a model of the illumination, ionisation and hydrodynamics of the molecular halo, demonstrating the power of JWST to investigate complex stellar outflows. Further, new measurements of the A-type visual companion allow us to derive the value for the mass of the progenitor of a central star to date with excellent precision: 2.86+/-0.06 Mo. These results serve as path finders for future JWST observations of PNe providing unique insight into fundamental astrophysical processes including colliding winds, and binary star interactions, with implications for supernovae and gravitational wave systems.

preprint2022arXiv

Detectability of Chlorofluorocarbons in the Atmospheres of Habitable M-dwarf Planets

The presence of chlorofluorocarbons (CFCs) in Earth's atmosphere is a direct result of technology. Ozone-depleting CFCs have been banned by most countries, but some CFCs have persistent in elevated concentrations due to their long stratospheric lifetimes. CFCs are effective greenhouse gases and could serve as a remotely detectable spectral signature of technology. Here we use a three-dimensional climate model and a synthetic spectrum generator to assess the detectability of CFC-11 and CFC-12 as a technosignature on exoplanets. We consider the case of TRAPPIST-1e as well as a habitable Earth-like planet around a 3300 K M-dwarf star, with CFC abundances ranging from one to five times present-day levels. Assuming an optimistic James Webb Space Telescope (JWST) Mid Infrared Instrument (MIRI) low resolution spectrometer (LRS) noise floor level of 10 ppm to multiple co-added observations, we find that spectral features potentially attributable to present or historic Earth-level CFC features could be detected with a SNR $\ge 3-5$ on TRAPPIST-1e, if present, in $\sim 100$ hours of in-transit time. However, applying a very conservative 50 ppm noise floor to co-added observations, even a 5x Earth-level CFC would not be detectable no matter the observation time. Such observations could be carried out simultaneously and at no additional cost with searches for biosignature gases. Non-detection would place upper limits on the CFC concentration. We find that with the launch of JWST, humanity may be approaching the cusp of being able to detect passive atmospheric technosignatures equal in strength to its own around the nearest stars.

preprint2022arXiv

Effects of Charge Exchange on the Evaporative Wind of HD 209458b

The role of charge exchange in shaping exoplanet photoevaporation remains a topic of contention. Exchange of electrons between stellar wind protons from the exoplanet's host star and neutral hydrogen from the planet's wind has been proposed as a mechanism to create "energetic neutral atoms" (ENAs), which could explain the high absorption line velocities observed in systems where mass loss is occurring. In this paper we present results from 3D hydrodynamic simulations of the mass loss of a planet similar to HD 209458b. We self-consistently launch a planetary wind by calculating the ionization and heating resulting from incident high-energy radiation, inject a stellar wind into the simulation, and allow electron exchange between the stellar and planetary winds. We predict the potential production of ENAs by the wind-wind interaction analytically, then present the results of our simulations, which confirm the analytic limits. Within the limits of our hydrodynamic simulation, we find that charge exchange with the stellar wind properties examined here is unable to explain the absorption observed at high Doppler velocities.

preprint2022arXiv

Excitation & Excavation of the Claws of the Southern Crab

We show that the Southern Crab (aka Hen2-104) presents an auspicious opportunity to study the form and speed of the invisible winds that excavate and shock the lobes of various types of bipolar nebulae associated with close and highly evolved binary stars. A deep three-color image overlay of Hen2-104 reveals that the ionization state of its lobe edges, or "claws", increases steadily from singly to doubly ionized values with increasing wall latitude. This "reverse" ionization pattern is unique among planetary nebulae (and similar objects) and incompatible with UV photoionization from a central source. We show that the most self-consistent explanation for the ionization pattern is shock ionization by a fast (~600 km s^-1) "tapered" stellar wind in which the speed and momentum flux of the wind increase with equatorial latitude. We present a hydrodynamic simulation that places the latitude-dependent form, the knotty walls, and the reverse ionization of the outer lobes of Hen2-104 into a unified context.

preprint2022arXiv

Jets from main sequence and white dwarf companions during common envelope evolution

It has long been speculated that jet feedback from accretion onto the companion during a common envelope (CE) event could affect the orbital evolution and envelope unbinding process, but this conjecture has heretofore remained largely untested. We present global 3D hydrodynamical simulations of CE evolution (CEE) that include a jet subgrid model and compare them with an otherwise identical model without a jet. Our binary consists of a $2M_\odot$ red giant branch primary and a $1M_\odot$ or $0.5M_\odot$ main sequence or white dwarf secondary companion modeled as a point particle. We run the simulations for 10 orbits (40 days). Our jet model adds mass at a constant rate $\dot{M}_\mathrm{j}$ of order the Eddington rate, with maximum velocity $v_\mathrm{j}$ of order the escape speed, to two spherical sectors with the jet axis perpendicular to the orbital plane, and supplies kinetic energy at the rate $\sim\dot{M}_\mathrm{j} v_\mathrm{j}^2/40$. We explore the influence of the jet on orbital evolution, envelope morphology and envelope unbinding, and assess the dependence of the results on jet mass-loss rate, launch speed, companion mass, opening angle, and whether or not subgrid accretion is turned on. In line with our theoretical estimates, we find that in all cases the jet becomes choked around the time of first periastron passage. We also find that jets lead to increases in unbound mass of up to $\sim10\%$, as compared to simulations which do not include a jet.

preprint2022arXiv

New Clues to the Evolution of Dwarf Carbon Stars From Their Variability and X-ray Emission

As main-sequence stars with C$>$O, dwarf carbon (dC) stars are never born alone but inherit carbon-enriched material from a former asymptotic giant branch (AGB) companion. In contrast to M dwarfs in post-mass transfer binaries, C$_2$ and/or CN molecular bands allow dCs to be identified with modest-resolution optical spectroscopy, even after the AGB remnant has cooled beyond detectability. Accretion of substantial material from the AGB stars should spin up the dCs, potentially causing a rejuvenation of activity detectable in X-rays. Indeed, a few dozen dCs have recently been found to have photometric variability with periods under a day. However, most of those are likely post-common-envelope binaries (PCEBs), spin-orbit locked by tidal forces, rather than solely spun-up by accretion. Here, we study the X-ray properties of a sample of the five nearest known dCs with $Chandra$. Two are detected in X-rays, the only two for which we also detected short-period photometric variability. We suggest that the coronal activity detected so far in dCs is attributable to rapid rotation due to tidal locking in short binary orbits after a common-envelope phase, late in the thermally pulsing (TP) phase of the former C-AGB primary (TP-AGB).

preprint2022arXiv

Opportunities for Technosignature Science in the Astro2020 Report

The Astro2020 report outlines numerous recommendations that could significantly advance technosignature science. Technosignatures refer to any observable manifestations of extraterrestrial technology, and the search for technosignatures is part of the continuum of the astrobiological search for biosignatures. The search for technosignatures is directly relevant to the "World and Suns in Context" theme and "Pathways to Habitable Worlds" program in the Astro2020 report. The relevance of technosignatures was explicitly mentioned in "E1 Report of the Panel on Exoplanets, Astrobiology, and the Solar System," which stated that "life's global impacts on a planet's atmosphere, surface, and temporal behavior may therefore manifest as potentially detectable exoplanet biosignatures, or technosignatures" and that potential technosignatures, much like biosignatures, must be carefully analyzed to mitigate false positives. The connection of technosignatures to this high-level theme and program can be emphasized, as the report makes clear the purpose is to address the question "Are we alone?" This question is also presented in the Explore Science 2020-2024 plan as a driver of NASA's mission. This white paper summarizes the potential technosignature opportunities within the recommendations of the Astro2020 report, should they be implemented by funding agencies. The objective of this paper is to demonstrate the relevance of technosignature science to a wide range of missions and urge the scientific community to include the search for technosignatures as part of the stated science justifications for the large and medium programs that include the Infrared/Optical/Ultraviolet space telescope, Extremely Large Telescopes, probe-class far-infrared and X-ray missions, and various facilities in radio astronomy.

preprint2022arXiv

The Case for Technosignatures: Why They May Be Abundant, Long-lived, Highly Detectable, and Unambiguous

The intuition suggested by the Drake equation implies that technology should be less prevalent than biology in the galaxy. However, it has been appreciated for decades in the SETI community that technosignatures could be more abundant, longer-lived, more detectable, and less ambiguous than biosignatures. We collect the arguments for and against technosignatures' ubiquity and discuss the implications of some properties of technological life that fundamentally differ from nontechnological life in the context of modern astrobiology: It can spread among the stars to many sites, it can be more easily detected at large distances, and it can produce signs that are unambiguously technological. As an illustration in terms of the Drake equation, we consider two Drake-like equations, for technosignatures (calculating N(tech)) and biosignatures (calculating N(bio)). We argue that Earth and humanity may be poor guides to the longevity term L and that its maximum value could be very large, in that technology can outlive its creators and even its host star. We conclude that while the Drake equation implies that N(bio)>>N(tech), it is also plausible that N(tech)>>N(bio). As a consequence, as we seek possible indicators of extraterrestrial life, for instance, via characterization of the atmospheres of habitable exoplanets, we should search for both biosignatures and technosignatures. This exercise also illustrates ways in which biosignature and technosignature searches can complement and supplement each other and how methods of technosignature search, including old ideas from SETI, can inform the search for biosignatures and life generally.

preprint2021arXiv

Contact Inequality -- First Contact Will Likely Be With An Older Civilization

First contact with another civilization, or simply another intelligence of some kind, will likely be quite different depending on whether that intelligence is more or less advanced than ourselves. If we assume that the lifetime distribution of intelligences follows an approximately exponential distribution, one might naively assume that the pile-up of short-lived entities dominates any detection or contact scenario. However, it is argued here that the probability of contact is proportional to the age of said intelligence (or possibly stronger), which introduces a selection effect. We demonstrate that detected intelligences will have a mean age twice that of the underlying (detected + undetected) population, using the exponential model. We find that our first contact will most likely be with an older intelligence, provided that the maximum allowed mean lifetime of the intelligence population, $τ_{\rm max}$, is $>e$ times larger than our own. Older intelligences may be rare but they disproportionality contribute to first contacts, introducing what we call a 'contact inequality', analogous to wealth inequality. This reasoning formalizes intuitional arguments and highlights that first contact would likely be one-sided, with ramifications for how we approach SETI.

preprint2021arXiv

Panchromatic HST/WFC3 Imaging Studies of Young, Rapidly Evolving Planetary Nebulae. I. NGC 6302

We present the results of a comprehensive, near-UV-to-near-IR Hubble Space Telescope WFC3 imaging study of the young planetary nebula (PN) NGC 6302, the archetype of the class of extreme bi-lobed, pinched-waist PNe that are rich in dust and molecular gas. The new WFC3 emission-line image suite clearly defines the dusty toroidal equatorial structure that bisects NGC 6302's polar lobes, and the fine structures (clumps, knots, and filaments) within the lobes. The most striking aspect of the new WFC3 image suite is the bright, S-shaped 1.64 micron [Fe II] emission that traces the southern interior of the east lobe rim and the northern interior of the west lobe rim, in point-symmetric fashion. We interpret this [Fe II] emitting region as a zone of shocks caused by ongoing, fast (~100 km/s), collimated, off-axis winds from NGC 6302's central star(s). The [Fe II] emission and a zone of dusty, N- and S-rich clumps near the nebular symmetry axis form wedge-shaped structures on opposite sides of the core, with boundaries marked by sharp azimuthal ionization gradients. Comparison of our new images with earlier HST/WFC3 imaging reveals that the object previously identified as NGC 6302's central star is a foreground field star. Shell-like inner lobe features may instead pinpoint the obscured central star's actual position within the nebula's dusty central torus. The juxtaposition of structures revealed in this HST/WFC3 imaging study of NGC 6302 presents a daunting challenge for models of the origin and evolution of bipolar PNe.

preprint2021arXiv

The Formation of Discs in the Interior of AGB Stars from the Tidal Disruption of Planets and Brown Dwarfs

A significant fraction of isolated white dwarfs host magnetic fields in excess of a MegaGauss. Observations suggest that these fields originate in interacting binary systems where the companion is destroyed thus leaving a singular, highly-magnetized white dwarf. In post-main-sequence evolution, radial expansion of the parent star may cause orbiting companions to become engulfed. During the common envelope phase, as the orbital separation rapidly decreases, low-mass companions will tidally disrupt as they approach the giant's core. We hydrodynamically simulate the tidal disruption of planets and brown dwarfs, and the subsequent accretion disc formation, in the interior of an asymptotic giant branch star. These dynamically formed discs are commensurate with previous estimates, suggesting strong magnetic fields may originate from these tidal disruption events.

preprint2020arXiv

Bipolar Planetary Nebulae from Outflow Collimation by Common Envelope Evolution

The morphology of bipolar planetary nebulae (PNe) can be attributed to interactions between a fast wind from the central engine and dense toroidal shaped ejecta left over from common envelope (CE) evolution. Here we use the 3-D hydrodynamic AMR code AstroBEAR to study the possibility that bipolar PN outflows can emerge collimated even from an uncollimated spherical wind in the aftermath of a CE event. The output of a single CE simulation via the SPH code PHANTOM serves as the initial conditions. Four cases of winds, all with high enough momenta to account for observed high momenta preplanetary nebula outflows, are injected spherically from the region of the CE binary remnant into the ejecta. We compare cases with two different momenta and cases with no radiative cooling versus application of optically thin emission via a cooling curve to the outflow. Our simulations show that in all cases highly collimated bipolar outflows result from deflection of the spherical wind via the interaction with the CE ejecta. Significant asymmetries between the top and bottom lobes are seen in all cases. The asymmetry is strongest for the lower momentum case with radiative cooling. While real post CE winds may be aspherical, our models show that collimation via "inertial confinement" will be strong enough to create jet-like outflows even beginning with maximally uncollimated drivers. Our simulations reveal detailed shock structures in the shock focused inertial confinement (SFIC) model and develop a lens-shaped inner shock that is a new feature of SFIC driven bipolar lobes.

preprint2020arXiv

Common Envelope Evolution on the Asymptotic Giant Branch: Unbinding within a Decade?

Common envelope (CE) evolution is a critical but still poorly understood progenitor phase of many high-energy astrophysical phenomena. Although 3D global hydrodynamic CE simulations have become more common in recent years, those involving an asymptotic giant branch (AGB) primary are scarce, due to the high computational cost from the larger dynamical range compared to red giant branch (RGB) primaries. But CE evolution with AGB progenitors is desirable to simulate because such events are the likely progenitors of most bi-polar planetary nebulae (PNe), and prominent observational testing grounds for CE physics. Here we present a high resolution global simulation of CE evolution involving an AGB primary and $1\,\mathrm{M}_\odot$ secondary, evolved for $20$ orbital revolutions. During the last $16$ of these orbits, the envelope unbinds at an almost constant rate of about $0.1$-$0.2\,\mathrm{M}_\odot\,\mathrm{yr}^{-1}$. If this rate were maintained, the envelope would be unbound in less than $10\,\mathrm{yr}$. The dominant source of this unbinding is consistent with inspiral; we assess the influence of the ambient medium to be subdominant. We compare this run with a previous run that used an RGB phase primary evolved from the same $2\,\mathrm{M}_\odot$ main sequence star to assess the influence of the evolutionary state of the primary. When scaled appropriately, the two runs are quite similar, but with some important differences.

preprint2020arXiv

Effects of Radiation Pressure on the Evaporative Wind of HD 209458b

The role of radiation pressure in shaping exoplanet photoevaporation remains a topic of contention. Radiation pressure from the exoplanet's host star has been proposed as a mechanism to drive the escaping atmosphere into a "cometary" tail and explain the high velocities observed in systems where mass loss is occurring. In this paper we present results from high-resolution 3-D hydrodynamic simulations of a planet similar to HD 209458b. We self-consistently launch a wind flowing outward from the planet by calculating the ionization and heating resulting from incident high-energy radiation, and account for radiation pressure. We first present a simplified calculation, setting a limit on the Lyman-$α$ flux required to drive the photo-evaporated planetary material to larger radii and line-of-sight velocities. We then present the results of our simulations, which confirm the limits determined by our analytic calculation. We thus demonstrate that, within the limits of our hydrodynamic simulation and for the Lyman-$α$ fluxes expected for HD 209458, radiation pressure is unlikely to significantly affect photoevaporative winds or to explain the high velocities at which wind material is observed, though further possibilities remain to be investigated.

preprint2020arXiv

The Fermi Paradox and the Aurora Effect: Exo-civilization Settlement, Expansion and Steady States

We model the settlement of the galaxy by space-faring civilizations in order to address issues related to the Fermi Paradox. We explore the problem in a way that avoids assumptions about the intent and motivation of any exo-civilization seeking to settle other planetary systems. We first consider the speed of an advancing settlement via probes of finite velocity and range to determine if the galaxy can become inhabited with space-faring civilizations on timescales shorter than its age. We also include the effect of stellar motions on the long term behavior of the settlement front which adds a diffusive component to its advance. The results of these models demonstrate that the Milky Way can be readily 'filled-in' with settled stellar systems under conservative assumptions about interstellar spacecraft velocities and launch rates. We then consider the question of the galactic steady-state achieved in terms of the fraction of settled planets. We do this by considering the effect of finite settlement civilization lifetimes on the steady states. We find a range of parameters for which the galaxy supports a population of interstellar space-faring civilizations even though some settleable systems are uninhabited. Both results point to ways in which Earth might remain unvisited in the midst of an inhabited galaxy. Finally we consider how our results can be combined with the finite horizon for evidence of previous settlements in Earth's geologic record. Our steady-state model can constrain the probabilities for an Earth visit by a settling civilization before a given time horizon. These results break the link between Hart's famous "Fact A" (no interstellar visitors on Earth now) and the conclusion that humans must, therefore, be the only technological civilization in the galaxy.

preprint2019arXiv

Models of the Mass-Ejection Histories of pre Planetary Nebulae, IV. Magnetized Winds and the Origins of Jets, Bullets, and FLIERs

The influences and consequences of toroidal magnetic fields in shaping the visible lobes of pre planetary nebulae ("prePNe") are explored in this, the last of a series of papers of parameter studies of prePN evolution. To probe these influences we start with the steady, diverging, and field-free wind model of our previous papers and add weak to moderate toroidal fields to the winds in order to generate arrays of outcomes after 500 y, after which the structures grow almost homologously. As expected, toroidal fields in the stellar winds invariably form very thin and dense axial features whose structure is best described as a thin cold jet with an ultra-dense and neutral leading knot, or "bullet", at its tip. The speed of the leading knot depends only on the density contrast (the ratio of injected to ambient gas densities at the nozzle) and wind injection speed, but not on the field strength or opening angle. The lobes formed by the ram pressure of the winds take a variety of forms and sizes that depend primarily on the geometric structure of the injected gas and the density contrast. About 20% of the HST images of prePNe show the unique signatures of shaping by toroidal fields. Pairs of low-ionization knots seen along the major axis of fully ionized PNe, often called "FLIERs" are easily explained as the very dense, cold , and neutral remnants of magnetically formed knots.

preprint2018arXiv

Accretion in Common Envelope Evolution

Common envelope evolution (CEE) occurs in some binary systems involving asymptotic giant branch (AGB) or red giant branch (RGB) stars, and understanding this process is crucial for understanding the origins of various transient phenomena. CEE has been shown to be highly asymmetrical and global 3D simulations are needed to help understand the dynamics. We perform and analyze hydrodynamic CEE simulations with the adaptive mesh refinement (AMR) code AstroBEAR, and focus on the role of accretion onto the companion star. We bracket the range of accretion rates by comparing a model that removes mass and pressure using a subgrid accretion prescription with one that does not. Provided a pressure-release valve, such as a bipolar jet, is available, super-Eddington accretion could be common. Finally, we summarize new results pertaining to the energy budget, and discuss the overall implications relating to the feasibility of unbinding the envelope in CEE simulations.

preprint2016arXiv

Molecular cloud formation in high-shear, magnetized colliding flows

The colliding flows (CF) model is a well-supported mechanism for generating molecular clouds. However, to-date most CF simulations have focused on the formation of clouds in the normal-shock layer between head-on colliding flows. We performed simulations of magnetized colliding flows that instead meet at an oblique-shock layer. Oblique shocks generate shear in the post-shock environment, and this shear creates inhospitable environments for star formation. As the degree of shear increases (i.e. the obliquity of the shock increases), we find that it takes longer for sink particles to form, they form in lower numbers, and they tend to be less massive. With regard to magnetic fields, we find that even a weak field stalls gravitational collapse within forming clouds. Additionally, an initially oblique collision interface tends to reorient over time in the presence of a magnetic field, so that it becomes normal to the oncoming flows. This was demonstrated by our most oblique shock interface, which became fully normal by the end of the simulation.

preprint2016arXiv

The Creation of AGB Fallback Shells

The possibility that mass ejected during Asymptotic Giant Branch (AGB) stellar evolution phases falls back towards the star has been suggested in applications ranging from the formation of accretion disks to the powering of late-thermal pulses. In this paper, we seek to explicate the properties of fallback flow trajectories from mass-loss events. We focus on a transient phase of mass ejection with sub-escape speeds, followed by a phase of a typical AGB wind. We solve the problem using both hydrodynamic simulations and a simplified one-dimensional analytic model that matches the simulations. For a given set of initial wind characteristics, we find a critical shell velocity that distinguishes between "shell fallback" and "shell escape". We discuss the relevance of our results for both single and binary AGB stars. In particular, we discuss how our results help to frame further studies of fallback as a mechanism for forming the substantial population of observed post-AGB stars with dusty disks.

preprint2016arXiv

The inner cavity of the circumnuclear disc

The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and second, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here.

preprint2016arXiv

Uncertain for A Century: Quantum Mechanics and The Dilemma of Interpretation

Quantum Mechanics, the physical theory describing the microworld, represents one of science's greatest triumphs. It lies at the root of all modern digital technologies and offers unparalleled correspondence between prediction and experiments. Remarkably, however, after more than 100 years it is still unclear what quantum mechanics means in terms of basic philosophical questions about the nature of reality. While there are many interpretations of the mathematical machinery of quantum physics, there remains no experimental means to distinguish between most of them. In this contribution, (based on a discussion at the NYAS), I wish to consider the ways in which the enduring lack of an agreed upon interpretation of quantum physics influences a number of critical philosophical debates about physics and reality. I briefly review two problems effected by quantum interpretations: the meaning of the term 'Universe' and the nature of consciousness. In what follows I am explicitly not advocating for any particular quantum interpretation. Instead, I am interested in how the explicit inability of modern physics to experimentally distinguish between interpretations with wildly divergent ontological/epistemological implications plays into discussions of physics and its description of the world.

preprint2015arXiv

A New Empirical Constraint on the Prevalence of Technological Species in the Universe

In this paper we address the cosmic frequency of technological species. Recent advances in exoplanet studies provide strong constraints on all astrophysical terms in the Drake Equation. Using these and modifying the form and intent of the Drake equation we show that we can set a firm lower bound on the probability that one or more additional technological species have evolved anywhere and at any time in the history of the observable Universe. We find that as long as the probability that a habitable zone planet develops a technological species is larger than ~$10^{-24}$, then humanity is not the only time technological intelligence has evolved. This constraint has important scientific and philosophical consequences.

preprint2014arXiv

MHD Effects on Pulsed YSO Jets. I. 2.5-D Simulations

In this paper we explore the dynamics of radiative axisymmetric MHD jets at high resolution using AMR methods. The goal of the study is to determine both the dynamics and emission properties of such jets. To that end we have implemented microphysics enabling us to produce synthetic maps of H$α$ and [S II]. The jets are pulsed either sinusoidally or randomly via a time-dependent ejection velocity which leads to a complicated structure of internal shocks and rarefactions as has been seen in previous simulations. The high resolution of our simulations allows us to explore in great detail the effect of pinch forces (due to the jet's toroidal magnetic field) within the "working surfaces" where pulses interact. We map the strong H$α$ emission marking shock fronts and the strong [S II] emission inside cooling regions behind shocks as observed with high-resolution images of jets. We find that pinch forces in the stronger field cases produce additional emission regions along the axis as compared with purely hydrodynamic runs. These simulations are a first step to understanding the full 3-D emission properties of radiative MHD jets.

preprint2014arXiv

Numerical Simulations of Mach Stem Formation via Intersecting Bow Shocks

Hubble Space Telescope observations show bright knots of H$α$ emission within outflowing young stellar jets. Velocity variations in the flow create secondary bow shocks that may intersect and lead to enhanced emission. When the bow shocks intersect at or above a certain critical angle, a planar shock called a Mach stem is formed. These shocks could produce brighter H$α$ emission since the incoming flow to the Mach stem is parallel to the shock normal. In this paper we report first results of a study using 2-D numerical simulations designed to explore Mach stem formation at the intersection of bow shocks formed by hypersonic "bullets" or "clumps". Our 2-D simulations show how the bow shock shapes and intersection angles change as the adiabatic index $γ$ changes. We show that the formation or lack of a Mach stem in our simulations is consistent with the steady-state Mach stem formation theory. Our ultimate goal, which is part of an ongoing research effort, is to characterize the physical and observational consequences of bow shock intersections including the formation of Mach stems.

preprint2014arXiv

On the role of ambient environments in the collapse of Bonnor-Ebert spheres

We consider the interaction between a marginally stable Bonnor-Ebert (BE) sphere and the surrounding ambient medium. In particular, we explore how the infall from an evolving ambient medium can trigger collapse of the sphere using 3-D adaptive mesh refinement simulations. We find the resulting collapse dynamics to vary considerably with ambient density. In the highest ambient density cases, infalling material drives a strong compression wave into the cloud. It is the propagation of this wave through the cloud interior that triggers the subsequent collapse. For lower ambient densities, we find the main trigger of collapse to be a quasistatic adjustment of the BE sphere to gravitational settling of the ambient gas. In all cases, we find that the classic "outside-in" collapse mode for super-critical BE spheres is recovered before a protostar (i.e. sink particle) forms. Our work supports scenarios in which BE dynamics naturally begins with either a compression-wave or infall dominated phase, and only later assumes the usual outside-in collapse behavior.

preprint2014arXiv

Triggered Star Formation and Its Consequences

Star formation can be triggered by compression from wind or supernova driven shock waves that sweep over molecular clouds. Because these shocks will likely contain processed elements, triggered star formation has been proposed as an explanation for short lived radioactive isotopes (SLRI) in the Solar System. Previous studies have tracked the triggering event to the earliest phases of collapse and have focused on the shock properties required for both successful star formation and mixing of SLRI's. In this paper, we use Adaptive Mesh Refinement (AMR) simulation methods, including sink particles, to simulate the full collapse and subsequent evolution of a stable Bonnor-Ebert sphere subjected to a shock and post-shock wind. We track the flow of the cloud material after a star (a sink particle) has formed. For non-rotating clouds we find robust triggered collapse and little bound circumstellar material remaining around the post-shock collapsed core. When we add initial cloud rotation we observe the formation of disks around the collapsed core which then interact with the post-shock flow. Our results indicate that these circumstellar disks are massive enough to form planets and are long-lived, in spite of the ablation driven by post-shock flow ram pressure. As a function of the initial conditions, we also track the time evolution of the accretion rates and particle mixing between between the ambient wind and cloud material. The latter is maximized for cases of highest mach number.

preprint2013arXiv

Efficient Parallelization for AMR MHD Multiphysics Calculations; Implementation in AstroBEAR

Current Adaptive Mesh Refinement (AMR) simulations require algorithms that are highly parallelized and manage memory efficiently. As compute engines grow larger, AMR simulations will require algorithms that achieve new levels of efficient parallelization and memory management. We have attempted to employ new techniques to achieve both of these goals. Patch or grid based AMR often employs ghost cells to decouple the hyperbolic advances of each grid on a given refinement level. This decoupling allows each grid to be advanced independently. In AstroBEAR we utilize this independence by threading the grid advances on each level with preference going to the finer level grids. This allows for global load balancing instead of level by level load balancing and allows for greater parallelization across both physical space and AMR level. Threading of level advances can also improve performance by interleaving communication with computation, especially in deep simulations with many levels of refinement. While we see improvements of up to 30% on deep simulations run on a few cores, the speedup is typically more modest (5-20%) for larger scale simulations. To improve memory management we have employed a distributed tree algorithm that requires processors to only store and communicate local sections of the AMR tree structure with neighboring processors. Using this distributed approach we are able to get reasonable scaling efficiency (> 80%) out to 12288 cores and up to 8 levels of AMR - independent of the use of threading.

preprint2013arXiv

Mass loading and knot formation in AGN jets by stellar winds

Jets from active galaxies propagate from the central black hole out to the radio lobes on scales of hundreds of kiloparsecs. The jets may encounter giant stars with strong stellar winds and produce observable signatures. For strong winds and weak jets, the interaction may truncate the jet flow during its transit via the mass loading. For weaker jets, the interaction can produce knots in the jet. We present recent 3DMHD numerical simulations to model the evolution of this jet-wind interaction and its observational consequences. We explore (i) the relative mechanical luminosity of the radio jets and the stellar winds (ii) the impact parameter between the jets' axis and the stellar orbital path (iii) the relative magnetic field strength of the jets and the stellar winds.

preprint2013arXiv

MHD Shock-Clump Evolution with Self-Contained Magnetic Fields

We study the interaction of strong shock waves with magnetized clumps. Previous numerical work focused on the simplified scenario in which shocked clumps are immersed in a globally uniform magnetic field that extends through both the clump and the ambient medium. Here we consider the complementary circumstance in which the field is completely self-contained within the clumps. This could arise naturally during clump formation via dynamical or thermal instabilities for example as magnetic field pinches off from the ambient medium. Using our AMR MHD code AstroBEAR, we carry out a series of simulations with magnetized clumps that have different self-contained magnetic field configurations. We find that the clump and magnetic evolution are sensitive to the fraction of magnetic field aligned with versus perpendicular to the shock normal. The relative strength of magnetic pressure and tension in the different field configurations allows us to analytically understand the different cases of post-shock evolution. We also show how turbulence and the mixing it implies depends of the initial field configuration and suggest ways in which observed shock-clump morphology may be used as a proxy for identifying internal field topologies a posteriori.

preprint2013arXiv

Outflows from Evolved Stars: The Rapidly Changing Fingers of CRL618

Our ultimate goal is to probe the nature of the collimator of the outflows in the pre PN CRL618. CRL618 is uniquely suited for this purpose owing to its multiple, bright, and carefully studied finger-shaped outflows east and west of its nucleus. We compare new HST images to images in the same filters observed as much as 11 y previously to uncover large proper motions and surface brightness changes in its multiple finger-shaped outflows. The expansion age of the ensemble of fingers is close to 100y. We find strong brightness variations at the fingertips during the past decade. Deep IR images reveal a multiple ring- like structure of the surrounding medium into which the outflows propagate and interact. Tightly constrained three-dimensional ("3D") hydrodynamic models link the properties of the fingers to their possible formation histories. We incorporate previously published complementary information to discern whether each of the fingers of CRL618 are the results of steady, collimated outflows or a brief ejection event that launched a set of bullets about a century ago. Finally, we argue on various physical grounds that fingers of CRL618 are likely to be the result of a spray of clumps ejected at the nucleus of CRL618 since any mechanism that form a sustained set of unaligned jets is unprecedented.

preprint2013arXiv

Size of discs formed by wind accretion in binaries can be underestimated if the role of wind-driving force is ignored

Binary systems consisting of a secondary accreting form a wind-emitting primary are ubiquitous in astrophysics. The phenomenology of such Bondi-Hoyle-Lyttleton (BHL) accretors is particularly rich when an accretion disc forms around the secondary. The outer radius of such discs is commonly estimated from the net angular momentum produced by a density variation of material across the BHL or Bondi accretion cylinder, as the latter is tilted with respect to the direction to the primary due to orbital motion. But this approach has ignored the fact that the wind experiences an outward driving force that the secondary does not. In actuality, the accretion stream falls toward a retarded point in the secondary's orbit as the secondary is pulled toward the primary relative to the stream. The result is a finite separation or "accretion stream impact parameter" (ASIP) separating the secondary and stream. When the orbital radius a_o exceeds the BHL radius r_b, the ratio of outer disc radius estimated as the ASIP to the conventional estimate a_o^{1/2}/r_b^{1/2}>1. We therefore predict that discs will form at larger radii from the secondary than traditional estimates. This agrees with the importance of the ASIP emphasized by Huarte-Espinosa et al. (2013) and the practical consequence that resolving the initial outer radius of such an accretion disc in numerical simulations can be less demanding than what earlier estimates would suggest.

preprint2013arXiv

Sustainability and the Astrobiological Perspective: Framing Human Futures in a Planetary Context

We explore how questions related to developing a sustainable human civilization can be cast in terms of astrobiology. In particular we show how ongoing astrobiological studies of the coupled relationship between life, planets and their co-evolution can inform new perspectives and direct new studies in sustainability science. Using the Drake Equation as a vehicle to explore the gamut of astrobiology, we focus on its most import factor for sustainability: the mean lifetime <L> of an ensemble of Species with Energy-Intensive Technology (SWEIT). We then cast the problem into the language of dynamical system theory and introduce the concept of a trajectory bundle for SWEIT evolution and discuss how astrobiological results usefully inform the creation of dynamical equations, their constraints and initial conditions. Three specific examples of how astrobiological considerations can be folded into discussions of sustainability are discussed: (1) concepts of planetary habitability, (2) mass extinctions and their possible relation to the current, so-called Anthropocene epoch, and (3) today's changes in atmospheric chemisty (and the climate change it entails) in the context of pervious epochs of biosphere-driven atmospheric and climate alteration (i.e. the Great Oxidation Event).

preprint2013arXiv

The Effects of Inhomogeneities within Colliding Flows on the Formation and Evolution of Molecular Clouds

Observational evidence from local star-forming regions mandates that star formation occurs shortly after, or even during, molecular cloud formation. Models of the formation of molecular clouds in large-scale converging flows have identified the physical mechanisms driving the necessary rapid fragmentation. They also point to global gravitational collapse driving supersonic turbulence in molecular clouds. Previous cloud formation models have focused on turbulence generation, gravitational collapse, magnetic fields, and feedback. Here, we explore the effect of structure in the flow on the resulting clouds and the ensuing gravitational collapse. We compare two extreme cases, one with a collision between two smooth streams, and one with streams containing small clumps. We find that structured converging flows lead to a delay of local gravitational collapse ("star formation"). Thus, more gas has time to accumulate, eventually leading to a strong global collapse, and thus to a high star formation rate. Uniform converging flows fragment hydrodynamically early on, leading to the rapid onset of local gravitational collapse and an overall low sink formation rate.

preprint2013arXiv

The Formation and Evolution of Wind-Capture Disks In Binary Systems

We study the formation, evolution and physical properties of accretion disks formed via wind capture in binary systems. Using the AMR code AstroBEAR, we have carried out high resolution 3D simulations that follow a stellar mass secondary in the co-rotating frame as it orbits a wind producing AGB primary. We first derive a resolution criteria, based on considerations of Bondi-Hoyle flows, that must be met in order to properly resolve the formation of accretion disks around the secondary. We then compare simulations of binaries with three different orbital radii (10, 15, 20 AU). Disks are formed in all three cases, however the size of the disk and, most importantly, its accretion rate decreases with orbital radii. In addition, the shape of the orbital motions of material within the disk becomes increasingly elliptical with increasing binary separation. The flow is mildly unsteady with "fluttering" around the bow shock observed. The disks are generally well aligned with the orbital plane after a few binary orbits. We do not observe the presence of any large scale, violent instabilities (such as the flip-flop mode). For the first time, moreover, it is observed that the wind component that is accreted towards the secondary has a vortex tube-like structure, rather than a column-like one as it was previously thought. In the context of AGB binary systems that might be precursors to Pre-Planetary and Planetary Nebula, we find that the wind accretion rates at the chosen orbital separations are generally too small to produce the most powerful outflows observed in these systems if the companions are main sequence stars but marginally capable if the companions are white dwarfs. It is likely that many of the more powerful PPN and PN involve closer binaries than the ones considered here. The results also demonstrate principles of broad relevance to all wind-capture binary systems.

preprint2012arXiv

Comparing Poynting flux dominated magnetic tower jets with kinetic-energy dominated jets

Magnetic Towers represent one of two fundamental forms of MHD outflows. Driven by magnetic pressure gradients, these flows have been less well studied than magneto-centrifugally launched jets even though magnetic towers may well be as common. Here we present new results exploring the behavior and evolution of magnetic tower outflows and demonstrate their connection with pulsed power experimental studies and purely hydrodynamic jets which might represent the asymptotic propagation regimes of magneto-centrifugally launched jets. High-resolution AMR MHD simulations (using the AstroBEAR code) provide insights into the underlying physics of magnetic towers and help us constrain models of their propagation. Our simulations have been designed to explore the effects of thermal energy losses and rotation on both tower flows and their hydro counterparts. We find these parameters have significant effects on the stability of magnetic towers, but mild effects on the stability of hydro jets. Current-driven perturbations in the Poynting Flux Dominated (PDF) towers are shown to be amplified in both the cooling and rotating cases. Our studies of the long term evolution of the towers show that the formation of weakly magnetized central jets within the tower are broken up by these instabilities becoming a series of collimated clumps which magnetization properties vary over time. In addition to discussing these results in light of laboratory experiments, we address their relevance to astrophysical observations of young star jets and outflow from highly evolved solar type stars.

preprint2012arXiv

From Bipolar to Elliptical: Simulating the Morphological Evolution of Planetary Nebulae

The majority of Proto-planetary nebulae (PPN) are observed to have bipolar morphologies. The majority of mature PN are observed to have elliptical shapes. In this paper we address the evolution of PPN/PN morphologies attempting to understand if a transition from strongly bipolar to elliptical shape can be driven by changes in the parameters of the mass loss process. To this end we present 2.5D hydrodynamical simulations of mass loss at the end stages of stellar evolution for intermediate mass stars. We track changes in wind velocity, mass loss rate and mass loss geometry. In particular we focus on the transition from mass loss dominated by a short duration jet flow (driven during the PPN phase) to mass loss driven by a spherical fast wind (produced by the central star of the PN). We address how such changes in outflow characteristics can change the nebula from a bipolar to an elliptical morphology. Our results show that including a period of jet formation in the temporal sequence of PPN to PN produces realistic nebular synthetic emission geometries. More importantly such a sequence provides insight, in principle, into the apparent difference in morphology statistics characterizing PPN and PN systems. In particular we find that while jet driven PPN can be expected to be dominated by bipolar morphologies, systems that begin with a jet but are followed by a spherical fast wind will evolve into elliptical nebulae. Furthermore, we find that spherical nebulae are highly unlikely to ever derive from either bipolar PPN or elliptical PN.

preprint2012arXiv

Magnetic Towers and Binary-Formed Disks: New Results for PN Evolution

We present new results of 3-D AMR MHD simulations focusing on two distinct aspects of PPN evolution. We first report new simulations of collimated outflows driven entirely by magnetic fields. These Poynting flux dominated "magnetic towers" hold promise for explaining key properties of PPN flows. Our simulations address magnetic tower evolution and stability. We also present results of a campaign of simulations to explore the development of accretion disks formed via wind capture. Our result focus on the limits of disk formation and the range of disk properties.

preprint2012arXiv

On the structure and stability of magnetic tower jets

Modern theoretical models of astrophysical jets combine accretion, rotation, and magnetic fields to launch and collimate supersonic flows from a central source. Near the source, magnetic field strengths must be large enough to collimate the jet requiring that the Poynting flux exceeds the kinetic-energy flux. The extent to which the Poynting flux dominates kinetic energy flux at large distances from the engine distinguishes two classes of models. In magneto-centrifugal launch (MCL) models, magnetic fields dominate only at scales $\lesssim 100$ engine radii, after which the jets become hydrodynamically dominated (HD). By contrast, in Poynting flux dominated (PFD) magnetic tower models, the field dominates even out to much larger scales. To compare the large distance propagation differences of these two paradigms, we perform 3-D ideal MHD AMR simulations of both HD and PFD stellar jets formed via the same energy flux. We also compare how thermal energy losses and rotation of the jet base affects the stability in these jets. For the conditions described, we show that PFD and HD exhibit observationally distinguishable features: PFD jets are lighter, slower, and less stable than HD jets. Unlike HD jets, PFD jets develop current-driven instabilities that are exacerbated as cooling and rotation increase, resulting in jets that are clumpier than those in the HD limit. Our PFD jet simulations also resemble the magnetic towers that have been recently created in laboratory astrophysical jet experiments.

preprint2011arXiv

Consequences of Magnetic Field Structure for Heat Transport in Magnetohydrodynamics

Interfaces between hot and cold magnetized plasmas exist in various astrophysical contexts, for example where hot outflows impinge on an ambient interstellar medium (ISM). It is of interest to understand how the structure of the magnetic field spanning the interface affects the temporal evolution of the temperature gradient. Here we explore the relation between the magnetic field topology and the heat transfer rate by adding various fractions of tangled vs. ordered field across a hot-cold interface allow the system to evolve to a steady state. We find a simple mathematical relation for the rate of heat conduction as a function of the initial ratio of ordered to tangled field across the interface. We discuss potential implications for the astrophysical context of magnetized wind blown bubbles (WBB) around evolved stars.

preprint2011arXiv

Efficient Parallelization for AMR MHD Multiphysics Calculations; Implementation in AstroBEAR

Current AMR simulations require algorithms that are highly parallelized and manage memory efficiently. As compute engines grow larger, AMR simulations will require algorithms that achieve new levels of efficient parallelization and memory management. We have attempted to employ new techniques to achieve both of these goals. Patch or grid based AMR often employs ghost cells to decouple the hyperbolic advances of each grid on a given refinement level. This decoupling allows each grid to be advanced independently. In AstroBEAR we utilize this independence by threading the grid advances on each level with preference going to the finer level grids. This allows for global load balancing instead of level by level load balancing and allows for greater parallelization across both physical space and AMR level. Threading of level advances can also improve performance by interleaving communication with computation, especially in deep simulations with many levels of refinement. To improve memory management we have employed a distributed tree algorithm that requires processors to only store and communicate local sections of the AMR tree structure with neighboring processors.

preprint2011arXiv

Laboratory experiments and simulations on jets

Astrophysical jets have been studied with observations, theoretical models and numerical simulations for decades. Recently, supersonic magnetized jets have been formed in laboratory experiments of high-energy density plasmas. I will review these studies and discuss the experimental setup that has been used to form millimeter-scale jets driven by strong toroidal magnetic fields in a MAGPIE generator. The physical conditions of these experiments are such that they can be scaled to astrophysical scenarios. These laboratory jets provide insights on the underlying physics of magnetic tower jets and help constrain some models of astrophysical jets. In this context, we also discuss the connection between the laboratory jets and recent 3D-MHD numerical simulations of Poynting flux dominated jets. The simulations allow us to investigate the effects of thermal energy losses and base rotation on the growth rate of kink mode perturbations, and to compare the evolution of PFD jets with a hydrodynamic counterpart of the same energy flux.

preprint2011arXiv

Modeling Poynting flux vs. kinetic-energy dominated jets

We present 3D-MHD AMR simulations of Poynting flux dominated (PFD) jets formed by injection of magnetic energy. We compare their evolution with a hydrodynamic jet which is formed by injecting kinetic energy with the same energy flux than the PFD jets. We predict characteristic emission distributions for each of these jets. Current-driven perturbations in PFD jets are amplified by both cooling and rotation for the regimes studied: Shocks and thermal pressure support are weakened by cooling, making the jets more susceptible to kinking. Rotation amplifies the toroidal magnetic field which also exacerbates the kink instability.

preprint2011arXiv

The Evolution of Hetergeneous "Clumpy Jets": A Parameter Study

We investigate the role discrete clumps embedded in an astrophysical jet play on the jet's morphology and line emission characteristics. By varying clumps' size, density, position, and velocity, we cover a range of parameter space motivated by observations of objects such as the Herbig Haro object HH~34. We here extend the results presented in Yirak et al. 2009, including how analysis of individual observations may lead to spurious sinusoidal variation whose parameters vary widely over time, owing chiefly to interacts between clumps. The goodness of the fits, while poor in all simulations, are best when clump-clump collisions are minimal. Our results indicate that a large velocity dispersion leads to a clump-clump collision-dominated flow which disrupts the jet beam. Finally, we present synthetic emission images of H-$α$ and [SII] and note an excess of [SII] emission along the jet length as compared to observations. This suggests that observed beams undergo earlier processing, if they are present at all.

preprint2011arXiv

The Illumination and Growth of CRL 2688: An Analysis of New & Archival HST Observations

We present four-color images of CRL 2688 obtained in 2009 using the Wide-Field Camera 3 on HST. The F606W image is compared with archival images in very similar filters to monitor the proper motions of nebular structure. We find that the bright N-S lobes have expanded uniformly by 2.5% and that the ensemble of rings has translated radially by 0.07 in 6.65 y. The rings were ejected every 100y for ~4 millennia until the lobes formed 250y ago. Starlight scattered from the edges of the dark E-W dust lane is coincident with extant H2 images and leading tips of eight pairs of CO outflows. We interpret this as evidence that fingers lie within geometrically opposite cones of opening angles {\approx} 30{\circ} like those in CRL618. By combining our results of the rings with 12CO absorption from the extended AGB wind we ascertain that the rings were ejected at ~18 km s-1 with very little variation and that the distance to CRL2688, v_{exp}$ / ${\dotθ}_exp$, is 300 - 350 pc. Our 2009 imaging program included filters that span 0.6 to 1.6μm. We constructed a two-dimensional dust scattering model of stellar radiation through CRL2688 that successfully reproduces the details of the nebular geometry, its integrated spectral energy distribution, and nearly all of its color variations. The model implies that the optical opacity of the lobes >~ 1, the dust particle density in the rings decreases as radius^{-3} and that the mass and momentum of the AGB winds and their rings have increased over time.

preprint2010arXiv

Explosive Outflows Powered by the Decay of Non-Hierarchical Multiple Systems of Massive Stars: Orion BN/KL

The explosive BN/KL outflow emerging from OMC1 behind the Orion Nebula may have been powered by the dynamical decay of a non-hierarchical multiple system $\sim$500 years ago that ejected the massive stars I, BN, and source n, with velocities of about 10 to 30 km s$^{-1}$. New proper motion measurements of H$_2$ features show that within the errors of measurement, the outflow originated from the site of stellar ejection. Combined with published data, these measurements indicate an outflow age of $\sim$500 years, similar to the time since stellar ejection. The total kinetic energy of the ejected stars and the outflow is about 2 to $6 \times 10^{47}$ ergs. It is proposed that the gravitational potential energy released by the formation of a short-period binary, most likely source I, resulted in stellar ejection and powered the outflow. A scenario is presented for the formation of a compact, non-hierarchical multiple star system, its decay into an ejected binary and two high-velocity stars, and launch of the outflow. Three mechanisms may have contributed to the explosion in the gas: (i) Unbinding of the circum-cluster envelope following stellar ejection, (ii) disruption of circumstellar disks and high-speed expulsion of the resulting debris during the final stellar encounter, and (iii) the release of stored magnetic energy. Plausible proto-stellar disk end envelope properties can produce the observed outflow mass, velocity, and kinetic energy distributions. The ejected stars may have acquired new disks by fall-back or Bondi-Hoyle accretion with axes roughly orthogonal to their velocities. The expulsion of gas and stars from OMC1 may have been driven by stellar interactions.

preprint2010arXiv

From Bipolar to Elliptical: Morphological Changes in the Temporal Evolution of PN

Proto-planetary nebulae (pPN) and planetary nebulae (PN) seem to be formed by interacting winds from asymptotic giant branch (AGB) stars. The observational issue that most pPN are bipolar but most older PN are elliptical is addressed. We present 2.5D hydrodynamical numerical simulations of episodic cooling interacting winds to investigate the long term evolution of PN morphologies. We track wind acceleration, decrease in mass-loss and episodic change in wind geometry from spherical (AGB) to collimated (pPN) and back to spherical again (PN). This outflow sequence is found to produce realistic PN dynamics and morphological histories. Effects from different AGB distributions and jet duty cycles are also investigated.

preprint2010arXiv

Isotropically Driven versus Outflow Driven Turbulence: Observational Consequences for Molecular Clouds

Feedback from protostellar outflows can influence the nature of turbulence in star forming regions even if they are not the primary source of velocity dispersion for all scales of molecular clouds. For the rate and power expected in star forming regions, we previously (Carroll et al. 2009) demonstrated that outflows could drive supersonic turbulence at levels consistent with the scaling relations from Matzner 2007 although with a steeper velocity power spectrum than expected for an isotropically driven supersonic turbulent cascade. Here we perform higher resolution simulations and combine simulations of outflow driven turbulence with those of isotropically forced turbulence. We find that the presence of outflows within an ambient isotropically driven turbulent environment produces a knee in the velocity power spectrum at the outflow scale and a steeper slope at sub-outflow scales than for a purely isotropically forced case. We also find that the presence of outflows flattens the density spectrum at large scales effectively reducing the formation of large scale turbulent density structures. These effects are qualitatively independent of resolution. We have also carried out Principal Component Analysis (PCA) for synthetic data from our simulations. We find that PCA as a tool for identifying the driving scale of turbulence has a misleading bias toward low amplitude large scale velocity structures even when they are not necessarily the dominant energy containing scales. This bias is absent for isotropically forced turbulence but manifests strongly for collimated outflow driven turbulence.

preprint2010arXiv

Modelling magnetically dominated and radiatively cooling jets

Using 3D-MHD Eulerian-grid numerical simulations, we study the formation and evolution of rising magnetic towers propagating into an ambient medium. The towers are generated from a localized injection of pure magnetic energy. No rotation is imposed on the plasma. We compare the evolution of a radiatively cooling tower with an adiabatic one, and find that both bend due to pinch instabilities. Collimation is stronger in the radiative cooling case; the adiabatic tower tends to expand radially. Structural similarities are found between these towers and the millimeter scale magnetic towers produced in laboratory experiments.

preprint2009arXiv

Laboratory Astrophysics and the State of Astronomy and Astrophysics

Laboratory astrophysics and complementary theoretical calculations are the foundations of astronomy and astrophysics and will remain so into the foreseeable future. The impact of laboratory astrophysics ranges from the scientific conception stage for ground-based, airborne, and space-based observatories, all the way through to the scientific return of these projects and missions. It is our understanding of the under-lying physical processes and the measurements of critical physical parameters that allows us to address fundamental questions in astronomy and astrophysics. In this regard, laboratory astrophysics is much like detector and instrument development at NASA, NSF, and DOE. These efforts are necessary for the success of astronomical research being funded by the agencies. Without concomitant efforts in all three directions (observational facilities, detector/instrument development, and laboratory astrophysics) the future progress of astronomy and astrophysics is imperiled. In addition, new developments in experimental technologies have allowed laboratory studies to take on a new role as some questions which previously could only be studied theoretically can now be addressed directly in the lab. With this in mind we, the members of the AAS Working Group on Laboratory Astrophysics, have prepared this State of the Profession Position Paper on the laboratory astrophysics infrastructure needed to ensure the advancement of astronomy and astrophysics in the next decade.

preprint2009arXiv

Magnetic Nested-wind Scenarios for Bipolar Outflows: Pre-planetary and YSO nebular shaping

We present results of a series of magnetohydrodynamic (MHD) and hydro- dynamic (HD) 2.5D simulations of the morphology of outflows driven by nested wide-angle winds - i.e. winds which eminate from a central star as well as from an orbiting accretion disk. While our results are broadly relevent to nested wind systems we have tuned the parameters of the simulations to touch on issues in both Young Stellar Objects and Planetary Nebula studies. In particular our studies connect to open issues in the early evolution of Planetary Nebulae. We find that nested MHD winds exhibit marked morphological differences from the single MHD wind case along both dimensions of the flow. Nested HD winds on the other hand give rise mainly to geometric distortions of an outflow that is topologically similar to the flow arising from a single stellar HD wind. Our MHD results are insensitive to changes in ambient temperature between ionized and un-ionized circumstellar environments. The results are sensitive to the relative mass-loss rates, and to the relative speeds of the stellar and disk winds. We also present synthetic emission maps of both nested MHD and HD simulations. We find that nested MHD winds show knots of emission appearing on-axis that do not appear in the HD case.

preprint2009arXiv

New Discoveries in Cosmology and Fundamental Physics through Advances in Laboratory Astrophysics

As the Cosmology and Fundamental Physics (CFP) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the early universe, the microwave background, the reionization and galaxy formation up to virialization of protogalaxies, large scale structure, the intergalactic medium, the determination of cosmological parameters, dark matter, dark energy, tests of gravity, astronomically determined physical constants, and high energy physics using astronomical messengers. Central to the progress in these areas are the corresponding advances in laboratory astrophysics which are required for fully realizing the CFP scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics which produce the observed astrophysical processes. The 5 areas of laboratory astrophysics which we have identified as relevant to the CFP panel are atomic, molecular, plasma, nuclear, and particle physics. Here, Section 2 describes some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the CFP scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

preprint2009arXiv

New Discoveries in Galaxies across Cosmic Time through Advances in Laboratory Astrophysics

As the Galaxies across Cosmic Time (GCT) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the formation, evolution, and global properties of galaxies and galaxy clusters, as well as active galactic nuclei and QSOs, mergers, star formation rate, gas accretion, and supermassive black holes. Central to the progress in these areas are the corresponding advances in laboratory astrophysics that are required for fully realizing the GCT scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produce the observed astrophysical processes. The 5 areas of laboratory astrophysics that we have identified as relevant to the CFP panel are atomic, molecular, solid matter, plasma, nuclear, and particle physics. In this white paper, we describe in Section 2 some of the new scientific opportunities and compelling scientific themes that will be enabled by advances in laboratory astrophysics. In Section 3, we provide the scientific context for these opportunities. Section 4 briefly discusses some of the experimental and theoretical advances in laboratory astrophysics required to realize the GCT scientific opportunities of the next decade. As requested in the Call for White Papers, Section 5 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 6.

preprint2009arXiv

New Discoveries in Stars and Stellar Evolution through Advances in Laboratory Astrophysics

As the Stars and Stellar Evolution (SSE) panel is fully aware, the next decade will see major advances in our understanding of these areas of research. To quote from their charge, these advances will occur in studies of the Sun as a star, stellar astrophysics, the structure and evolution of single and multiple stars, compact objects, SNe, gamma-ray bursts, solar neutrinos, and extreme physics on stellar scales. Central to the progress in these areas are the corresponding advances in laboratory astrophysics, required to fully realize the SSE scientific opportunities within the decade 2010-2020. Laboratory astrophysics comprises both theoretical and experimental studies of the underlying physics that produces the observed astrophysical processes. The 6 areas of laboratory astrophysics, which we have identified as relevant to the CFP panel, are atomic, molecular, solid matter, plasma, nuclear physics, and particle physics. In this white paper, we describe in Section 2 the scientific context and some of the new scientific opportunities and compelling scientific themes which will be enabled by advances in laboratory astrophysics. In Section 3, we discuss some of the experimental and theoretical advances in laboratory astrophysics required to realize the SSE scientific opportunities of the next decade. As requested in the Call for White Papers, Section 4 presents four central questions and one area with unusual discovery potential. Lastly, we give a short postlude in Section 5.

preprint2009arXiv

Self-Convergence of Radiatively Cooling Clumps

Numeric convergence studies demonstrate that the evolution of an adiabatic clump is well-captured by roughly 100 cells per clump radius. The presence of radiative cooling, however, imposes limits on the problem due to the removal of thermal energy. Numerical studies which include radiative cooling typically adopt the 100--200 cells per clump radius resolution. In this paper we present the results of a convergence study for radiatively cooling clumps undertaken over a broad range of resolutions, from 12 to 1,536 cells per clump radius, employing adaptive mesh refinement (AMR) in a 2D axisymmetric geometry ("2.5D"). We also provide a fully 3D simulation, at 192 cells per clump radius, which supports our 2.5D results. We find no appreciable self-convergence at ~100 cells per clump radius as small-scale differences owing to increasingly resolving the "cooling length" have global effects. We therefore conclude that self-convergence is an insufficient criterion to apply on its own when addressing the question of sufficient resolution for radiatively cooled shocked clump simulations. We suggest the adoption of alternate criteria to support a statement of sufficient resolution, such as the demonstration of adequate resolution of the cooling layers behind shocks. We discuss an associated refinement criteria for AMR codes.

preprint2009arXiv

Serendipitous XMM-Newton Detection of X-ray Emission from the Bipolar Planetary Nebula Hb 5

We report the serendipitous detection by the XMM-Newton X-ray Observatory of an X-ray source at the position of the Type I (He- and N-rich) bipolar planetary nebula Hb 5. The Hb 5 X-ray source appears marginally resolved. While the small number of total counts (~170) and significant off-axis angle of the X-ray source (~7.8') precludes a definitive spatial analysis, the morphology of the X-ray emission appears to trace the brightest features seen in optical images of Hb 5. The X-ray spectrum is indicative of a thermal plasma at a temperature between 2.4 and 3.7 MK and appears to display strong Neon emission. The inferred X-ray luminosity is L_X = 1.5 x 10^32 ergs/s. These results suggest that the detected X-ray emission is dominated by shock-heated gas in the bipolar nebula, although we cannot rule out the presence of a point-like component at the position of the central star. The implications for and correspondence with current models of shock-heated gas in planetary nebulae is discussed.

preprint2008arXiv

Hypersonic Buckshot: Astrophysical Jets as Heterogeneous Collimated Plasmoids

Herbig-Haro (HH) jets are commonly thought of as homogeneous beams of plasma traveling at hypersonic velocities. Structure within jet beams is often attributed to periodic or ``pulsed'' variations of conditions at the jet source. Simulations based on this scenario result in knots extending across the jet diameter. Observations and recent high energy density laboratory experiments shed new light on structures below this scale and indicate they may be important for understanding the fundamentals of jet dynamics. In this paper we offer an alternative to ``pulsed'' models of protostellar jets. Using direct numerical simulations we explore the possibility that jets are chains of sub-radial clumps propagating through a moving inter-clump medium. Our models explore an idealization of this scenario by injecting small ($r<r_{jet}$), dense ($ρ>ρ_{jet}$) spheres embedded in an otherwise smooth inter-clump jet flow. The spheres are initialized with velocities differing from the jet velocity by $\sim15$%. We find the consequences of shifting from homogeneous to heterogeneous flows are significant as clumps interact with each other and with the inter-clump medium in a variety of ways. Structures which mimic what is expected from pulsed-jet models can form, as can previously unseen ``sub-radial'' behaviors including backward facing bow shocks and off-axis working surfaces. While these small-scale structures have not been seen before in simulation studies, they are found in high resolution jet observations. We discuss implications of our simulations for the interpretation of protostellar jets with regard to characterization of knots by a ``lifetime'' or ``velocity history'' approach as well as linking observed structures with central engines which produce the jets.

preprint2008arXiv

Outflow Driven Turbulence in Molecular Clouds

In this paper we explore the relationship between protostellar outflows and turbulence in molecular clouds. Using 3-D numerical simulations we focus on the hydrodynamics of multiple outflows interacting within a parsec scale volume. We explore the extent to which transient outflows injecting directed energy and momentum into a sub-volume of a molecular cloud can be converted into random turbulent motions. We show that turbulence can readily be sustained by these interactions and show that it is possible to broadly characterize an effective driving scale of the outflows. We compare the velocity spectrum obtained in our studies to that of isotropically forced hydrodynamic turbulence finding that in outflow driven turbulence a power law is indeed achieved. However we find a steeper spectrum (beta ~ 3) is obtained in outflow driven turbulence models than in isotropically forced simulations (beta ~ 2). We discuss possible physical mechanisms responsible for these results as well and their implications for turbulence in molecular clouds where outflows will act in concert with other processes such as gravitational collapse.

preprint2001arXiv

Dynamos in Asymptotic-Giant-Branch Stars As the Origin of Magnetic Fields Shaping Planetary Nebulale

Planetary nebulae are thought to be formed when a slow wind from the progenitor giant star is overtaken by a subsequent fast wind generated as the star enters its white dwarf stage$^{1}$. A shock forms near the boundary between the winds, which creates a relatively dense shell that provides the characteristic appearance of a planetary nebula. A spherically symmetric wind will produce a spherically symmetric shell, yet over half of known planetary nebulae are not spherical; rather, they are elliptical or bipolar in shape$^{2}$. While a magnetic field could launch and collimate a bipolar outflow, the origin of such a field has hitherto been unclear, as previous work suggested that a field could not be generated${^3}$. Here we show that an asymptotic-giant-branch (AGB) star can indeed generate a strong magnetic field, in a dynamo at the interface between a rapidly rotating core and the more slowly rotating envelope of the star. The field is strong enough to shape the bipolar outflows that produce the observed bipolar planetary nebulae. Magnetic braking of the stellar core during this process may also explain the puzzlingly$^{4}$ slow rotation of most white dwarf stars.

preprint2000arXiv

The Magnetic Geometry of Pulsed Astrophysical Jets

Hypersonic, highly collimated, mass outflows ({\it jets}) are a ubiquitous phenomena in astrophysics. While the character of the jets differ, many exhibit some form of quasi-periodic clumping indicating the jet source is episodic or {\it pulsed}. The presence of pulsed jets in so many astrophysical contexts suggests a common formation mechanism. Such a process seems to have been found in {\it Magneto-centrifugal launching}, the combination of magnetic and centrifugal forces that occurs when a magnetized gaseous accretion disk orbits a central gravitating source. Observations of strong magnetic fields in jets are, however, rare or indirect. Thus the presence and effects of magnetic fields in YSO jets remains an unresolved issue of the highest importance. In this letter we focus on what should be expected of the structure of the fields in pulsed YSO jets. We show that combining velocity variability with an initial field configuration consistent with collimated, Magneto-centrifugally launched jets leads to a clear set of predictions concerning the geometry and relative strength of the magnetic field components in evolving YSO (and perhaps other) jets.

preprint1999arXiv

The MHD Kelvin-Helmholtz Instability III: The Role of Sheared Magnetic Field in Planar Flows

We have carried out simulations of the nonlinear evolution of the magnetohydrodynamic (MHD) Kelvin-Helmholtz (KH) instability for compressible fluids in $2\frac{1}{2}$-dimensions, extending our previous work by Frank et al (1996) and Jones \etal (1997). In the present work we have simulated flows in the x-y plane in which a ``sheared'' magnetic field of uniform strength ``smoothly'' rotates across a thin velocity shear layer from the z direction to the x direction, aligned with the flow field. We focus on dynamical evolution of fluid features, kinetic energy dissipation, and mixing of the fluid between the two layers, considering their dependence on magnetic field strength for this geometry. The introduction of magnetic shear can allow a Cat's Eye-like vortex to form, even when the field is stronger than the nominal linear instability limit given above. For strong fields that vortex is asymmetric with respect to the preliminary shear layer, however, so the subsequent dissipation is enhanced over the uniform field cases of comparable field strength. In fact, so long as the magnetic field achieves some level of dynamical importance during an eddy turnover time, the asymmetries introduced through the magnetic shear will increase flow complexity, and, with that, dissipation and mixing. The degree of the fluid mixing between the two layers is strongly influenced by the magnetic field strength. Mixing of the fluid is most effective when the vortex is disrupted by magnetic tension during transient reconnection, through local chaotic behavior that follows.

preprint1997arXiv

Outflow collimation in young stellar objects

In this paper we explore the effect of radiative losses on purely hydrodynamic jet collimation models applicable to Young Stellar Objects (YSOs). In our models aspherical bubbles form from the interaction of a central YSO wind with an aspherical circum-protostellar density distribution. Building on a previous non-radiative study (Frank & Mellema 1996) we demonstrate that supersonic jets are a natural and robust consequence of aspherical wind-blown bubble evolution. The simulations show that the addition of radiative cooling makes the hydrodynamic collimation mechanisms studied by Frank & Mellema (1996) more effective. An unexpected result is the production of cool, dense jets forming through conical converging flows at the poles of the bubbles. For steady winds the formation of these jets occurs early in the bubble evolution. At later times we find that the dynamical and cooling time scales for the jet material become similar. Since observations of YSO jets show considerable variability in the jet beam we present a simple one-dimensional (1-D) model for the evolution of a variable wind interacting with an accreting environment. These models indicate that, in the presence of a varying protostellar wind, the hydrodynamic collimation processes studied in our simulations can produce cool jets with sizes and time scales consistent with observations.

preprint1995arXiv

Numerical Magnetohydrodynamics in Astrophysics: Algorithm and Tests for Multi-Dimensional Flow

We present for astrophysical use a multi-dimensional numerical code to solve the equations for ideal magnetohydrodynamics (MHD). It is based on an explicit finite difference method on an Eulerian grid, called the Total Variation Diminishing (TVD) scheme, which is a second-order-accurate extension of the Roe-type upwind scheme. Multiple spatial dimensions are treated through a Strang-type operator splitting. The constraint of a divergence-free field is enforced exactly by calculating a correction via a gauge transformation in each time step. Results from two-dimensional shock tube tests show that the code captures correctly discontinuities in all three MHD waves families as well as contact discontinuities. The numerical viscosities and resistivity in the code, which are useful in order to understand simulations involving turbulent flows, are estimated through the decay of two-dimensional linear waves. Finally, the robustness of the code in two-dimensions is demonstrated through calculations of the Kelvin-Helmholtz instability and the Orszag-Tang vortex.