Search NASA⌕ Search

SEARCH · Search NASA

Results for “MAGNETIC STARS”

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,189 records · Page 66

RADIO ASTROMETRY OF THE CLOSE ACTIVE BINARY HR 5110

The close active binary HR 5110 was observed at six epochs over 26 days using a global very long baseline interferometry array at 15.4 GHz. We used phase referencing to determine the position of the radio centroid at each epoch with an uncertainty significantly smaller than the component separation. After correcting for proper motion and parallax, we find that the centroid locations of all six epochs have barycenter separations consistent with an emission source located on the KIV secondary, and not in an interaction region between the stars or on the F primary. We used a homogeneous power-law gyrosynchrotron emission model to reproduce the observed flux densities and fractional circular polarization. The resulting ranges of mean magnetic field strength and relativistic electron densities are of the order of 10 G and 10{sup 5} cm{sup −3}, respectively, in the source region.

79 ASTRONOMY AND ASTROPHYSICS↗

Radiation magnetohydrodynamic simulations of protostellar collapse: Low-metallicity environments

Among many physical processes involved in star formation, radiation transfer is one of the key processes because it dominantly controls the thermodynamics. Because metallicities control opacities, they are one of the important environmental parameters that affect star formation processes. In this work, I investigate protostellar collapse in solar-metallicity and low-metallicity (Z = 0.1 Z {sub ☉}) environments using three-dimensional radiation hydrodynamic and magnetohydrodynamic simulations. Because radiation cooling in high-density gas is more effective in low-metallicity environments, first cores are colder and have lower entropies. As a result, first cores are smaller, less massive, and have shorter lifetimes in low-metallicity clouds. Therefore, first cores would be less likely to be found in low-metallicity star forming clouds. This also implies that first cores tend to be more gravitationally unstable and susceptible to fragmentation. The evolution and structure of protostellar cores formed after the second collapse weakly depend on metallicities in the spherical and magnetized models, despite the large difference in the metallicities. Because this is due to the change of the heat capacity by dissociation and ionization of hydrogen, it is a general consequence of the second collapse as long as the effects of radiation cooling are not very large during the second collapse. On the other hand, the effects of different metallicities are more significant in the rotating models without magnetic fields, because they evolve slower than other models and therefore are more affected by radiation cooling.

79 ASTRONOMY AND ASTROPHYSICS↗

Where Are the r-modes? Chandra Observations of Millisecond Pulsars

We present the results of Chandra observations of two non-accreting millisecond pulsars, PSRs J1640+2224 (J1640) and J1709+2313 (J1709), with low inferred magnetic fields and spin-down rates in order to constrain their surface temperatures, obtain limits on the amplitude of unstable r -modes in them, and make comparisons with similar limits obtained for a sample of accreting low-mass X-ray binary (LMXB) neutron stars. We detect both pulsars in the X-ray band for the first time. They are faint, with inferred soft X-ray fluxes (0.3–3 keV) of ≈6 × 10{sup −15} and 3 × 10{sup −15} erg cm{sup −2} s{sup −1} for J1640 and J1709, respectively. Spectral analysis assuming hydrogen atmosphere emission gives global effective temperature upper limits (90% confidence) of 3.3–4.3 × 10{sup 5} K for J1640 and 3.6–4.7 × 10{sup 5} K for J1709, where the low end of the range corresponds to canonical neutron stars ( M = 1.4 M {sub ⊙}), and the upper end corresponds to higher-mass stars ( M = 2.21 M {sub ⊙}). Under the assumption that r -mode heating provides the thermal support, we obtain dimensionless r -mode amplitude upper limits of 3.2–4.8 × 10{sup −8} and 1.8–2.8 × 10{sup −7} for J1640 and J1709, respectively, where again the low end of the range corresponds to lower-mass, canonical neutron stars ( M = 1.4 M {sub ⊙}). These limits are about an order of magnitude lower than those we derived previously for a sample of LMXBs, except for the accreting millisecond X-ray pulsar SAX J1808.4–3658, which has a comparable amplitude limit to J1640 and J1709.

79 ASTRONOMY AND ASTROPHYSICS↗

The 2014–2017 outburst of the young star ASASSN-13db: A time-resolved picture of a very-low-mass star between EXors and FUors*

Context. Accretion outbursts are key elements in star formation. ASASSN-13db is a M5-type star with a protoplanetary disk, the lowest-mass star known to experience accretion outbursts. Since its discovery in 2013, it has experienced two outbursts, the second of which started in November 2014 and lasted until February 2017. Aims. We observe the photometric and spectroscopic behavior of ASASSN-13db during the 2014–2017 outburst. Methods. We use high- and low-resolution spectroscopy and time-resolved photometry from the ASAS-SN survey, the LCOGT and the Beacon Observatory to study the light curve of ASASSN-13db and the dynamical and physical properties of the accretion flow. Results. The 2014–2017 outburst lasted for nearly 800 days. A 4.15 d period in the light curve likely corresponds to rotational modulation of a star with hot spot(s). The spectra demonstrate multiple emission lines with variable inverse P-Cygni profiles and a highly variable blue-shifted absorption below the continuum. Line ratios from metallic emission lines (Fe I/Fe II, Ti I/Ti II) suggest temperatures of ~5800–6000 K in the accretion flow. Conclusions. Photometrically and spectroscopically, the 2014–2017 event displays an intermediate behavior between EXors and FUors. The accretion rate ([$\dot{M}$]= 1–3 × 10 -7 M ⊙ /yr), about two orders of magnitude higher than the accretion rate in quiescence, is not vastly different from the accretion rate observed in 2013. The absorption features in the spectra suggest that the system is viewed at a high angle and drives a powerful, non-axisymmetric wind, maybe related to magnetic reconnection. The properties of ASASSN-13db suggest that temperatures lower than those for solar-type stars are needed for modeling accretion in very-low-mass systems. Finally, the rotational modulation during the outburst reveals that accretion-related structures settle after the beginning of the outburst and can be relatively stable and long-lived. Our work also demonstrates the power of time-resolved photometry and spectroscopy to explore the properties of variable and outbursting stars.

79 ASTRONOMY AND ASTROPHYSICS↗

Near 3:2 and 2:1 mean motion resonance formation in the systems observed by Kepler

The Kepler mission has released ∼4229 transiting planet candidates. There are approximately 222 candidate systems with three planets. Among them, the period ratios of planet pairs near 1.5 and 2.0 reveal that two peaks exist for which the proportions of the candidate systems are ∼7.0% and 18.0%, respectively. In this work, we study the formation of mean motion resonance (MMR) systems, particularly for the planetary configurations near 3:2 and 2:1 MMRs, and we concentrate on the interplay between the resonant configuration and the combination of stellar accretion rate, stellar magnetic field, speed of migration, and additional planets. We perform more than 1000 runs by assuming a system with a solar-like star and three surrounding planets. From the statistical results, we find that under the formation scenario, the proportions near 1.5 and 2.0 can reach 14.5% and 26.0%, respectively. In addition, M-dot =0.1×10{sup −8} M{sub ⊙} yr{sup −1} is propitious toward the formation of 3:2 resonance, whereas M-dot =2×10{sup −8} M{sub ⊙} yr{sup −1} contributes to the formation of 2:1 resonance. The speed-reduction factor of type I migration f {sub 1} ≥ 0.3 facilitates 3:2 MMRs, whereas f {sub 1} ≥ 0.1 facilitates 2:1 MMRs. If additional planets are present in orbits within the innermost or beyond the outermost planet in a three-planet system, 3:2:1 MMRs can be formed, but the original systems trapped in 4:2:1 MMRs are not affected by the supposed planets. In summary, we conclude that this formation scenario will provide a likely explanation for Kepler candidates involved in 2:1 and 3:2 MMRs.

79 ASTRONOMY AND ASTROPHYSICS↗

Stellar Winds and Dust Avalanches in the AU Mic Debris Disk

We explain the fast-moving, ripple-like features in the edge-on debris disk orbiting the young M dwarf AU Mic. The bright features are clouds of submicron dust repelled by the host star’s wind. The clouds are produced by avalanches: radial outflows of dust that gain exponentially more mass as they shatter background disk particles in collisional chain reactions. The avalanches are triggered from a region a few au across—the “avalanche zone”—located on AU Mic’s primary “birth” ring at a true distance of ∼35 au from the star but at a projected distance more than a factor of 10 smaller: the avalanche zone sits directly along the line of sight to the star, on the side of the ring nearest Earth, launching clouds that disk rotation sends wholly to the southeast, as observed. The avalanche zone marks where the primary ring intersects a secondary ring of debris left by the catastrophic disruption of a progenitor up to Varuna in size, less than tens of thousands of years ago. Only where the rings intersect are particle collisions sufficiently violent to spawn the submicron dust needed to seed the avalanches. We show that this picture works quantitatively, reproducing the masses, sizes, and velocities of the observed escaping clouds. The Lorentz force exerted by the wind’s magnetic field, whose polarity reverses periodically according to the stellar magnetic cycle, promises to explain the observed vertical undulations. The timescale between avalanches, about 10 yr, might be set by time variability of the wind mass loss rate or, more speculatively, by some self-regulating limit cycle.

79 ASTRONOMY AND ASTROPHYSICS↗

An experimental platform for investigating astrophysically relevant magnetized plasma jets on the COBRA facility

A new platform has been developed for the 1-MA COBRA generator to investigate the physical processes affecting the formation, collimation, and stability of high-speed outflows in magnetically driven laboratory plasma jets. Such experiments serve as diagnostically accessible surrogates for astrophysical jets under the assumption that the underlying dynamics are scale invariant. In contrast to previous current driven high energy density laboratory jet experiments that use radial/conical wire arrays or foils, the platform described here uses azimuthally symmetric gas-puff injection. This avoids the ablation phase from a solid target, allowing the jets to develop earlier and be driven longer without depleting their mass source and disrupting. A permanent magnet provides an initial poloidal magnetic field, which links the two concentric electrodes and mimics the boundary conditions of a star-accretion disk system. Extended magnetohydrodynamic effects can be assessed using a polarity convolute, which allows for reversal of the electrode bias. The resulting plasma jets exhibit remarkable stability, persisting for hundreds of nanoseconds and achieving aspect ratios ≳30:1.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

PHASE-RESOLVED INFRARED SPECTROSCOPY AND PHOTOMETRY OF V1500 CYGNI, AND A SEARCH FOR SIMILAR OLD CLASSICAL NOVAE

We present phase-resolved near-infrared photometry and spectroscopy of the classical nova (CN) V1500 Cyg to explore whether cyclotron emission is present in this system. While the spectroscopy do not indicate the presence of discrete cyclotron harmonic emission, the light curves suggest that a sizable fraction of its near-infrared fluxes are due to this component. The light curves of V1500 Cyg appear to remain dominated by emission from the heated face of the secondary star in this system. We have used infrared spectroscopy and photometry to search for other potential magnetic systems among old CNe. We have found that the infrared light curves of V1974 Cyg superficially resemble those of V1500 Cyg, suggesting a highly irradiated companion. The old novae V446 Her and QV Vul have light curves with large amplitude variations like those seen in polars, suggesting they might have magnetic primaries. We extract photometry for 79 old novae from the Two Micron All Sky Survey Point Source Catalog and use those data to derive the mean, un-reddened infrared colors of quiescent novae. We also extract WISE data for these objects and find that 45 of them were detected. Surprisingly, a number of these systems were detected in the WISE 22 {mu}m band. While two of those objects produced significant dust shells (V705 Cas and V445 Pup), the others did not. It appears that line emission from their ionized ejected shells is the most likely explanation for those detections.

79 ASTRONOMY AND ASTROPHYSICS↗

A COMPREHENSIVE PHOTOMETRIC STUDY OF THE CONTACT BINARY GN BOO WITH POSSIBLE MAGNETIC ACTIVITIES

Multi-color photometric data of GN Boo observed from 2010 to 2013 are presented. The intrinsic variations of the light curves are remarkable, and their phenomena are probable evidence of stronger magnetic activities on the surfaces of the components. Based on all CCD times of minimum light, a secular increase superimposed on a cyclic oscillation is found. The orbital period increases at a rate of dP/dt = +1.74 × 10{sup −7} days yr{sup −1}, which can be explained by mass transfer from the less massive component to the more massive one. The period and amplitude of the cyclic variation are P = 9.5632 yr and A = 0.0046 days, respectively, which correspond to the previous published results. Using the 2010 version of the W–D code, five sets of photometric solutions were derived from our new data. The results imply that the stellar spot, the degree of fill-out, and the temperature difference ΔT between the components of GN Boo are variable. It is inferred that the magnetic activities perhaps influence the outer radius of the component, causing the temperature of the component star and the level of contact to change.

79 ASTRONOMY AND ASTROPHYSICS↗

EARLY STAGES OF CLUSTER FORMATION: FRAGMENTATION OF MASSIVE DENSE CORES DOWN TO {approx}< 1000 AU

In order to study the fragmentation of massive dense cores, which constitute the cluster cradles, we observed the continuum at 1.3 mm and the CO (2-1) emission of four massive cores with the Plateau de Bure Interferometer in the most extended configuration. We detected dust condensations down to {approx}0.3 M {sub Sun} and separate millimeter sources down to 0.''4 or {approx}< 1000 AU, comparable to the sensitivities and separations reached in optical/infrared studies of clusters. The CO (2-1) high angular resolution images reveal high-velocity knots usually aligned with previously known outflow directions. This, in combination with additional cores from the literature observed at similar mass sensitivity and spatial resolution, allowed us to build a sample of 18 protoclusters with luminosities spanning three orders of magnitude. Among the 18 regions, {approx}30% show no signs of fragmentation, while 50% split up into {approx}> 4 millimeter sources. We compiled a list of properties for the 18 massive dense cores, such as bolometric luminosity, total mass, and mean density, and found no correlation of any of these parameters with the fragmentation level. In order to investigate the combined effects of the magnetic field, radiative feedback, and turbulence in the fragmentation process, we compared our observations to radiation magnetohydrodynamic simulations and found that the low-fragmented regions are reproduced well in the magnetized core case, while the highly fragmented regions are consistent with cores where turbulence dominates over the magnetic field. Overall, our study suggests that the fragmentation in massive dense cores could be determined by the initial magnetic field/turbulence balance in each particular core.

79 ASTRONOMY AND ASTROPHYSICS↗

DELO-BEZIER FORMAL SOLUTIONS OF THE POLARIZED RADIATIVE TRANSFER EQUATION

We present two new accurate and efficient methods to compute the formal solution of the polarized radiative transfer equation. In this work, the source function and the absorption matrix are approximated using quadratic and cubic Bezier spline interpolants. These schemes provide second- and third-order approximations, respectively, and do not suffer from erratic behavior of the polynomial approximation (overshooting). The accuracy and the convergence of the new method are studied along with other popular solutions of the radiative transfer equation, using stellar atmospheres with strong gradients in the line-of-sight velocity and in the magnetic-field vector.

79 ASTRONOMY AND ASTROPHYSICS↗

LOW-MASS PLANETS IN PROTOPLANETARY DISKS WITH NET VERTICAL MAGNETIC FIELDS: THE PLANETARY WAKE AND GAP OPENING

Some regions in protoplanetary disks are turbulent, while some regions are quiescent (e.g. the dead zone). In order to study how planets open gaps in both inviscid hydrodynamic disk (e.g. the dead zone) and the disk subject to magnetorotational instability (MRI), we carried out both shearing box two-dimensional inviscid hydrodynamical simulations and three-dimensional unstratified magnetohydrodynamical (MHD) simulations (having net vertical magnetic fields) with a planet at the box center. We found that, due to the nonlinear wave steepening, even a low mass planet can open gaps in both cases, in contradiction to the ''thermal criterion'' for gap opening. In order to understand if we can represent the MRI turbulent stress with the viscous {alpha} prescription for studying gap opening, we compare gap properties in MRI-turbulent disks to those in viscous HD disks having the same stress, and found that the same mass planet opens a significantly deeper and wider gap in net vertical flux MHD disks than in viscous HD disks. This difference arises due to the efficient magnetic field transport into the gap region in MRI disks, leading to a larger effective {alpha} within the gap. Thus, across the gap, the Maxwell stress profile is smoother than the gap density profile, and a deeper gap is needed for the Maxwell stress gradient to balance the planetary torque density. Comparison with previous results from net toroidal flux/zero flux MHD simulations indicates that the magnetic field geometry plays an important role in the gap opening process. We also found that long-lived density features (termed zonal flows) produced by the MRI can affect planet migration. Overall, our results suggest that gaps can be commonly produced by low mass planets in realistic protoplanetary disks, and caution the use of a constant {alpha}-viscosity to model gaps in protoplanetary disks.

79 ASTRONOMY AND ASTROPHYSICS↗

MISALIGNMENT OF MAGNETIC FIELDS AND OUTFLOWS IN PROTOSTELLAR CORES

We present results of {lambda}1.3 mm dust-polarization observations toward 16 nearby, low-mass protostars, mapped with {approx}2.''5 resolution at CARMA. The results show that magnetic fields in protostellar cores on scales of {approx}1000 AU are not tightly aligned with outflows from the protostars. Rather, the data are consistent with scenarios where outflows and magnetic fields are preferentially misaligned (perpendicular), or where they are randomly aligned. If one assumes that outflows emerge along the rotation axes of circumstellar disks, and that the outflows have not disrupted the fields in the surrounding material, then our results imply that the disks are not aligned with the fields in the cores from which they formed.

79 ASTRONOMY AND ASTROPHYSICS↗

SURVIVAL OF INTERSTELLAR MOLECULES TO PRESTELLAR DENSE CORE COLLAPSE AND EARLY PHASES OF DISK FORMATION

An outstanding question of astrobiology is the link between the chemical composition of planets, comets, and other solar system bodies and the molecules formed in the interstellar medium. Understanding the chemical and physical evolution of the matter leading to the formation of protoplanetary disks is an important step for this. We provide some new clues to this long-standing problem using three-dimensional chemical simulations of the early phases of disk formation: we interfaced the full gas-grain chemical model Nautilus with the radiation-magnetohydrodynamic model RAMSES, for different configurations and intensities of the magnetic field. Our results show that the chemical content (gas and ices) is globally conserved during the collapsing process, from the parent molecular cloud to the young disk surrounding the first Larson core. A qualitative comparison with cometary composition suggests that comets are constituted of different phases, some molecules being direct tracers of interstellar chemistry, while others, including complex molecules, seem to have been formed in disks, where higher densities and temperatures allow for an active grain surface chemistry. The latter phase, and its connection with the formation of the first Larson core, remains to be modeled.

79 ASTRONOMY AND ASTROPHYSICS↗

On the relaxation behaviors of slow and classical glitches: Observational biases and their opposite recovery trends

We study the pulsar timing properties and the data analysis methods of glitch recoveries. In some cases one first fits the times of arrival (TOAs) to obtain the 'time-averaged' frequency ν and its first derivative ν-dot , and then fits models to them. However, our simulations show that ν and ν-dot obtained in this way are systematically biased unless the time intervals between the nearby data points of TOAs are smaller than about 10{sup 4} s, which is much shorter than typical observation intervals. Alternatively, glitch parameters can be obtained by fitting the phases directly with relatively smaller biases, but the initial recovery time scale is usually chosen by eye, which may introduce a strong bias. We also construct a phenomenological model by assuming a pulsar spin-down law of ν-dot ν{sup −3}=−H{sub 0}G(t) with G(t) = 1 + κe {sup –t/τ} for a glitch recovery, where H {sub 0} is a constant and κ and τ are the glitch parameters to be found. This model can reproduce the observed data of slow glitches from B1822–09 and the giant classical glitch of B2334+61, with κ < 0 or κ > 0, respectively. We then use this model to simulate TOA data and test several fitting procedures for a glitch recovery. The best procedure is (1) to use a very high order polynomial (e.g., to 50th order) to precisely describe the phase, (2) obtain ν(t) and ν-dot (t) from the polynomial, then (3) obtain the glitch parameters from ν(t) or ν-dot (t). Finally, the uncertainty in the starting time t {sub 0} of a classical glitch causes uncertainties in some glitch parameters, but less so for a slow glitch, t {sub 0}, of which can be determined from data.

79 ASTRONOMY AND ASTROPHYSICS↗

Physical conditions in the reconnection layer in pulsar magnetospheres

The magnetosphere of a rotating pulsar naturally develops a current sheet (CS) beyond the light cylinder (LC). Magnetic reconnection in this CS inevitably dissipates a nontrivial fraction of the pulsar spin-down power within a few LC radii. We develop a basic physical picture of reconnection in this environment and discuss its implications for the observed pulsed gamma-ray emission. We argue that reconnection proceeds in the plasmoid-dominated regime, via a hierarchical chain of multiple secondary islands/flux ropes. The inter-plasmoid reconnection layers are subject to strong synchrotron cooling, leading to significant plasma compression. Using the conditions of pressure balance across these current layers, the balance between the heating by magnetic energy dissipation and synchrotron cooling, and Ampere's law, we obtain simple estimates for key parameters of the layers—temperature, density, and layer thickness. In the comoving frame of the relativistic pulsar wind just outside of the equatorial CS, these basic parameters are uniquely determined by the strength of the reconnecting upstream magnetic field. For the case of the Crab pulsar, we find them to be of order 10 GeV, 10{sup 13} cm{sup –3}, and 10 cm, respectively. After accounting for the bulk Doppler boosting due to the pulsar wind, the synchrotron and inverse-Compton emission from the reconnecting CS can explain the observed pulsed high-energy (GeV) and very high energy (∼100 GeV) radiation, respectively. Also, we suggest that the rapid relative motions of the secondary plasmoids in the hierarchical chain may contribute to the production of the pulsar radio emission.

79 ASTRONOMY AND ASTROPHYSICS↗

Ultra high-energy neutrinos via heavy-meson synchrotron emission in strong magnetic fields

We explore the generation and possibility for the detection of heavy-meson synchrotron emission due to the acceleration of ultra-relativistic protons (and possibly nuclei) in the presence of strong magnetic fields (H ≳ 10{sup 15} G) in transient astrophysical environments such as magnetar flares. We show that, in addition to the well-known pion synchrotron emission, heavy vector mesons like ρ, D{sub S} , J/Ψ, and Y could be generated. For high enough energies and magnetic field strengths, such heavy vector mesons can be formed with high intensity (∼10{sup 3} times the photon intensity) through strong couplings to the ultra-relativistic nucleons. We examine in particular the synchrotron emission and subsequent cooling and decay of the heavy ρ{sup 0} and Y(1S) mesons, e.g., via p → p' + Y(1S), Y(1S) → τ{sup +} + τ{sup –}, τ{sup −}→μ{sup −}+ ν-bar {sub μ}+ν{sub τ} and e{sup −}+ ν-bar {sub e}+ν{sub τ}. We evaluate the spectra of escaping ν {sub e}, ν{sub μ}, and ν{sub τ} due to the decay of short-lived τ mesons. We deduce the possible event rate in a terrestrial TeV neutrino detector. We estimate that neutrinos produced from the heavy vector-meson synchrotron radiation from a strong magnetar soft gamma repeater burst will only be detectable with the current generation of detectors if the source is very nearby (<30 pc). Nevertheless, if ever detected, the existence of heavy meson synchrotron emission might be identifiable by the unique signature of energetic tau neutrinos emanating from the source.

79 ASTRONOMY AND ASTROPHYSICS↗

Fragmentation of massive dense cores down to ≲ 1000 AU: Relation between fragmentation and density structure

In order to shed light on the main physical processes controlling fragmentation of massive dense cores, we present a uniform study of the density structure of 19 massive dense cores, selected to be at similar evolutionary stages, for which their relative fragmentation level was assessed in a previous work. We inferred the density structure of the 19 cores through a simultaneous fit of the radial intensity profiles at 450 and 850 μm (or 1.2 mm in two cases) and the spectral energy distribution, assuming spherical symmetry and that the density and temperature of the cores decrease with radius following power-laws. Even though the estimated fragmentation level is strictly speaking a lower limit, its relative value is significant and several trends could be explored with our data. We find a weak (inverse) trend of fragmentation level and density power-law index, with steeper density profiles tending to show lower fragmentation, and vice versa. In addition, we find a trend of fragmentation increasing with density within a given radius, which arises from a combination of flat density profile and high central density and is consistent with Jeans fragmentation. We considered the effects of rotational-to-gravitational energy ratio, non-thermal velocity dispersion, and turbulence mode on the density structure of the cores, and found that compressive turbulence seems to yield higher central densities. Finally, a possible explanation for the origin of cores with concentrated density profiles, which are the cores showing no fragmentation, could be related with a strong magnetic field, consistent with the outcome of radiation magnetohydrodynamic simulations.

79 ASTRONOMY AND ASTROPHYSICS↗