B-field Strength Research Articles

Probing the Dust Grain Alignment Mechanisms in Spiral Galaxies with M51 as the Case Study

Magnetic fields (B fields) in galaxies have recently been traced using far-infrared and submillimeter polarimetric observations with the Stratospheric Observatory for Infrared Astronomy, James Clerk Maxwell Telescope, and Atacama Large Millimeter/submillimeter Array. The main assumption is that dust grains are magnetically aligned with the local B-field in the interstellar medium (ISM). However, the range of physical conditions of the ISM, dust grain sizes, and B-field strengths in galaxies where this assumption is valid has not been characterized yet. Here, we use the well-studied spiral galaxy M51 as a case study. We find that the timescale for the alignment mechanism arising from magnetically aligned dust grains (B-RAT) dominates over other alignment mechanisms, including radiative precession (k-RAT) and mechanical alignment (v-MAT), as well as the randomization effect (gas damping). We estimate the sizes of the aligned dust grain to be in the range of 0.009–0.182 μm and 0.019–0.452 μm for arms and interarms, respectively. We show that the difference in the polarization fraction between arms and interarms may arise from the enhancement of small dust grain sizes in the arms resulting from large grains being broken into small grains as an effect of the grain alignment disruption (RAT-D) mechanism. We argue that the RAT-D mechanism needs to have additional effects, e.g., intrinsic variations of the B-field structure and turbulence, in the galaxy’s components to fully explain the polarization fraction variations within the arms and interarms.

Mapping and characterizing magnetic fields in the Rho Ophiuchus-A molecular cloud with SOFIA/HAWC+

Context. Together with gravity, turbulence, and stellar feedback, magnetic fields (B-fields) are thought to play a critical role in the evolution of molecular clouds and star formation processes. The polarization of thermal dust emission is a popular tracer of B-fields in star-forming regions. Aims. We aim to map the morphology and measure the strength of B-fields of the nearby molecular cloud, rho Ophiuchus-A (ρ Oph-A), to understand the role of B-fields in regulating star formation and in shaping the cloud. Methods. We analyzed the far-infrared (FIR) polarization of thermal dust emission observed by SOFIA/HAWC+ at 89 and 154 μm toward the densest part of ρ Oph-A, which is irradiated by the nearby B3/4 star, Oph-S1. These FIR polarimetric maps cover an area of ~4.5′ × 4.5′ (corresponding to 0″.18 × 0″.18 pc2) with an angular resolution of 7.8″ and 13.6″ respectively. Results. The ρ Oph-A cloud exhibits well-ordered B-fields with magnetic orientations that are mainly perpendicular to the ridge of the cloud toward the densest region. We obtained a map of B-field strengths in the range of 0.2–2.5 mG, using the Davis-Chandrasekhar-Fermi (DCF) method. The B-fields are strongest at the densest part of the cloud, which is associated with the starless core SM1, and then decrease toward the outskirts of the cloud. By calculating the map of the mass-to-flux ratio, Alfvén Mach number, and plasma β parameter in ρ Oph-A, we find that the cloud is predominantly magnetically sub-critical, sub-Alfvénic, which implies that the cloud is supported by strong B-fields that dominate over gravity, turbulence, and thermal gas energy. The measured B-field strengths at the two densest subsregions using other methods that account for the compressible mode are relatively lower than that measured with the DCF method. However, these results do not significantly change our conclusions on the roles of B-fields relative to gravity and turbulence on star formation. Our virial analysis suggests that the cloud is gravitationally unbound, which is consistent with the previous detection of numerous starless cores in the cloud. By comparing the magnetic pressure with the radiation pressure from the Oph-S1 star, we find that B-fields are sufficiently strong to support the cloud against radiative feedback and to regulate the shape of the cloud.

Massive Star Formation Starts in Subvirial Dense Clumps Unless Resisted by Strong Magnetic Fields

Knowledge of the initial conditions of high-mass star formation is critical for theoretical models, but are not well observed. Built on our previous characterization of a Galaxy-wide sample of 463 candidate high-mass starless clumps (HMSCs), here we investigate the dynamical state of a representative subsample of 44 HMSCs (radii 0.13–1.12 pc) using Green Bank Telescope NH3 (1,1) and (2,2) data from the Radio Ammonia Mid-Plane Survey pilot data release. By fitting the two NH3 lines simultaneously, we obtain velocity dispersion, gas kinetic temperature, NH3 column density and abundance, Mach number, and virial parameter. Thermodynamic analysis reveals that most HMSCs have Mach number <5, inconsistent with what have been considered in theoretical models. All but one (43 out of 44) of the HMSCs are gravitationally bound with virial parameter α vir < 2. Either these massive clumps are collapsing or magnetic field strengths of 0.10–2.65 mG (average 0.51 mG) would be needed to prevent them from collapsing. The estimated B-field strength correlates tightly with density, Best/mG=0.269(nH2/104cm−3)0.61 , with a similar power-law index as found in observations but a factor of 4.6 higher in strength. For the first time, the initial dynamical state of high-mass formation regions has been statistically constrained to be subvirial, in contradiction to theoretical models in virial equilibrium and in agreement with the lack of observed massive starless cores. The findings urge future observations to quantify the magnetic field support in the prestellar stage of massive clumps, which has rarely been explored so far, toward a full understanding of the physical conditions that initiate massive star formation.

Magnetic field at the Galactic centre from multiwavelength dust polarization

ABSTRACT We have mapped the magnetic field (B-field) for a region of about 30 pc around the centre of our Galaxy, which encompasses the circumnuclear disc (CND), the minispiral, and the 20 and 50 km s−1 molecular clouds, using thermal dust polarization observations obtained from SOFIA/HAWC+ and JCMT/SCUPOL. We decompose the spectra of 12CO (J = 3 → 2) transition from this region into individual cloud components and find the polarization observed at different wavelengths might be tracing completely different layers of dust along the line of sight. We use modified Davis–Chandrasekhar–Fermi methods to measure the strength of B-field projected in the plane of the sky ($B_{{}_{\mathrm{POS}}}$). The mean $B_{{}_{\mathrm{POS}}}$ of the CND and the minispiral, probed at 53 μm is of the order of ∼2 mG. $B_{{}_{\mathrm{POS}}}\!\!\!\lt \!1$ mG close to the Galactic Centre, in the region of the ionized mini-cavity within the CND, and increases outwards. However, the longer wavelength polarization at 216 μm appears to come from a dust layer that is cooler and behind the CND and has a stronger B-field of about 7 mG. The B-field strength is lowest along the Eastern Arm of the minispiral, which is also the only region with Alfvén Mach number, $\mathcal {M}_{\mathrm{A}}\gt 1$ and mass-to-flux ratio, λ ≳ 1. Such an observed weak B-field could be a result of the low resolution of the observation, where the tangled B-fields due to the strong turbulence in the high density clumps of the CND are lost within the beam size of the observation.

Polarized Light from Massive Protoclusters (POLIMAP). I. Dissecting the Role of Magnetic Fields in the Massive Infrared Dark Cloud G28.37+0.07

Magnetic fields may play a crucial role in setting the initial conditions of massive star and star cluster formation. To investigate this, we report SOFIA-HAWC+ 214 μm observations of polarized thermal dust emission and high-resolution GBT-Argus C18O(1-0) observations toward the massive Infrared Dark Cloud (IRDC) G28.37+0.07. Considering the local dispersion of B-field orientations, we produce a map of the B-field strength of the IRDC, which exhibits values between ∼0.03 and 1 mG based on a refined Davis–Chandrasekhar–Fermi method proposed by Skalidis & Tassis. Comparing to a map of inferred density, the IRDC exhibits a B–n relation with a power-law index of 0.51 ± 0.02, which is consistent with a scenario of magnetically regulated anisotropic collapse. Consideration of the mass-to-flux ratio map indicates that magnetic fields are dynamically important in most regions of the IRDC. A virial analysis of a sample of massive, dense cores in the IRDC, including evaluation of magnetic and kinetic internal and surface terms, indicates consistency with virial equilibrium, sub-Alfvénic conditions, and a dominant role for B-fields in regulating collapse. A clear alignment of magnetic field morphology with the direction of the steepest column density gradient is also detected. However, there is no preferred orientation of protostellar outflow directions with the B-field. Overall, these results indicate that magnetic fields play a crucial role in regulating massive star and star cluster formation, and therefore they need to be accounted for in theoretical models of these processes.

A Study of the Radio Spectrum of Mrk 421

We present the results of a spectral analysis using simultaneous multifrequency (22, 43, 86, and 129 GHz) very long baseline interferometry (VLBI) observations of the Korean VLBI Network on BL Lac object, Markarian 421. The data we used were obtained from 2013 January to 2018 June. The light curves showed several flux enhancements with global decreases. To separate the variable and quiescent components in the multifrequency light curves for milliarcsecond-scale emission regions, we assumed that the quiescent radiation comes from the emission regions radiating constant optically thin synchrotron emissions (i.e., a minimum flux density with an optically thin spectral index). The quiescent spectrum determined from the multifrequency light curves was subtracted from the total CLEAN flux density, yielding a variable component in the flux that produces the time-dependent spectrum. We found that the observed spectra were flat at 22–43 GHz, and relatively steep at 43–86 GHz, whereas the quiescent-corrected spectra are sometimes quite different from the observed spectra (e.g., sometimes inverted at 22–43 GHz). The quiescent-corrected spectral indices were much more variable than the observed spectral indices. This spectral investigation implies that the quiescent-spectrum correction can significantly affect the multifrequency spectral index of variable compact radio sources such as blazars. Therefore, the synchrotron self-absorption B-field strength (B SSA) can be significantly affected because B SSA is proportional to the fifth power of turnover frequency.

Effects of grain magnetic properties and grain growth on synthetic dust polarization of MHD simulations of low-mass Class 0/I YSOs

ABSTRACT Thermal dust polarization is a powerful tool to probe magnetic fields ($\boldsymbol{B}$) and grain properties. However, a systematic study of the dependence of dust polarization on grain properties in protostellar environments is not yet available. In this paper, we post-process a non-ideal MHD simulation of a collapsing protostellar core with our updated POLARIS code to study in detail the effects of iron inclusions and grain growth on thermal dust polarization. We found that superparamagnetic (SPM) grains can produce high polarization degree of $p \sim 10\!-\!40~{{\ \rm per\ cent}}$ beyond ∼500 au from the protostar because of their efficient alignment by magnetically enhanced radiative torque mechanism. The magnetic field turbulence in the envelope causes the decrease in p with increasing emission intensity I as p ∝ Iα with the slope α ∼ −0.3. But within 500 au, SPM grains tend to have inefficient internal alignment and be aligned with $\boldsymbol{B}$ by RATs only, producing lower $p \sim 1~{{\ \rm per\ cent}}$ and a steeper slope of α ∼ −0.6. For paramagnetic (PM) grains, the alignment loss of grains above $1\, {\mu \rm {m}}$ in the inner ∼200 au produces $p \lt \lt 1~{{\ \rm per\ cent}}$ and the polarization hole with α ∼ −0.9. Grain growth can increase p in the envelope for SPM grains, but cause stronger depolarization for SPM grains in the inner ∼500 au and for PM grains in the entire protostellar core. Finally, we found the increase of polarization angle dispersion function S with iron inclusions and grain growth, implying the dependence of B-field strength measured using the David–Chandrasekhar–Fermi technique on grain alignment and grain properties.

Prospects for detecting proto-neutron star rotation and spin-down using supernova neutrinos

ABSTRACT After a successful supernova, a proto-neutron star (PNS) cools by emitting neutrinos on ∼1–100 s time-scales. Provided that there are neutrino emission ‘hotspots’ or ‘cold-spots’ on the surface of the rotating PNS, we can expect a periodic modulation in the number of neutrinos observable by detectors. We show that Fourier transform techniques can be used to determine the PNS rotation rate from the neutrino arrival times. Provided there is no spin-down, a 1-parameter Discrete Fourier Transform (DFT) is sufficient to determine the spin period of the PNS. If the PNS is born as a magnetar with polar magnetic field strength B0 ≳ 1015 G and is ‘slowly’ rotating with an initial spin period ≳100 ms, then it can spin-down to periods of the order of seconds during the cooling phase. We propose a modified DFT technique with three frequency parameters to detect spin-down. Due to lack of neutrino data from a nearby supernova except the ∼20 neutrinos detected from SN1987A, we use toy models and one physically motivated modulating function to generate neutrino arrival times. We use the false alarm rate (FAR) to quantify the significance of the Fourier power spectrum peaks. We show that PNS rotation and spin-down are detected with $\rm FAR\,\lt\, 2~{{\ \rm per\ cent}}$ (2σ) for periodic signal content $\rm M\gtrsim 13-15~{{\ \rm per\ cent}}$ if 5 × 103 neutrinos are detected in ∼3 s and with $\rm FAR\,\lt\, 1{{\ \rm per\ cent}}$ for $\rm M\,\ge 5{{\ \rm per\ cent}}$ if 5 × 104 neutrinos are detected in ∼3 s. Since we can expect ∼104−105 neutrino detections from a supernova at 10 kpc, detection of PNS rotation and spin-down is possible using the neutrinos from the next Galactic supernova.

Comptonization by reconnection plasmoids in black hole coronae – III. Dependence on the guide field in pair plasma

ABSTRACT We perform non-radiative two-dimensional particle-in-cell simulations of magnetic reconnection for various strengths of the guide field (perpendicular to the reversing field), in magnetically dominated electron–positron plasmas. Magnetic reconnection under such conditions could operate in accretion disc coronae around black holes. There, it has been suggested that the transrelativistic bulk motions of reconnection plasmoids containing inverse-Compton-cooled electrons could Compton-upscatter soft photons to produce the observed non-thermal hard X-rays. Our simulations are performed for magnetizations 3 ≤ σ ≤ 40 (defined as the ratio of enthalpy density of the reversing field to plasma enthalpy density) and guide field strengths 0 ≤ Bg/B0 ≤ 1 (normalized to the reversing field strength B0). We find that the mean bulk energy of the reconnected plasma depends only weakly on the flow magnetization but strongly on the guide field strength – with Bg/B0 = 1 yielding a mean bulk energy twice smaller than Bg/B0 = 0. Similarly, the dispersion of bulk motions around the mean – a signature of stochasticity in the plasmoid chain’s motions – is weakly dependent on magnetization (for σ ≳ 10) but strongly dependent on the guide field strength – dropping by more than a factor of two from Bg/B0 = 0 to Bg/B0 = 1. In short, reconnection in strong guide fields (Bg/B0 ∼ 1) leads to slower and more ordered plasmoid bulk motions than its weak guide field (Bg/B0 ∼ 0) counterpart.

The early evolution of magnetar rotation – II. Rapidly rotating magnetars: implications for gamma-ray bursts and superluminous supernovae

ABSTRACT Rapidly rotating magnetars have been associated with gamma-ray bursts (GRBs) and superluminous supernovae (SLSNe). Using a suite of two-dimensional magnetohydrodynamic simulations at fixed neutrino luminosity and a couple of evolutionary models with evolving neutrino luminosity and magnetar spin period, we show that magnetars are viable central engines for powering GRBs and SLSNe. We also present analytical estimates of the energy outflow rate from the proto-neutron star (PNS) as a function of polar magnetic field strength B0, PNS angular velocity Ω⋆, PNS radius R⋆, and mass outflow rate $\dot{M}$. We show that rapidly rotating magnetars with spin periods P⋆ ≲ 4 ms and polar magnetic field strength B0 ≳ 1015 G can release 1050 to 5 × 1051 erg of energy during the first ∼2 s of the cooling phase. Based on this result, it is plausible that sustained energy injection by magnetars through the relativistic wind phase can power GRBs. We also show that magnetars with moderate field strengths of B0 ≲ 5 × 1014 G do not release a large fraction of their rotational kinetic energy during the cooling phase and, hence, are not likely to power GRBs. Although we cannot simulate to times greater than ∼3–5 s after a supernova, we can hypothesize that moderate field strength magnetars can brighten the supernova light curves by releasing their rotational kinetic energy via magnetic dipole radiation on time-scales of days to weeks, since these do not expend most of their rotational kinetic energy during the early cooling phase.

B-fields and Dust in Interstellar Filaments Using Dust Polarization (BALLAD-POL). I. The Massive Filament G11.11–0.12 Observed by SOFIA/HAWC+

We report the first measurement of polarized thermal dust emission toward the entire infrared dark cloud G11.11−0.12 taken by the polarimeter SOFIA/HAWC+ at 214 μm. The obtained magnetic fields (B-fields) from the polarized emission of the early-stage and massive filament tend to be perpendicular to its spine. We produce a map of B-field strengths for the center region of the filament. The strengths vary in the range of 100–600 μG and are strongest along the filament's spine. The central region is sub-Alfvénic and mostly subcritical, meaning that B-fields dominate over turbulence and are strong enough to resist gravitational collapse. The alignment and properties of dust grains in the filament are studied using radiative torque (RAT) theory. We find the decrease of polarization degree P with emission intensity I, i.e., depolarization effect, of the form P ∝ I −α with α ∼ 0.8–0.9, implying a significant loss of grain alignment in the filament's spine. The depolarization can be explained by the decrease in RAT alignment efficiency toward the denser regions with weaker radiation field, which cannot be explained by B-field tangling. We study the effect of the enhanced magnetic relaxation by embedded iron inclusions on RAT alignment and find that the high polarization fraction P ∼ 20%–30% in the outer layer of the filament is potential evidence for the magnetically enhanced RAT alignment mechanism. This is the first time this effect is evaluated in a filament. Based on the polarization fraction and RAT alignment theory, we also find evidence for grain growth in the filament.

The Strength of the Sheared Magnetic Field in the Galactic’s Circumnuclear Disk

Recent high-resolution 53 μm polarimetric observations from SOFIA/HAWC+ have revealed the inferred plane-of-the-sky magnetic field (B-field) orientation in the Galactic center’s circumnuclear disk (CND). The B-field is mostly aligned with the steamers of ionized material falling onto Sgr A* at large, differential velocities (shear). In such conditions, estimating the B-field strength with the “classical” Davis–Chandrasekhar–Fermi (DCF) method does not provide accurate results. We derive a “modified” DCF method by solving the ideal-MHD equations from first principles considering the effects of a large-scale, shear flow on the propagation of a fast magnetosonic wave. In the context of the DCF approximation, both the value of the shear and its Laplacian affect the inferred B-field strength. Using synthetic polarization data from MHD simulations for a medium dominated by shear flows, we find that the “classical” DCF determines B-field strengths only within >50% of the true value where the “modified” DCF results are improved significantly (∼3%–22%). Applying our “modified” DCF method to the CND revealed B-field strengths of 1–16 mG in the northern arm, 1–13 mG in the eastern arm, and 3–27 mG in the western arc at spatial scales ≲1 pc, with median values of 5.1 ± 0.8, 4.0 ± 1.2, and 8.5 ± 2.3 mG, respectively. The balance between turbulent gas energy (kinetic plus hydrostatic) and turbulent magnetic energy densities suggest that, along the magnetic-field-flow direction, magnetic effects become less dominant as the shear flow increases and weakens the B-field via magnetic convection. Our results indicate that the transition from magnetically to gravitationally dominated accretion of material onto Sgr A* starts at distances ∼1 pc.

Double SSA spectrum and magnetic field strength of the FSRQ 3C 454.3

ABSTRACT We present the results of a radio multifrequency ($\rm 3{-}340~GHz$) study of the blazar 3C 454.3. After subtracting the quiescent spectrum corresponding to optically thin emission, we found two individual synchrotron self-absorption (SSA) features in the wide-band spectrum. The one SSA had a relatively low turnover frequency (νm) in the range of $\rm 3{-}37~GHz$ (lower νm SSA spectrum, LSS), and the other one had a relatively high νm of $\rm 55{-}124~GHz$ (higher νm SSA spectrum, HSS). Using the SSA parameters, we estimated B-field strengths at the surface where optical depth τ = 1. The estimated B-field strengths were $\rm \gt 7$ and $\rm \gt 0.2~mG$ for the LSS and HSS, respectively. The LSS-emitting region was magnetically dominated before the 2014 June γ-ray flare. The quasi-stationary component (C), ∼0.6 mas apart from the 43 -GHz radio core, became brighter than the core with decreasing observing frequency, and we found that component C was related to the LSS. A decrease in jet width was found near component C. As a moving component, K14 approached component C, and the flux density of the component was enhanced while the angular size decreased. The high intrinsic brightness temperature in the fluid frame was obtained as TB, int ≈ (7.0 ± 1.0) × 1011 K from the jet component after the 2015 August γ-ray flare, suggesting that component C is a high-energy emitting region. The observed local minimum of jet width and re-brightening behaviour suggest a possible recollimation shock in component C.

SOFIA Observations of 30 Doradus. II. Magnetic Fields and Large-scale Gas Kinematics

The heart of the Large Magellanic Cloud, 30 Doradus, is a complex region with a clear core-halo structure. Feedback from the stellar cluster R136 has been shown to be the main source of energy creating multiple parsec-scale expanding-shells in the outer region, and carving a nebula core in the proximity of the ionization source. We present the morphology and strength of the magnetic fields (B-fields) of 30 Doradus inferred from the far-infrared polarimetric observations by SOFIA/HAWC+ at 89, 154, and 214 μm. The B-field morphology is complex, showing bending structures around R136. In addition, we use high spectral and angular resolution [C ii] observations from SOFIA/GREAT and CO(2-1) from APEX. The kinematic structure of the region correlates with the B-field morphology and shows evidence of multiple expanding-shells. Our B-field strength maps, estimated using the Davis–Chandrasekhar–Fermi method and structure-function, show variations across the cloud within a maximum of 600, 450, and 350 μG at 89, 154, and 214 μm, respectively. We estimated that the majority of the 30 Doradus clouds are subcritical and sub-Alfvénic. The probability distribution function of the gas density shows that the turbulence is mainly compressively driven, while the plasma beta parameter indicates supersonic turbulence. We show that the B-field is sufficient to hold the cloud structure integrity under feedback from R136. We suggest that supersonic compressive turbulence enables the local gravitational collapse and triggers a new generation of stars to form. The velocity gradient technique using [C ii] and CO(2-1) is likely to confirm these suggestions.

Radio-loud Exoplanet-exomoon Survey: GMRT Search for Electron Cyclotron Maser Emission

We conducted the first dedicated search for signatures of exoplanet–exomoon interactions using the Giant Metrewave Radio Telescope (GMRT) as part of the radio-loud exoplanet-exomoon survey. Due to stellar tidal heating, irradiation, and subsequent atmospheric escape, candidate “exo-Io” systems are expected to emit up to 106 times more plasma flux than the Jupiter-Io DC circuit. This can induce detectable radio emission from the exoplanet-exomoon system. We analyze three “exo-Io” candidate stars: WASP-49, HAT-P 12, and HD 189733. We perform 12 hr phase-curve observations of WASP-49b at 400 MHz during primary & secondary transit, as well as first & third quadratures achieving a 3σ upper limit of 0.18 mJy beam−1 averaged over four days. HAT-P 12 was observed with GMRT at 150 and 325 MHz. We further analyzed the archival data of HD 189733 at 325 MHz. No emission was detected from the three systems. However, we place strong upper limits on radio flux density. Given that most exo-Io candidates orbit hot Saturns, we encourage more multiwavelength searches (in particular low frequencies) to span the lower range of exoplanet B-field strengths constrained here.

Experimental evidence of early-time saturation of the ion-Weibel instability in counterstreaming plasmas of CH, Al, and Cu.

The collisionless ion-Weibel instability is a leading candidate mechanism for the formation of collisionless shocks in many astrophysical systems, where the typical distance between particle collisions is much larger than the system size. Multiple laboratory experiments aimed at studying this process utilize laser-driven (I≳10^{15} W/cm^{2}), counterstreaming plasma flows (V≲2000km/s) to create conditions unstable to Weibel-filamentation and growth. This technique intrinsically produces temporally varying plasma conditions at the midplane of the interaction where Weibel-driven B fields are generated and studied. Experiments discussed herein demonstrate robust formation of Weibel-driven B fields under multiple plasma conditions using CH, Al, and Cu plasmas. Linear theory based on benchmarked radiation-hydrodynamic FLASH calculations is compared with Fourier analyses of proton images taken ∼5-6 linear growth times into the evolution. The new analyses presented here indicate that the low-density, high-velocity plasma-conditions present during the first linear-growth time (∼300-500ps) sets the spectral characteristics of Weibel filaments during the entire evolution. It is shown that the dominant wavelength (∼300μm) at saturation persists well into the nonlinear phase, consistent with theory under these experimental conditions. However, estimates of B-field strength, while difficult to determine accurately due to the path-integrated nature of proton imaging, are shown to be in the ∼10-30T range, an order of magnitude above the expected saturation limit in homogenous plamas but consistent with enhanced B fields in the midplane due to temporally varying plasma conditions in experiments.

Wannier Diagrams for Semiconductor Artificial Graphene

Quantum transport has been simulated in hexagonal semiconductor lattices of antidots with a period of 80 nm and short-range disorder. Wannier diagrams, i.e., DoS(n, B) maps of the density of states, where n is the electron density and B is the magnetic field strength, have been calculated for several potential modulation amplitudes comparable to or much larger than the Fermi energy. Deep dips in the maps of the density of states have the form of rays with positive, zero, and negative slopes. In addition to the fan of the rays separating the first and second, as well as the second and third Landau levels, the maps include rays that are parallel to them and are shifted in n and B by integers of the characteristic electron density n0 and the characteristic magnetic field strength B0, respectively. It has been shown that the sign and magnitude of the slope of the rays in the density of states correspond to the centers of the plateaus of quantized Hall resistances Rxy. The lattice is brightly manifested in the Rxy(n, B) maps as the replicas of the first and second plateaus in Rxy and as oscillations of Rxy between negative and positive values at a fixed magnetic field or a fixed electron density, which indicates the interchange between the hole and electron charge carriers.

The impact of ion chamber components on k B,Q for reference dosimetry in MRgRT

For reference dosimetry in MRgRT, k B,Q is used to correct for the impact of the magnetic field on the chamber calibration coefficient. It has been demonstrated that for accurate simulation of k B,Q the dead volume (DV) must be considered. This work goes one step further by analysing the contribution of secondary electrons generated in the various chamber components to k B,Q . The Farmer-type chamber PTW 30013 geometry was modelled for two different DVs. Monte Carlo simulations were performed for a 60Co source and a 7 MV MRI-linac and the model was validated against measurements. Both parallel (α = 0° or 180°) and perpendicular (α = 90° or 270°) orientations of the chamber and the magnetic (B) field were considered, and several B-field strengths between 0 T and 1.5 T. To study the dose contribution to the reduced volume (RV = cavity — DV) from the secondary electrons produced in certain components of the chamber the labelling of the particles was implemented in the PENELOPE user code PENMAIN. A separate model with each solid component of the chamber modelled as liquid water was used to investigate the impact of material choice on k B,Q . Results show that simulated k B,Q values agree better with the measured k B,Q when the DV is considered. It is demonstrated that small components of the chamber impact k B,Q considerably, since the contribution to the RV-dose from the bodies closer to the RV is higher than without B. Moreover, it is seen that the impact to the dose in the RV is reduced when the material of each component is modelled as liquid water. Therefore, chamber design and, to a lesser extent, choice of material affect k B,Q , and an accurate geometrical model of the chamber components and its further validation are important for correct calculations of k B,Q .

Studying Magnetic Fields and Dust in M17 Using Polarized Thermal Dust Emission Observed by SOFIA/HAWC+

Abstract We report on the highest spatial resolution measurement to date of magnetic fields (B-fields) in M17 using thermal dust polarization measurements taken by SOFIA/HAWC+ centered at a wavelength of 154 μm. Using the Davis–Chandrasekhar–Fermi method, in which the polarization angle dispersion calculated using the structure function technique is the quantity directly observed by SOFIA/HAWC+, we found the presence of strong B-fields of 980 ± 230 and 1665 ± 885 μG in the lower-density M17-N and higher-density M17-S regions, respectively. The B-field morphology in M17-N possibly mimics the fields in gravitationally collapsing molecular cores, while in M17-S the fields run perpendicular to the density structure. M17-S also displays a pillar feature and an asymmetric large-scale hourglass-shaped field. We use the mean B-field strengths to determine Alfvénic Mach numbers for both regions, finding that B-fields dominate over turbulence. We calculate the mass-to-flux ratio, λ, finding λ = 0.07 for M17-N and 0.28 for M17-S. These subcritical λ values are consistent with the lack of massive stars formed in M17. To study dust physics, we analyze the relationship between dust polarization fraction, p, emission intensity, I, gas column density, N(H2), polarization angle dispersion function, S, and dust temperature, T d. p decreases with intensity as I −α with α = 0.51. p tends to first increase with T d, but then decreases at higher T d. The latter feature, seen in M17-N at high T d when N(H2) and S decrease, is evidence of the radiative torque disruption effect.

Development and application of OpenFOAM based magnetohydrodynamic solver

We develop a compressible magnetohydrodynamic solver to simulate the transonic flows based on an open-source computational fluid dynamics platform OpenFOAM. The solver is achieved by modifying the density-based Riemann solver &lt;i&gt;rhoCentralFoam&lt;/i&gt; which adopts a central scheme and is available in OpenFOAM. To improve simulation accuracy and avoid non-physical oscillations, a specialized pressure-implicit algorithm with the splitting of operators is implemented to guarantee the incompressibility of magnetic field. The solver is benchmarked and the convergence rate is between the first and the second order. After benchmark, we apply this solver to magnetohydrodynamic simulations of intense-laser-produced plasma. The influences of uniform axial magnetic field and nonuniform coil-current-induced magnetic field on laser-produced plasma jets are investigated. With the uniform axial magnetic field, the positions of nozzle and the distance between knots are linearly related to square root of thermal over magnetic pressure. With the nonuniform magnetic field generated in the coil, knots are nonlinearly distributed in space and the nozzle position is modulated according to preliminary simulations. In the two kinds of magnetic fields, when the B-field strength is the same at coil center, the magnetic field of relatively small coils can shorten the times of forming nozzles and knots, suggesting that the coil magnetic field is equivalent to a higher uniform one. The simulations can be used as a reference for our future experiment on magnetized laser-produced plasma jet. Meanwhile, our simulation investigation shows that this magnetohydrodynamic solver is suitable for engineering calculation for laser plasma experiments and can deal with the situation with relatively complex configurations.