Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak localization”

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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

EVIDENCE OF THE DYNAMICS OF RELATIVISTIC JET LAUNCHING IN QUASARS

Hubble Space Telescope (HST) spectra of the EUV, the optically thick emission from the innermost accretion flow onto the central supermassive black hole, indicate that radio loud quasars (RLQs) tend to be EUV weak compared to the radio-quiet quasars; yet the remainder of the optically thick thermal continuum is indistinguishable. The deficit of EUV emission in RLQs has a straightforward interpretation as a missing or a suppressed innermost region of local energy dissipation in the accretion flow. This article is an examination of the evidence for a distribution of magnetic flux tubes in the innermost accretion flow that results in magnetically arrested accretion (MAA) and creates the EUV deficit. These same flux tubes and possibly the interior magnetic flux that they encircle are the sources of the jet power as well. In the MAA scenario, islands of large-scale vertical magnetic flux perforate the innermost accretion flow of RLQs. The first prediction of the theory that is supported by the HST data is that the strength of the (large-scale poloidal magnetic fields) jets in the MAA region is regulated by the ram pressure of the accretion flow in the quasar environment. The second prediction that is supported by the HST data is that the rotating magnetic islands remove energy from the accretion flow as a Poynting flux dominated jet in proportion to the square of the fraction of the EUV emitting gas that is displaced by these islands.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Effect of anisotropic fast ions on internal kink stability in DIII-D negative and positive triangularity plasmas

Recent DIII-D experiments show that sawtooth stability is strongly affected by anisotropic fast ions from neutral beam injection (NBI) in both negative and positive triangularity plasmas. Fast ions from co-current NBI are stabilizing for the sawtooth stability, resulting in longer sawtooth periods. On the other hand, fast ions from counter-current NBI are destabilizing, leading to small and frequent sawteeth. The relative change of sawtooth period and amplitude is more than a factor of two. These observations appear to hold in both plasma shapes. Non-perturbative toroidal modeling, utilizing the magnetohydrodynamic-kinetic hybrid stability code MARS-K (Liu et al 2008 Phys. Plasmas 15 112503), reveals an asymmetric dependence of the stability of the n = 1 (n is the toroidal mode number) internal kink mode on the injection direction of NBI, being qualitatively consistent with the experimentally observed sawtooth behavior. The MARS-K modeling results suggest that anisotropic fast ions affect the mode growth rate and frequency through both adiabatic and non-adiabatic contributions. Here, the asymmetry of the internal kink mode instability relative to the NBI direction is mainly due to the non-adiabatic contribution of passing fast ions, which stabilize (destabilize) the internal kink with the co-(counter-) current NBI as compared to the fluid counterpart. However, finite orbit width (FOW) correction to passing particles partially cancels the asymmetry. Trapped particles are always stabilizing due to precessional drift resonance. Modeling also shows that fast ions affect the internal kink in a similar manner in both negative and positive triangularity plasmas, although being slightly more unstable in the negative triangularity configuration already in the fluid limit. The similarity is mainly attributed to the fact that the mode is localized in the plasma core region, with very similar eigenmode structures in both negative and positive configurations. Furthermore, MARS-K modeling indicates that other factors, such as the plasma rotation and the drift kinetic effects of thermal plasmas, weakly modify the mode stability as compared to the drift kinetic resonance effects and FOW correction of fast ions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Magnetism and quantum melting in moiré-material Wigner crystals

Recent experiments have established that semiconductor-based moir´e materials can host incompressible states at a series of fractional moir´e-miniband fillings. These states have been identified as generalized Wigner crystals in which electrons localize on a subset of the available triangular-lattice moir´e superlattice sites. In this article, we use momentum-space exact diagonalization to investigate the many-body ground state evolution at rational fillings from the weak-hopping classical lattice gas limit, in which only spin degrees-of-freedom are active at low energies, to the strong-hopping metallic regime where the Wigner crystals melt. We specifically address the nature of the magnetic ground states of the generalized Wigner crystals at fillings ν = 1/3 and ν = 2/3.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Effects of strain and defects on the topological properties of HfTe 5

Manipulating topological states is a topic of vigorous research. However, the impact of such topological phase transitions on the quasiparticle dynamics remains elusive. Here, in this work, we present the effects of a transition from a strong to weak topological insulator in HfTe 5 as a function of Te vacancy concentration. We observed a significant transition from distinct sharp surface states and Dirac crossing to a Fermi-liquid-like quasiparticle state in which these surface-localized features are heavily suppressed. Additionally, by inducing the same effect through applied uniaxial stress, we demonstrate that changes in the lattice constants play the foremost role in determining the electronic structure, self-energy, and topological states of HfTe 5 . Our results demonstrate the possibility of using both defect chemistry and strain as control parameters for topological phase transitions and associated many-body physics.

36 MATERIALS SCIENCE↗
DOE OSTI2017

Quantified Gamow shell model interaction for p s d -shell nuclei

Background: The structure of weakly bound and unbound nuclei close to particle drip lines is one of the major science drivers of nuclear physics. A comprehensive understanding of these systems goes beyond the traditional configuration interaction approach formulated in the Hilbert space of localized states (nuclear shell model) and requires an open quantum system description. The complex-energy Gamow shell model (GSM) provides such a framework as it is capable of describing resonant and nonresonant many-body states on equal footing. Purpose: To make reliable predictions, quality input is needed that allows for the full uncertainty quantification of theoretical results. In this study, we carry out the optimization of an effective GSM (one-body and two-body) interaction in the psdf-shell-model space. The resulting interaction is expected to describe nuclei with 5 ≤ A ≲ 12 at the p-sd-shell interface. Method: The one-body potential of the 4 He core is modeled by a Woods-Saxon + spin-orbit + Coulomb potential, and the finite-range nucleon-nucleon interaction between the valence nucleons consists of central, spin-orbit, tensor, and Coulomb terms. The GSM is used to compute key fit observables. The χ 2 optimization is performed using the Gauss-Newton algorithm augmented by the singular value decomposition technique. The resulting covariance matrix enables quantification of statistical errors within the linear regression approach. Results: The optimized one-body potential reproduces nucleon- 4 He scattering phase shifts up to an excitation energy of 20 MeV. The two-body interaction built on top of the optimized one-body field is adjusted to the bound and unbound ground-state binding energies and selected excited states of the helium, lithium, and beryllium isotopes up to A = 9 . A very good agreement with experimental results was obtained for binding energies. First applications of the optimized interaction include predictions for two-nucleon correlation densities and excitation spectra of light nuclei with quantified uncertainties. In conclusion: The new interaction will enable comprehensive and fully quantified studies of structure and reactions aspects of nuclei from the psd region of the nuclear chart.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2018

Empirical Determination of Dark Matter Velocities Using Metal-Poor Stars

The Milky Way dark matter halo is formed from the accretion of smaller subhalos. These sub-units also harbor stars—typically old and metal-poor—that are deposited in the Galactic inner regions by disruption events. In this report, we show that the dark matter and metal-poor stars in the Solar neighborhood share similar kinematics due to their common origin. Using the high-resolution eris simulation, which traces the evolution of both the dark matter and baryons in a realistic Milky Way analog galaxy, we demonstrate that metal-poor stars are indeed effective tracers for the local, virialized dark matter velocity distribution. The local dark matter velocities can therefore be inferred from observations of the stellar halo made by the Sloan Digital Sky Survey within 4 kpc of the Sun. This empirical distribution differs from the standard halo model in important ways and suggests that the bounds on the spin-independent scattering cross section may be weakened for dark matter masses below ~10 GeV. Data from Gaia will allow us to further refine the expected distribution for the smooth dark matter component, and to test for the presence of local substructure.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

Formation of gapless Z{sub 2} spin liquid phase manganites in the (Sm{sub 1–y}Gd{sub y}){sub 0.55}Sr{sub 0.45}MnO{sub 3} system in zero magnetic field: Topological phase transitions to states with low and high density of 2D-vortex pairs induced by the magnetic field

The evolution of the ground state of the manganese spin ensemble in the (Sm{sub 1–y}Gd{sub y}){sub 0.55}Sr{sub 0.45}MnO{sub 3} in the case of isovalent substitution of rare-earth samarium ions with large radii with gadolinium ions with significantly smaller radii is studied. The measured temperature dependences of the ac magnetic susceptibility and the field dependences of the dc magnetizations are analyzed using the Heisenberg–Kitaev model describing the transition from the ordered spin state with classical isotropic AFM exchange to the frustrated spin state with quantum highly anisotropic FM exchange. A continuous transition from the 3D ferromagnetic state of manganese spins in the initial sample with y = 0 to zigzag AFM ordering of CE-type spins in ab planes for y = 0.5, coexisting in samples with y = 0.5, 0.6, and 0.7 at temperatures below T{sub N} ≅ 48.5 K with a disordered phase such as a quantum Griffiths phase is identified. As the gadolinium concentration further increases, the CE-type zigzag AFM structure is molten, which leads to the appearance of an unusual phase in Gd{sub 0.55}Sr{sub 0.45}MnO{sub 3} in the temperature range close to the absolute zero. This phase has characteristic features of a gapless Z{sub 2} quantum spin liquid in zero external magnetic field. The step changes in the magnetization isotherms measured at 4.2 K in the field range of ±75 kOe are explained by quantum phase transitions of the Z{sub 2} spin liquid to a phase with topological order in weak magnetic fields and a polarized phase in strong fields. The significant difference between critical fields and magnetization jumps in isotherms indicates the existence of hysteretic phenomena in quantum spin liquid magnetization–demagnetization processes caused by the difference between localization–delocalization of 2D vortex pairs induced by a magnetic field in a quantum spin liquid with disorder.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2016

STUDYING THE INTERSTELLAR MEDIUM AND THE INNER REGION OF NPS/LOOP 1 WITH SHADOW OBSERVATIONS TOWARD MBM36

We analyzed data from a shadow observation of the high density molecular cloud MBM36 (l ∼ 4°, b ∼ 35°) with Suzaku. MBM36 is located in a region that emits relatively weakly in the 3/4 keV band compared to the surrounding North Polar Spur (NPS)/Loop 1 structure and the Galactic Bulge (GB). The contrast between high and low density targets in the MBM36 area allows one to separate the local and distant contributors to the soft diffuse X-ray background, providing a much better characterization of the individual components compared to single pointing observations. We identify two non-local thermal components, one at kT ≈ 0.12 keV and one at kT ≈ 0.29 keV. The colder component matches well with models of emission from the higher latitude region of the GB. The emission of the warmer component is in agreement with models predicting that the NPS is due to a hypershell from the center of the Milky Way. Geometrical and pressure calculations rule out a nearby bubble as responsible for the emission associated with the NPS. Any Galactic Halo/circumgalactic halo emission, if present, is outshined by the other components. We also report an excess emission around 0.9 keV, likely due to an overabundance of Ne ix.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1994

Effect of high-latitude ionospheric convection on Sun-aligned polar caps

A coupled magnetospheric-ionospheric (M-I) magnetohydrodynamic (MHD) model has been used to simulate the formation of Sun-aligned polar cap arcs for a variety of interplanetary magnetic field (IMF) dependent polar cap convection fields. The formation process involves launching an Alfven shear wave from the magnetosphere to the ionosphere where the ionospheric conductance can react self-consistently to changes in the upward currents. We assume that the initial Alfven shear wave is the result of solar wind-magnetosphere interactions. The simulations show how the E region density is affected by the changes in the electron precipitation that are associated with the upward currents. These changes in conductance lead to both a modified Alfven wave reflection at the ionosphere and the generation of secondary Alfven waves in the ionosphere. The ensuing bouncing of the Alfven waves between the ionosphere and magnetosphere is followed until an asymptotic solution is obtained. At the magnetosphere the Alfven waves reflect at a fixed boundary. The coupled M-I Sun-aligned polar cap arc model of Zhu et al.(1993a) is used to carry out the simulations. This study focuses on the dependence of the polar cap arc formation on the background (global) convection pattern. Since the polar cap arcs occur for northward and strong B(sub y) IMF conditions, a variety of background convection patterns can exist when the arcs are present. The study shows that polar cap arcs can be formed for all these convection patterns; however, the arc features are dramatically different for the different patterns. For weak sunward convection a relatively confined single pair of current sheets is associated with the imposed Alfven shear wave structure. However, when the electric field exceeds a threshold, the arc structure intensifies, and the conductance increases as does the local Joule heating rate. These increases are faster than a linear dependence on the background electric field strength. Furthermore, above the threshold, the single current sheet pair splits into multiple current sheet pairs. For the fixed initial ionospheric and magnetospheric conditions used in this study, the separation distance between the current pairs was found to be almost independent of the background electric field strength. For either three-cell or distorted two-cell background convection patterns the arc formation favored the positive B(sub y) case in the northern hemisphere.

Sojka, J. J.↗
DOE OSTI2015

A high-transparency, micro-patternable chip for X-ray diffraction analysis of microcrystals under native growth conditions

A highly X-ray-transparent, silicon nitride-based device has been designed and fabricated to harvest protein microcrystals for high-resolution X-ray diffraction data collection using microfocus beamlines and XFELs. Microcrystals present a significant impediment to the determination of macromolecular structures by X-ray diffraction methods. Although microfocus synchrotron beamlines and X-ray free-electron lasers (XFELs) can enable the collection of interpretable diffraction data from microcrystals, there is a need for efficient methods of harvesting small volumes (<2 µl) of microcrystals grown under common laboratory formats and delivering them to an X-ray beam source under native growth conditions. One approach that shows promise in overcoming the challenges intrinsic to microcrystal analysis is to pair so-called ‘fixed-target’ sample-delivery devices with microbeam-based X-ray diffraction methods. However, to record weak diffraction patterns it is necessary to fabricate devices from X-ray-transparent materials that minimize background scattering. Presented here is the design of a new micro-diffraction device consisting of three layers fabricated from silicon nitride, photoresist and polyimide film. The chip features low X-ray scattering and X-ray absorption properties, and uses a customizable blend of hydrophobic and hydrophilic surface patterns to help localize microcrystals to defined regions. Microcrystals in their native growth conditions can be loaded into the chips with a standard pipette, allowing data collection at room temperature. Diffraction data collected from hen egg-white lysozyme microcrystals (10–15 µm) loaded into the chips yielded a complete, high-resolution (<1.6 Å) data set sufficient to determine a high-quality structure by molecular replacement. The features of the chip allow the rapid and user-friendly analysis of microcrystals grown under virtually any laboratory format at microfocus synchrotron beamlines and XFELs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1988

Particle-laden weakly swirling free jets: Measurements and predictions

A theoretical and experimental investigation of particle-laden, weakly swirling, turbulent free jets was conducted. Glass particles, having a Sauter mean diameter of 39 microns, with a standard deviation of 15 microns, were used. A single loading ratio (the mass flow rate of particles per unit mass flow rate of air) of 0.2 was used in the experiments. Measurements are reported for three swirl numbers, ranging from 0 to 0.33. The measurements included mean and fluctuating velocities of both phases, and particle mass flux distributions. Measurements were also completed for single-phase non-swirling and swirling jets, as baselines. Measurements were compared with predictions from three types of multiphase flow analysis, as follows: (1) locally homogeneous flow (LHF) where slip between the phases was neglected; (2) deterministic separated flow (DSF), where slip was considered but effects of turbulence/particle interactions were neglected; and (3) stochastic separated flow (SSF), where effects of both interphase slip and turbulence/particle interactions were considered using random sampling for turbulence properties in conjunction with random-walk computations for particle motion. Single-phase weakly swirling jets were considered first. Predictions using a standard k-epsilon turbulence model, as well as two versions modified to account for effects of streamline curvature, were compared with measurements. Predictions using a streamline curvature modification based on the flux Richardson number gave better agreement with measurements for the single-phase swirling jets than the standard k-epsilon model. For the particle-laden jets, the LHF and DSF models did not provide very satisfactory predictions. The LHF model generally overestimated the rate of decay of particle mean axial and angular velocities with streamwise distance, and predicted particle mass fluxes also showed poor agreement with measurements, due to the assumption of no-slip between phases. The DSF model also performed quite poorly for predictions of particle mass flux because turbulent dispersion of the particles was neglected. The SSF model, which accounts for both particle inertia and turbulent dispersion of the particles, yielded reasonably good predictions throughout the flow field for the particle-laden jets.

Bulzan, Daniel L.↗
DOE OSTI2019

On the Mn 4+ R-line emission intensity and its tunability in solids

The spectroscopic properties of Mn 4+ ions in solids are analyzed to establish key parameters that are responsible for the energy and the intensity of the 2 E g → 4 A 2g zero-phonon emission transition (R-line). These two factors of the emission spectrum are critical in determining the phosphor luminosity or brightness. The emission energy is chiefly determined by the “Mn 4+ -ligand” bonding covalence. In fluorides, the Mn-F bonding covalence is weak, and the R-line occur at higher energy. In oxides, the stronger Mn-O bonding covalence reduces the energy of the 2 E g state. In oxides, the R-line shifts to higher energy with increasing local distortion of the MnO 6 octahedral moiety. Those host crystals that locate the Mn 4+ ion at a noncentrosymmetric site (absence of inversion symmetry) relax the parity section rule, which enhances the R-line intensity. Increased R-line intensity in the emission spectrum produces color changes that increases phosphor luminosity. The influence of second coordination sphere (next-nearest neighbors) on the R-line intensity is considered. This article establishes guidelines for the development of host crystals that meet the requirements of high luminosity phosphors for general illumination applications.

36 MATERIALS SCIENCE↗
DOE OSTI2015

Imaging ballistic carrier trajectories in graphene using scanning gate microscopy

We use scanning gate microscopy to map out the trajectories of ballistic carriers in high-mobility graphene encapsulated by hexagonal boron nitride and subject to a weak magnetic field. We employ a magnetic focusing geometry to image carriers that emerge ballistically from an injector, follow a cyclotron path due to the Lorentz force from an applied magnetic field, and land on an adjacent collector probe. The local electric field generated by the scanning tip in the vicinity of the carriers deflects their trajectories, modifying the proportion of carriers focused into the collector. By measuring the voltage at the collector while scanning the tip, we are able to obtain images with arcs that are consistent with the expected cyclotron motion. We also demonstrate that the tip can be used to redirect misaligned carriers back to the collector.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

Simulation study of the effects of polymer network dynamics and mesh confinement on the diffusion and structural relaxation of penetrants

The diffusion of small molecular penetrants through polymeric materials represents an important fundamental problem, relevant to the design of materials for applications such as coatings and membranes. Polymer networks hold promise in these applications because dramatic differences in molecular diffusion can result from subtle changes in the network structure. In this paper, we use molecular simulation to understand the role that cross-linked network polymers have in governing the molecular motion of penetrants. By considering the local, activated alpha relaxation time of the penetrant and its long-time diffusive dynamics, we can determine the relative importance of activated glassy dynamics on penetrants at the segmental scale vs entropic mesh confinement on penetrant diffusion. We vary several parameters, such as the cross-linking density, temperature, and penetrant size, to show that cross-links primarily affect molecular diffusion through the modification of the matrix glass transition, with local penetrant hopping at least partially coupled to the segmental relaxation of the polymer network. This coupling is very sensitive to the local activated segmental dynamics of the surrounding matrix, and we also show that penetrant transport is affected by dynamic heterogeneity at low temperatures. To contrast, only at high temperatures and for large penetrants or when the dynamic heterogeneity effect is weak, does the effect of mesh confinement become significant, even though penetrant diffusion more broadly empirically follows similar trends as established models of mesh confinement-based transport.

Chemistry↗
DOE OSTI2021

Network analysis of earthquake ground motion spatial correlation: a case study with the San Jacinto seismic nodal array

SUMMARY The spatial correlation of earthquake ground motion intensity can be measured from strong motion data; however, the data used in past studies is sparsely sampled in space, and only the interstation distance was considered as a correlation variable. These limitations mean that we have only weak constraints on the true correlation structure of ground motion and that potentially important aspects of spatial correlation are unconstrained. In this study, we combine a large-N seismic array and graph analytics to explore this issue at a local scale using small local and regional earthquakes. Our result suggests site conditions, and how they interact with the incident seismic wavefield, strongly condition the spatial correlation of ground motion. Future progress in characterizing ground motion spatial variability will require dense wavefield measurements, either through nodal deployments, or perhaps distributed acoustic sensing measurements, of seismic wavefields. Aftershock sequences of major earthquakes would provide particularly data-rich targets of opportunity.

58 GEOSCIENCES↗
DOE OSTI2018

Passive state preparation in the Gaussian-modulated coherent-states quantum key distribution

In the Gaussian-modulated coherent-states (GMCS) quantum key distribution (QKD) protocol, Alice prepares quantum states actively: For each transmission, Alice generates a pair of Gaussian-distributed random numbers, encodes them on a weak coherent pulse using optical amplitude and phase modulators, and then transmits the Gaussian-modulated weak coherent pulse to Bob. Here we propose a passive state preparation scheme using a thermal source. In our scheme, Alice splits the output of a thermal source into two spatial modes using a beam splitter. She measures one mode locally using conjugate optical homodyne detectors, and transmits the other mode to Bob after applying appropriate optical attenuation. Under normal conditions, Alice's measurement results are correlated to Bob's, and they can work out a secure key, as in the active state preparation scheme. Given the initial thermal state generated by the source is strong enough, this scheme can tolerate high detector noise at Alice's side. Furthermore, the output of the source does not need to be single mode, since an optical homodyne detector can selectively measure a single mode determined by the local oscillator. Preliminary experimental results suggest that the proposed scheme could be implemented using an off-the-shelf amplified spontaneous emission source.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Spin-orbit coupled systems in the atomic limit: rhenates, osmates, iridates

Motivated by RIXS experiments on a wide range of complex heavy oxides, including rhenates, osmates, and iridates, we discuss the theory of RIXS for site-localized t 2 g orbital systems with strong spin-orbit coupling. For such systems, we present exact diagonalization results for the spectrum at different electron fillings, showing that it accesses “single-particle” and “multiparticle” excitations. This leads to a simple picture for the energies and intensities of the RIXS spectra in Mott insulators such as double perovskites which feature highly localized electrons, and yields estimates of the spin-orbit coupling and Hund's coupling in correlated 5d oxides. We present new higher resolution RIXS data at the Re L 3 edge in Ba 2 YReO 6 which finds a previously unresolved peak splitting, providing further confirmation of our theoretical predictions. Using ab initio electronic structure calculations on Ba 2 MReO 6 (with M=Re, Os, Ir) we show that while the atomic limit yields a reasonable effective Hamiltonian description of the experimental observations, effects such as t 2 g –e g interactions and hybridization with oxygen are important. Our ab initio estimate for the strength of the intersite exchange coupling shows that, compared to the d 3 systems, the exchange is one or two orders of magnitude weaker in the d 2 and d 4 materials, which may partly explain the suppression of long-range magnetic order in the latter compounds. As a way to interpolate between the site-localized picture and our electronic structure band calculations, we discuss the spin-orbital levels of the MO 6 cluster. In conclusion, this suggests a possible role for intracluster excitons in Ba 2 YIrO 6 which may lead to a weak breakdown of the atomic J eff = 0 picture and to small magnetic moments.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

The Role of Strong Magnetic Fields in Stabilizing Highly Luminous Thin Disks

Abstract We present and analyze a set of three-dimensional, global, general relativistic radiation magnetohydrodynamic simulations of thin, radiation-pressure-dominated accretion disks surrounding a nonrotating, stellar-mass black hole. The simulations are initialized using the Shakura–Sunyaev model with a mass accretion rate of M ̇ = 3 L Edd / c 2 (corresponding to L = 0.17 L Edd ). Our previous work demonstrated that such disks are thermally unstable when accretion is driven by an α -viscosity. In the present work, we test the hypothesis that strong magnetic fields can both drive accretion through magnetorotational instability and restore stability to such disks. We test four initial magnetic field configurations: (1) a zero-net-flux case with a single, radially extended set of magnetic field loops (dipole), (2) a zero-net-flux case with two radially extended sets of magnetic field loops of opposite polarity stacked vertically (quadrupole), (3) a zero-net-flux case with multiple radially concentric rings of alternating polarity (multiloop), and (4) a net-flux, vertical magnetic field configuration (vertical). In all cases, the fields are initially weak, with a gas-to-magnetic pressure ratio ≳100. Based on the results of these simulations, we find that the dipole and multiloop configurations remain thermally unstable like their α -viscosity counterpart, in our case collapsing vertically on the local thermal timescale and never fully recovering. The vertical case, on the other hand, stabilizes and remains so for the duration of our tests (many thermal timescales). The quadrupole case is intermediate, showing signs of both stability and instability. The key stabilizing factor is the ability of specific field configurations to build up and sustain strong, P mag ≳ 0.5 P tot , toroidal fields near the midplane of the disk. We discuss the reasons why certain configurations are able to do this effectively and others are not. We then compare our stable simulations to the standard Shakura–Sunyaev disk.

79 ASTRONOMY AND ASTROPHYSICS↗