Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “finding”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,639 records · Page 91
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Star cluster formation in cosmological simulations. I. Properties of young clusters

We present a new implementation of star formation in cosmological simulations by considering star clusters as a unit of star formation. Cluster particles grow in mass over several million years at the rate determined by local gas properties, with high time resolution. The particle growth is terminated by its own energy and momentum feedback on the interstellar medium. We test this implementation for Milky Way-sized galaxies at high redshift by comparing the properties of model clusters with observations of young star clusters. We find that the cluster initial mass function is best described by a Schechter function rather than a single power law. In agreement with observations, at low masses the logarithmic slope is $\alpha \approx 1.8\mbox{–}2$, while the cutoff at high mass scales with the star formation rate (SFR). A related trend is a positive correlation between the surface density of the SFR and fraction of stars contained in massive clusters. Both trends indicate that the formation of massive star clusters is preferred during bursts of star formation. These bursts are often associated with major-merger events. We also find that the median timescale for cluster formation ranges from 0.5 to 4 Myr and decreases systematically with increasing star formation efficiency. Local variations in the gas density and cluster accretion rate naturally lead to the scatter of the overall formation efficiency by an order of magnitude, even when the instantaneous efficiency is kept constant. As a result, comparison of the formation timescale with the observed age spread of young star clusters provides an additional important constraint on the modeling of star formation and feedback schemes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

ORIGINS OF SCATTER IN THE RELATIONSHIP BETWEEN HCN 1-0 AND DENSE GAS MASS IN THE GALACTIC CENTER

We investigate the correlation of HCN 1-0 with gas mass in the central 300 pc of the Galaxy. We find that on the ∼10 pc size scale of individual cloud cores, HCN 1-0 is well correlated with dense gas mass when plotted as a log–log relationship. There is ∼0.75 dex of scatter in this relationship from clouds like Sgr B2, which has an integrated HCN 1-0 intensity of a cloud less than half its mass, and others that have HCN 1-0 enhanced by a factor of 2–3 relative to clouds of comparable mass. We identify the two primary sources of scatter to be self-absorption and variations in HCN abundance. We also find that the extended HCN 1-0 emission is more intense per unit mass than in individual cloud cores. In fact the majority (80%) of HCN 1-0 emission comes from extended gas with column densities below 7 × 10{sup 22} cm{sup −2}, accounting for 68% of the total mass. We find variations in the brightness of HCN 1-0 would only yield a ∼10% error in the dense gas mass inferred from this line in the Galactic center. However, the observed order of magnitude HCN abundance variations, and the systematic nature of these variations, warn of potential biases in the use of HCN as dense gas mass tracer in more extreme environments such as an active galactic nucleus and shock-dominated regions. We also investigate other 3 mm tracers, finding that HNCO is better correlated with mass than HCN, and might be a better tracer of cloud mass in this environment.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Long-lived Dust Asymmetries at Dead Zone Edges in Protoplanetary Disks

A number of transition disks exhibit significant azimuthal asymmetries in thermal dust emission. One possible origin for these asymmetries is dust trapping in vortices formed at the edges of dead zones. We carry out high-resolution, two-dimensional hydrodynamic simulations of this scenario, including the effects of dust feedback. We find that, although feedback weakens the vortices and slows down the process of dust accumulation, the dust distribution in the disk can nonetheless remain asymmetric for many thousands of orbits. We show that even after 10{sup 4} orbits, or 2.5 Myr when scaled to the parameters of Oph IRS 48 (a significant fraction of its age), the dust is not dispersed into an axisymmetric ring, in contrast to the case of a vortex formed by a planet. This is because accumulation of mass at the dead zone edge constantly replenishes the vortex, preventing it from being fully destroyed. We produce synthetic dust emission images using our simulation results. We find that multiple small clumps of dust may be distributed azimuthally. These clumps, if not resolved from one another, appear as a single large feature. A defining characteristic of a disk with a dead zone edge is that an asymmetric feature is accompanied by a ring of dust located about twice as far from the central star.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Polarization Signatures of Kink Instabilities in the Blazar Emission Region from Relativistic Magnetohydrodynamic Simulations

Kink instabilities are likely to occur in the current-carrying magnetized plasma jets. Recent observations of the blazar radiation and polarization signatures suggest that the blazar emission region may be considerably magnetized. While the kink instability has been studied with first-principle magnetohydrodynamic (MHD) simulations, the corresponding time-dependent radiation and polarization signatures have not been investigated. In this paper, we perform comprehensive polarization-dependent radiation modeling of the kink instability in the blazar emission region based on relativistic MHD (RMHD) simulations. We find that the kink instability may give rise to strong flares with polarization angle (PA) swings or weak flares with polarization fluctuations, depending on the initial magnetic topology and magnetization. These findings are consistent with observations. Compared with the shock model, the kink model generates polarization signatures that are in better agreement with the general polarization observations. Therefore, we suggest that kink instabilities may widely exist in the jet environment and provide an efficient way to convert the magnetic energy and produce multiwavelength flares and polarization variations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Thermal pressure in the cold neutral medium of nearby galaxies

Dynamic and thermal processes regulate the structure of the multiphase interstellar medium (ISM), and ultimately establish how galaxies evolve through star formation. Thus, to constrain ISM models and better understand the interplay of these processes, it is of great interest to measure the thermal pressure (P{sub th}) of the diffuse, neutral gas. By combining [C ii] 158 μm, H i, and CO data from 31 galaxies selected from the Herschel KINGFISH sample, we have measured thermal pressures in 534 predominantly atomic regions with typical sizes of ∼1 kiloparsec. We find a distribution of thermal pressures in the P{sub th}/k∼10{sup 3}--10{sup 5} K cm{sup −3} range. For a sub-sample of regions with conditions similar to those of the diffuse, neutral gas in the Galactic plane, we find thermal pressures that follow a log-normal distribution with a median value of P{sub th}/k ≈ 3600 K cm{sup −3}. These results are consistent with thermal pressure measurements using other observational methods. We find that P{sub th} increases with radiation field strength and star formation activity, as expected from the close link between the heating of the gas and the star formation rate. Our thermal pressure measurements fall in the regime where a two-phase ISM with cold and warm neutral media could exist in pressure equilibrium. Finally, we find that the midplane thermal pressure of the diffuse gas is about ∼30% of the vertical weight of the overlying ISM, consistent with results from hydrodynamical simulations of self-regulated star formation in galactic disks.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Evolution Models of Helium White Dwarf–Main-sequence Star Merger Remnants

It is predicted that orbital decay by gravitational-wave radiation and tidal interaction will cause some close binary stars to merge within a Hubble time. The merger of a helium-core white dwarf with a main-sequence (MS) star can produce a red giant branch star that has a low-mass hydrogen envelope when helium is ignited and thus become a hot subdwarf. Because detailed calculations have not been made, we compute post-merger models with a stellar evolution code. We find the evolutionary paths available to merger remnants and find the pre-merger conditions that lead to the formation of hot subdwarfs. We find that some such mergers result in the formation of stars with intermediate helium-rich surfaces. These stars later develop helium-poor surfaces owing to diffusion. Combining our results with a model population and comparing to observed stars, we find that some observed intermediate helium-rich hot subdwarfs can be explained as the remnants of the mergers of helium-core white dwarfs with low-mass MS stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Molecular Gas Kinematics and Star Formation Properties of the Strongly-lensed Quasar Host Galaxy RXS J1131–1231

We report observations of CO(J = 2 → 1) and CO(J=3→2) line emission toward the quadruply-lensed quasar RXS J1131−1231 at z = 0.654 obtained using the Plateau de Bure Interferometer (PdBI) and the Combined Array for Research in Millimeter-wave Astronomy (CARMA). Our lens modeling shows that the asymmetry in the double-horned CO(J = 2 → 1) line profile is mainly a result of differential lensing, where the magnification factor varies from ∼3 to ∼9 across different kinematic components. The intrinsically symmetric line profile and a smooth source-plane velocity gradient suggest that the host galaxy is an extended rotating disk, with a CO size of R{sub CO}∼6 kpc and a dynamical mass of M{sub dyn}∼8×10{sup 10} M {sub ⊙}. We also find a secondary CO-emitting source near RXS J1131−1231, the location of which is consistent with the optically-faint companion reported in previous studies. The lensing-corrected molecular gas masses are M {sub gas} = (1.4 ± 0.3) × 10{sup 10} M {sub ⊙} and (2.0 ± 0.1) × 10{sup 9} M {sub ⊙} for RXS J1131−1231 and the companion, respectively. We find a lensing-corrected stellar mass of M {sub *} = (3 ± 1) × 10{sup 10} M {sub ⊙} and a star formation rate of SFR{sub FIR} = (120 ± 63) M {sub ⊙} yr{sup −1}, corresponding to a specific SFR and star formation efficiency comparable to z ∼ 1 disk galaxies not hosting quasars. The implied gas mass fraction of ∼18 ± 4% is consistent with the previously observed cosmic decline since z ∼ 2. We thus find no evidence for quenching of star formation in RXS J1131−1231. This agrees with our finding of an elevated M{sub BH}/M{sub bulge} ratio of >0.27{sub −0.08}{sup +0.11}% compared to the local value, suggesting that the bulk of its black hole mass is largely in place while its stellar bulge is still assembling.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Velocity segregation and systematic biases in velocity dispersion estimates with the SPT-GMOS spectroscopic survey

The velocity distribution of galaxies in clusters is not universal; rather, galaxies are segregated according to their spectral type and relative luminosity. We examine the velocity distributions of different populations of galaxies within 89 Sunyaev Zel'dovich (SZ) selected galaxy clusters spanning $ 0.28 < z < 1.08$. Our sample is primarily draw from the SPT-GMOS spectroscopic survey, supplemented by additional published spectroscopy, resulting in a final spectroscopic sample of 4148 galaxy spectra---2868 cluster members. The velocity dispersion of star-forming cluster galaxies is $17\pm4$% greater than that of passive cluster galaxies, and the velocity dispersion of bright ($m < m^{*}-0.5$) cluster galaxies is $11\pm4$% lower than the velocity dispersion of our total member population. We find good agreement with simulations regarding the shape of the relationship between the measured velocity dispersion and the fraction of passive vs. star-forming galaxies used to measure it, but we find a small offset between this relationship as measured in data and simulations in which suggests that our dispersions are systematically low by as much as 3\% relative to simulations. We argue that this offset could be interpreted as a measurement of the effective velocity bias that describes the ratio of our observed velocity dispersions and the intrinsic velocity dispersion of dark matter particles in a published simulation result. Here, by measuring velocity bias in this way suggests that large spectroscopic surveys can improve dispersion-based mass-observable scaling relations for cosmology even in the face of velocity biases, by quantifying and ultimately calibrating them out.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Mass, Color, and Structural Evolution of Today’s Massive Galaxies Since z ∼ 5

In this paper, we use stacking analysis to trace the mass growth, color evolution, and structural evolution of present-day massive galaxies (log(M{sub ∗}/M{sub ⊙})=11.5) out to z = 5. We utilize the exceptional depth and area of the latest UltraVISTA data release, combined with the depth and unparalleled seeing of CANDELS to gather a large, mass-selected sample of galaxies in the NIR (rest-frame optical to UV). Progenitors of present-day massive galaxies are identified via an evolving cumulative number density selection, which accounts for the effects of merging to correct for the systematic biases introduced using a fixed cumulative number density selection, and find progenitors grow in stellar mass by ≈1.5 dex since z = 5. Using stacking, we analyze the structural parameters of the progenitors and find that most of the stellar mass content in the central regions was in place by z∼2, and while galaxies continue to assemble mass at all radii, the outskirts experience the largest fractional increase in stellar mass. However, we find evidence of significant stellar mass build-up at r<3 kpc beyond z>4 probing an era of significant mass assembly in the interiors of present-day massive galaxies. We also compare mass assembly from progenitors in this study to the EAGLE simulation and find qualitatively similar assembly with z at r<3 kpc. We identify z∼1.5 as a distinct epoch in the evolution of massive galaxies where progenitors transitioned from growing in mass and size primarily through in situ star formation in disks to a period of efficient growth in r {sub e} consistent with the minor merger scenario.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Light curve and SED modeling of the gamma-ray binary 1FGL J1018.6–5856: Constraints on the orbital geometry and relativistic flow

We present broadband spectral energy distributions and light curves of the gamma-ray binary 1FGL J1018.6–5856 measured in the X-ray and the gamma-ray bands. We find that the orbital modulation in the low-energy gamma-ray band is similar to that in the X-ray band, suggesting a common spectral component. However, above a GeV the orbital light curve changes significantly. We suggest that the GeV band contains significant flux from a pulsar magnetosphere, while the X-ray to TeV light curves are dominated by synchrotron and Compton emission from an intrabinary shock (IBS). We find that a simple one-zone model is inadequate to explain the IBS emission, but that beamed Synchrotron-self Compton radiation from adiabatically accelerated plasma in the shocked pulsar wind can reproduce the complex multiband light curves, including the variable X-ray spike coincident with the gamma-ray maximum. Furthermore, the model requires an inclination of ~50° and an orbital eccentricity of ~0.35, consistent with the limited constraints from existing optical observations. This picture motivates searches for pulsations from the energetic young pulsar powering the wind shock.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

First Discoveries of z > 6 Quasars with the DECam Legacy Survey and UKIRT Hemisphere Survey

In this paper, we present the first discoveries from a survey of z ≳ 6 quasars using imaging data from the DECam Legacy Survey (DECaLS) in the optical, the UKIRT Deep Infrared Sky Survey (UKIDSS) and a preliminary version of the UKIRT Hemisphere Survey (UHS) in the near-IR, and ALLWISE in the mid-IR. DECaLS will image 9000 deg 2 of sky down to z AB ~ 23.0, and UKIDSS and UHS will map the northern sky at 0 < decl. < +60°, reaching J VEGA ~ 19.6 (5-σ). The combination of these data sets allows us to discover quasars at redshift z ≳ 7 and to conduct a complete census of the faint quasar population at z ≳ 6. In this paper, we report on the selection method of our search, and on the initial discoveries of two new, faint z ≳ 6 quasars and one new z = 6.63 quasar in our pilot spectroscopic observations. The two new z ~ 6 quasars are at z = 6.07 and z = 6.17 with absolute magnitudes at rest-frame wavelength 1450 Å being M 1450 = -25.83 and M 1450 = -25.76, respectively. These discoveries suggest that we can find quasars close to or fainter than the break magnitude of the Quasar Luminosity Function (QLF) at z ≳ 6. The new z = 6.63 quasar has an absolute magnitude of M 1450 = -25.95. This demonstrates the potential of using the combined DECaLS and UKIDSS/UHS data sets to find z ≳ 7 quasars. Finally, extrapolating from previous QLF measurements, we predict that these combined data sets will yield ~200 z ~ 6 quasars to z AB < 21.5, ~1000 z ~ 6 quasars to z AB < 23, and ~30 quasars at z > 6.5 to J VEGA < 19.5.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Extended Gamma-Ray Emission from the G25.0+0.0 Region: A Star-forming Region Powered by the Newly Found OB Association?

We report a study of extended γ -ray emission with the Large Area Telescope (LAT) on board the Fermi Gamma-ray Space Telescope , which is likely to be the second case of a γ -ray detection from a star-forming region (SFR) in our Galaxy. The LAT source is located in the G25 region, 1.°7 × 2.°1 around ( l , b ) = (25.°0, 0.°0). The γ -ray emission is found to be composed of two extended sources and one pointlike source. The extended sources have similar sizes of about 1.°4 × 0.°6. An ∼0.°4 diameter subregion of one has a photon index of Γ = 1.53 ± 0.15, and is spatially coincident with HESS J1837−069, likely a pulsar wind nebula. The other parts of the extended sources have a photon index of Γ = 2.1 ± 0.2 without significant spectral curvature. Given their spatial and spectral properties, they have no clear associations with sources at other wavelengths. Their γ -ray properties are similar to those of the Cygnus cocoon SFR, the only firmly established γ -ray detection of an SFR in the Galaxy. Indeed, we find bubble-like structures of atomic and molecular gas in G25, which may be created by a putative OB association/cluster. The γ -ray emitting regions appear confined in the bubble-like structure; similar properties are also found in the Cygnus cocoon. In addition, using observations with the XMM-Newton , we find a candidate young massive OB association/cluster G25.18+0.26 in the G25 region. We propose that the extended γ -ray emission in G25 is associated with an SFR driven by G25.18+0.26. Based on this scenario, we discuss possible acceleration processes in the SFR and compare them with the Cygnus cocoon.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Implications of z ~ 6 Quasar Proximity Zones for the Epoch of Reionization and Quasar Lifetimes

In this paper, we study quasar proximity zones in the redshift range $5.77\leqslant z\leqslant 6.54$ by homogeneously analyzing 34 medium-resolution spectra, encompassing both archival and newly obtained data, and exploiting recently updated systemic redshift and magnitude measurements. Whereas previous studies found strong evolution of proximity zone sizes with redshift and argued that this provides evidence for a rapidly evolving intergalactic medium (IGM) neutral fraction during reionization, we measure a much shallower trend $\propto {(1+z)}^{-1.44}$. We compare our measured proximity zone sizes to predictions from hydrodynamical simulations post-processed with one-dimensional radiative transfer and find good agreement between observations and theory, irrespective of the ionization state of the ambient IGM. This insensitivity to IGM ionization state has been previously noted, and results from the fact that the definition of proximity zone size as the first drop of the smoothed quasar spectrum below the 10% flux transmission level probes locations where the ionizing radiation from the quasar is an order of magnitude larger than the expected ultraviolet ionizing background that sets the neutral fraction of the IGM. Our analysis also uncovered three objects with exceptionally small proximity zones (two have ${R}_{p}\lt 1$ proper Mpc), which constitute outliers from the observed distribution and are challenging to explain with our radiative transfer simulations. Finally, we consider various explanations for their origin, such as strong absorption line systems associated with the quasar or patchy reionization, but find that the most compelling scenario is that these quasars have been shining for ≲10 5 years.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Is There a Disk of Satellites around the Milky Way?

The “disk of satellites” (DoS) around the Milky Way is a highly debated topic with conflicting interpretations of observations and their theoretical models. We perform a comprehensive analysis of all of the dwarfs detected in the Milky Way and find that the DoS structure depends strongly on the plane identification method and the sample size. In particular, we demonstrate that a small sample size can artificially produce a highly anisotropic spatial distribution and a strong clustering of the angular momentum of the satellites. Moreover, we calculate the evolution of the 11 classical satellites with proper motion measurements and find that the thin DoS in which they currently reside is transient. Furthermore, we analyze two cosmological simulations using the same initial conditions of a Milky-Way-sized galaxy, an N -body run with dark matter only, and a hydrodynamic one with both baryonic and dark matter, and find that the hydrodynamic simulation produces a more anisotropic distribution of satellites than the N -body one. Our results suggest that an anisotropic distribution of satellites in galaxies can originate from baryonic processes in the hierarchical structure formation model, but the claimed highly flattened, coherently rotating DoS of the Milky Way may be biased by the small-number selection effect. These findings may help resolve the contradictory claims of DoS in galaxies and the discrepancy among numerical simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

The Viewing Geometry of Brown Dwarfs Influences Their Observed Colors and Variability Amplitudes

In this paper we study the full sample of known Spitzer [3.6 μm] and J-band variable brown dwarfs. We calculate the rotational velocities, vsin i, of 16 variable brown dwarfs using archival Keck NIRSPEC data and compute the inclination angles of 19 variable brown dwarfs. The results obtained show that all objects in the sample with mid-IR variability detections are inclined at an angle >20{sup ∘}, while all objects in the sample displaying J-band variability have an inclination angle >35{sup ∘}. J-band variability appears to be more affected by inclination than Spitzer [3.6 μm] variability, and is strongly attenuated at lower inclinations. Since J-band observations probe deeper into the atmosphere than mid-IR observations, this effect may be due to the increased atmospheric path length of J-band flux at lower inclinations. We find a statistically significant correlation between the color anomaly and inclination of our sample, where field objects viewed equator-on appear redder than objects viewed at lower inclinations. Considering the full sample of known variable L, T, and Y spectral type objects in the literature, we find that the variability properties of the two bands display notably different trends that are due to both intrinsic differences between bands and the sensitivity of ground-based versus space-based searches. However, in both bands we find that variability amplitude may reach a maximum at ∼7–9 hr periods. Finally, we find a strong correlation between color anomaly and variability amplitude for both the J-band and mid-IR variability detections, where redder objects display higher variability amplitudes.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

(Almost) Dark Galaxies in the ALFALFA Survey: Isolated H i-bearing Ultra-diffuse Galaxies

Here, we present a sample of 115 very low optical surface brightness, highly extended, H i-rich galaxies carefully selected from the ALFALFA survey that have similar optical absolute magnitudes, surface brightnesses, and radii to recently discovered "ultra-diffuse" galaxies (UDGs). However, these systems are bluer and have more irregular morphologies than other UDGs, are isolated, and contain significant reservoirs of H i. We find that while these sources have normal star formation rates for H i-selected galaxies of similar stellar mass, they have very low star formation efficiencies. We further present deep optical and H i-synthesis follow-up imaging of three of these H i-bearing ultra-diffuse sources. We measure H i diameters extending to ~40 kpc, but note that while all three sources have large H i diameters for their stellar mass, they are consistent with the H i mass–H i radius relation. We further analyze the H i velocity widths and rotation velocities for the unresolved and resolved sources, respectively, and find that the sources appear to inhabit halos of dwarf galaxies. Finally, we estimate spin parameters, and suggest that these sources may exist in high spin parameter halos, and as such may be potential H i-rich progenitors to the ultra-diffuse galaxies observed in cluster environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

PTF11kx: A Type Ia Supernova with Hydrogen Emission Persisting after 3.5 Years

The optical transient PTF11kx exhibited both the characteristic spectral features of Type Ia supernovae (SNe Ia) and the signature of ejecta interacting with circumstellar material (CSM) containing hydrogen, indicating the presence of a nondegenerate companion. In this paper, we present an optical spectrum at 1342 days after peak from Keck Observatory, in which the broad component of Hα emission persists with a similar profile as in early-time observations. We also present Spitzer IRAC detections obtained 1237 and 1818 days after peak, and an upper limit from Hubble Space Telescope ultraviolet imaging at 2133 days. We interpret our late-time observations in the context of published results—and reinterpret the early-time observations—in order to constrain the CSM's physical parameters and to compare to theoretical predictions for recurrent-nova systems. We find that the CSM's radial extent may be several times the distance between the star and the CSM's inner edge, and that the CSM column density may be two orders of magnitude lower than previous estimates. We show that the Hα luminosity decline is similar to other SNe with CSM interaction and demonstrate how our infrared photometry is evidence for newly formed, collisionally heated dust. We create a model for PTF11kx's late-time CSM interaction and find that X-ray reprocessing by photoionization and recombination cannot reproduce the observed Hα luminosity, suggesting that the X-rays are thermalized and that Hα radiates from collisional excitation. In conclusion, we discuss the implications of our results regarding the progenitor scenario and the geometric properties of the CSM for the PTF11kx system.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Ab Initio Simulations of a Supernova-driven Galactic Dynamo in an Isolated Disk Galaxy

We study the magnetic field evolution of an isolated spiral galaxy, using isolated Milky Way–mass galaxy formation simulations and a novel prescription for magnetohydrodynamic (MHD) supernova feedback. Our main result is that a galactic dynamo can be seeded and driven by supernova explosions, resulting in magnetic fields whose strength and morphology are consistent with observations. In our model, supernovae supply thermal energy and a low-level magnetic field along with their ejecta. The thermal expansion drives turbulence, which serves a dual role by efficiently mixing the magnetic field into the interstellar medium and amplifying it by means of a turbulent dynamo. The computational prescription for MHD supernova feedback has been implemented within the publicly available ENZO code and is fully described in this paper. This improves upon ENZO 's existing modules for hydrodynamic feedback from stars and active galaxies. We find that the field attains microgauss levels over gigayear timescales throughout the disk. The field also develops a large-scale structure, which appears to be correlated with the disk’s spiral arm density structure. We find that seeding of the galactic dynamo by supernova ejecta predicts a persistent correlation between gas metallicity and magnetic field strength. We also generate all-sky maps of the Faraday rotation measure from the simulation-predicted magnetic field, and we present a direct comparison with observations.

79 ASTRONOMY AND ASTROPHYSICS↗