Showing 1–50 of 94 results for author: Qiu, J

Kinetic temperature of massive star-forming molecular clumps measured with formaldehyde V. The massive filament DR21

Authors: X. Zhao, X. D. Tang, C. Henkel, Y. Gong, Y. Lin, D. L. Li, Y. X. He, Y. P. Ao, X. Lu, T. Liu, Y. Sun, K. Wang, X. P. Chen, J. Esimbek, J. J. Zhou, J. W. Wu, J. J. Qiu, X. W. Zheng, J. S. Li, C. S. Luo, Q. Zhao

Abstract: The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense g… ▽ More The kinetic temperature structure of the massive filament DR21 has been mapped using the IRAM 30 m telescope. This mapping employed the para-H$_2$CO triplet ($J_{\rm K_aK_c}$ = 3$_{03}$--2$_{02}$, 3$_{22}$--2$_{21}$, and 3$_{21}$--2$_{20}$) on a scale of $\sim$0.1 pc. By modeling the averaged line ratios of para-H$_{2}$CO with RADEX under non-LTE assumptions, the kinetic temperature of the dense gas was derived at a density of $n$(H$_{2}$) = 10$^{5}$ cm$^{-3}$. The para-H$_2$CO lines reveal significantly higher temperatures than NH$_3$ (1,1)/(2,2) and FIR wavelengths. The dense clumps appear to correlate with the notable kinetic temperature. Among the four dense cores (N44, N46, N48, and N54), temperature gradients are observed on a scale of $\sim$0.1-0.3 pc. This suggests that the warm dense gas is influenced by internal star formation activity. With the exception of N54, the temperature profiles of these cores were fitted with power-law indices ranging from $-$0.3 to $-$0.5. This indicates that the warm dense gas is heated by radiation emitted from internally embedded protostar(s) and/or clusters. While there is no direct evidence supporting the idea that the dense gas is heated by shocks resulting from a past explosive event in the DR21 region, our measurements toward the DR21W1 region provide compelling evidence that the dense gas is indeed heated by shocks originating from the western DR21 flow. Higher temperatures appear to be associated with turbulence. The physical parameters of the dense gas in the DR21 filament exhibit a remarkable similarity to the results obtained in OMC-1 and N113. This may imply that the physical mechanisms governing the dynamics and thermodynamics of dense gas traced by H$_{2}$CO in diverse star formation regions may be dominated by common underlying principles despite variations in specific environmental conditions. (abbreviated) △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 16 pages, 8 figures, 3 tabels. Accepted for publication by Astronomy & Astrophysics

Journal ref: A&A 687, A207 (2024)

Provably Convergent and Robust Newton-Raphson Method: A New Dawn in Primitive Variable Recovery for Relativistic MHD

Authors: Chaoyi Cai, Jianxian Qiu, Kailiang Wu

Abstract: A long-standing and formidable challenge faced by all conservative schemes for relativistic magnetohydrodynamics (RMHD) is the recovery of primitive variables from conservative ones. This process involves solving highly nonlinear equations subject to physical constraints. An ideal solver should be "robust, accurate, and fast -- it is at the heart of all conservative RMHD schemes," as emphasized in… ▽ More A long-standing and formidable challenge faced by all conservative schemes for relativistic magnetohydrodynamics (RMHD) is the recovery of primitive variables from conservative ones. This process involves solving highly nonlinear equations subject to physical constraints. An ideal solver should be "robust, accurate, and fast -- it is at the heart of all conservative RMHD schemes," as emphasized in [S.C. Noble et al., ApJ, 641:626-637, 2006]. Despite over three decades of research, seeking efficient solvers that can provably guarantee stability and convergence remains an open problem. This paper presents the first theoretical analysis for designing a robust, physical-constraint-preserving (PCP), and provably (quadratically) convergent Newton-Raphson (NR) method for primitive variable recovery in RMHD. Our key innovation is a unified approach for the initial guess, devised based on sophisticated analysis. It ensures that the NR iteration consistently converges and adheres to physical constraints. Given the extreme nonlinearity and complexity of the iterative function, the theoretical analysis is highly nontrivial and technical. We discover a pivotal inequality for delineating the convexity and concavity of the iterative function and establish theories to guarantee the PCP property and convergence. We also develop theories to determine a computable initial guess within a theoretical "safe" interval. Intriguingly, we find that the unique positive root of a cubic polynomial always falls within this interval. Our PCP NR method is versatile and can be seamlessly integrated into any RMHD scheme that requires the recovery of primitive variables, potentially leading to a broad impact in this field. As an application, we incorporate it into a discontinuous Galerkin method, resulting in fully PCP schemes. Several numerical experiments demonstrate the efficiency and robustness of the PCP NR method. △ Less

Submitted 8 April, 2024; originally announced April 2024.

Comments: 26 pages, 7 figures

Internal magnetic field structures observed by PSP/WISPR in a filament related coronal mass ejection

Authors: G. M. Cappello, M. Temmer, A. Vourlidas, C. Braga, P. C. Liewer, J. Qiu, G. Stenborg, A. Kouloumvakos, A. M. Veronig, V. Bothmer

Abstract: We track and investigate from white-light data taken with the Wide-field Instrument for Solar PRobe (WISPR) aboard Parker Solar Probe (PSP), localized density enhancements, reflecting small-scale magnetic structures belonging to a filament-related coronal mass ejection (CME). We aim to investigate the 3D location, morphology, and evolution of the internal magnetic fine structures of CMEs. Specific… ▽ More We track and investigate from white-light data taken with the Wide-field Instrument for Solar PRobe (WISPR) aboard Parker Solar Probe (PSP), localized density enhancements, reflecting small-scale magnetic structures belonging to a filament-related coronal mass ejection (CME). We aim to investigate the 3D location, morphology, and evolution of the internal magnetic fine structures of CMEs. Specifically, we ask: what is their relationship with the filament/source region and the flux rope? The fast tangential motion of the PSP spacecraft during its perihelion permits viewing the same event from multiple angles in short times relative to the event's evolution. Hence, we can derive the three-dimensional information of selected CME features from a single spacecraft using triangulation techniques. We group small-scale structures with roughly similar speeds, longitude and latitude, into three distinct morphological groups. We find twisted magnetic field patterns close to the eastern leg of the CME that may be related to 'horns' outlining the edges of the flux-rope cavity. Aligned thread-like bundles are identified close to the western leg. They may be related to confined density enhancements evolving during the filament eruption. High density blob-like features (magnetic islands) are widely spread in longitude ($\sim$40°) close to the flanks and rear part of the CME. We demonstrate that CME flux ropes may comprise different morphological groups with a cluster behavior, apart from the blobs which instead span a wide range of longitudes. This may hint either to the three-dimensionality of the post-CME current sheet (CS) or to the influence of the ambient corona in the evolutionary behavior of the CS. Importantly, we show that the global appearance of the CME can be very different in WISPR (0.11--0.16~AU) and instruments near 1~AU because of shorter line-of-sight integration of WISPR. △ Less

Submitted 22 February, 2024; originally announced February 2024.

A $λ$ 3 mm line survey towards the circumstellar envelope of the carbon-rich AGB star IRC+10216 (CW Leo)

Authors: Juan Tuo, Xiaohu Li, Jixian Sun, Tom J. Millar, Yong Zhang, Jianjie Qiu, Donghui Quan, Jarken Esimbek, Jianjun Zhou, Yu Gao, Qiang Chang, Lin Xiao, Yanan Feng, Zhenzhen Miao, Rong Ma, Ryszard Szczerba, Xuan Fang

Abstract: We present an unbiased $λ$ 3 mm spectral line survey (between 84.5 and 115.8 GHz), conducted by the Purple Mountain Observatory 13.7 meter radio telescope, together with updated modeling results, towards the carbon-rich Asymptotic Giant Branch star, IRC+10216 (CW Leo). A total of 75 spectral lines (96 transitions) are detected, and identified to arise from 19 molecules: C$_2$H, $l$-C$_3$H, C$_4$H,… ▽ More We present an unbiased $λ$ 3 mm spectral line survey (between 84.5 and 115.8 GHz), conducted by the Purple Mountain Observatory 13.7 meter radio telescope, together with updated modeling results, towards the carbon-rich Asymptotic Giant Branch star, IRC+10216 (CW Leo). A total of 75 spectral lines (96 transitions) are detected, and identified to arise from 19 molecules: C$_2$H, $l$-C$_3$H, C$_4$H, CN, C$_3$N, HC$_3$N, HC$_5$N, HCN, HNC, CH$_3$CN, MgNC, CO, $c$-C$_3$H$_2$, SiC$_2$, SiO, SiS, CS, C$_2$S, C$_3$S, and their isotopologues. Among them, one molecular emission line (H$^{13}$CCCN $J=13-12$) is discovered in IRC+10216 for the first time. The excitation temperature, column density, and fractional abundance of the detected species are deduced by assuming they are in local thermodynamic equilibrium. In addition, the isotopic ratios of [$^{12}$C]/[$^{13}$C], [$^{32}$S]/[$^{34}$S], [$^{28}$Si]/[$^{29}$Si], and [$^{12}$C$^{34}$S]/[$^{13}$C$^{32}$S] are obtained and found to be consistent with previous studies. Finally, we summarize all of the 106 species detected in IRC+10216 to date with their observed and modeled column densities for the convenience of future studies. △ Less

Submitted 30 January, 2024; v1 submitted 29 January, 2024; originally announced January 2024.

Comments: 71 pages, 39 figures, 10 tables. Accepted for publication in ApJS

Physical Parameters of 11,100 Short-Period ASAS-SN Eclipsing Contact Binaries

Authors: Xu-Zhi Li, Qing-Feng Zhu, Xu Ding, Xiao-Hui Xu, Hang Zheng, Jin-Sheng Qiu, Ming-Chao Liu

Abstract: Starting from more than 11,200 short-period (less than 0.5 days) EW-type eclipsing binary candidates with the All-Sky Automated Survey for Supernovae (ASAS-SN) V-band light curves, we use MCMC and neural networks (NNs) to obtain the mass ratio ($q$), orbital inclination ($incl$), fill-out factor ($f$) and temperature ratio ($T_s/T_p$). After cross-matching with the Gaia DR3 database, the final sam… ▽ More Starting from more than 11,200 short-period (less than 0.5 days) EW-type eclipsing binary candidates with the All-Sky Automated Survey for Supernovae (ASAS-SN) V-band light curves, we use MCMC and neural networks (NNs) to obtain the mass ratio ($q$), orbital inclination ($incl$), fill-out factor ($f$) and temperature ratio ($T_s/T_p$). After cross-matching with the Gaia DR3 database, the final sample contains parameters of 2,399 A-type and 8,712 W-type contact binaries (CBs). We present the distributions of parameters of these 11,111 short-period CBs. The mass ratio ($q$) and fill-out factor ($f$) are found to obey log-normal distributions, and the remaining parameters obey normal distributions. There is a significant period-temperature correlation of these CBs. Additionally, the temperature ratio (${T_s}$/${T_p}$) tends to increase as the orbital period decreases for W-type CBs. There is no significant correlation between them for A-type CBs. The mass ratio and fill-out factor ($q-f$) diagram suggest there is no significant correlation between these two parameters. A clear correlation exists between the mass ratio and radius ratio. The radius ratio increases with the mass ratio. Moreover, the deep fill-out CBs tend to fall on the upper boundary of the $q$$-$${R_s}$/${R_p}$ distribution, while the shallow fill-out CBs fall on the lower boundary. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 15 pages, 8 figures, 1 table.Accepted for publication in ApJS. The full dataset will be available with the ApJS article

Gas-phase molecules in protoplanetary nebulae with the 21 μm emission feature

Authors: Jian-Jie Qiu, Yong Zhang, Jun-ichi Nakashima, Jiang-Shui Zhang, Fei Li, Deng-Rong Lu, Xin-Di Tang, Xiao-Ling Yu, Lan-Wei Jia

Abstract: It has been more than 30 years since the enigmatic 21 μm emission feature was first discovered in protoplanetary nebulae (PPNs). Although dozens of different dust carrier candidates have been proposed, there is as yet no widely accepted one. We present the results of molecular observations toward 21μm objects using the 10m Submillimeter Telescope of Arizona Radio Observatory at the 1.3 mm band and… ▽ More It has been more than 30 years since the enigmatic 21 μm emission feature was first discovered in protoplanetary nebulae (PPNs). Although dozens of different dust carrier candidates have been proposed, there is as yet no widely accepted one. We present the results of molecular observations toward 21μm objects using the 10m Submillimeter Telescope of Arizona Radio Observatory at the 1.3 mm band and the 13.7 m telescope of Purple Mountain Observatory at the 3mm band, aiming to investigate whether the gas-phase environments of these unusual sources have some peculiarities compared to normal PPNs. We detect 31 emission lines belonging to seven different molecular species, most of which are the first detection in 21 μm PPNs. The observations provide clues on the identification of the 21 μm feature. We report a correlation study between the fractional abundance of gas-phase molecules and the strengths of the 21 μm emission. Our study shows that given the small sample size, the 21 μm feature has weak or no correlations with the gas-phase molecules. Future radio observations of high spatial and spectral resolution toward a large sample are desirable to elucidate the 21 μm emission phenomena. △ Less

Submitted 30 December, 2023; originally announced January 2024.

Comments: 40 pages, 10 figures, 11 tables, accepted for publication in The Astronomical Journal

Rotating vector model and radius-to-frequency mapping in the presence of multipole magnetic field

Authors: J. L. Qiu, H. Tong, H. G. Wang

Abstract: The rotating vector model and radius-to-frequency mapping in the presence of multipole magnetic field in pulsars and magnetars are considered. An axisymmetric potential field is assumed. It is found that: (1) The radiation beam in the case of multipole field is wider than the dipole case. This may account the increasing pulse width at higher frequency of pulsars (anti-radius-to-frequency mapping).… ▽ More The rotating vector model and radius-to-frequency mapping in the presence of multipole magnetic field in pulsars and magnetars are considered. An axisymmetric potential field is assumed. It is found that: (1) The radiation beam in the case of multipole field is wider than the dipole case. This may account the increasing pulse width at higher frequency of pulsars (anti-radius-to-frequency mapping). (2) The expression for the polarization position angle is unchanged. Only the inclination angle α and phase constant φ_0 will change. The angle between the rotational axis and line of sight, and the position angle constant ψ_0 will not change. When fitting the varying position angle of magnetars, these constraints should be considered. The appearance and disappearance of multipole field may account for the changing slope of position angle in the radio emitting magnetar Swift J1818.0-1607. Similar but more active process in magnetar magnetospheres may account for the diverse position angle in fast radius bursts. △ Less

Submitted 10 October, 2023; originally announced October 2023.

Comments: accepted in ApJ. significantly rewritten compared with the original version

Systematic Search for Water Fountain Candidates using the Databases of Circumstellar Maser Sources

Authors: Haichen Fan, Jun-ichi Nakashima, D. Engels, Yong Zhang, Jian-Jie Qiu, Huan-Xue Feng, Jia-Yong Xie, Hiroshi Imai, Chih-Hao Hsia

Abstract: Water fountains (WFs) are thought to be objects in the morphological evolution of the circumstellar envelopes of low- and intermediate-mass evolved stars, transitioning from spherically symmetric to asymmetric shapes. We used databases of circumstellar 1612 MHz OH and 22.235 GHz H$_2$O maser sources to search for new WF candidates using the criterion of a larger velocity range of the H$_2$O maser… ▽ More Water fountains (WFs) are thought to be objects in the morphological evolution of the circumstellar envelopes of low- and intermediate-mass evolved stars, transitioning from spherically symmetric to asymmetric shapes. We used databases of circumstellar 1612 MHz OH and 22.235 GHz H$_2$O maser sources to search for new WF candidates using the criterion of a larger velocity range of the H$_2$O maser emission compared to that of the OH maser emission. Thus, it is in principle possible to identify WFs with an H$_2$O velocity ranges smaller than those for the previously known WFs. For the OH maser line, we analyzed database entries of 8,474 observations from 2,195 sources, and 6,085 observations from 3,642 sources for H$_2$O maser line. After a close examination of the velocity ranges and line profiles, we identified 11 sources that meet the criterion mentioned above. We examined the IRAS colors of the selected sources and found that two of them (IRAS 19069+0916 and IRAS 19319+2214) are in the color region for post-AGB stars. We find that the maser velocity criterion can discover other astrophysically interesting objects than just WFs. Such objects may include peculiar planetary nebulae with maser emissions and stellar merger remnants. △ Less

Submitted 8 October, 2023; originally announced October 2023.

Comments: 34 pages, 10 figures, Accepted for publication in ApJS

A 3 mm Molecular Line Survey toward The C-star Envelope CIT 6

Authors: Kui-Lu Yang, Yong Zhang, Jian-Jie Qiu, Yi-Na Ao, Xiao-Hu Li

Abstract: We present an unbiased molecular line survey toward the carbon-rich circumstellar envelope CIT 6 carried out between 90 and 116 GHz with the Arizona Radio Observatory 12 m telescope. A total of 42 lines assigned to 10 molecular species and 4 isotopologues are detected. Despite the absence of any newly identified circumstellar molecules, several transitions are freshly reported for this object. Thi… ▽ More We present an unbiased molecular line survey toward the carbon-rich circumstellar envelope CIT 6 carried out between 90 and 116 GHz with the Arizona Radio Observatory 12 m telescope. A total of 42 lines assigned to 10 molecular species and 4 isotopologues are detected. Despite the absence of any newly identified circumstellar molecules, several transitions are freshly reported for this object. This work is a complement to our previous line survey toward CIT 6 in the frequency ranges of 36-49 GHz and 131-268 GHz. Based on the measurements and in combination with previously published data, we perform a rotation-diagram analysis to determine the column densities, excitation temperatures, and fractional abundances of the molecules. The excitation temperature varies along the radius. The abundance pattern of CIT6 is broadly consistent with that of the prototype C-star IRC+10216 at the specified detection sensitivity, indicating that the molecular richness in IRC+10216 cannot be attributed to an interpretation of unusual chemistry. Nevertheless, subtle distinctions between the two C-stars are found. The higher abundance of carbon-chain molecules and lower 12C/13C ratio in CIT 6 compared to IRC+10216 imply that the former is more massive or in a more evolved phase than the latter. △ Less

Submitted 29 July, 2023; originally announced July 2023.

Comments: 34 pages, 14 figures, 5 tables

Journal ref: PASJ (2023)

The Role of Magnetic Shear in Reconnection-Driven Flare Energy Release

Authors: J. Qiu, M. Alaoui, S. K. Antiochos, J. T. Dahlin, M. Swisdak, J. F. Drake, A. Robison, C. R. DeVore, V. M. Uritsky

Abstract: Using observations from the Solar Dynamics Observatory's Atmosphere Imaging Assembly and the Ramaty High Energy Solar Spectroscopic Imager, we present novel measurements of the shear of post-reconnection flare loops (PRFLs) in SOL20141218T21:40 and study its evolution with respect to magnetic reconnection and flare emission. Two quasi-parallel ribbons form adjacent to the magnetic polarity inversi… ▽ More Using observations from the Solar Dynamics Observatory's Atmosphere Imaging Assembly and the Ramaty High Energy Solar Spectroscopic Imager, we present novel measurements of the shear of post-reconnection flare loops (PRFLs) in SOL20141218T21:40 and study its evolution with respect to magnetic reconnection and flare emission. Two quasi-parallel ribbons form adjacent to the magnetic polarity inversion line (PIL), spreading in time first parallel to the PIL and then mostly in a perpendicular direction. We measure magnetic reconnection rate from the ribbon evolution, and also the shear angle of a large number of PRFLs observed in extreme ultraviolet passbands ($\lesssim$1 MK). For the first time, the shear angle measurements are conducted using several complementary techniques allowing for a cross-validation of the results. In this flare, the total reconnection rate is much enhanced before a sharp increase of the hard X-ray emission, and the median shear decreases from 60$^\circ$-70$^\circ$ to 20$^\circ$, on a time scale of ten minutes. We find a correlation between the shear-modulated total reconnection rate and the non-thermal electron flux. These results confirm the strong-to-weak shear evolution suggested in previous observational studies and reproduced in numerical models, and also confirm that, in this flare, reconnection is not an efficient producer of energetic non-thermal electrons during the first ten minutes when the strongly sheared PRFLs are formed. We conclude that an intermediate shear angle, $\le 40^\circ$, is needed for efficient particle acceleration via reconnection, and we propose a theoretical interpretation. △ Less

Submitted 15 August, 2023; v1 submitted 26 June, 2023; originally announced June 2023.

Comments: 19 pages, 10 figures

Provably convergent Newton-Raphson methods for recovering primitive variables with applications to physical-constraint-preserving Hermite WENO schemes for relativistic hydrodynamics

Authors: Chaoyi Cai, Jianxian Qiu, Kailiang Wu

Abstract: The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as… ▽ More The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as the primitive variables cannot be explicitly reformulated using conservative variables due to relativistic effects. In this paper, we propose three efficient Newton--Raphson (NR) methods for robustly recovering primitive variables from conservative variables. Importantly, we rigorously prove that these NR methods are always convergent and PCP, meaning they preserve the physical constraints throughout the NR iterations. The discovery of these robust NR methods and their PCP convergence analyses are highly nontrivial and technical. As an application, we apply the proposed NR methods to design PCP finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the RHD equations. Our PCP HWENO schemes incorporate high-order HWENO reconstruction, a PCP limiter, and strong-stability-preserving time discretization. We rigorously prove the PCP property of the fully discrete schemes using convex decomposition techniques. Moreover, we suggest the characteristic decomposition with rescaled eigenvectors and scale-invariant nonlinear weights to enhance the performance of the HWENO schemes in simulating large-scale RHD problems. Several demanding numerical tests are conducted to demonstrate the robustness, accuracy, and high resolution of the proposed PCP HWENO schemes and to validate the efficiency of our NR methods. △ Less

Submitted 24 May, 2023; originally announced May 2023.

Comments: 49 pages

A New Period Determination Method for Periodic Variable Stars

Authors: Xiao-Hui Xu, Qing-Feng Zhu, Xu-Zhi Li, Bin Li, Hang Zheng, Jin-Sheng Qiu, Hai-Bin Zhao

Abstract: Variable stars play a key role in understanding the Milky Way and the universe. The era of astronomical big data presents new challenges for quick identification of interesting and important variable stars. Accurately estimating the periods is the most important step to distinguish different types of variable stars. Here, we propose a new method of determining the variability periods. By combining… ▽ More Variable stars play a key role in understanding the Milky Way and the universe. The era of astronomical big data presents new challenges for quick identification of interesting and important variable stars. Accurately estimating the periods is the most important step to distinguish different types of variable stars. Here, we propose a new method of determining the variability periods. By combining the statistical parameters of the light curves, the colors of the variables, the window function and the GLS algorithm, the aperiodic variables are excluded and the periodic variables are divided into eclipsing binaries and NEB variables (other types of periodic variable stars other than eclipsing binaries), the periods of the two main types of variables are derived. We construct a random forest classifier based on 241,154 periodic variables from the ASAS-SN and OGLE datasets of variables. The random forest classifier is trained on 17 features, among which 11 are extracted from the light curves and 6 are from the Gaia Early DR3, ALLWISE and 2MASS catalogs. The variables are classified into 7 superclasses and 17 subclasses. In comparison with the ASAS-SN and OGLE catalogs, the classification accuracy is generally above approximately 82% and the period accuracy is 70%-99%. To further test the reliability of the new method and classifier, we compare our results with the results of Chen et al. (2020) for ZTF DR2. The classification accuracy is generally above 70%. The period accuracy of the EW and SR variables is 50% and 53%, respectively. And the period accuracy of other types of variables is 65%-98%. △ Less

Submitted 20 December, 2022; originally announced December 2022.

Comments: 23 pages, 10 figures

Journal ref: PASP2022

Investigating pre-eruptive magnetic properties at the footprints of erupting magnetic flux ropes

Authors: Wensi Wang, Jiong Qiu, Rui Liu, Chunming Zhu, Kai E Yang, Qiang Hu, Yuming Wang

Abstract: It is well established that solar eruptions are powered by free magnetic energy stored in current-carrying magnetic field in the corona. It has also been generally accepted that magnetic flux ropes (MFRs) are a critical component of many coronal mass ejections (CMEs). What remains controversial is whether MFRs are present well before the eruption. Our aim is to identify progenitors of MFRs, and in… ▽ More It is well established that solar eruptions are powered by free magnetic energy stored in current-carrying magnetic field in the corona. It has also been generally accepted that magnetic flux ropes (MFRs) are a critical component of many coronal mass ejections (CMEs). What remains controversial is whether MFRs are present well before the eruption. Our aim is to identify progenitors of MFRs, and investigate pre-eruptive magnetic properties associated with these progenitors. Here we analyze 28 MFRs erupting within 45 deg from the disk center from 2010 to 2015. All MFRs'feet are well identified by conjugate coronal dimmings. We then calculate magnetic properties at the feet of the MFRs, prior to their eruptions, using Helioseismic and Magnetic Imager (HMI) vector magnetograms. Our results show that only 8 erupting MFRs are associated with significant non-neutralized electric currents, 4 of which also exhibit pre-eruptive dimmings at the foot-prints. Twist and current distributions are asymmetric at the two feet of these MFRs. The presence of pre-eruption dimmings associated with non-neutralized currents suggests the pre-existing MFRs. Furthermore, evolution of conjugate dimmings and electric currents within the foot-prints can provide clues about the internal structure of MFRs and their formation mechanism. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: 40 pages, 13 figures, 6 tables, Accepted for publication in ApJ

Molecules in the peculiar age-defying source IRAS 19312+1950

Authors: Jian-Jie Qiu, Yong Zhang, Jun-ichi Nakashima, Jiang-Shui Zhang, Nico Koning, Xin-Di Tang, Yao-Ting Yan, Huan-Xue Feng

Abstract: Context. IRAS 19312+1950 is an isolated infrared source that exhibits a characteristic quasi-point-symmetric morphology in the near- and mid-infrared images and is also very bright in molecular radio lines. Because of its unique observational characteristics, various observational studies have been conducted and several hypotheses have been proposed regarding its origin, which is still unclear. So… ▽ More Context. IRAS 19312+1950 is an isolated infrared source that exhibits a characteristic quasi-point-symmetric morphology in the near- and mid-infrared images and is also very bright in molecular radio lines. Because of its unique observational characteristics, various observational studies have been conducted and several hypotheses have been proposed regarding its origin, which is still unclear. So far, it has been suggested that it could be a peculiar evolved star, a young stellar object, or even a red nova remnant. Regardless of which type of object it is ultimately classified as, IRAS 19312+1950 is exceptionally bright in the infrared and molecular radio lines and therefore will undoubtedly be crucial as a prototype of this kind of object having a peculiar nature or unusual evolutionary phase. Aims. This study aims to reveal the molecular composition of the central part of IRAS 19312+1950 by performing an unbiased molecular radio line survey and discussing the origin of the object from a molecular chemical point of view. Methods. We carried out a spectral line survey with the IRAM 30 m telescope towards the center of IRAS 19312+1950 in the 3 and 1.3 mm windows. Results. In total, 28 transition lines of 22 molecular species and those isotopologues are detected towards IRAS 19312+1950, some of which exhibit a broad and a narrow components. Seventeen thermal lines and 1 maser line are newly detected. The molecular species of C$^{17}$O, $^{30}$SiO, HN$^{13}$C, HC$^{18}$O$^{+}$, H$_{2}$CO, and $c$-C$_{3}$H$_{2}$ are detected for the first time in this object. Conclusions. Our results, in combination with previous studies, favor the hypothesis that IRAS 19312+1950 might be a red nova remnant, in which the progenitors that merged to become a red nova may have contained at least two evolved stars with oxygen-rich and carbon-rich chemistry, respectively. △ Less

Submitted 24 November, 2022; originally announced November 2022.

Comments: 30 pages, 24 figures, accepted for publication in A&A

Journal ref: A&A 669, A121 (2023)

Extracting the Heliographic Coordinates of Coronal Rays using Images from WISPR/Parker Solar Probe

Authors: P. C. Liewer, J. Qiu, F. Ark, P. Penteado, G. Stenborg, A. Vourlidas, J. R. Hall, P. Riley

Abstract: The Wide-field Imager for Solar Probe (WISPR) onboard Parker Solar Probe (PSP), observing in white light, has a fixed angular field of view, extending from 13.5 degree to 108 degree from the Sun and approximately 50 degree in the transverse direction. In January 2021, on its seventh orbit, PSP crossed the heliospheric current sheet (HCS) near perihelion at a distance of 20 solar radii. At this tim… ▽ More The Wide-field Imager for Solar Probe (WISPR) onboard Parker Solar Probe (PSP), observing in white light, has a fixed angular field of view, extending from 13.5 degree to 108 degree from the Sun and approximately 50 degree in the transverse direction. In January 2021, on its seventh orbit, PSP crossed the heliospheric current sheet (HCS) near perihelion at a distance of 20 solar radii. At this time, WISPR observed a broad band of highly variable solar wind and multiple coronal rays. For six days around perihelion, PSP was moving with an angular velocity exceeding that of the Sun. During this period, WISPR was able to image coronal rays as PSP approached and then passed under or over them. We have developed a technique for using the multiple viewpoints of the coronal rays to determine their location (longitude and latitude) in a heliocentric coordinate system and used the technique to determine the coordinates of three coronal rays. The technique was validated by comparing the results to observations of the coronal rays from Solar and Heliophysics Observatory (SOHO) / Large Angle and Spectrometric COronagraph (LASCO)/C3 and Solar Terrestrial Relations Observatory (STEREO)-A/COR2. Comparison of the rays' locations were also made with the HCS predicted by a 3D MHD model. In the future, results from this technique can be used to validate dynamic models of the corona. △ Less

Submitted 6 September, 2022; originally announced September 2022.

Comments: 14 pages, 13 figures

Development of Integral Field Spectrographs to Revolutionize Spectroscopic Observations of Solar Flares and other Energetic Solar Eruptions

Authors: Haosheng Lin, Tetsu Anan, Gianna Cauzzi, Lyndsay Fletcher, Pei Huang, Adam Kowalski, Maxim Kramar, Jiong Qiu, Jenna Samra, Constance Spittler, Takashi Sukegawa, Gregory Wirth

Abstract: The Sun's proximity offers us a unique opportunity to study in detail the physical processes on a star's surface; however, the highly dynamic nature of the stellar surface -- in particular, energetic eruptions such as flares and coronal mass ejections -- presents tremendous observational challenges. Spectroscopy probes the physical state of the solar atmosphere, but conventional scanning spectrogr… ▽ More The Sun's proximity offers us a unique opportunity to study in detail the physical processes on a star's surface; however, the highly dynamic nature of the stellar surface -- in particular, energetic eruptions such as flares and coronal mass ejections -- presents tremendous observational challenges. Spectroscopy probes the physical state of the solar atmosphere, but conventional scanning spectrographs and spectrometers are unable to capture the full evolutionary history of these dynamic events with a sufficiently wide field of view and high spatial, spectral, and temporal resolution. Resolving the physics of the dynamic sun requires gathering simultaneous spectra across a contiguous area over the full duration of these events, a goal now tantalizingly close to achievable with continued investment in developing powerful new Integral Field Spectrographs to serve as the foundation of both future ground- and space-based missions. This technology promises to revolutionize our ability to study solar flares and CMEs, addressing NASA's strategic objective to "understand the Sun, solar system, and universe." Since such events generate electromagnetic radiation and high-energy particles that disrupt terrestrial electric infrastructure, this investment not only advances humanity's scientific endeavors but also enhances our space weather forecasting capability to protect against threats to our technology-reliant civilization. △ Less

Submitted 1 September, 2022; originally announced September 2022.

Comments: 14 pages, 4 figures, White Paper submitted for 2022 Solar and Space Physics (Heliophysics) Decadal Survey

Interstellar Nitrogen Isotope Ratios: New NH3 Data from the Galactic Center out to the Perseus Arm

Authors: J. L. Chen, J. S. Zhang, C. Henkel, Y. T. Yan, H. Z. Yu, J. J. Qiu, X. D. Tang, J. Wang, W. Liu, Y. X. Wang, Y. H. Zheng, J. Y. Zhao, Y. P. Zou

Abstract: Our aim is to measure the interstellar 14N/15N ratio across the Galaxy, to establish a standard data set on interstellar ammonia isotope ratios, and to provide new constraints on the Galactic chemical evolution. The (J, K ) = (1, 1), (2, 2), and (3, 3) lines of 14NH3 and 15NH3 were observed with the Shanghai Tianma 65 m radio telescope (TMRT) and the Effelsberg 100 m telescope toward a large sampl… ▽ More Our aim is to measure the interstellar 14N/15N ratio across the Galaxy, to establish a standard data set on interstellar ammonia isotope ratios, and to provide new constraints on the Galactic chemical evolution. The (J, K ) = (1, 1), (2, 2), and (3, 3) lines of 14NH3 and 15NH3 were observed with the Shanghai Tianma 65 m radio telescope (TMRT) and the Effelsberg 100 m telescope toward a large sample of 210 sources. One hundred fourty-one of these sources were detected by the TMRT in 14NH3. Eight of them were also detected in 15NH3. For 10 of the 36 sources with strong NH3 emission, the Effelsberg 100 m telescope successfully detected their 15NH3(1, 1) lines, including 3 sources (G081.7522, W51D, and Orion-KL) with detections by the TMRT telescope. Thus, a total of 15 sources are detected in both the 14NH3 and 15NH3 lines. Line and physical parameters for these 15 sources are derived, including optical depths, rotation and kinetic temperatures, and total column densities. 14N/15N isotope ratios were determined from the 14NH3/15NH3 abundance ratios. The isotope ratios obtained from both telescopes agree for a given source within the uncertainties, and no dependence on heliocentric distance and kinetic temperature is seen. 14N/15N ratios tend to increase with galactocentric distance, confirming a radial nitrogen isotope gradient. This is consistent with results from recent Galactic chemical model calculations, including the impact of superasymptotic giant branch stars and novae. △ Less

Submitted 8 August, 2022; originally announced August 2022.

Cyanopolyyne line survey towards high-mass star-forming regions with TMRT

Authors: Y. X. Wang, J. S. Zhang, Y. T. Yan, J. J. Qiu, J. L. Chen, J. Y. Zhao, Y. P. Zou, X. C. Wu, X. L. He, Y. B. Gong, J. H. Cai

Abstract: We carried out a cyanopolyyne line survey towards a large sample of HMSFRs using the Shanghai Tian Ma 65m Radio Telescope (TMRT). Our sample consisted of 123 targets taken from the TMRT C band line survey. It included three kinds of sources, namely those with detection of the 6.7 GHz CH3OH maser alone, with detection of the radio recombination line (RRL) alone, and with detection of both (hereafte… ▽ More We carried out a cyanopolyyne line survey towards a large sample of HMSFRs using the Shanghai Tian Ma 65m Radio Telescope (TMRT). Our sample consisted of 123 targets taken from the TMRT C band line survey. It included three kinds of sources, namely those with detection of the 6.7 GHz CH3OH maser alone, with detection of the radio recombination line (RRL) alone, and with detection of both (hereafter referred to as Maser-only, RRL-only, and Maser-RRL sources, respectively). We detected HC3N in 38 sources, HC5N in 11 sources, and HC7N in G24.790+0.084, with the highest detection rate being found for Maser-RRL sources and a very low detection rate found for RRL-only sources. Their column densities were derived using the rotational temperature measured from the NH3 lines. And we constructed and fitted the far-infrared (FIR) spectral energy distributions. Based on these, we derive their dust temperatures, H2 column densities, and abundances of cyanopolyynes relative to H2. The detection rate, the column density, and the relative abundance of HC3N increase from Maser-only to Maser-RRL sources and decrease from Maser-RRL to RRL-only sources. This trend is consistent with the proposed evolutionary trend of HC3N under the assumption that our Maser-only, Maser-RRL, and RRL-only sources correspond to massive young stellar objects, ultra-compact HII regions, and normal classical HII regions, respectively. Furthermore, a statistical analysis of the integrated line intensity and column density of HC3N and shock-tracing molecules (SiO, H2CO) enabled us to find positive correlations between them. This suggests that HC3N may be another tracer of shocks, and should therefore be the subject of further observations and corresponding chemical simulations. Our results indirectly support the idea that the neutral--neutral reaction between C2H2 and CN is the dominant formation pathway of HC3N. △ Less

Submitted 16 May, 2022; originally announced May 2022.

Comments: 23 pages, 5 figures, 4 tables, Accepted to A&A

Journal ref: A&A 663, A177 (2022)

Properties and Energetics of Magnetic Reconnection: I. Evolution of Flare Ribbons

Authors: Jiong Qiu, Jianxia Cheng

Abstract: In this article, we measure the mean magnetic shear from the morphological evolution of flare ribbons, and examine the evolution of flare thermal and non-thermal X-ray emissions during the progress of flare reconnection. We analyze three eruptive flares and three confined flares ranging from GOES class C8.0 to M7.0. They exhibit well-defined two ribbons along the magnetic polarity inversion line (… ▽ More In this article, we measure the mean magnetic shear from the morphological evolution of flare ribbons, and examine the evolution of flare thermal and non-thermal X-ray emissions during the progress of flare reconnection. We analyze three eruptive flares and three confined flares ranging from GOES class C8.0 to M7.0. They exhibit well-defined two ribbons along the magnetic polarity inversion line (PIL), and have been observed by the Atmospheric Imaging Assembly and the Ramaty High Energy Solar Spectroscopic Imager from the onset of the flare throughout the impulsive phase. The analysis confirms the strong-to-weak shear evolution in the core region of the flare, and the flare hard X-ray emission rises as the shear decreases. In eruptive flares in this sample, significant non-thermal hard X-ray emission lags the ultraviolet emission from flare ribbons, and rises rapidly when the shear is modest. In all flares, we observe that the plasma temperature rises in the early phase when the flare ribbons rapidly spread along the PIL and the shear is high. We compare these results with prior studies, and discuss their implications, as well as complications, related to physical mechanisms governing energy partition during flare reconnection. △ Less

Submitted 6 May, 2022; originally announced May 2022.

Comments: 28 pages, 8 figures; Solar Physics

Validation and interpretation of three-dimensional configuration of a magnetic cloud flux rope

Authors: Qiang Hu, Chunming Zhu, Wen He, Jiong Qiu, Lan K. Jian, Avijeet Prasad

Abstract: One "strong" magnetic cloud (MC) with the magnetic field magnitude reaching $\sim$ 40 nT at 1 au during 2012 June 16-17 is examined in association with a pre-existing magnetic flux rope (MFR) identified on the Sun. The MC is characterized by a quasi-three dimensional (3D) flux rope model based on in situ measurements from the Wind spacecraft. The magnetic flux contents and other parameters are qua… ▽ More One "strong" magnetic cloud (MC) with the magnetic field magnitude reaching $\sim$ 40 nT at 1 au during 2012 June 16-17 is examined in association with a pre-existing magnetic flux rope (MFR) identified on the Sun. The MC is characterized by a quasi-three dimensional (3D) flux rope model based on in situ measurements from the Wind spacecraft. The magnetic flux contents and other parameters are quantified. In addition, a correlative study with the corresponding measurements of the same structure crossed by the Venus Express (VEX) spacecraft at a heliocentric distance 0.7 au and with an angular separation $\sim 6^\circ$ in longitude is performed to validate the MC modeling results. The spatial variation between the Wind and VEX magnetic field measurements is attributed to the 3D configuration of the structure as featured by a knotted bundle of flux. The comparison of the magnetic flux contents between the MC and the source region on the Sun indicates that the 3D reconnection process accompanying an M1.9 flare may correspond to the magnetic reconnection between the field lines of the pre-existing MFR rooted in the opposite polarity footpoints. Such a process reduces the amount of the axial magnetic flux in the erupted flux rope, by approximately 50\%, in this case. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: Submitted to ApJ

Molecules in the carbon-rich protoplanetary nebula CRL 2688

Authors: Jian-Jie Qiu, Yong Zhang, Jiang-Shui Zhang, Jun-ichi Nakashima

Abstract: We present observations of the carbon-rich protoplanetary nebula (PPN) CRL 2688 made with the Institut de Radioastronomie Millimetrique (IRAM) 30 m telescope in the 3mm and 2mm bands. In total, 196 transition lines belonging to 38 molecular species and isotopologues are detected, among which, to our best knowledge, 153 transition lines and 13 species are the first report for this object. Additiona… ▽ More We present observations of the carbon-rich protoplanetary nebula (PPN) CRL 2688 made with the Institut de Radioastronomie Millimetrique (IRAM) 30 m telescope in the 3mm and 2mm bands. In total, 196 transition lines belonging to 38 molecular species and isotopologues are detected, among which, to our best knowledge, 153 transition lines and 13 species are the first report for this object. Additionally, in order to contribute to future research, we have collected observational data on the molecular lines of CRL 2688 from the literature and compiled them into a single unified catalog. We find that the molecular abundance of CRL 2688 cannot be explained by the standard model of a circumstellar envelope. The implications of metal-bearing molecules on circumstellar chemistry are discussed. △ Less

Submitted 6 April, 2022; originally announced April 2022.

Comments: 55 pages, 10 figures, Accepted for publication in ApJS

Journal ref: ApJS, 2022, 259, 56

Quantitative Characterization of Magnetic Flux Rope Properties for Two Solar Eruption Events

Authors: Wen He, Qiang Hu, Chaowei Jiang, Jiong Qiu, Avijeet Prasad

Abstract: In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR du… ▽ More In order to bridge the gap between heliospheric and solar observations of coronal mass ejections (CMEs), one of the key steps is to improve the understanding of their corresponding magnetic structures like the magnetic flux ropes (MFRs). But it remains a challenge to confirm the existence of a coherent MFR before or upon the CME eruption on the Sun and to quantitatively characterize the CME-MFR due to the lack of direct magnetic field measurement in the corona. In this study, we investigate the MFR structures, originating from two active regions (ARs), AR 11719 and AR 12158, and estimate their magnetic properties quantitatively. We perform the nonlinear force-free field extrapolations with preprocessed photospheric vector magnetograms. In addition, remote-sensing observations are employed to find indirect evidence of MFRs on the Sun and to analyze the time evolution of magnetic reconnection flux associated with the flare ribbons during the eruption. A coherent "pre-existing" MFR structure prior to the flare eruption is identified quantitatively for one event from the combined analysis of the extrapolation and observation. Then the characteristics of MFRs for two events on the Sun before and during the eruption, forming the CME-MFR, including the axial magnetic flux, field-line twist, and reconnection flux, are estimated and compared with the corresponding in situ modeling results. We find that the magnetic reconnection associated with the accompanying flares for both events injects significant amount of flux into the erupted CME-MFRs. △ Less

Submitted 9 January, 2022; originally announced January 2022.

Comments: 28 pages, 10 figures, 3 tables, submitted to ApJ

Connecting Chromospheric Condensation Signatures to Reconnection Driven Heating Rates in an Observed Flare

Authors: William H. Ashfield IV, Dana W. Longcope, Chunming Zhu, Jiong Qiu

Abstract: Observations of solar flare reconnection at very high spatial and temporal resolution can be made indirectly at the footpoints of reconnected loops into which flare energy is deposited. The response of the lower atmosphere to this energy input includes a downward-propagating shock called chromospheric condensation, which can be observed in the UV and visible. In order to characterize reconnection… ▽ More Observations of solar flare reconnection at very high spatial and temporal resolution can be made indirectly at the footpoints of reconnected loops into which flare energy is deposited. The response of the lower atmosphere to this energy input includes a downward-propagating shock called chromospheric condensation, which can be observed in the UV and visible. In order to characterize reconnection using high-resolution observations of this response, one must develop a quantitative relationship between the two. Such a relation was recently developed and here we test it on observations of chromospheric condensation in a single footpoint from a flare ribbon of the X1.0 flare on 25 Oct. 2014 (SOL2014-10-25T16:56:36). Measurements taken of Si iv 1402.77 Å emission spectra using the Interface Region Imaging Spectrograph (IRIS) in a single pixel show red-shifted component undergoing characteristic condensation evolution. We apply the technique called the Ultraviolet Footpoint Calorimeter (UFC) to infer energy deposition into the one footpoint. This energy profile, persisting much longer than the observed condensation, is input into a one-dimensional, hydrodynamic simulation to compute the chromospheric response, which contains a very brief condensation episode. From this simulation we synthesize Si iv spectra and compute the time-evolving Doppler velocity. The synthetic velocity evolution is found to compare reasonably well with the IRIS observation, thus corroborating our reconnection-condensation relationship. The exercise reveals that the chromospheric condensation characterizes a particular portion of the reconnection energy release rather than its entirety, and that the time scale of condensation does not necessarily reflect the time scale of energy input. △ Less

Submitted 3 December, 2021; originally announced December 2021.

The Nearby Evolved Stars Survey II: Constructing a volume-limited sample and first results from the James Clerk Maxwell Telescope

Authors: P. Scicluna, F. Kemper, I. McDonald, S. Srinivasan, A. Trejo, S. H. J. Wallström, J. G. A. Wouterloot, J. Cami, J. Greaves, Jinhua He, D. T. Hoai, Hyosun Kim, O. C. Jones, H. Shinnaga, C. J. R. Clark, T. Dharmawardena, W. Holland, H. Imai, J. Th. van Loon, K. M. Menten, R. Wesson, H. Chawner, S. Feng, S. Goldman, F. C. Liu , et al. (67 additional authors not shown)

Abstract: The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a n… ▽ More The Nearby Evolved Stars Survey (NESS) is a volume-complete sample of $\sim$850 Galactic evolved stars within 3\,kpc at (sub-)mm wavelengths, observed in the CO $J = $ (2$-$1) and (3$-$2) rotational lines, and the sub-mm continuum, using the James Clark Maxwell Telescope and Atacama Pathfinder Experiment. NESS consists of five tiers, based on distances and dust-production rate (DPR). We define a new metric for estimating the distances to evolved stars and compare its results to \emph{Gaia} EDR3. Replicating other studies, the most-evolved, highly enshrouded objects in the Galactic Plane dominate the dust returned by our sources, and we initially estimate a total DPR of $4.7\times 10^{-5}$ M$_\odot$ yr$^{-1}$ from our sample. Our sub-mm fluxes are systematically higher and spectral indices are typically shallower than dust models typically predict. The 450/850 $μ$m spectral indices are consistent with the blackbody Rayleigh--Jeans regime, suggesting a large fraction of evolved stars have unexpectedly large envelopes of cold dust. △ Less

Submitted 24 October, 2021; originally announced October 2021.

Comments: 20 pages, 13 figures, 5 tables, accepted for publications in MNRAS

Variability of the Reconnection Guide Field in Solar Flares

Authors: Joel T. Dahlin, Spiro K. Antiochos, Jiong Qiu, C. Richard DeVore

Abstract: Solar flares may be the best-known examples of the explosive conversion of magnetic energy into bulk motion, plasma heating, and particle acceleration via magnetic reconnection. The energy source for all flares is the highly sheared magnetic field of a filament channel above a polarity inversion line (PIL). During the flare, this shear field becomes the so-called reconnection guide field (i.e., th… ▽ More Solar flares may be the best-known examples of the explosive conversion of magnetic energy into bulk motion, plasma heating, and particle acceleration via magnetic reconnection. The energy source for all flares is the highly sheared magnetic field of a filament channel above a polarity inversion line (PIL). During the flare, this shear field becomes the so-called reconnection guide field (i.e., the non-reconnecting component), which has been shown to play a major role in determining key properties of the reconnection including the efficiency of particle acceleration. We present new high-resolution, three-dimensional, magnetohydrodynamics simulations that reveal the detailed evolution of the magnetic shear/guide field throughout an eruptive flare. The magnetic shear evolves in three distinct phases: shear first builds up in a narrow region about the PIL, then expands outward to form a thin vertical current sheet, and finally is transferred by flare reconnection into an arcade of sheared flare loops and an erupting flux rope. We demonstrate how the guide field may be inferred from observations of the sheared flare loops. Our results indicate that initially the guide field is larger by about a factor of 5 than the reconnecting component, but it weakens by more than an order of magnitude over the course of the flare. Instantaneously, the guide field also varies spatially over a similar range along the three-dimensional current sheet. We discuss the implications of the remarkable variability of the guide field for the timing and localization of efficient particle acceleration in flares. △ Less

Submitted 29 July, 2022; v1 submitted 8 October, 2021; originally announced October 2021.

Comments: 15 pages, published in the Astrophysical Journal

Correlated Spatiotemporal Evolution of Extreme-Ultraviolet Ribbons and Hard X-rays in a Solar Flare

Authors: Stephen J. Naus, Jiong Qiu, C. Richard DeVore, Spiro K. Antiochos, Joel T. Dahlin, James F. Drake, Marc Swisdak

Abstract: We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet (UV) images from IRIS and SDO/AIA, magnetic data from SDO/HMI, hard X-ray (HXR) images from RHESSI, and light curves from Fermi/GBM, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. T… ▽ More We analyze the structure and evolution of ribbons from the M7.3 SOL2014-04-18T13 flare using ultraviolet (UV) images from IRIS and SDO/AIA, magnetic data from SDO/HMI, hard X-ray (HXR) images from RHESSI, and light curves from Fermi/GBM, in order to infer properties of coronal magnetic reconnection. As the event progresses, two flare ribbons spread away from the magnetic polarity inversion line. The width of the newly brightened front along the extension of the ribbon is highly intermittent in both space and time, presumably reflecting non-uniformities in the structure and/or dynamics of the flare current sheet. Furthermore, the ribbon width grows most rapidly in regions exhibiting concentrated non-thermal HXR emission, with sharp increases slightly preceding the HXR bursts. The light curve of the ultraviolet emission matches the HXR light curve at photon energies above 25 keV. In other regions the ribbon-width evolution and light curves do not temporally correlate with the HXR emission. This indicates that the production of non-thermal electrons is highly non-uniform within the flare current sheet. Our results suggest a strong connection between the production of non-thermal electrons and the locally enhanced perpendicular extent of flare ribbon fronts, which in turn reflects inhomogeneous structure and/or reconnection dynamics of the current sheet. Despite this variability, the ribbon fronts remain nearly continuous, quasi-one-dimensional features. Thus, although the reconnecting coronal current sheets are highly structured, they remain quasi-two-dimensional and the magnetic energy release occurs systematically, rather than stochastically, through the volume of reconnecting magnetic flux. △ Less

Submitted 5 December, 2021; v1 submitted 30 September, 2021; originally announced September 2021.

The Neupert Effect of Flare UltraViolet and Soft X-ray Emissions

Authors: Jiong Qiu

Abstract: We model the Neupert effect that relates flare heating energy with the observed SXR emission. The traditional form of the Neupert effect refers to the correlation between the time-integrated HXR or microwave light curve and the SXR light curve. In this paper, instead, we use as the proxy for heating energy the ultraviolet (UV) emission at the foot-points of flare loops, and modify the model of the… ▽ More We model the Neupert effect that relates flare heating energy with the observed SXR emission. The traditional form of the Neupert effect refers to the correlation between the time-integrated HXR or microwave light curve and the SXR light curve. In this paper, instead, we use as the proxy for heating energy the ultraviolet (UV) emission at the foot-points of flare loops, and modify the model of the Neupert effect by taking into account the discrete nature of flare heating as well as cooling. In the modified empirical model, spatially resolved UV lightcurves from the transition region or upper chromosphere are each convolved with a kernel function characterizing the decay of the flare loop emission. Contributions by all loops are summed to compare with the observed total SXR emission. The model has successfully reproduced the observed SXR emission from its rise to decay. To estimate heating energies in flare loops, we also employ the UV Foot-point Calorimeter (UFC) method that infers heating rates in flare loops from these UV light curves and models evolution of flare loops with a zero-dimensional hydrodynamic code. The experiments show that a multitude of impulsive heating events do not well reproduce the observed flare SXR light curve, but a two-phase heating model leads to better agreement with observations. Comparison of the two models of the Neupert effect further allows us to calibrate the UFC method, and improve the estimate of heating rates in flare loops continuously formed by magnetic reconnection throughout the flare evolution. △ Less

Submitted 26 January, 2021; originally announced January 2021.

Comments: accepted by ApJ

Evolution of a Streamer-Blowout CME as Observed by Imagers on Parker Solar Probe and the Solar Terrestrial Relations Observatory

Authors: P. C. Liewer, J. Qiu, A. Vourlidas, J. R. Hall, P. Penteado

Abstract: Context: On 26-27 January 2020, the wide-field imager WISPR on Parker Solar Probe (PSP) observed a coronal mass ejection (CME) from a distance of approximately 30 solar radii as it passed through the instrument's 95 degree field-of-view, providing an unprecedented view of the flux rope morphology of the CME's internal structure. The same CME was seen by STEREO, beginning on 25 January. Aims: Our… ▽ More Context: On 26-27 January 2020, the wide-field imager WISPR on Parker Solar Probe (PSP) observed a coronal mass ejection (CME) from a distance of approximately 30 solar radii as it passed through the instrument's 95 degree field-of-view, providing an unprecedented view of the flux rope morphology of the CME's internal structure. The same CME was seen by STEREO, beginning on 25 January. Aims: Our goal was to understand the origin and determine the trajectory of this CME. Methods: We analyzed data from three well-placed spacecraft: Parker Solar Probe (PSP), Solar Terrestrial Relations Observatory-Ahead (STEREO-A), and Solar Dynamics Observatory (SDO). The CME trajectory was determined using the method described in Liewer et al. (2020) and verified using simultaneous images of the CME propagation from STEREO-A. The fortuitous alignment with STEREO-A also provided views of coronal activity leading up to the eruption. Observations from SDO, in conjunction with potential magnetic field models of the corona, were used to analyze the coronal magnetic evolution for the three days leading up to the flux rope ejection from the corona on 25 January. Results: We found that the 25 January CME is likely the end result of a slow magnetic flux rope eruption that began on 23 January and was observed by STEREO-A/Extreme Ultraviolet Imager (EUVI). Analysis of these observations suggest that the flux rope was apparently constrained in the corona for more than a day before its final ejection on 25 January. STEREO-A/COR2 observations of swelling and brightening of the overlying streamer for several hours prior to eruption on January 25 led us to classify this as a streamer-blowout CME. The analysis of the SDO data suggests that restructuring of the coronal magnetic fields caused by an emerging active region led to the final ejection of the flux rope. △ Less

Submitted 14 December, 2020; v1 submitted 9 December, 2020; originally announced December 2020.

Comments: 12 pages, 12 figures. Accepted in A&A

Journal ref: A&A 650, A32 (2021)

On the quasi-three dimensional configuration of magnetic clouds

Authors: Qiang Hu, Wen He, Jiong Qiu, A. Vourlidas, Chunming Zhu

Abstract: We develop an optimization approach to model the magnetic field configuration of magnetic clouds, based on a linear-force free formulation in three dimensions. Such a solution, dubbed the Freidberg solution, is kin to the axi-symmetric Lundquist solution, but with more general "helical symmetry". The merit of our approach is demonstrated via its application to two case studies of in-situ measured… ▽ More We develop an optimization approach to model the magnetic field configuration of magnetic clouds, based on a linear-force free formulation in three dimensions. Such a solution, dubbed the Freidberg solution, is kin to the axi-symmetric Lundquist solution, but with more general "helical symmetry". The merit of our approach is demonstrated via its application to two case studies of in-situ measured magnetic clouds. Both yield results of reduced $χ^2\approx 1$. Case 1 shows a winding flux rope configuration with one major polarity. Case 2 exhibits a double-helix configuration with two flux bundles winding around each other and rooted on regions of mixed polarities. This study demonstrates the three-dimensional complexity of the magnetic cloud structures. △ Less

Submitted 17 December, 2020; v1 submitted 22 October, 2020; originally announced October 2020.

Comments: Submitted to ApJL on Aug 24, 2020; rejected by Editor without review... Resubmitted to GRL, in press (see DOI; Open Access)

Journal ref: GRL, 2020

Trajectory Determination for Coronal Ejecta Observed by WISPR/Parker Solar Probe

Authors: P. C. Liewer, J. Qiu, P. Penteado, J. R. Hall, A. Vourlidas, R. A. Howard

Abstract: The {\it Wide-field Imager for Solar Probe} (WISPR) on {\it Parker Solar Probe} (PSP), observing in white light, has a fixed angular field of view, extending from 13.5$^{\circ}$ to 108$^{\circ}$ from the Sun and approximately 50$^{\circ}$ in the transverse directions. Because of the highly elliptical orbit of PSP, the physical extent of the imaged coronal region varies directly as the distance fro… ▽ More The {\it Wide-field Imager for Solar Probe} (WISPR) on {\it Parker Solar Probe} (PSP), observing in white light, has a fixed angular field of view, extending from 13.5$^{\circ}$ to 108$^{\circ}$ from the Sun and approximately 50$^{\circ}$ in the transverse directions. Because of the highly elliptical orbit of PSP, the physical extent of the imaged coronal region varies directly as the distance from the Sun, requiring new techniques for analysis of the motions of observed density features. Here, we present a technique for determining the 3D trajectory of CMEs and other coronal ejecta moving radially at a constant velocity by first tracking the motion in a sequence of images and then applying a curve-fitting procedure to determine the trajectory parameters (distance vs. time, velocity, longitude and latitude). To validate the technique, we have determined the trajectory of two CMEs observed by WISPR that were also observed by another white-light imager, either LASCO/C3 or STEREO-A/HI1. The second viewpoint was used to verify the trajectory results from this new technique and help determine its uncertainty. △ Less

Submitted 19 September, 2020; originally announced September 2020.

Comments: 28 pages, 11 figures, accepted by Solar Physics

Solar Flare Energy Partitioning and Transport -- the Gradual Phase (a Heliophysics 2050 White Paper)

Authors: Graham S. Kerr, Meriem Alaoui, Joel C. Allred, Nicholas H. Bian, Brian R. Dennis, A. Gordon Emslie, Lyndsay Fletcher, Silvina Guidoni, Laura A. Hayes, Gordon D. Holman, Hugh S. Hudson, Judith T. Karpen, Adam F. Kowalski, Ryan O. Milligan, Vanessa Polito, Jiong Qiu, Daniel F. Ryan

Abstract: Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous… ▽ More Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare gradual phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, the flare gradual phase persists much longer than predicted so, by 2050, we must identify the characteristics of the significant energy deposition sustaining the gradual phase, and address the fundamental processes of turbulence and non-local heat flux. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: White paper submitted to the Heliophysics 2050 Workshop

Solar Flare Energy Partitioning and Transport -- the Impulsive Phase (a Heliophysics 2050 White Paper)

Authors: Graham S. Kerr, Meriem Alaoui, Joel C. Allred, Nicholas H. Bian, Brian R. Dennis, A. Gordon Emslie, Lyndsay Fletcher, Silvina Guidoni, Laura A. Hayes, Gordon D. Holman, Hugh S. Hudson, Judith T. Karpen, Adam F. Kowalski, Ryan O. Milligan, Vanessa Polito, Jiong Qiu, Daniel F. Ryan

Abstract: Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous… ▽ More Solar flares are a fundamental component of solar eruptive events (SEEs; along with solar energetic particles, SEPs, and coronal mass ejections, CMEs). Flares are the first component of the SEE to impact our atmosphere, which can set the stage for the arrival of the associated SEPs and CME. Magnetic reconnection drives SEEs by restructuring the solar coronal magnetic field, liberating a tremendous amount of energy which is partitioned into various physical manifestations: particle acceleration, mass and magnetic-field eruption, atmospheric heating, and the subsequent emission of radiation as solar flares. To explain and ultimately predict these geoeffective events, the heliophysics community requires a comprehensive understanding of the processes that transform and distribute stored magnetic energy into other forms, including the broadband radiative enhancement that characterises flares. This white paper, submitted to the Heliophysics 2050 Workshop, discusses the flare impulsive phase part of SEEs, setting out the questions that need addressing via a combination of theoretical, modelling, and observational research. In short, by 2050 we must determine the mechanisms of particle acceleration and propagation, and must push beyond the paradigm of energy transport via nonthermal electron beams, to also account for accelerated protons & ions and downward directed Alfven waves. △ Less

Submitted 17 September, 2020; originally announced September 2020.

Comments: White paper submitted to the Heliophysics 2050 Workshop

Critical Science Plan for the Daniel K. Inouye Solar Telescope (DKIST)

Authors: Mark P. Rast, Nazaret Bello González, Luis Bellot Rubio, Wenda Cao, Gianna Cauzzi, Edward DeLuca, Bart De Pontieu, Lyndsay Fletcher, Sarah E. Gibson, Philip G. Judge, Yukio Katsukawa, Maria D. Kazachenko, Elena Khomenko, Enrico Landi, Valentin Martínez Pillet, Gordon J. D. Petrie, Jiong Qiu, Laurel A. Rachmeler, Matthias Rempel, Wolfgang Schmidt, Eamon Scullion, Xudong Sun, Brian T. Welsch, Vincenzo Andretta, Patrick Antolin , et al. (62 additional authors not shown)

Abstract: The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With… ▽ More The Daniel K. Inouye Solar Telescope (DKIST) will revolutionize our ability to measure, understand and model the basic physical processes that control the structure and dynamics of the Sun and its atmosphere. The first-light DKIST images, released publicly on 29 January 2020, only hint at the extraordinary capabilities which will accompany full commissioning of the five facility instruments. With this Critical Science Plan (CSP) we attempt to anticipate some of what those capabilities will enable, providing a snapshot of some of the scientific pursuits that the Daniel K. Inouye Solar Telescope hopes to engage as start-of-operations nears. The work builds on the combined contributions of the DKIST Science Working Group (SWG) and CSP Community members, who generously shared their experiences, plans, knowledge and dreams. Discussion is primarily focused on those issues to which DKIST will uniquely contribute. △ Less

Submitted 20 August, 2020; v1 submitted 18 August, 2020; originally announced August 2020.

How Does Magnetic Reconnection Drive the Early Stage Evolution of Coronal Mass Ejections?

Authors: Chunming Zhu, Jiong Qiu, Paulett Liewer, Angelos Vourlidas, Michael Spiegel, Qiang Hu

Abstract: Theoretically, CME kinematics are related to magnetic reconnection processes in the solar corona. However, the current quantitative understanding of this relationship is based on the analysis of only a handful of events. Here we report a statistical study of 60 CME-flare events from August 2010 to December 2013. We investigate kinematic properties of CMEs and magnetic reconnection in the low coron… ▽ More Theoretically, CME kinematics are related to magnetic reconnection processes in the solar corona. However, the current quantitative understanding of this relationship is based on the analysis of only a handful of events. Here we report a statistical study of 60 CME-flare events from August 2010 to December 2013. We investigate kinematic properties of CMEs and magnetic reconnection in the low corona during the early phase of the eruptions, by combining limb observations from STEREO with simultaneous on-disk views from SDO. For a subset of 42 events with reconnection rate evaluated by the magnetic fluxes swept by the flare ribbons on the solar disk observed from SDO, we find a strong correlation between the peak CME acceleration and the peak reconnection rate. Also, the maximum velocities of relatively fast CMEs (> 600 km/s) are positively correlated with the reconnection flux, but no such correlation is found for slow CMEs. A time-lagged correlation analysis suggests that the distribution of the time lag of CME acceleration relative to reconnection rate exhibits three peaks, approximately 10 minutes apart, and on average, acceleration-lead events have smaller reconnection rates. We further compare the CME total mechanical energy with the estimated energy in the current sheet. The comparison suggests that, for small-flare events, reconnection in the current sheet alone is insufficient to fuel CMEs. Results from this study suggest that flare reconnection may dominate the acceleration of fast CMEs, but for events of slow CMEs and weak reconnection, other mechanisms may be more important. △ Less

Submitted 24 March, 2020; originally announced March 2020.

Comments: 31 pages, 13 figures, 1 table. Accepted for publication by ApJ

λ = 2 mm spectroscopy observations toward the circumnuclear disk of NGC 1068

Authors: Jianjie Qiu, Jiangshui Zhang, Yong Zhang, Lanwei Jia, Xindi Tang

Abstract: Aims. We investigate the physical and chemical conditions of molecular gas in the circumnuclear disk (CND) region of NGC 1068. Methods. We carried out a spectral line survey with the IRAM 30m telescope toward the center of NGC 1068 and mainly focused on the 2 mm band with a frequency coverage of 160.7-168.6 GHz and 176.5-184.3 GHz. Results. Fifteen lines are detected in NGC 1068, eight of which ar… ▽ More Aims. We investigate the physical and chemical conditions of molecular gas in the circumnuclear disk (CND) region of NGC 1068. Methods. We carried out a spectral line survey with the IRAM 30m telescope toward the center of NGC 1068 and mainly focused on the 2 mm band with a frequency coverage of 160.7-168.6 GHz and 176.5-184.3 GHz. Results. Fifteen lines are detected in NGC 1068, eight of which are new detections for this galaxy. We derive the rotation temperatures and column densities of fourteen molecular species. Conclusions. Based on the [HCO+ (2 - 1)]/[HOC+ (2 - 1)] ratio, we obtain a high ionization degree in the CND of NGC 1068. It is found that HC3N is concentrated in the east knot, while 13CCH, CH3CN, SO, HOC+, CS, CH3CCH, and H2CO are concentrated in the west knot. Compared to the star-forming galaxies M 82 and NGC 253, the chemistry of NGC 1068 might be less strongly affected by the UV radiation field, and its kinetic temperature might be lower. △ Less

Submitted 23 January, 2020; originally announced January 2020.

Comments: 28 pages, 10 figures

Journal ref: A&A 634, A125 (2020)

Quantifying the Toroidal Flux of Pre-existing Flux Ropes of CMEs

Authors: C. Xing, X. Cheng, Jiong Qiu, Qiang Hu, E. R. Priest, M. D. Ding

Abstract: In past decades, much progress has been achieved on the origin and evolution of coronal mass ejections (CMEs). In-situ observations of the counterparts of CMEs, especially magnetic clouds (MCs) near the Earth, have provided measurements of the structure and total flux of CME flux ropes. However, it has been difficult to measure these properties in the erupting CME flux rope, in particular in the p… ▽ More In past decades, much progress has been achieved on the origin and evolution of coronal mass ejections (CMEs). In-situ observations of the counterparts of CMEs, especially magnetic clouds (MCs) near the Earth, have provided measurements of the structure and total flux of CME flux ropes. However, it has been difficult to measure these properties in the erupting CME flux rope, in particular in the pre-existing flux rope. In this work, we propose a model to estimate the toroidal flux of the pre-existing flux rope by subtracting the flux contributed by magnetic reconnection during the eruption from the flux measured in the MC. The flux by the reconnection is derived from geometric properties of two-ribbon flares based on a quasi-2D reconnection model. We then apply the model to four CME/flare events and find that the ratio of toroidal flux in the pre-existing flux rope to that of the associated MC lies in the range of 0.40--0.88. It indicates that the toroidal flux of the pre-existing flux rope has an important contribution to that of the CME flux rope and is usually at least as large as the flux arising from the eruption process for the selected events. △ Less

Submitted 23 December, 2019; originally announced December 2019.

Comments: 10 pages, 9 figures, 1 table; accepted for publication in ApJ

Evolution of A Magnetic Flux Rope toward Eruption

Authors: Wensi Wang, Chunming Zhu, Jiong Qiu, Rui Liu, Kai E. Yang, Qiang Hu

Abstract: It is well accepted that a magnetic flux rope (MFR) is a critical component of many coronal mass ejections (CMEs), yet how it evolves toward eruption remains unclear. Here we investigate the continuous evolution of a pre-existing MFR, which is rooted in strong photospheric magnetic fields and electric currents. The evolution of the MFR is observed by the Solar Terrestrial Relations Observatory (ST… ▽ More It is well accepted that a magnetic flux rope (MFR) is a critical component of many coronal mass ejections (CMEs), yet how it evolves toward eruption remains unclear. Here we investigate the continuous evolution of a pre-existing MFR, which is rooted in strong photospheric magnetic fields and electric currents. The evolution of the MFR is observed by the Solar Terrestrial Relations Observatory (STEREO) and the Solar Dynamics Observatory (SDO) from multiple viewpoints. From STEREO's perspective, the MFR starts to rise slowly above the limb five hours before it erupts as a halo CME on 2012 June 14. In SDO observations, conjugate dimmings develop on the disk, simultaneously with the gradual expansion of the MFR, suggesting that the dimmings map the MFR's feet. The evolution comprises a two-stage gradual expansion followed by another stage of rapid acceleration/eruption. Quantitative measurements indicate that magnetic twist of the MFR increases from 1.0 +/- 0.5 to 2.0 +/- 0.5 turns during the five-hour expansion, and further increases to about 4.0 turns per AU when detected as a magnetic cloud at 1 AU two day later. In addition, each stage is preceded by flare(s), implying reconnection is actively involved in the evolution and eruption of the MFR. The implications of these measurements on the CME initiation mechanisms are discussed. △ Less

Submitted 9 December, 2018; originally announced December 2018.

Comments: Accepted by The Astrophysical Journal

Supermassive Black Holes with High Accretion Rates in Active Galactic Nuclei. VIII. Structure of the Broad-Line Region and Mass of the Central Black Hole in Mrk 142

Authors: Yan-Rong Li, Yu-Yang Songsheng, Jie Qiu, Chen Hu, Pu Du, Kai-Xing Lu, Ying-Ke Huang, Jin-Ming Bai, Wei-Hao Bian, Ye-Fei Yuan, Luis C. Ho, Jian-Min Wang

Abstract: This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring se… ▽ More This is the eighth in a series of papers reporting on a large reverberation mapping campaign to measure black hole (BH) mass in high accretion rate active galactic nuclei (AGNs). We employ the recently developed dynamical modeling approach for broad-line regions (BLRs) based on the method of Pancoast et al. to analyze the reverberation mapping dataset of Mrk 142 observed in the first monitoring season. In this approach, continuum variations are reconstructed using a damped random walk process, and BLR structure is delineated using a flexible disk-like geometry, in which BLR clouds move around the central BH with Keplerian orbits or inflow/outflow motion. The approach also includes the possibilities of anisotropic emission of BLR clouds, non-linear response of the line emission to the continuum, and different long-term trends in the continuum and emission-line variations. We implement the approach in a Bayesian framework that is apt for parallel computation and use a Markov Chain Monte Carlo technique to recover the parameters and uncertainties for the modeling, including mass of the central BH. We apply three BLR models with different prescriptions of BLR clouds distributions and find that the best model for fitting the data of Mrk 142 is a two-zone BLR model, consistent with the theoretical BLR model surrounding slim accretion disks. The best model yields a BH mass of $\log (M_\bullet/M_\odot)=6.23_{-0.45}^{+0.26}$, resulting in a virial factor of $\log f=-0.36_{-0.54}^{+0.33}$ for the full width at half maximum of the H$β$ line measured from the mean spectrum. The virial factors for the other measures of the H$β$ line width are also presented. △ Less

Submitted 15 November, 2018; originally announced November 2018.

Comments: 21 pages, 12 figures, 5 tables; ApJ accepted; The software BRAINS available at https://github.com/LiyrAstroph/BRAINS

Supermassive black holes with high accretion rates in active galactic nuclei. VII. Reconstruction of velocity-delay maps by maximum entropy method

Authors: Ming Xiao, Pu Du, Keith D. Horne, Chen Hu, Yan-Rong Li, Ying-Ke Huang, Kai-Xing Lu, Jie Qiu, Fang Wang, Jin-Ming Bai, Weihao Bian, Luis C. Ho, Ye-Fei Yuan, Jian-Min Wang

Abstract: As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their Hβ and Hγ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of Hβ and Hγ emission lines of Mrk 335 indi… ▽ More As one of the series of papers reporting on a large reverberation mapping campaign, we apply the maximum entropy method (MEM) to 9 narrow-line Seyfert 1 galaxies with super-Eddington accretion rates observed during 2012-2013 for the velocity-delay maps of their Hβ and Hγ emission lines. The maps of 6 objects are reliably reconstructed using MEM. The maps of Hβ and Hγ emission lines of Mrk 335 indicate that the gas of its broad-line region (BLR) is infalling. For Mrk 142, its Hβ and Hγ lines show signatures of outflow. The Hβ and Hγ maps of Mrk 1044 demonstrate complex kinematics -- a virialized motion accompanied by an outflow signature, and the Hβ map of IRAS F12397+3333 is consistent with a disk or a spherical shell. The Hβ maps of Mrk 486 and MCG +06-26-012 suggest the presence of an inflow and outflow, respectively. These super-Eddington accretors show diverse geometry and kinematics. Brief discussions of their BLRs are provided for each individual object. △ Less

Submitted 2 August, 2018; originally announced August 2018.

Comments: 15 pages, 18 figures, 2 tables, accepted for publication in The Astrophysical Journal

The MALATANG Survey: the L_gas-L_IR correlation on sub-kiloparsec scale in six nearby star-forming galaxies as traced by HCN J=4-3 and HCO^+ J=4-3

Authors: Qing-Hua Tan, Yu Gao, Zhi-Yu Zhang, Thomas R. Greve, Xue-Jian Jiang, Christine D. Wilson, Chen-Tao Yang, Ashley Bemis, Aeree Chung, Satoki Matsushita, Yong Shi, Yi-Ping Ao, Elias Brinks, Malcolm J. Currie, Timothy A. Davis, Richard de Grijs, Luis C. Ho, Masatoshi Imanishi, Kotaro Kohno, Bumhyun Lee, Harriet Parsons, Mark G. Rawlings, Dimitra Rigopoulou, Erik Rosolowsky, Joanna Bulger , et al. (21 additional authors not shown)

Abstract: We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where… ▽ More We present HCN J=4-3 and HCO^+ J=4-3 maps of six nearby star-forming galaxies, NGC 253, NGC 1068, IC 342, M82, M83, and NGC 6946, obtained with the James Clerk Maxwell Telescope as part of the MALATANG survey. All galaxies were mapped in the central 2 arcmin $\times$ 2 arcmin region at 14 arcsec (FWHM) resolution (corresponding to linear scales of ~ 0.2-1.0 kpc). The L_IR-L'_dense relation, where the dense gas is traced by the HCN J=4-3 and the HCO^+ J=4-3 emission, measured in our sample of spatially-resolved galaxies is found to follow the linear correlation established globally in galaxies within the scatter. We find that the luminosity ratio, L_IR/L'_dense, shows systematic variations with L_IR within individual spatially resolved galaxies, whereas the galaxy-integrated ratios vary little. A rising trend is also found between L_IR/L'_dense ratio and the warm-dust temperature gauged by the 70 μm/100 μm flux ratio. We find the luminosity ratios of IR/HCN(4-3) and IR/HCO^+(4-3), which can be taken as a proxy for the efficiency of star formation in the dense molecular gas (SFE_dense), appears to be nearly independent of the dense-gas fraction (f_dense) for our sample of galaxies. The SFE of the total molecular gas (SFE_mol) is found to increase substantially with f_dense when combining our data with that on local (ultra)luminous infrared galaxies and high-z quasars. The mean L'_HCN(4-3)/L'_HCO^+(4-3) line ratio measured for the six targeted galaxies is 0.9+/-0.6. No significant correlation is found for the L'_HCN(4-3)/L'_HCO^+(4-3) ratio with the SFR as traced by L_IR, nor with the warm-dust temperature, for the different populations of galaxies. △ Less

Submitted 15 May, 2018; originally announced May 2018.

Comments: 29 pages, 12 figures, accepted for publication in ApJ

Two-Phase Heating in Flaring Loops

Authors: Chunming Zhu, Jiong Qiu, Dana W. Longcope

Abstract: We analyze and model a C5.7 two-ribbon solar flare observed by SDO, Hinode and GOES on 2011 December 26. The flare is made of many loops formed and heated successively over one and half hours, and their footpoints are brightened in the UV 1600 A before enhanced soft X-ray and EUV missions are observed in flare loops. Assuming that anchored at each brightened UV pixel is a half flaring loop, we ide… ▽ More We analyze and model a C5.7 two-ribbon solar flare observed by SDO, Hinode and GOES on 2011 December 26. The flare is made of many loops formed and heated successively over one and half hours, and their footpoints are brightened in the UV 1600 A before enhanced soft X-ray and EUV missions are observed in flare loops. Assuming that anchored at each brightened UV pixel is a half flaring loop, we identify more than 6,700 half flaring loops, and infer the heating rate of each loop from the UV light curve at the foot-point. In each half loop, the heating rate consists of two phases, an intense impulsive heating followed by a low-rate heating persistent for more than 20 minutes. Using these heating rates, we simulate the evolution of their coronal temperatures and densities with the model of "enthalpy-based thermal evolution of loops" (EBTEL). In the model, suppression of thermal conduction is also considered. This model successfully reproduces total soft X-ray and EUV light curves observed in fifteen pass-bands by four instruments GOES, AIA, XRT, and EVE. In this flare, a total energy of 4.9x10^30 ergs is required to heat the corona, around 40% of this energy is in the slow heating phase. About two fifth of the total energy used to heat the corona is radiated by the coronal plasmas, and the other three fifth transported to the lower atmosphere by thermal conduction. △ Less

Submitted 2 February, 2018; originally announced February 2018.

Comments: 18 pages, 10 figures, accepted by ApJ

Adaptive inversion algorithm for 1.5 um visibility lidar incorporating in situ Angstrom wavelength exponent

Authors: Xiang Shang, Haiyun Xia, Xiankang Dou, Mingjia Shangguan, Manyi Li, Chong Wang, Jiawei Qiu, Lijie Zhao, Shengfu Lin

Abstract: As one of the most popular applications of lidar systems, the atmospheric visibility is defined to be inversely proportional to the atmospheric extinction coefficient at 0.55 um. Since the laser at 1.5 um shows the highest maximum permissible exposure in the wavelength ranging from 0.3 um to 10 um, the eye-safe 1.5 um lidar can be deployed in urban areas. In such a case, the measured extinction co… ▽ More As one of the most popular applications of lidar systems, the atmospheric visibility is defined to be inversely proportional to the atmospheric extinction coefficient at 0.55 um. Since the laser at 1.5 um shows the highest maximum permissible exposure in the wavelength ranging from 0.3 um to 10 um, the eye-safe 1.5 um lidar can be deployed in urban areas. In such a case, the measured extinction coefficient at 1.5 um should be converted to that at 0.55 um for visibility retrieval. Although several models have been established since 1962, the accurate wavelength conversion remains a challenge. An adaptive inversion algorithm for 1.5 um visibility lidar is proposed and demonstrated by using the in situ Angstrom wavelength exponent, which is derived from either a sun photometer or an aerosol spectrometer. The impact of the particle size distribution of atmospheric aerosols and the Rayleigh backscattering of atmospheric molecules are taken into account. In comparison, the Angstrom wavelength exponent derived from the sun photometer is 7.7% higher than that derived from the aerosol spectrometer. Then, using 1.5 um visibility lidar, the visibility with a temporal resolution of 5 min is detected over 48 hours in Hefei (31.83 N, 117.25 E). The average visibility error between the new method and a visibility sensor (Vaisala, PWD52) is 5.2% with the R-square value of 0.96, while the relative error between another reference visibility lidar at 532 nm and the visibility sensor is 6.7% with the R-square value of 0.91. All results agree with each other well, demonstrating the accuracy and stability of the algorithm. △ Less

Submitted 19 November, 2017; originally announced November 2017.

Comments: 9 pages, 7 figures

MSC Class: 86-05

Comparison of Helioseismic Far-side Active Region Detections with STEREO Far-Side EUV Observations of Solar Activity

Authors: P. C. Liewer, J. Qiu, C. Lindsey

Abstract: Seismic maps of the Sun's far hemisphere, computed from Doppler data from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) are now being used routinely to detect strong magnetic regions on the far side of the Sun (http:/jsoc.stanford.edu/data/farside/). To test the reliability of this technique, the helioseismically inferred active region detections are comp… ▽ More Seismic maps of the Sun's far hemisphere, computed from Doppler data from the Helioseismic and Magnetic Imager (HMI) on board the Solar Dynamics Observatory (SDO) are now being used routinely to detect strong magnetic regions on the far side of the Sun (http:/jsoc.stanford.edu/data/farside/). To test the reliability of this technique, the helioseismically inferred active region detections are compared with far-side observation of solar activity from the Solar TErrestrial RElations Observatory (STEREO), using brightness in extreme ultraviolet light (EUV) as a proxy for magnetic fields. Two approaches are used to analyze nine months of STEREO and HMI data. In the first approach, we determine whether or not new large east-limb active regions are detected seismically on the far side before they appear Earth side and study how the detectability of these regions relates to their EUV intensity. We find that, while there is a range of EUV intensities for which far-side regions may or may not be detected seismically, there appears to be an intensity level above which they are almost always detected and an intensity level below which they are never detected. In the second approach, we analyze concurrent extreme ultraviolet and helioseismic far-side observations. We find that 100% (22) of the far-side seismic regions correspond to an extreme ultraviolet plage; 95% of these either became a NOAA-designated magnetic region when reaching the east limb or were one before crossing to the far side. A low but significant correlation is found between the seismic signature strength and the EUV intensity of a farside region. △ Less

Submitted 22 September, 2017; originally announced September 2017.

Quasi-periodic oscillations in flares and coronal mass ejections associated with magnetic reconnection

Authors: Takuya Takahashi, Jiong Qiu, Kazunari Shibata

Abstract: We propose a mechanism for quasi-periodic oscillations of both coronal mass ejections (CMEs) and flare loops as related to magnetic reconnection in eruptive solar flares. We perform two-dimensional numerical MHD simulations of magnetic flux rope eruption, with three different values of the global Lundquist number. In the low Lundquist number run, no oscillatory behavior is found. In the moderate L… ▽ More We propose a mechanism for quasi-periodic oscillations of both coronal mass ejections (CMEs) and flare loops as related to magnetic reconnection in eruptive solar flares. We perform two-dimensional numerical MHD simulations of magnetic flux rope eruption, with three different values of the global Lundquist number. In the low Lundquist number run, no oscillatory behavior is found. In the moderate Lunquist number run, on the other hand, quasi-periodic oscillations are excited both at the bottom of the flux rope and at the flare loop-top. In the high Lundquist number run, quasi-periodic oscillations are also excited; in the meanwhile, the dynamics become turbulent due to the formation of multiple plasmoids in the reconnection current sheet. In high and moderate Lundquist number runs, thin reconnection jet collide with the flux rope bottom or flare loop-top and dig them deeply. Steep oblique shocks are formed as termination shocks where reconnection jet is bent (rather than decelerated) in horizontal direction, resulting in supersonic back-flows. The structure becomes unstable, and quasi-periodic oscillation of supersonic back-flows appear at locally confined high-beta region at both the flux rope bottom and flare loop-top. We compare the observational characteristics of quasi-periodic oscillations in erupting flux ropes, post-CME current sheets, flare ribbons and light curves, with corresponding dynamical structures found in our simulation. △ Less

Submitted 15 September, 2017; originally announced September 2017.

Comments: 30 pages, 20 figures, accepted for publication in ApJ

Spectroscopic Observations of Magnetic Reconnection and Chromospheric Evaporation in an X-shaped Solar Flare

Authors: Y. Li, M. Kelly, M. D. Ding, J. Qiu, X. S. Zhu, W. Q. Gan

Abstract: We present observations of distinct UV spectral properties at different locations during an atypical X-shaped flare (SOL2014-11-09T15:32) observed by the Interface Region Imaging Spectrograph (IRIS). In this flare, four chromospheric ribbons appear and converge at an X-point where a separator is anchored. Above the X-point, two sets of non-coplanar coronal loops approach laterally and reconnect at… ▽ More We present observations of distinct UV spectral properties at different locations during an atypical X-shaped flare (SOL2014-11-09T15:32) observed by the Interface Region Imaging Spectrograph (IRIS). In this flare, four chromospheric ribbons appear and converge at an X-point where a separator is anchored. Above the X-point, two sets of non-coplanar coronal loops approach laterally and reconnect at the separator. The IRIS slit was located close to the X-point, cutting across some of the flare ribbons and loops. Near the location of the separator, the Si IV 1402.77 A line exhibits significantly broadened line wings extending to 200 km/s but an unshifted line core. These spectral features suggest the presence of bidirectional flows possibly related to the separator reconnection. While at the flare ribbons, the hot Fe XXI 1354.08 A line shows blueshifts and the cool Si IV 1402.77 A, C II 1335.71 A, and Mg II 2803.52 A lines show evident redshifts up to a velocity of 80 km/s, which are consistent with the scenario of chromospheric evaporation/condensation. △ Less

Submitted 29 August, 2017; originally announced August 2017.

Comments: ApJ accepted. Comments are welcome

The dependence of the IR-radio correlation on the metallicity

Authors: Jianjie Qiu, Yong Shi, Junzhi Wang, Zhi-Yu Zhang, Luwenjia Zhou

Abstract: We have compiled a sample of 26 metal-poor galaxies with 12 + log(O/H) < 8.1 with both infrared continuum and 1.4 GHz radio continuum data. By comparing to galaxies at higher metallicity, we have investigated the dependence on the metallicity of the IR-radio relationship at 24 um, 70 um, 100 um and 160 um bands as well as the integrated FIR luminosity. It is found that metal-poor galaxies have on… ▽ More We have compiled a sample of 26 metal-poor galaxies with 12 + log(O/H) < 8.1 with both infrared continuum and 1.4 GHz radio continuum data. By comparing to galaxies at higher metallicity, we have investigated the dependence on the metallicity of the IR-radio relationship at 24 um, 70 um, 100 um and 160 um bands as well as the integrated FIR luminosity. It is found that metal-poor galaxies have on average lower qIR than metal-rich ones with larger offsets at longer IR wavelengths, from -0.06 dex in q24um to -0.6 dex in q160um. The qIR of all galaxies as a whole at 160 um show positive trends with the metallicity and IR-to-FUV ratio, and negative trends with the IR color, while those at lower IR wavelengths show weaker correlations. We proposed a mechanism that invokes combined effects of low obscured-SFR/total-SFR fraction and warm dust temperature at low metallicity to interpret the above behavior of qIR, with the former reducing the IR radiation and the latter further reducing the IR emission at longer IR wavelength. Other mechanisms that are related to the radio emission including the enhanced magnetic field strength and increased thermal radio contribution are unable to reconcile the IR-wavelength-dependent differences of qIR between metal-poor and metal- rich galaxies. In contrast to qIR, the mean total-SFR/radio ratio of metal-poor galaxies is the same as the metal-rich one, indicating the 1.4 GHz radio emission is still an effective tracer of SFRs at low metallicity. △ Less

Submitted 8 August, 2017; originally announced August 2017.

Comments: 25 pages, 11 figures, 4 tables. ApJ in press

Gradual Solar Coronal Dimming and Evolution of Coronal Mass Ejection in the Early Phase

Authors: Jiong Qiu, Jianxia Cheng

Abstract: We report observations of a two-stage coronal dimming in an eruptive event of a two-ribbon flare and a fast coronal mass ejection (CME). Weak gradual dimming persists for more than half an hour before the onset of the two-ribbon flare and the fast rise of the CME. It is followed by abrupt rapid dimming. The two-stage dimming occurs in a pair of conjugate dimming regions adjacent to the two flare r… ▽ More We report observations of a two-stage coronal dimming in an eruptive event of a two-ribbon flare and a fast coronal mass ejection (CME). Weak gradual dimming persists for more than half an hour before the onset of the two-ribbon flare and the fast rise of the CME. It is followed by abrupt rapid dimming. The two-stage dimming occurs in a pair of conjugate dimming regions adjacent to the two flare ribbons, and the flare onset marks the transition between the two stages of dimming. At the onset of the two-ribbon flare, transient brightenings are also observed inside the dimming regions, before rapid dimming occurs at the same places. These observations suggest that the CME structure, most probably anchored at the twin dimming regions, undergoes a slow rise before the flare onset, and its kinematic evolution has significantly changed at the onset of flare reconnection. We explore diagnostics of the CME evolution in the early phase with analysis of the gradual dimming signatures prior to the CME eruption. △ Less

Submitted 8 July, 2017; originally announced July 2017.

Journal ref: The Astrophysical Journal Letters, Volume 838, Issue 1, article id. L6, 6 pp. (2017)

Elongation of Flare Ribbons

Authors: Jiong Qiu, Dana W. Longcope, Paul A. Cassak, Eric R. Priest

Abstract: We present an analysis of the apparent elongation motion of flare ribbons along the polarity inversion line (PIL) as well as the shear of flare loops in several two-ribbon flares. Flare ribbons and loops spread along the PIL at a speed ranging from a few to a hundred km/s. The shear measured from conjugate foot-points is consistent with the measurement from flare loops, and both show the decrease… ▽ More We present an analysis of the apparent elongation motion of flare ribbons along the polarity inversion line (PIL) as well as the shear of flare loops in several two-ribbon flares. Flare ribbons and loops spread along the PIL at a speed ranging from a few to a hundred km/s. The shear measured from conjugate foot-points is consistent with the measurement from flare loops, and both show the decrease of shear toward a potential field as a flare evolves and ribbons and loops spread along the PIL. Flares exhibiting fast bi-directional elongation appear to have a strong shear, which may indicate a large magnetic guide field relative to the reconnection field in the coronal current sheet. We discuss how the analysis of ribbon motion could help infer properties in the corona where reconnection takes place. △ Less

Submitted 8 July, 2017; originally announced July 2017.

Journal ref: The Astrophysical Journal, Volume 838, Issue 1, article id. 17, 17 pp. (2017)

Imaging Observations of Magnetic Reconnection in a Solar Eruptive Flare

Authors: Y. Li, X. Sun, M. D. Ding, J. Qiu, E. R. Priest

Abstract: Solar flares are one of the most energetic events in the solar atmosphere. It is widely accepted that flares are powered by magnetic reconnection in the corona. An eruptive flare is usually accompanied by a coronal mass ejection, both of which are probably driven by the eruption of a magnetic flux rope (MFR). Here we report an eruptive flare on 2016 March 23 observed by the Atmospheric Imaging Ass… ▽ More Solar flares are one of the most energetic events in the solar atmosphere. It is widely accepted that flares are powered by magnetic reconnection in the corona. An eruptive flare is usually accompanied by a coronal mass ejection, both of which are probably driven by the eruption of a magnetic flux rope (MFR). Here we report an eruptive flare on 2016 March 23 observed by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory. The extreme-ultraviolet imaging observations exhibit the clear rise and eruption of an MFR. In particular, the observations reveal solid evidence for magnetic reconnection from both the corona and chromosphere during the flare. Moreover, weak reconnection is observed before the start of the flare. We find that the preflare weak reconnection is of tether-cutting type and helps the MFR to rise slowly. Induced by a further rise of the MFR, strong reconnection occurs in the rise phases of the flare, which is temporally related to the MFR eruption. We also find that the magnetic reconnection is more of 3D-type in the early phase, as manifested in a strong-to-weak shear transition in flare loops, and becomes more 2D-like in the later phase, as shown by the apparent rising motion of an arcade of flare loops. △ Less

Submitted 30 December, 2016; originally announced December 2016.

Comments: ApJ accepted. Comments are welcome

Coronal Holes and Open Magnetic Flux over Cycles 23 and 24

Authors: Chris Lowder, Jiong Qiu, Robert Leamon

Abstract: As the observational signature of the footprints of solar magnetic field lines open into the heliosphere, coronal holes provide a critical measure of the structure and evolution of these lines. Using a combination of Solar and Heliospheric Observatory / Extreme ultraviolet Imaging Telescope (SOHO/EIT), Solar Dynamics Observatory / Atmospheric Imaging Assembly (SDO/AIA) and Solar Terrestrial Relati… ▽ More As the observational signature of the footprints of solar magnetic field lines open into the heliosphere, coronal holes provide a critical measure of the structure and evolution of these lines. Using a combination of Solar and Heliospheric Observatory / Extreme ultraviolet Imaging Telescope (SOHO/EIT), Solar Dynamics Observatory / Atmospheric Imaging Assembly (SDO/AIA) and Solar Terrestrial Relations Observatory / Extreme Ultraviolet Imager (STEREO/EUVI A/B) extreme ultraviolet (EUV) observations spanning 1996-2015 (nearly two solar cycles), coronal holes are automatically detected and characterized. Coronal hole area distributions show distinct behavior in latitude, defining the domain of polar and low-latitude coronal holes. The northern and southern polar regions show a clear asymmetry, with a lag between hemispheres in the appearance and disappearance of polar coronal holes. △ Less

Submitted 22 December, 2016; originally announced December 2016.

Comments: Accepted for publication in Solar Physics