Researcher profile

Alexander Hubbard

Alexander Hubbard contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
22works
0followers
9topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

22 published item(s)

preprint2016arXiv

Ferromagnetism and particle collisions: applications to protoplanetary disks and the meteoritical record

The meteoritical record shows both iron partitioning and tungsten isotopic partitioning between matrix and chondrules. Tungsten is not abundant enough to have driven its own isotopic partitioning, but if tungsten were correlated with iron, then ferromagnetic interactions grains could help explain both observations. We derive a practical parameterization for the increase in particle-particle collision rates caused by mutually attracting particle magnetic dipole moments. While the appropriate magnetic parameters remain uncertain, we show that ambient magnetic fields in protoplanetary disks are expected to be strong enough to magnetize iron metal bearing dust grains sufficiently to drive large increases in their collision rates. Such increased collision rates between iron metal rich grains could help preserve primordial iron and W isotopic inhomogeneities; and would help explain why the meteoritical record shows their partitioning in the solar nebula. The importance of magnetic interactions for larger grains whose growth is balanced by fragmentation is less clear, and will require future laboratory or numerical studies.

preprint2016arXiv

Generating potassium abundance variations in the Solar Nebula

An intriguing aspect of chondritic meteorites is that they are complementary: while their separate components have wildly varying abundances, bulk chondrites have nearly solar composition. This implies that the nearly-solar reservoirs in which chondrites were born were in turn assembled from sub-reservoirs of differing compositions that birthed the different components. We focus on explaining the potassium abundance variations between chondrules even within a single chondrite, while maintaining the observed CI $^{41}$K to $^{39}$K ratios. This requires physically separating potassium and chondrules while the temperature is high enough for K to be in the gas phase. We examine several mechanisms which could drive the dust through gas and show that to do so locally would have required long (sub-orbital to many orbits) time scales; with shortest potassium depletion time scales occurring in a scenario where chondrules formed high above the midplane and settled out of the evaporated potassium. While orbital time scales are at odds with laboratory chondrule cooling rate estimates, any other model for the origin for the potassium abundance variation has to wrestle with the severe logistical difficulty of generating a plethora of correlated reservoirs which varied strongly in their potassium abundances, but not in their potassium isotope ratios.

preprint2016arXiv

Partitioning tungsten between matrix precursors and chondrule precursors through relative settling

Recent studies of chondrites have found a tungsten isotopic anomaly between chondrules and matrix. Given the refractory nature of tungsten, this implies that W was carried into the solar nebula by at least two distinct families of pre-solar grains. The observed chondrule/matrix split requires that the distinct families were kept separate during the dust coagulation process, and that the two families of grain interacted with the chondrule formation mechanism differently. We take the co-existence of different families of solids in the same general orbital region at the chondrule-precursor size as given, and explore the requirements for them to have interacted with the chondrule formation process at significantly different rates. We show that this sorting of families of solids into chondrule and matrix destined dust had to have been at least as powerful a sorting mechanism as the relative settling of aerodynamically distinct grains at at least two scale heights above the midplane. The requirement that the chondrule formation mechanism was correlated in some fashion with a dust grain sorting mechanism argues strongly for spatially localized chondrule formation mechanisms such as turbulent dissipation in non-thermally ionized disk surface layers, and argues against volume filling mechanisms such as planetesimal bow shocks.

preprint2015arXiv

Compound Chondrules fused Cold

About 4-5% of chondrules are compound: two separate chondrules stuck together. This is commonly believed to be the result of the two component chondrules having collided shortly after forming, while still molten. This allows high velocity impacts to result in sticking. However, at T ~ 1100K, the temperature below which chondrules collide as solids (and hence usually bounce), coalescence times for droplets of appropriate composition are measured in tens of seconds. Even at 1025K, at which temperature theory predicts that the chondrules must have collided extremely slowly to have stuck together, the coalescence time scale is still less than an hour. These coalescence time scales are too short for the collision of molten chondrules to explain the observed frequency of compound chondrules. We suggest instead a scenario where chondrules stuck together in slow collisions while fully solid; and the resulting chondrule pair was subsequently briefly heated to a temperature in the range of 900-1025K. In that temperature window the coalescence time is finite but long, covering a span of hours to a decade. This is particularly interesting because those temperatures are precisely the critical window for thermally ionized MRI activity, so compound chondrules provide a possible probe into that vital regime.

preprint2015arXiv

Photophoresis in a Dilute, Optically Thick Medium and Dust Motion in Protoplanetary Disks

We derive expressions for the photophoretic force on opaque spherical particles in a dilute gas in the optically thick regime where the radiation field is in local thermal equilibrium. Under those conditions, the radiation field has a simple form, leading to well defined analytical approximations for the photophoretic force that also consider both the internal thermal conduction within the particle, and the effects of heat conduction and radiation to the surrounding gas. We derive these results for homogeneous spherical particles; and for the double layered spheres appropriate for modeling solid grains with porous aggregate mantles. Then, as a specific astrophysical application of these general physical results, we explore the parameter space relevant to the photophoresis driven drift of dust in protoplanetary disks. We show that highly porous silicate grains have sufficiently low thermal conductivities that photophoretic effects, such as significant relative velocities between particles with differing porosity or levitation above the midplane, are expected to occur.

preprint2015arXiv

Turbulent Thermal Diffusion: A Way to Concentrate Dust in Protoplanetary Discs

Turbulence acting on mixes of gas and particles generally evenly diffuses the latter through the former. However, in the presence of background gas temperature gradients a phenomenon known as turbulent thermal diffusion appears as a particle drift velocity (rather than a diffusive term). This process moves particles from hot regions to cold ones. We rederive turbulent thermal diffusion using astrophysical language and demonstrate that it could play a major role in protoplanetary discs by concentrating particles by factors of tens. Such a concentration would set the stage for collective behavior such as the streaming instability and hence planetesimal formation.

preprint2014arXiv

Convective Overstability in radially stratified accretion disks under thermal relaxation

This letter expands the stability criterion for radially stratified, vertically {unstratified} accretion disks incorporating thermal relaxation. We find a linear amplification of epicyclic oscillations in these disks that depends on the effective cooling time, i.e. an overstability. The growth rates of the overstability vanish for both extreme cases, e.g. infinite cooling time and instantaneous cooling, i.e. the adiabatic and fully isothermal cases. However, for thermal relaxation times $τ$ on the order of the orbital frequency, $τΩ\sim 1$, modes grow at a rate proportional to the square of the Brunt-Väisälä frequency. The overstability is based on epicyclic motions, with the thermal relaxation causing gas to heat while radially displaced inwards, and cool while radially displaced outwards. This causes the gas to have a lower density when moving outwards compared to when it moves inwards, so it feels the outwards directed pressure force more strongly on that leg of the journey. We suggest the term ``Convective Overstability" for the phenomenon that has already been numerically studied in the non-linear regime in the context of amplifying vortices in disks, under the name ``Subcritical Baroclinic Instability". The point of the present paper is to make clear that vortex formation in three-dimensional disks is neither subcritical, i.e. does not need a finite perturbation, nor is it baroclinic in the sense of geophysical fluid dynamics, which requires on vertical shear. We find that Convective Overstability is a linear instability that will operate under a wide range of physical conditions for circumstellar disks.

preprint2014arXiv

Explaining Mercury's Density through Magnetic Erosion

In protoplanetary disks, dust grains rich in metallic iron can attract each other magnetically. If they are magnetized to values near saturation, the magnetically induced collision speeds are high enough to knock off the non-magnetized, loosely bound silicates. This process enriches the surviving portions of the dust grains in metallic iron, which further enhances the magnetically mediated collisions. The magnetic enhancement to the collisional cross-section between the iron rich dust results in rapid grain growth leading to planetesimal formation. While this process of knocking off silicates, which we term "magnetic erosion", occurs only in a very limited portion of a protoplanetary disk, it is a possible explanation for Mercury's disproportionately large iron core.

preprint2014arXiv

Protoplanetary dust porosity and FU Orionis Outbursts: Solving the mystery of Earth's missing volatiles

The Earth is known to be depleted in volatile lithophile elements in a fashion that defies easy explanation. We resolve this anomaly with a model that combines the porosity of collisionally grown dust grains in protoplanetary disks with heating from FU Orionis events that dramatically raise protoplanetary disk temperatures. The heating from an FU Orionis event alters the aerodynamical properties of the dust while evaporating the volatiles. This causes the dust to settle, abandoning those volatiles. The success of this model in explaining the elemental composition of the Earth is a strong argument in favor of highly porous collisionally grown dust grains in protoplanetary disks outside our Solar System. Further, it demonstrates how thermal (or condensation based) alterations of dust porosity, and hence aerodynamics, can be a strong factor in planet formation, leading to the onset of rapid gravitational instabilities in the dust disk and the subsequent collapse that forms planetesimals.

preprint2014arXiv

Radial Stresses and Energy Transport in Accretion Disks

Early in the study of viscous accretion disks it was realized that energy transfers from distant sources must be important, not least because the flow at the disk midplane in the bulk of the disk is likely outwards, out of the gravitational potential well. If the source of the viscosity is powered by accretion, such as in the case of the magneto-rotational instability, such distant energy sources must lie in the innermost regions of the disk, where accretion occurs even at the midplane. We argue here that modulations in this energy supply can alter the accretion rate on dynamical, rather than far longer viscous, time scales. This means that both the steady state value of and fluctuations in the inner disk's accretion rate, depending on the details of the inner boundary condition and occurring on the inner disk's rapid evolution time, can affect the outer disk. This is particularly interesting because observations have shown that disk accretion is not steady (e.g.~EX Lupi type objects). We also note that the power supplied to shearing boxes is set by the boxes themselves rather than the physical energy fluxes in a global disk. That is, their saturated magnetic field is not subject to the full set of energy constraints present in an actual disk. Our analysis suggests that large scale radial transport of energy has a critical impact on the evolution and variability of accretion disks.

preprint2014arXiv

Ribbons characterize magnetohydrodynamic magnetic fields better than lines: a lesson from dynamo theory

Blackman & Brandenburg argued that magnetic helicity conservation in dynamo theory can in principle be captured by diagrams of mean field dynamos when the magnetic fields are represented by ribbons or tubes, but not by lines. Here we present such a schematic ribbon diagram for the $α^2$ dynamo that tracks magnetic helicity and provides distinct scales of large scale magnetic helicity, small scale magnetic helicity, and kinetic helicity involved in the process. This also motivates our construction of a new "2.5 scale" minimalist generalization of the helicity-evolving equations for the α^2 dynamo that separately allows for these three distinct length scales while keeping only two dynamical equations. We solve these equations and, as in previous studies, find that the large scale field first grows at a rate independent of the magnetic Reynolds number R_M before quenching to an R_M dependent regime. But we also show that the larger the ratio of the wavenumber where the small scale current helicity resides to that of the forcing scale, the earlier the non-linear dynamo quenching occurs, and the weaker the large scale field is at the turnoff from linear growth. The harmony between the theory and the schematic diagram exemplifies a general lesson that magnetic fields in MHD are better visualized as two-dimensional ribbons (or pairs of lines) rather than single lines.

preprint2014arXiv

Semarkona: Lessons for chondrule and chondrite formation

We consider the evidence presented by the LL3.0 chondrite Semarkona, including its chondrule fraction, chondrule size distribution and matrix thermal history. We show that no more than a modest fraction of the ambient matrix material in the Solar Nebula could have been melted into chondrules; and that much of the unprocessed matrix material must have been filtered out at some stage of Semarkona's parent body formation process. We conclude that agglomerations of many chondrules must have formed in the Solar Nebula, which implies that chondrules and matrix grains had quite different collisional sticking parameters. Further, we note that the absence of large melted objects in Semarkona means that chondrules must have exited the melting zone rapidly, before the chondrule agglomerations could form. The simplest explanation for this rapid exit is that chondrule melting occurred in surface layers of the disk. The newly formed, compact, chondrules then settled out of those layers on short time scales.

preprint2014arXiv

Temperature Fluctuations driven by Magnetorotational Instability in Protoplanetary Disks

The magnetorotational instability (MRI) drives magnetized turbulence in sufficiently ionized regions of protoplanetary disks, leading to mass accretion. The dissipation of the potential energy associated with this accretion determines the thermal structure of accreting regions. Until recently, the heating from the turbulence has only been treated in an azimuthally averaged sense, neglecting local fluctuations. However, magnetized turbulence dissipates its energy intermittently in current sheet structures. We study this intermittent energy dissipation using high resolution numerical models including a treatment of radiative thermal diffusion in an optically thick regime. Our models predict that these turbulent current sheets drive order unity temperature variations even where the MRI is damped strongly by Ohmic resistivity. This implies that the current sheet structures where energy dissipation occurs must be well resolved to correctly capture the flow structure in numerical models. Higher resolutions are required to resolve energy dissipation than to resolve the magnetic field strength or accretion stresses. The temperature variations are large enough to have major consequences for mineral formation in disks, including melting chondrules, remelting calcium-aluminum rich inclusions, and annealing silicates; and may drive hysteresis: current sheets in MRI active regions could be significantly more conductive than the remainder of the disk.

preprint2013arXiv

Mineral Processing by Short Circuits in Protoplanetary Disks

Meteoritic chondrules were formed in the early solar system by brief heating of silicate dust to melting temperatures. Some highly refractory grains (Type B calcium-aluminum-rich inclusions, CAIs) also show signs of transient heating. A similar process may occur in other protoplanetary disks, as evidenced by observations of spectra characteristic of crystalline silicates. One possible environment for this process is the turbulent magnetohydrodynamic flow thought to drive accretion in these disks. Such flows generally form thin current sheets, which are sites of magnetic reconnection, and dissipate the magnetic fields amplified by a disk dynamo. We suggest that it is possible to heat precursor grains for chondrules and other high-temperature minerals in current sheets that have been concentrated by our recently described short-circuit instability. We extend our work on this process by including the effects of radiative cooling, taking into account the temperature dependence of the opacity; and by examining current sheet geometry in three-dimensional, global models of magnetorotational instability. We find that temperatures above 1600 K can be reached for favorable parameters that match the ideal global models. This mechanism could provide an efficient means of tapping the gravitational potential energy of the protoplanetary disk to heat grains strongly enough to form high-temperature minerals. The volume-filling nature of turbulent magnetic reconnection is compatible with constraints from chondrule-matrix complementarity, chondrule-chondrule complementarity, the occurrence of igneous rims, and compound chondrules. The same short-circuit mechanism may perform other high-temperature mineral processing in protoplanetary disks such as the production of crystalline silicates and CAIs.

preprint2013arXiv

Turbulence Induced Collision Velocities and Rates between Different Sized Dust Grains

We study the collision rates and velocities for point-particles of different sizes in turbulent flows. We construct fits for the collision rates at specified velocities (effectively a collisional velocity probability distribution) for particle stopping time ratios up to four; already by that point the collisional partners are very poorly correlated and so the results should be robust for even larger stopping time ratios. Significantly, we find that while particles of very different masses have approximately Maxwellian collisional statistics, as the mass ratio shrinks the distribution changes dramatically. At small stopping time ratios, the collisional partners are highly correlated and we find a population of high number density (clustered), low relative-velocity particle pairs. Unlike in the case of identical stopping time collisional partners, this low relative-velocity clustered population is collisional, but the clustering is barely adequate to trigger bulk effects such as the streaming instability. We conclude our analysis by constructing a master fit to the collisional statistics as a function only of the stopping time ratio. Together with our previous work for identical stopping time particle pairs, this provides a recipe for including collisional velocity probability distributions in dust coagulation models for protoplanetary disks. We also include our recipe for determining particle collisional diagnostics from numerical simulations.

preprint2012arXiv

Catastrophic quenching in alpha-Omega dynamos revisited

At large magnetic Reynolds numbers, magnetic helicity evolution plays an important role in astrophysical large-scale dynamos. The recognition of this fact led to the development of the dynamical alpha quenching formalism, which predicts catastrophically low mean fields in open systems. Here we show that in oscillatory alpha-Omega dynamos this formalism predicts an unphysical magnetic helicity transfer between scales. An alternative technique is proposed where this artifact is removed by using the evolution equation for the magnetic helicity of the total field in the shearing--advective gauge. In the traditional dynamical alpha quenching formalism, this can be described as an additional magnetic helicity flux of small-scale fields that does not appear in homogeneous alpha-squared dynamos. In alpha-Omega dynamos, the alternative formalism is shown to lead to larger saturation fields than previously obtained with the traditional formalism.

preprint2012arXiv

Short Circuits in Thermally Ionized Plasmas: A Mechanism for Intermittent Heating of Protoplanetary Disks

Many astrophysical systems of interest, including protoplanetary accretion disks, are made of turbu- lent magnetized gas with near solar metallicity. Thermal ionization of alkali metals in such gas exceeds non-thermal ionization when temperatures climb above roughly 1000 K. As a result, the conductiv- ity, proportional to the ionization fraction, gains a strong, positive dependence on temperature. In this paper, we demonstrate that this relation between the temperature and the conductivity triggers an exponential instability that acts similarly to an electrical short, where the increased conductivity concentrates the current and locally increases the Ohmic heating. This contrasts with the resistiv- ity increase expected in an ideal magnetic reconnection region. The instability acts to focus narrow current sheets into even narrower sheets with far higher currents and temparatures. We lay out the basic principles of this behavior in this paper using protoplanetary disks as our example host system, motivated by observations of chondritic meteorites and their ancestors, dust grains in protoplanetary disks, that reveal the existence of strong, frequent heating events that this instability could explain.

preprint2012arXiv

Turbulence induced collisional velocities and density enhancements: large inertial range results from shell models

To understand the earliest stages of planet formation, it is crucial to be able to predict the rate and the outcome of dust grains collisions, be it sticking and growth, bouncing, or fragmentation. The outcome of such collisions depends on the collision speed, so we need a solid understanding of the rate and velocity distribution of turbulence-induced dust grain collisions. The rate of the collisions depends both on the speed of the collisions and the degree of clustering experienced by the dust grains, which is a known outcome of turbulence. We evolve the motion of dust grains in simulated turbulence, an approach that allows a large turbulent inertial range making it possible to investigate the effect of turbulence on meso-scale grains (millimeter and centimeter). We find three populations of dust grains: one highly clustered, cold and collisionless; one warm; and the third "hot". Our results can be fit by a simple formula, and predict both significantly slower typical collisional velocities for a given turbulent strength than previously considered, and modest effective clustering of the collisional populations, easing difficulties associated with bouncing and fragmentation barriers to dust grain growth. Nonetheless, the rate of high velocity collisions falls off merely exponentially with relative velocity so some mid- or high-velocity collisions will still occur, promising some fragmentation.

preprint2011arXiv

Magnetic helicity transport in the advective gauge family

Magnetic helicity fluxes are investigated in a family of gauges in which the contribution from ideal magnetohydrodynamics takes the form of a purely advective flux. Numerical simulations of magnetohydrodynamic turbulence in this advective gauge family exhibit instabilities triggered by the build-up of unphysical irrotational contributions to the magnetic vector potential. As a remedy, the vector potential is evolved in a numerically well behaved gauge, from which the advective vector potential is obtained by a gauge transformation. In the kinematic regime, the magnetic helicity density evolves similarly to a passive scalar when resistivity is small and turbulent mixing is mild, i.e. when the fluid Reynolds number is not too large. In the dynamical regime, resistive contributions to the magnetic helicity flux in the advective gauge are found to be significant owing to the development of small length scales in the irrotational part of the magnetic vector potential.

preprint2011arXiv

The fratricide of alpha-Omega dynamos by their alpha-squared siblings

Context. Helically forced magneto-hydrodynamic shearing-sheet turbulence can support different large-scale dynamo modes, although the α - Ω mode is generally expected to dominate because it is the fastest growing. In an α - Ω dynamo, most of the field amplification is produced by the shear. As differential rotation is an ubiquitous source of shear in astrophysics, such dynamos are believed to be the source of most astrophysical large-scale magnetic fields. Aims. We study the stability of oscillatory migratory α - Ω type dynamos in turbulence simulations. Methods. We use shearing-sheet simulations of hydromagnetic turbulence that is helically forced at a wavenumber that is about three times larger than the lowest wavenumber in the domain so that both α - Ω and α2 dynamo action is possible. Results. After initial dominance and saturation, the α - Ω mode is found to be destroyed by an orthogonal α2 mode sustained by the helical turbulence alone. We show that there are at least two processes through which this transition can occur. Conclusions. The fratricide of α - Ω dynamos by its α2 sibling is discussed in the context of grand minima of solar and stellar activity. However, the genesis of α - Ω dynamos from an α2 dynamo has not yet been found.

preprint2010arXiv

Magnetic helicity flux in the presence of shear

Magnetic helicity has risen to be a major player in dynamo theory, with the helicity of the small-scale field being linked to the dynamo saturation process for the large-scale field. It is a nearly conserved quantity, which allows its evolution equation to be written in terms of production and flux terms. The flux term can be decomposed in a variety of fashions. One particular contribution that has been expected to play a significant role in dynamos in the presence of mean shear was isolated by Vishniac & Cho (2001, ApJ 550, 752). Magnetic helicity fluxes are explicitly gauge dependent however, and the correlations that have come to be called the Vishniac-Cho flux were determined in the Coulomb gauge, which turns out to be fraught with complications in shearing systems. While the fluxes of small-scale helicity are explicitly gauge dependent, their divergences can be gauge independent. We use this property to investigate magnetic helicity fluxes of small-scale field through direct numerical simulations in a shearing-box system and find that in a numerically usable gauge the divergence of the small-scale helicity flux vanishes, while the divergence of the Vishniac-Cho flux remains finite. We attribute this seeming contradiction to the existence of horizontal fluxes of small-scale magnetic helicity with finite divergences even in our shearing-periodic domain.

preprint2008arXiv

Identifying Deficiencies of Standard Accretion Disk Theory: Lessons from a Mean-Field Approach

Turbulent viscosity is frequently used in accretion disk theory to replace the microphysical viscosity in order to accomodate the observational need for in- stabilities in disks that lead to enhanced transport. However, simply replacing the microphysical transport coefficient by a single turbulent transport coeffi- cient hides the fact that the procedure should formally arise as part of a closure in which the hydrodynamic or magnetohydrodynamic equations are averaged, and correlations of turbulent fluctuations are replaced by transport coefficients. Here we show how a mean field approach leads quite naturally two transport coefficients, not one, that govern mass and angular momentum transport. In particular, we highlight that the conventional approach suffers from a seemingly inconsistent neglect of turbulent diffusion in the surface density equation. We constrain these new transport coefficients for specific cases of inward, outward, and zero net mass transport. In addition, we find that one of the new transport terms can lead to oscillations in the mean surface density which then requires a constant or small inverse Rossby number for disks to maintain a monotonic power-law surface density.