Source author record

T. W. Hartquist

T. W. Hartquist 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

21works
6topics
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

21 published item(s)

preprint2016arXiv

Ambipolar diffusion regulated collapse of filaments threaded by perpendicular magnetic fields

We numerically reproduce the density profiles for filaments that are in magnetohydrostatic and pressure equilibrium with their surroundings obtained in Tomisaka (2014) and show that these equilibria are dynamically stable. If the effect of ambipolar diffusion is considered, these filaments lose magnetic support initiating cloud collapse. The filaments do not lose magnetic flux. Rather the magnetic flux is redistributed within the filament from the centre towards the envelope. The rate of the collapse is inversely proportional to the fractional ionisation and two gravitationally-driven ambipolar diffusion regimes for the collapse are observed as predicted in Mouschovias & Morton (1991). For high values of the ionisation coefficient, that is $X \geq 10^{-7}$, the gas is strongly coupled to the magnetic field and the Jeans length is larger than the ambipolar diffusion length scale. Then the collapse is governed by magnetically-regulated ambipolar diffusion. For $X \lesssim 10^{-8}$, the gas is weakly coupled to the magnetic field and the magnetic support is removed by gravitationally-dominated ambipolar diffusion. Here, neutrals and ions only collide sporadically, that is the ambipolar diffusion length scale is larger than the Jeans length. When decaying turbulence is included, additional support is provided to the filament. This slows down the collapse of the filament even in the absence of a magnetic field. When a magnetic field is present, the collapse rate increases by a ratio smaller than for the non-magnetic case. This is because of a speed-up of the ambipolar diffusion due to larger magnetic field gradients generated by the turbulence and because the ambipolar diffusion aids the dissipation of turbulence below the ambipolar diffusion length scale. The highest increase in the rate is observed for the lowest ionisation coefficient and the highest turbulent intensity.

preprint2016arXiv

The interaction of hydrodynamic shocks with self-gravitating clouds

We describe the results of 3D simulations of the interaction of hydrodynamic shocks with Bonnor-Ebert spheres performed with an Adaptive Mesh Refinement code. The calculations are isothermal and the clouds are embedded in a medium in which the sound speed is either four or ten times that in the cloud. The strengths of the shocks are such that they induce gravitational collapse in some cases and not in others and we derive a simple estimate for the shock strength required for this to occur. These results are relevant to dense cores and Bok globules in star forming regions subjected to shocks produced by stellar feedback.

preprint2015arXiv

Gravitational instabilities in a protosolar-like disc I: dynamics and chemistry

To date, most simulations of the chemistry in protoplanetary discs have used 1+1D or 2D axisymmetric $α$-disc models to determine chemical compositions within young systems. This assumption is inappropriate for non-axisymmetric, gravitationally unstable discs, which may be a significant stage in early protoplanetary disc evolution. Using 3D radiative hydrodynamics, we have modelled the physical and chemical evolution of a 0.17 M$_{\odot}$ self-gravitating disc over a period of 2000 yr. The 0.8 M$_{\odot}$ central protostar is likely to evolve into a solar-like star, and hence this Class 0 or early Class I young stellar object may be analogous to our early Solar System. Shocks driven by gravitational instabilities enhance the desorption rates, which dominate the changes in gas-phase fractional abundances for most species. We find that at the end of the simulation, a number of species distinctly trace the spiral structure of our relatively low-mass disc, particularly CN. We compare our simulation to that of a more massive disc, and conclude that mass differences between gravitationally unstable discs may not have a strong impact on the chemical composition. We find that over the duration of our simulation, successive shock heating has a permanent effect on the abundances of HNO, CN and NH$_3$, which may have significant implications for both simulations and observations. We also find that HCO$^+$ may be a useful tracer of disc mass. We conclude that gravitational instabilities induced in lower mass discs can significantly, and permanently, affect the chemical evolution, and that observations with high-resolution instruments such as ALMA offer a promising means of characterising gravitational instabilities in protosolar discs.

preprint2014arXiv

Numerical Simulations of a Shock Interacting with Multiple Magnetized Clouds

We present 2D adiabatic magnetohydrodynamic (MHD) simulations of a shock interacting with groups of two or three cylindrical clouds. We study how the presence of a nearby cloud influences the dynamics of this interaction, and explore the resulting differences and similarities in the evolution of each cloud. The understanding gained from this small-scale study will help to interpret the behaviour of systems with many 10's or 100's of clouds. We observe a wide variety of behaviour in the interactions studied, which is dependent on the initial positions of the clouds and the orientation and strength of the magnetic field. We find: i) some clouds are stretched along their field-lines, whereas others are confined by their field-lines; ii) upstream clouds may accelerate past downstream clouds (though magnetic tension can prevent this); iii) clouds may also change their relative positions transverse to the direction of shock propagation as they "slingshot" past each other; iv) downstream clouds may be offered some protection from the oncoming flow as a result of being in the lee of an upstream cloud; v) the cycle of cloud compression and re-expansion is generally weaker when there are nearby neighbouring clouds; vi) the plasma $β$ in cloud material can vary rapidly as clouds collide with one another, but low values of $β$ are always transitory. This work is relevant to studies of multi-phase regions, where fast, low-density gas interacts with dense clouds, such as in circumstellar bubbles, supernova remnants, superbubbles and galactic winds.

preprint2013arXiv

Ram pressure stripping of the hot gaseous haloes of galaxies using the k-epsilon sub-grid turbulence model

We perform three dimensional hydrodynamic simulations of the ram pressure stripping of the hot extended gaseous halo of a massive galaxy using the k-epsilon sub-grid turbulence model at Mach numbers 0.9, 1.1 and 1.9. The k-epsilon model is used to simulate high Reynolds number flows by increasing the transport coefficients in regions of high turbulence. We find that the initial, instantaneous stripping is the same whether or not the k-epsilon model is implemented and is in agreement with the results of other studies. However the use of the k-epsilon model leads to five times less gas remaining after stripping by a supersonic flow has proceeded for 10 Gyr, which is more consistent with what simple analytic calculations indicate. Hence the continual Kelvin-Helmholtz stripping plays a significant role in the ram pressure stripping of the haloes of massive galaxies. To properly account for this, simulations of galaxy clusters will require the use of sub-grid turbulence models

preprint2013arXiv

Simulated Observations of Young Gravitationally Unstable Protoplanetary Discs

The formation and earliest stages of protoplanetary discs remain poorly constrained by observations. ALMA will soon revolutionise this field. Therefore, it is important to provide predictions which will be valuable for the interpretation of future high sensitivity and high angular resolution observations. Here we present simulated ALMA observations based on radiative transfer modelling of a relatively massive (0.39 M_solar) self-gravitating disc embedded in a 10 M_solar dense core, with structure similar to the pre-stellar core L1544. We focus on simple species and conclude that C17O 3-2, HCO+ 3-2, OCS 26-25 and H2CO 404-303 lines can be used to probe the disc structure and kinematics at all scales.

preprint2013arXiv

The Responses of Magnetically Sub-Critical Cores to Shocks

An ideal magnetohydrodynamics (MHD) code with adaptive mesh refinement (AMR) was used to investigate the interactions of fast-mode shocks with self-gravitating, isothermal cores with mass-to-flux ratios that are somewhat below the minimum value required for gravitational collapse. We find that shock focussing produces colliding flows along the field lines that generate very high densities, even for relatively weak shocks. Self-gravity plays only a minor role in determining the highest density that is reached, but it does play a role in the subsequent evolution. The densities at comparable times differ by a factor of a few for shocks initially propagating perpendicularly or obliquely to the magnetic field in the ambient medium.

preprint2012arXiv

Magnetic Fields and Star Formation

Research performed in the 1950s and 1960s by Leon Mestel on the roles of magnetic fields in star formation established the framework within which he and other key figures have conducted subsequent investigations on the subject. This short tribute to Leon contains a brief summary of some, but not all, of his ground breaking contributions in the area. It also mentions of some of the relevant problems that have received attention in the last few years. The coverage is not comprehensive, and the authors have drawn on their own results more and touched more briefly on those of others than they would in a normal review. Theirs is a personal contribution to the issue honouring Leon, one of the truly great gentlemen, wits, and most insightful of astrophysicists.

preprint2012arXiv

Numerical simulations of shocks encountering clumpy regions

We present numerical simulations of the adiabatic interaction of a shock with a clumpy region containing many individual clouds. Our work incorporates a sub-grid turbulence model which for the first time makes this investigation feasible. We vary the Mach number of the shock, the density contrast of the clouds, and the ratio of total cloud mass to inter-cloud mass within the clumpy region. Cloud material becomes incorporated into the flow. This "mass-loading" reduces the Mach number of the shock, and leads to the formation of a dense shell. In cases in which the mass-loading is sufficient the flow slows enough that the shock degenerates into a wave. The interaction evolves through up to four stages: initially the shock decelerates; then its speed is nearly constant; next the shock accelerates as it leaves the clumpy region; finally it moves at a constant speed close to its initial speed. Turbulence is generated in the post-shock flow as the shock sweeps through the clumpy region. Clouds exposed to turbulence can be destroyed more rapidly than a similar cloud in an "isolated" environment. The lifetime of a downstream cloud decreases with increasing cloud-to-intercloud mass ratio. We briefly discuss the significance of these results for starburst superwinds and galaxy evolution.

preprint2012arXiv

Sputtering in oblique C-type shocks

We present the first results for the sputtering of grain mantles and cores obtained with self-consistent multifluid hydromagnetic models of C-type shocks propagating through dusty media. The threshold shock speed for mantle sputtering is about 10 km/s and is independent of density. The mantles are completely vapourised in shocks with speeds of 20-25 km/s. At such shock speeds core sputtering commences and gas-phase SiO forms. Core destruction is not total in any C-type shock because grains are not completely destroyed in shocks with speeds near the minimum speeds at which J-type shocks appear. Due to the density-dependence of the critical shock speed for this transition, higher gas-phase SiO fractional abundances are produced behind shocks propagating in lower density gas. For shock speeds near the threshold speeds for both core and mantle sputtering, sputtering is much greater for shock velocities at smaller angles relative to the upstream magnetic field. At higher shock speeds, the angular variation is still present but less pronounced.

preprint2011arXiv

Chemistry in a gravitationally unstable protoplanetary disc

Until now, axisymmetric, alpha-disc models have been adopted for calculations of the chemical composition of protoplanetary discs. While this approach is reasonable for many discs, it is not appropriate when self-gravity is important. In this case, spiral waves and shocks cause temperature and density variations that affect the chemistry. We have adopted a dynamical model of a solar-mass star surrounded by a massive (0.39 Msun), self-gravitating disc, similar to those that may be found around Class 0 and early Class I protostars, in a study of disc chemistry. We find that for each of a number of species, e.g. H2O, adsorption and desorption dominate the changes in the gas-phase fractional abundance; because the desorption rates are very sensitive to temperature, maps of the emissions from such species should reveal the locations of shocks of varying strengths. The gas-phase fractional abundances of some other species, e.g. CS, are also affected by gas-phase reactions, particularly in warm shocked regions. We conclude that the dynamics of massive discs have a strong impact on how they appear when imaged in the emission lines of various molecular species.

preprint2011arXiv

Chemistry in Cosmic-Ray Dominated Regions (CRDRs)

Molecular line observations may serve as diagnostics of the degree to which the number density of cosmic ray protons, having energies of 10s to 100s of MeVs each, is enhanced in starburst galaxies and galaxies with active nuclei. Results, obtained with the UCL\_PDR code, for the fractional abundances of molecules as functions of the cosmic-ray induced ionisation rate, $ζ$, are presented. The aim is not to model any particular external galaxies. Rather, it is to identify characteristics of the dependencies of molecular abundances on $ζ$, in part to enable the development of suitable observational programmes for cosmic ray dominated regions (CRDRs) which will then stimulate detailed modelling. For a number density of hydrogen nuclei of of $10^4$ cm$^{-3}$, and high visual extinction, the fractional abundances of some species increase as $ζ$ increases to $10^{-14}$ s$^{-1}$, but for much higher values of $ζ$ the fractional abundances of all molecular species are significantly below their peak values. We show in particular that OH, H$_{2}$O, H$_{3}^{+}$, H$_{3}$O$^{+}$ and OH$^{+}$ attain large fractional abundances ($\geqslant 10^{-8}$) for $ζ$ as large as $10^{-12}$ s$^{-1}$. HCO$^{+}$ is a poor tracer of CRDRs when $ζ> 10^{-13}$ s$^{-1}$. Sulphur-bearing species may be useful tracers of CRDRs gas in which $ζ\sim 10^{-16}$ s$^{-1}$. Ammonia has a large fractional abundance for $ζ\leqslant 10^{-16}$ s$^{-1}$ and nitrogen appears in CN-bearing species at significant levels as $ζ$ increases, even up to $\sim 10^{-14}$ s$^{-1}$. In this paper, we also discuss our model predictions, comparing them to recent detections in both galactic and extragalactic sources. We show that they agree well, to a first approximation, with the observational constraints.

preprint2011arXiv

Cosmic Rays and Radiative Instabilities

In the absence of magnetic fields and cosmic rays, radiative cooling laws with a range of dependences on temperature affect the stability of interstellar gas. For about four and a half decades, astrophysicists have recognised the importance of the thermal instablity for the formation of clouds in the interstellar medium. Even in the past several years, many papers have concerned the role of the thermal instability in the production of molecular clouds. About three and a half decades ago, astrophysicists investigating radiative shocks noticed that for many cooling laws such shocks are unstable. Attempts to address the effects of cosmic rays on the stablity of radiative media that are initially uniform or that have just passed through shocks have been made. The simplest approach to such studies involves the assumption that the cosmic rays behave as a fluid. Work based on such an approach is described. Cosmic rays have no effect on the stability of initially uniform, static media with respect to isobaric perturbations, though they do affect the stability of such media with respect to isentropic perturbations. The effect of cosmic rays on the stability of radiative shocked media depends greatly on the efficiency of the conversion of energy in accelerated cosmic rays into thermal energy in the thermalized fluid. If that efficiency is low, radiative cooling makes weak shocks propagating into upstream media with low cosmic-ray pressures more likely to be cosmic-ray dominated than adiabatic shocks of comparable strength. The cosmic-ray dominated shocks do not display radiative overstability. Highly efficient conversion of cosmic-ray energy into thermal energy leads shocked media to behave as they do when cosmic rays are absent.

preprint2011arXiv

First observations of CN(2-1), HCO$^{+}$(3-2) and C$_{2}$H(3-2) emission lines in the Perseus cluster: constraints on heating mechanisms in the cluster gas

We present the first observations of emission lines of CN(2-1), HCO$^{+}$(3-2) and C$_{2}$H(3-2) in the Perseus cluster. We observed at two positions: directly at the central galaxy, NGC 1275 and also at a position about 20$''$ to the east where associated filamentary structure has been shown to have strong CO emission. Clear detections in CN and HCO$^{+}$ transitions and a weak detection of the C$_{2}$H transition were made towards NGC 1275, while weak detections of CN and HCO$^{+}$ were made towards the eastern filamentary structure. Crude estimates of the column densities and fractional abundances (mostly upper limits) as functions of an unknown rotational temperature were made to both sources. These observational data were compared with the outputs of thermal/chemical models previously published by \citet{Baye10c} in an attempt to constrain the heating mechanisms in cluster gas. We find that models in which heating is dominated by cosmic rays can account for the molecular observations. This conclusion is consistent with that of \citet{Ferl09} in their study of gas traced by optical and infrared radiation. The cosmic ray heating rate in the regions probed by molecular emissions is required to be at least two orders of magnitude larger than that in the Milky Way.

preprint2010arXiv

Dusty Magnetohydrodynamics in Star Forming Regions

Star formation occurs in dark molecular regions where the number density of hydrogen nuclei, nH, exceeds 10^4 cc and the fractional ionization is 10^-7 or less. Dust grains with sizes ranging up to tenths of microns and perhaps down to tens of nanometers contain just under one percent of the mass. Recombination on grains is important for the removal of gas phase ions, which are produced by cosmic rays penetrating the dark regions. Collisions of neutrals with charged grains contribute significantly to the coupling of the magnetic field to the neutral gas. Consequently, the dynamics of the grains must be included in the magnetohydrodynamic models of large scale collapse, the evolution of waves and the structures of shocks important in star formation.

preprint2010arXiv

Shock-triggered formation of magnetically-dominated clouds. II. Weak shock-cloud interaction in three dimensions

To understand the formation of a magnetically dominated molecular cloud from an atomic cloud, we study the interaction of a weak, radiative shock with a magnetised cloud. The thermally stable warm atomic cloud is initially in static equilibrium with the surrounding hot ionised gas. A shock propagating through the hot medium then interacts with the cloud. We follow the dynamical evolution of the shocked cloud with a time-dependent ideal magnetohydrodynamic code. By performing the simulations in 3D, we investigate the effect of different magnetic field orientations including parallel, perpendicular and oblique to the shock normal. We find that the angle between the shock normal and the magnetic field must be small to produce clouds with properties similar to observed molecular clouds.

preprint2010arXiv

The Turbulent Destruction of Clouds - II. Mach Number Dependence, Mass-loss Rates, and Tail Formation

The turbulent destruction of a cloud subject to the passage of an adiabatic shock is studied. We find large discrepancies between the lifetime of the cloud and the analytical result of Hartquist et al. (1986). These differences appear to be due to the assumption in Hartquist et al. that mass-loss occurs largely as a result of lower pressure regions on the surface of the cloud away from the stagnation point, whereas in reality Kelvin-Helmholtz (KH) instabilities play a dominant role in the cloud destruction. We find that the true lifetime of the cloud (defined as when all of the material from the core of the cloud is well mixed with the intercloud material in the hydrodynamic cells) is about 6 times t_KHD, where t_KHD is the growth timescale for the most disruptive, long-wavelength, KH instabilities. These findings have wide implications for diffuse sources where there is transfer of material between hot and cool phases. The properties of the interaction as a function of Mach number and cloud density contrast are also studied. The interaction is milder at lower Mach numbers with the most marked differences occuring at low shock Mach numbers when the postshock gas is subsonic with respect to the cloud (i.e. M < 2.76). Material stripped off the cloud only forms a long "tail-like" feature if the density contrast of the cloud to the ambient medium, chi > 1e3.

preprint2009arXiv

Transient evolution of C-type shocks in dusty regions of varying density

Outflows of young stars drive shocks into dusty, molecular regions. Most models of such shocks assume that they are steady and propagating perpendicular to the magnetic field. Real shocks often violate both of these assumptions and the media through which they propagate are inhomogeneous. We use the code employed previously to produce the first time-dependent simulations of fast-mode, oblique C-type shocks interacting with density perturbations. We include a self-consistent calculation of the thermal and ionisation balances and a fluid treatment of grains. We identify features that develop when a multifluid shock encounters a density inhomogeneity to investigate whether any part of the precursor region ever behaves in a quasi-steady fashion. If it does the shock may be modelled approximately without solving the time-dependent hydromagnetic equations. Simulations were made for initially steady oblique C-type shocks encountering density inhomogeneities. For a semi-finite inhomogeneity with a density larger than the surrounding medium, a transmitted shock evolves from being J-type to a steady C-type shock on a timescale comparable to the ion-flow time through it. A sufficiently upstream part of the precursor of an evolving J-type shock is quasi-steady. The ion-flow timescale is also relevant for the evolution of a shock moving into a region of decreasing density. The models for shocks propagating into regions in which the density increases and then decreases to its initial value cannot be entirely described in terms of the results obtained for monotonically increasing and decreasing densities. For the latter model, the long-term evolution to a C-type shock cannot be approximated by quasi-steady models.

preprint2008arXiv

A cosmic-ray precursor model for a Balmer-dominated shock in Tycho's supernova remnant

We present a time-dependent cosmic-ray modified shock model for which the calculated H-alpha emissivity profile agrees well with the H-alpha flux increase ahead of the Balmer-dominated shock at knot g in Tycho's supernova remnant, observed by Lee et al (2007). The backreaction of the cosmic ray component on the thermal component is treated in the two-fluid approximation, and we include thermal particle injection and energy transfer due to the acoustic instability in the precursor. The transient state of our model that describes the current state of the shock at knot g, occurs during the evolution from a thermal gas dominated shock to a smooth cosmic-ray dominated shock. Assuming a distance of 2.3 kpc to Tycho's remnant we obtain values for the cosmic ray diffusion coefficient, the injection parameter, and the time scale for the energy transfer of 10^{24} cm^{2} s^{-1}, 4.2x10^{-3}, and 426 y, respectively. We have also studied the parameter space for fast (300 km s^{-1} - 3000 km s^{-1}), time-asymptotically steady shocks and have identified a branch of solutions, for which the temperature in the cosmic ray precursor typically reaches 2-6x10^{4} K and the bulk acceleration of the flow through the precursor is less than 10 km s^{-1}. These solutions fall into the low cosmic ray acceleration efficiency regime and are relatively insensitive to shock parameters. This low cosmic ray acceleration efficiency branch of solutions may provide a natural explanation for the line broadening of the H-alpha narrow component observed in non-radiative shocks in many supernova remnants.

preprint2008arXiv

Dust in Interstellar Clouds, Evolved Stars and Supernovae

Outflows of pre-main-sequence stars drive shocks into molecular material within 0.01 - 1 pc of the young stars. The shock-heated gas emits infrared, millimeter and submillimeter lines of many species including. Dust grains are important charge carriers and play a large role in coupling the magnetic field and flow of neutral gas. Some effects of the dust on the dynamics of oblique shocks began to emerge in the 1990s. However, detailed models of these shocks are required for the calculation of the grain sputtering contribution to gas phase abundances of species producing observed emissions. We are developing such models. Some of the molecular species introduced into the gas phase by sputtering in shocks or by thermally driven desorption in hot cores form on grain surfaces. Recently laboratory studies have begun to contribute to the understanding of surface reactions and thermally driven desorption important for the chemistry of star forming clouds. Dusty plasmas are prevalent in many evolved stars just as well as in star forming regions. Radiation pressure on dust plays a significant role in mass loss from some post-main-sequence stars. The mechanisms leading to the formation of carbonaceous dust in the stellar outflows are similar to those important for soot formation in flames. However, nucleation in oxygen-rich outflows is less well understood and remains a challenging research area. Dust is observed in supernova ejecta that have not passed through the reverse shocks that develop in the interaction of ejecta with ambient media. Dust is detected in high redshift galaxies that are sufficiently young that the only stars that could have produced the dust were so massive that they became supernovae. Consequently, the issue of the survival of dust in strong supernova shocks is of considerable interest.

preprint2003arXiv

The Evolution of Mass Loaded Supernova Remnants. II. Temperature Dependent Mass Injection Rates

We investigate the evolution of spherically symmetric supernova remnants in which mass loading takes place due to conductively driven evaporation of embedded clouds. Numerical simulations reveal significant differences between the evolution of conductively mass loaded and the ablatively mass loaded remnants studied in Paper I. A main difference is the way in which conductive mass loading is extinguished at fairly early times, once the interior temperature of the remnant falls below ~10 million K. Thus, at late times remnants that ablatively mass load are dominated by loaded mass and thermal energy, while those that conductively mass load are dominated by swept-up mass and kinetic energy. Simple approximations to the remnant evolution, complementary to those in Paper I, are given.