sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Zhu, Ling

Normalized to: Zhu, L.

86 article(s) in total. 574 co-authors, from 1 to 29 common article(s). Median position in authors list is 4,0.

[1]  oai:arXiv.org:2007.00176  [pdf] - 2126337
The JCMT BISTRO Survey: Magnetic fields associated with a network of filaments in NGC 1333
Doi, Yasuo; Hasegawa, Tetsuo; Furuya, Ray S.; Coud'e, Simon; Hull, Charles L.; Arzoumanian, Doris; Bastien, Pierre; Chen, Mr. Michael C.; Francesco, James; Friesen, Rachel K.; Houde, Martin; Inutsuka, Shu-ichiro; Mairs, Steve; Matsumura, Masafumi; Onaka, Takashi; Sadavoy, Sarah I.; Shimajiri, Yoshito; Tahani, Mehrnoosh; Tomisaka, Kohji; Eswaraiah, Chakali; Koch, Patrick M.; Pattle, Kate; Lee, Chang; Tamura, Motohide; Berry, David; Ching, Tao-Chung; Hwang, Ms. Jihye; Kwon, Woojin; Soam, Archana; Wang, Mr. Jia-Wei; Lai, Shih-Ping; Qiu, Keping; Ward-Thompson, Derek; Byun, Do-Young; Chen, Huei-Ru V.; Chen, Wen; Chen, Zhiwei; Cho, Jungyeon; Choi, Minho; Choi, Yunhee; Chrysostomou, Antonio; Chung, Eun; Diep, Pham; Duan, Hao-Yuan; Fanciullo, Lapo; Fiege, Jason; Franzmann, Ms. Erica; Friberg, Per; Fuller, Gary A.; Gledhill, Tim; Graves, Sarah F.; Greaves, Jane; Griffin, Matt J.; Gu, Mr. Qilao; Han, Ilseung; Hatchell, Jennifer; Hayashi, Saeko; Hoang, Thiem; Inoue, Tsuyoshi; Iwasaki, Kazunari; Jeong, Il-Gyo; Johnstone, Doug; Kanamori, Yoshihiro; Kang, Ji-hyun; Kang, Miju; Kang, Sung-Ju; Kataoka, Akimasa; Kawabata, Koji S.; Kemper, Francisca; Kim, Gwanjeong; Kim, Jongsoo; Kim, Kee-Tae; Kim, Kyoung; Kim, Mi-Ryang; Kim, Mr. Shinyoung; Kirk, Jason; Kobayashi, Mr. Masato I.; Konyves, Vera; Kusune, Mr. Takayoshi; Kwon, Jungmi; Lacaille, Mr. Kevin M.; Law, Chi-Yan; Lee, Chin-Fei; Lee, Miss Hyeseung; Lee, Jeong-Eun; Lee, Sang-Sung; Lee, Yong-Hee; Li, Dalei; Li, Di; Li, Hua-Bai; Liu, Hong-Li; Liu, Mr. Junhao; Liu, Sheng-Yuan; Liu, Tie; Looze, Ilse D.; Lyo, A-Ran; Matthews, Brenda; Moriarty-Schieven, Gerald H.; Nagata, Tetsuya; Nakamura, Fumitaka; Nakanishi, Hiroyuki; Ohashi, Nagayoshi; Park, Geumsook; Parsons, Harriet A.; Peretto, Nicolas; Pyo, Tae-Soo; Qian, Lei; Rao, Ramprasad; Rawlings, Mark; Retter, Brendan; Richer, John; Rigby, Andrew; Saito, Hiro; Savini, Giorgio; Scaife, Anna M.; Seta, Masumichi; Shinnaga, Hiroko; Tang, Ya-Wen; Tsukamoto, Yusuke; Viti, Serena; Wang, Hongchi; Whitworth, Anthony P.; Yen, Hsi-Wei; Yoo, Hyunju; Yuan, Jinghua; Yun, Hyeong-Sik; Zenko, Tetsuya; Zhang, Chuan-Peng; Zhang, Guoyin; Zhang, Ya-Peng; Zhou, Jianjun; Zhu, Lei; Andr'e, Philippe; Dowell, C. D.; Eyres, Stewart P.; Falle, Sam; Loo, Sven; Robitaille, Jean-Franois
Comments: 33 pages, 21 figures, accepted for publication in ApJ
Submitted: 2020-06-30
We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.
[2]  oai:arXiv.org:2006.13410  [pdf] - 2121547
Filament Intersections and Cold Dense Cores in Orion A North
Comments: 20 pages, 10 figures, 2 tables. Accepted to MNRAS
Submitted: 2020-06-23
We studied the filament structures and dense cores in OMC-2,3 region in Orion A North molecular cloud using the high-resolution N2H+ (1-0) spectral cube observed with the Atacama Large Millimeter/Submillimeter Array (ALMA). The filament network over a total length of 2 pc is found to contain 170 intersections and 128 candidate dense cores. The dense cores are all displaced from the infrared point sources (possible young stars), and the major fraction of cores (103) are located around the intersections. Towards the intersections, there is also an increasing trend for the total column density Ntot as well as the the power-law index of the column-density Probability Distribution Function (N-PDF), suggesting that the intersections would in general have more significant gas assembly than the other part of the filament paths. The virial analysis shows that the dense cores mostly have virial mass ratio of alpha_vir=M_vir/M_gas<1.0, suggesting them to be bounded by the self gravity. In the mean time, only about 23 percent of the cores have critical mass ratio of alpha_crit=M_crit/M_gas<1.0, suggesting them to be unstable against core collapse. Combining these results, it shows that the major fraction of the cold starless and possible prestellar cores in OMC-2,3 are being assembled around the intersections, and currently in a gravitationally bound state. But more extensive core collapse and star formation may still require continuous core-mass growth or other perturbatio
[3]  oai:arXiv.org:2006.04168  [pdf] - 2112023
Resolution-dependent Subsonic Non-thermal Line Dispersion Revealed by ALMA
Comments: 13 pages, 6 figures, accepted by RAA
Submitted: 2020-06-07, last modified: 2020-06-09
We report here Atacama Large Millimeter/submillimeter Array (ALMA) N$_2$H$^+$ (1-0) images of the Orion Molecular Cloud 2 and 3 (OMC-2/3) with high angular resolution (3'' or 1200 au) and high spatial dynamic range. Combining dataset from the ALMA main array, ALMA Compact Array (ACA), the Nobeyama 45m Telescope, and the JVLA (providing temperature measurement on matching scales), we find that most of the dense gas in OMC-2/3 is subsonic ($\rm \sigma_{NT}/c_{s}$ = 0.62) with a mean line width ($\Delta\upsilon$) of 0.39 km s$^{-1}$ FWHM. This is markedly different from the majority of previous observations of massive star-forming regions. In contrast, line widths from the Nobeyama Telescope are transonic at 0.69 km s$^{-1}$ ($\rm \sigma_{NT}/c_{s}$ = 1.08). We demonstrated that the larger line widths obtained by the single-dish telescope arose from unresolved sub-structures within their respective beams. The dispersions from larger scales $\sigma_{ls}$ (as traced by the Nobeyama Telescope) can be decomposed into three components $\rm \sigma_{ls}^2 = \sigma_{ss}^2+ \sigma_{bm}^2+ \sigma_{rd}^2$, where small-scale $\sigma_{ss}$ is the line dispersion of each ALMA beam, bulk motion $\sigma_{bm}$ is dispersion between peak velocity of each ALMA beam, and $\sigma_{rd}$ is the residual dispersion. Such decomposition, though purely empirical, appears to be robust throughout our data cubes. Apparent supersonic line widths, commonly found in massive molecular clouds, are thus likely due to the effect of poor spatial resolution. The observed non-thermal line dispersion (sometimes referred to as 'turbulence') transits from supersonic to subsonic at $\sim 0.05$ pc scales in OMC-2/3 region. Such transition could be commonly found with sufficient spatial (not just angular) resolution, even in regions with massive young clusters, such as Orion molecular clouds studied here.
[4]  oai:arXiv.org:2006.01147  [pdf] - 2105809
ALMA [N \i\i ] 205 \mu m Imaging Spectroscopy of the Lensed Submillimeter galaxy ID 141 at redshift 4.24
Comments: 13 pages, 6 figures, 2 tables, accepted by ApJ, lensing model code can be found here https://gitlab.com/cxylzlx/tiny_lens
Submitted: 2020-06-01
We present the Atacama Large Millimeter/submillimeter Array (ALMA) observation of the Sub-millimeter galaxy (SMG) ID 141 at z=4.24 in the [N II] 205 $\mu$m line (hereafter [N II]) and the underlying continuum at (rest-frame) 197.6 $\mu$m. Benefiting from lensing magnification by a galaxy pair at z=0.595, ID 141 is one of the brightest z$>4$ SMGs. At the angular resolutions of $\sim1.2''$ to $1.5''$ ($1'' \sim6.9$ kpc), our observation clearly separates, and moderately resolves the two lensed images in both continuum and line emission at $\rm S/N>5$ . Our continuum-based lensing model implies an averaged amplification factor of $\sim5.8$ and reveals that the de-lensed continuum image has the S\'ersic index $\simeq 0.95$ and the S\'ersic radius of $\sim0.18'' (\sim 1.24$ kpc). Furthermore, the reconstructed [N II] velocity field in the source plane is dominated by a rotation component with a maximum velocity of $\sim 300$ km/s at large radii, indicating a dark matter halo mass of $\sim 10^{12}M_{\odot}$. This, together with the reconstructed velocity dispersion field being smooth and modest in value ($<100$ km/s) over much of the outer parts of the galaxy, favours the interpretation of ID 141 being a disk galaxy dynamically supported by rotation. The observed [N II]/CO (7-6) and [N II]/[C II] 158 $\mu$m line luminosity ratios, which are consistent with the corresponding line ratio vs. far-infrared color correlation from local luminous infrared galaxies, imply a de-lensed star formation rate of ($1.8\pm 0.6)\times10^3M_\odot$/yr and provide an independent estimate on the size of the star-forming region $0.7^{+0.3}_{-0.3}$ kpc in radius.
[5]  oai:arXiv.org:2005.09410  [pdf] - 2097533
A discrete chemo-dynamical model of M87's globular clusters: Kinematics extending to ~ 400 kpc
Comments:
Submitted: 2020-05-19
We study the mass distribution and kinematics of the giant elliptical galaxy M87 (NGC 4486) using discrete chemo-dynamical, axisymmetric Jeans equation modelling. Our catalogue comprises 894 globular clusters (GCs) extending to a projected radius of $\sim 430$ kpc with line-of-sight velocities and colours, and Multi Unit Spectroscopic Explorer (MUSE) integral field unit data within the central $2.4$ kpc of the main galaxy. The gravitational potential for our models is a combination of a luminous matter potential with a varying mass-to-light ratio for the main galaxy, a supermassive black hole and a dark matter (DM) potential with a cusped or cored DM halo. The best-fitting models with either a cusped or a cored DM halo show no significant differences and both are acceptable. We obtain a total mass of $(2.16 \pm 0.38) \times 10^{13} M_{\odot}$ within $\sim$ 400 kpc. By including the stellar mass-to-light ratio gradient, the DM fraction increases from $\sim$ 26 percent (with no gradient) to $\sim$ 73 percent within $1\,R_e^{\rm maj}$ (major axis of half-light isophote, 14.2 kpc), and from $\sim$ 84 percent to $\sim$ 94 percent within $5\,R_e^{\rm maj}$ (71.2 kpc). Red GCs have moderate rotation with $V_{\rm max}/\sigma \sim$ 0.4, and blue GCs have weak rotation with $V_{\rm max}/\sigma \sim$ 0.1. Red GCs have tangential velocity dispersion anisotropy, while blue GCs are consistent with being nearly isotropic. Our results suggest that red GCs are more likely to be born in-situ, while blue GCs are more likely to be accreted.
[6]  oai:arXiv.org:2005.09149  [pdf] - 2097514
A Universal fundamental plane and the $M_{dyn}-M_{\star}$ relation for galaxies with CALIFA and MaNGA
Comments: Submitted to ApJ. 21 pages, 12 figures
Submitted: 2020-05-18
We use the stellar kinematics for $2458$ galaxies from the MaNGA survey to explore dynamical scaling relations between the stellar mass $M_{\star}$ and the total velocity parameter at the effective radius, $R_e$, defined as $S_{K}^{2}=KV_{R_e}^{2}+\sigma_{\star_e}^{2}$, which combines rotation velocity $V_{R_e}$, and velocity dispersion $\sigma_{\star_e}$. We confirm that spheroidal and spiral galaxies follow the same $M_{\star}-S_{0.5}$ scaling relation with lower scatter than the $M_{\star}-V_{R_e}$ and $M_{\star}-\sigma_{\star_e}$ ones. We also explore a more general Universal Fundamental Plane described by the equation $log(\Upsilon_{e}) = log (S_{0.5}^{2}) - log (I_{e}) - log (R_{e}) + C$, which in addition to kinematics, $S_{0.5}$, and effective radius, $R_e$, it includes surface brightness, $I_e$, and dynamical mass-to-light ratio, $\Upsilon_e$. We use sophisticated Schwarzschild dynamical models for a sub-sample of 300 galaxies from the CALIFA survey to calibrate the so called Universal Fundamental Plane. That calibration allows us to propose both: (i) a parametrization to estimate the difficult-to-measure dynamical mass-to-light ratio at the effective radius; and (ii) a new dynamical mass proxy consistent with dynamical models within $0.09\ dex$. We reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. We use the estimated dynamical mass-to-light ratio from our analysis, $\Upsilon_{e}^{fit}$, to explore the Universal Fundamental Plane with the MaNGA data set. We find that all classes of galaxies, from spheroids to disks, follow this Universal Fundamental Plane with a scatter significantly smaller $(0.05\ dex)$ than the one reported for the $M_{\star}-S_{0.5}$ relation $(0.1\ dex)$, the Fundamental Plane $(\sim 0.09\ dex)$ and comparable with Tully-Fisher studies $(\sim 0.05\ dex)$, but for a wider range of galaxy types.
[7]  oai:arXiv.org:2003.13705  [pdf] - 2097381
The Fornax 3D project: Globular clusters tracing kinematics and metallicities
Comments: 13 pages + appendix, accepted for publication in A&A, updated to published version
Submitted: 2020-03-30, last modified: 2020-05-13
Globular clusters (GCs) are found ubiquitously in massive galaxies and due to their old ages, they are regarded as fossil records of galaxy evolution. Spectroscopic studies of GC systems are often limited to the outskirts of galaxies, where GCs stand out against the galaxy background and serve as bright tracers of galaxy assembly. In this work, we use the capabilities of the Multi Unit Explorer Spectrograph (MUSE) to extract a spectroscopic sample of 722 GCs in the inner regions ($\lesssim 3 R_\text{eff}$) of 32 galaxies in the Fornax cluster. These galaxies were observed as part of the Fornax 3D project, a MUSE survey that targets early and late-type galaxies within the virial radius of Fornax. After accounting for the galaxy background in the GC spectra, we extracted line-of-sight velocities and determined metallicities of a sub-sample of 238 GCs. We found signatures of rotation within GC systems, and comparing the GC kinematics and that of the stellar body shows that the GCs trace the spheroid of the galaxies. While the red GCs prove to closely follow the metallicity profile of the host galaxy, the blue GCs show a large spread of metallicities but they are generally more metal-poor than the host.
[8]  oai:arXiv.org:2003.13707  [pdf] - 2097382
The Fornax 3D project: Non-linear colour-metallicity relation of globular clusters
Comments: 12 pages + appendix, updated to published version
Submitted: 2020-03-30, last modified: 2020-05-13
Globular cluster (GC) systems of massive galaxies often show a bimodal colour distribution. This has been interpreted as a metallicity bimodality, created by a two-stage galaxy formation where the red, metal-rich GCs were formed in the parent halo and the blue metal-poor GCs were accreted. This interpretation, however, crucially depends on the assumption that GCs are exclusively old stellar systems with a linear colour-metallicity relation (CZR). The shape of the CZR and range of GC ages are currently under debate, because their study requires high quality spectra to derive reliable stellar population properties. We determined metallicities with full spectral fitting from a sample of 187 GCs with high spectral signal-to-noise ratio in 23 galaxies of the Fornax cluster that were observed as part of the Fornax 3D project. The derived CZR from this sample is non-linear and can be described by a piecewise linear function with a break point at ($g - z$) $\sim$ 1.1 mag. The less massive galaxies in our sample ($M_\ast < 10^{10} M_\odot$) appear to have slightly younger GCs, but the shape of the CZR is insensitive to the GC ages. Although the least massive galaxies lack red, metal-rich GCs, a non-linear CZR is found irrespective of the galaxy mass, even in the most massive galaxies ($M_\ast \geq 10^{11} M_\odot$). Our CZR predicts narrow unimodal GC metallicity distributions for low mass and broad unimodal distributions for very massive galaxies, dominated by a metal-poor and metal-rich peak, respectively, and bimodal distributions for galaxies with intermediate masses (10$^{10}$ $\leq$ $M_\ast < 10^{11} M_\odot$) as a consequence of the relative fraction of red and blue GCs. The diverse metallicity distributions challenge the simple differentiation of GC populations solely based on their colour.
[9]  oai:arXiv.org:2005.04052  [pdf] - 2091825
Unravelling stellar populations in the Andromeda Galaxy
Comments: 14 pages, 9 figures, 1.1MB, submitted to A&A
Submitted: 2020-05-08
To understand the history and formation mechanisms of galaxies it is crucial to determine their current multidimensional structure. Here we focus on stellar population properties, such as metallicity and [$\alpha$/Fe] enhancement. We devise a new technique to recover the distribution of these parameters using spatially resolved, line-of-sight averaged data. Our chemodynamical method is based on the made-to-measure (M2M) framework and results in an $N$-body model for the abundance distribution. We test our method on a mock dataset and find that the radial and azimuthal profiles are well-recovered, however only the overall shape of the vertical profile matches the true profile. We apply our procedure to spatially resolved maps of mean [Z/H] and [$\alpha$/Fe] for the Andromeda Galaxy, using an earlier barred dynamical model of M31. We found that the metallicity is enhanced along the bar, with possible maxima at the ansae. In the edge-on view the [Z/H] distribution has an X shape due to the boxy/peanut bulge; the average vertical metallicity gradient is equal to $-0.132\pm0.006$ dex/kpc. We identify a metallicity-enhanced ring around the bar which also has relatively lower [$\alpha$/Fe]. The highest [$\alpha$/Fe] is found in the centre, due to the classical bulge. Away from the centre, the $\alpha$-overabundance in the bar region increases with height, which could be an indication of a thick disc. We argue that the galaxy assembly resulted in a sharp peak of metallicity in the central few hundred parsecs and a more gentle negative gradient in the remaining disc, but no [$\alpha$/Fe] gradient. The formation of the bar lead to the re-arrangement of the [Z/H] distribution, causing a flat gradient along the bar. Subsequent star formation close to the bar ends may have produced the metallicity enhancements at the ansae and the [Z/H] enhanced lower-$\alpha$ ring.
[10]  oai:arXiv.org:2004.14029  [pdf] - 2101499
A Fast Radio Burst discovered in FAST drift scan survey
Comments: 7 pages, 3 figures, accepted for publication in ApJL
Submitted: 2020-04-29
We report the discovery of a highly dispersed fast radio burst, FRB~181123, from an analysis of $\sim$1500~hr of drift-scan survey data taken using the Five-hundred-meter Aperture Spherical radio Telescope (FAST). The pulse has three distinct emission components, which vary with frequency across our 1.0--1.5~GHz observing band. We measure the peak flux density to be $>0.065$~Jy and the corresponding fluence $>0.2$~Jy~ms. Based on the observed dispersion measure of 1812~cm$^{-3}$~pc, we infer a redshift of $\sim 1.9$. From this, we estimate the peak luminosity and isotropic energy to be $\lesssim 2\times10^{43}$~erg~s$^{-1}$ and $\lesssim 2\times10^{40}$~erg, respectively. With only one FRB from the survey detected so far, our constraints on the event rate are limited. We derive a 95\% confidence lower limit for the event rate of 900 FRBs per day for FRBs with fluences $>0.025$~Jy~ms. We performed follow-up observations of the source with FAST for four hours and have not found a repeated burst. We discuss the implications of this discovery for our understanding of the physical mechanisms of FRBs.
[11]  oai:arXiv.org:2003.10456  [pdf] - 2096352
Fornax 3D project: automated detection of planetary nebulae in the centres of early-type galaxies and first results
Comments: 20 pages, 17 figures, 5 tables. Accepted for publication in Astronomy and Astrophysics Journal. Replaced with language edited version
Submitted: 2020-03-23, last modified: 2020-04-23
Extragalactic planetary nebulae (PNe) are detectable through relatively strong nebulous [OIII] emission and act as direct probes into the local stellar population. Because they have an apparently universal invariant magnitude cut-off, PNe are also considered to be a remarkable standard candle for distance estimation. Through detecting PNe within the galaxies, we aim to connect the relative abundances of PNe to the properties of their host galaxy stellar population. By removing the stellar background components from FCC 167 and FCC 219, we aim to produce PN luminosity functions (PNLF) of these galaxies, and thereby also estimate the distance modulus to these two systems. Finally, we test the reliability and robustness of our novel detection and analysis method. It detects unresolved point sources by their [OIII] 5007{\AA} emission within regions that have previously been unexplored. We model the [OIII] emissions in the spatial and spectral dimensions together, as afforded to us by the Multi Unit Spectroscopic Explorer (MUSE), and we draw on data gathered as part of the Fornax3D survey. For each source, we inspect the properties of the nebular emission lines to remove other sources that might hinder the safe construction of the PNLF, such as supernova remnants and HII regions. As a further step, we characterise any potential limitations and draw conclusions about the reliability of our modelling approach through a set of simulations. By applying this novel detection and modelling approach to integral field unit observations, we report for the distance estimates and luminosity-specific PNe frequency values for the two galaxies. Furthermore, we include an overview of source contamination, galaxy differences, and possible effects on the PNe populations in the dense stellar environments.
[12]  oai:arXiv.org:2002.02130  [pdf] - 2070116
First SETI Observations with China's Five-hundred-meter Aperture Spherical radio Telescope (FAST)
Comments: 21 pages, 19 figures, 3 tables, published in ApJ
Submitted: 2020-02-06, last modified: 2020-03-25
The Search for Extraterrestrial Intelligence (SETI) attempts to address the possibility of the presence of technological civilizations beyond the Earth. Benefiting from high sensitivity, large sky coverage, an innovative feed cabin for China's Five-hundred-meter Aperture Spherical radio Telescope (FAST), we performed the SETI first observations with FAST's newly commisioned 19-beam receiver; we report preliminary results in this paper. Using the data stream produced by the SERENDIP VI realtime multibeam SETI spectrometer installed at FAST, as well as its off-line data processing pipelines, we identify and remove four kinds of radio frequency interference(RFI): zone, broadband, multi-beam, and drifting, utilizing the Nebula SETI software pipeline combined with machine learning algorithms. After RFI mitigation, the Nebula pipeline identifies and ranks interesting narrow band candidate ET signals, scoring candidates by the number of times candidate signals have been seen at roughly the same sky position and same frequency, signal strength, proximity to a nearby star or object of interest, along with several other scoring criteria. We show four example candidates groups that demonstrate these RFI mitigation and candidate selection. This preliminary testing on FAST data helps to validate our SETI instrumentation techniques as well as our data processing pipeline.
[13]  oai:arXiv.org:2003.05561  [pdf] - 2115226
Disentangling the formation history of galaxies via population-orbit superposition: method validation
Comments: 20 pages, submitted to MNRAS
Submitted: 2020-03-11
We present population-orbit superposition models for external galaxies based on Schwarzschild's orbit-superposition method, by tagging the orbits with age and metallicity. The models fit the density distributions, as well as kinematic, age and metallicity maps from Integral Field Unit (IFU) spectroscopy observations. We validate the method and demonstrate its power by applying it to mock data, similar to those obtained by the Multi-Unit Spectroscopic Explorer (MUSE) IFU on the Very Large Telescope (VLT). These mock data are created from Auriga galaxy simulations, viewed at three different inclination angles ($\vartheta=40^o, 60^o, 80^o$). Constrained by MUSE-like mock data, our model can recover the galaxy's stellar orbit distribution projected in orbital circularity $\lambda_z$ vs. radius $r$, the intrinsic stellar population distribution in age $t$ vs. metallicity $Z$, and the correlation between orbits' circularity $\lambda_z$ and stellar age $t$. A physically motivated age-metallicity relation improves recovering the intrinsic stellar population distributions. We decompose galaxies into cold, warm and hot + counter-rotating components based on their orbit circularity distribution, and find that the surface density, mean velocity, velocity dispersion, age and metallicity maps of each component from our models well reproduce those from simulation, especially for projections close to edge-on. These galaxies exhibit strong global age vs. $\sigma_z$ relation, which is well recovered by our model. The method has the power to reveal the detailed build-up of stellar structures in galaxies, and offers a complement to local resolved, and high-redshift studies of galaxy evolution.
[14]  oai:arXiv.org:1909.06004  [pdf] - 2061570
Three dimensional projection effects on chemistry in a Planck galactic cold clump
Comments: Accepted for publication in ApJ. 17 pages, 11 figures, 3 tables
Submitted: 2019-09-12, last modified: 2020-01-28
Offsets of molecular line emission peaks from continuum peaks are very common but frequently difficult to explain with a single spherical cloud chemical model. We propose that the spatial projection effects of an irregular three dimensional (3D) cloud structure can be a solution. This work shows that the idea can be successfully applied to the Planck cold clump G224.4-0.6 by approximating it with four individual spherically symmetric cloud cores whose chemical patterns overlap with each other to produce observable line maps. With the empirical physical structures inferred from the observation data of this clump and a gas-grain chemical model, the four cores can satisfactorily reproduce its 850 $\mu$m continuum map and the diverse peak offsets of CCS, HC$_3$N and N$_2$H$^+$ simultaneously at chemical ages of about $8\times 10^5\sim 3\times 10^6$ yrs. The 3D projection effects on chemistry has the potential to explain such asymmetrical distributions of chemicals in many other molecular clouds.
[15]  oai:arXiv.org:1911.00777  [pdf] - 2021253
SDSS-IV MaNGA: Internal mass distributions and orbital structures of early-type galaxies and their dependence on environment
Comments: 19 pages, 18 figures
Submitted: 2019-11-02, last modified: 2019-12-27
In our earlier 2019 paper, we evaluated the reliability of Schwarzschild's orbit-superposition dynamical modelling method in estimating the internal mass distribution, intrinsic stellar shapes and orbit distributions of early-type galaxies (ETGs) taken from the Illustris cosmological simulation. We now apply the same techniques to galaxies taken from the integral-field survey Mapping Nearby Galaxies with APO (MaNGA), using a sample of 149 ETGs in the mass range of $10^{9.90}\sim10^{11.80} M_{\odot}$ and made up of 105 central and 44 satellite galaxies. We find that low-mass ETGs with $\log(M_*/M_{\odot})<11.1$ have an average dark matter fraction of $\sim0.2$ within one effective radius $R_{\rm e}$, tend to be oblate-like, and are dominated by rotation about their minor axis. High-mass ETGs with $\log(M_*/M_{\odot})>11.1$ have an average dark matter fraction of $\sim0.4$ within one effective radius $R_{\rm e}$, tend to be prolate-like, and are dominated by rotation about their major axis and by centrophilic orbits. The changes of internal structures within one $R_{\rm e}$ are dominated by the total stellar mass of the individual galaxies. We find no differences of internal structures between central and satellite ETGs for the same stellar masses. However, for similar stellar mass and colour distributions, we find that ETGs more prolate-like, or with more hot orbits, tend to have higher close neighbour counts at $r_p\sim40$ kpc.
[16]  oai:arXiv.org:1911.11058  [pdf] - 2003201
Mapping the dark matter halo of early-type galaxy NGC 2974 through orbit-based models with combined stellar and cold gas kinematics
Comments: 13 pages, 11 figures, accepted for publication in MNRAS
Submitted: 2019-11-25
We present an orbit-based method of combining stellar and cold gas kinematics to constrain the dark matter profile of early-type galaxies. We apply this method to early-type galaxy NGC 2974, using Pan-STARRS imaging and SAURON stellar kinematics to model the stellar orbits, and introducing HI kinematics from VLA observation as a tracer of the gravitational potential. The introduction of the cold gas kinematics shows a significant effect on the confidence limits of especially the dark halo properties: we exclude more than $95\%$ of models within the $1-\sigma$ confidence level of Schwarzschild modelling with only stellar kinematics, and reduce the relative uncertainty of the dark matter fraction significantly to $10\%$ within $5 R_\mathrm{e}$. Adopting a generalized-NFW dark matter profile, we measure a shallow cuspy inner slope of $0.6^{+0.2}_{-0.3}$ when including the cold gas kinematics in our model. We cannot constrain the inner slope with the stellar kinematics alone.
[17]  oai:arXiv.org:1910.06236  [pdf] - 2008006
The CALIFA view on stellar angular momentum across the Hubble sequence
Comments: 17 pages, 7 figures, A&A accepted
Submitted: 2019-10-14
[Abridged] We present the apparent stellar angular momentum of 300 galaxies across the Hubble sequence, using integral-field spectroscopic data from the CALIFA survey. Adopting the same $\lambda_\mathrm{R}$ parameter previously used to distinguish between slow and fast rotating early-type (elliptical and lenticular) galaxies, we show that spiral galaxies as expected are almost all fast rotators. Given the extent of our data, we provide relations for $\lambda_\mathrm{R}$ measured in different apertures, including conversions to long-slit 1D apertures. Our sample displays a wide range of $\lambda_\mathrm{Re}$ values, consistent with previous IFS studies. The fastest rotators are dominated by relatively massive and highly star-forming Sb galaxies, which preferentially reside in the main star-forming sequence. These galaxies reach $\lambda_\mathrm{Re}$ values of $\sim$0.85, are the largest galaxies at a given mass, and display some of the strongest stellar population gradients. Compared to the population of S0 galaxies, our findings suggest that fading may not be the dominant mechanism transforming spirals into lenticulars. Interestingly, we find that $\lambda_\mathrm{Re}$ decreases for late-type Sc and Sd spiral galaxies, with values than in occasions puts them in the slow-rotator regime. While for some of them this can be explained by their irregular morphologies and/or face-on configurations, others are edge-on systems with no signs of significant dust obscuration. The latter are typically at the low-mass end, but this does not explain their location in the classical ($V/\sigma$,$\varepsilon$) and ($\lambda_\mathrm{Re}$,$\varepsilon$) diagrams. Our initial investigations, based on dynamical models, suggest that these are dynamically hot disks, probably influenced by the observed important fraction of dark matter within R$_\mathrm{e}$.
[18]  oai:arXiv.org:1909.11806  [pdf] - 1969076
Unit panel nodes detection by CNN on FAST reflector
Comments: 13 pages, 12 figures, 2 tables, CNN applied on FAST's reflector measurement; matches the published version in RAA
Submitted: 2019-09-25
The 500-meter Aperture Spherical Radio Telescope(FAST) has an active reflector. During the observation, the reflector will be deformed into a paraboloid of 300-meters. To improve its surface accuracy, we propose a scheme for photogrammetry to measure the positions of 2226 nodes on the reflector. And the way to detect the nodes in the photos is the key problem in photogrammetry. This paper applies Convolutional Neural Network(CNN) with candidate regions to detect the nodes in the photos. The experiment results show a high recognition rate of 91.5%, which is much higher than the recognition rate of traditional edge detection.
[19]  oai:arXiv.org:1907.03770  [pdf] - 1913021
NGC 7457: Evidence for merger-driven cylindrical rotation in disc galaxies
Comments: 15 pages, 12 figures, 2 tables. Accepted by MNRAS
Submitted: 2019-07-08
We construct Schwarzschild orbit-based models of NGC 7457, known as a peculiar low-mass lenticular galaxy. Our best-fitting model successfully retrieves most of the unusual kinematics behaviours of this galaxy, in which, the orbital distribution of stars is dominated by warm and hot orbits. The reconstructed surface brightness of the hot component matches fairly well the photometric bulge and the reconstructed LOSVD map of this component shows clear rotation around the major photometric axis of the galaxy. In the absence of a dominant cold component, the outer part of our model is dominated by warm orbits, representing an exponential thick disc. Our orbital analysis also confirm the existence of a counter-rotating orbital substructure in the very centre, reported in previous observational studies. By comparing our model with a variety of simulation studies, and considering the stellar kinematics and populations properties of this galaxy, we suggest that the thick disc is most likely a dynamically heated structure, formed through the interactions and accretion of satellite(s) with near-polar initial inclination. We also suggest a merger-driven process as the most plausible scenario to explain the observed and dynamically-modelled properties of the bulge of NGC 7457. We conclude that both the high level of cylindrical rotation and unusually low velocity dispersion reported for the NGC 7457 have most-likely external origins. Therefore, NGC 7457 could be considered as a candidate for merger-driven cylindrical rotation in the absence of a strong bar in disc galaxies.
[20]  oai:arXiv.org:1907.01007  [pdf] - 1938516
Constraining nuclear star cluster formation using MUSE-AO observations of the early-type galaxy FCC47
Comments: 21 pages, accepted in A&A
Submitted: 2019-07-01
Nuclear star clusters (NSCs) are found in at least 70% of all galaxies, but their formation path is still unclear. In the most common scenarios, NSCs form in-situ from the galaxy's central gas reservoir, through merging of globular clusters (GCs), or through a combination of the two. As the scenarios pose different expectations for angular momentum and stellar population properties of the NSC in comparison to the host galaxy and the GC system, it is necessary to characterise the stellar light, NSC and GCs simultaneously. The large NSC (r$_\rm{eff} = 66$ pc) and rich GC system of the early-type Fornax cluster galaxy FCC47 (NGC1336) render this galaxy an ideal laboratory to constrain NSC formation. Using MUSE science verification data assisted by adaptive optics, we obtained maps for the stellar kinematics and for stellar-population properties of FCC47. We extracted the spectra of the central NSC and determined line-of-sight velocities of 24 GCs and metallicities of five. FCC47 shows two decoupled components (KDCs): a rotating disk and the NSC. Our orbit-based dynamical Schwarzschild model revealed that the NSC is a distinct kinematic feature and it constitutes the peak of metallicity and old ages in the galaxy. The main body consists of two counter-rotating populations and is dominated by a more metal-poor population. The GC system is bimodal with a dominant metal-poor population and the total GC system mass is $\sim 17\%$ of the NSC mass ($\sim$ 7 $\times$ $10^8 M_\odot$). The rotation, high metallicity and high mass of the NSC cannot be uniquely explained by GC-inspiral and most likely requires additional, but quickly quenched, in-situ formation. The presence of two KDCs most probably are evidence of a major merger that has altered the structure of FCC47 significantly, indicating the important role of galaxy mergers in forming the complex kinematics in the galaxy-NSC system.
[21]  oai:arXiv.org:1804.06635  [pdf] - 1901726
NIHAO XVI: The properties and evolution of kinematically selected discs, bulges and stellar haloes
Comments: Accepted by MNRAS. First application of the methods described in arXiv:1804.05576
Submitted: 2018-04-18, last modified: 2019-06-04
We use 25 simulated galaxies from the NIHAO project to define and characterize a variety of kinematic stellar structures: thin and thick discs, large scale single discs, classical and pseudo bulges, spheroids, inner discs, and stellar haloes. These structures have masses, spins, shapes and rotational support in good agreement with theoretical expectations and observational data. Above a dark matter halo mass of $2.5\times10^{\rm~11}M_{\rm\odot}$, all galaxies have a classical bulge and 70\% have a thin and thick disc. The kinematic (thin) discs follow a power-law relation between angular momentum and stellar mass $J_{\rm *}=3.4M_{\rm *}^{\rm1.26\pm0.06}$, in very good agreement with the prediction based on the empirical stellar-to-halo mass relation in the same mass range, and show a strong correlation between maximum `observed' rotation velocity and dark matter halo circular velocity $v_{\rm c}=6.4v_{\rm max}^{0.64\pm0.04}$. Tracing back in time these structures' progenitors, we find all to lose a fraction $1-f_j$ of their maximum angular momentum. Thin discs are significantly better at retaining their high-redshift spins ($f_j\sim0.70$) than thick ones ($f_j\sim0.40$). Stellar haloes have their progenitor baryons assembled the latest ($z_{\rm~1/2}\sim1.1$) and over the longest timescales ($\tau\sim6.2$~Gyr), and have the smallest fraction of stars born in-situ ($f_{\rm in-situ}=0.35\pm0.14$). All other structures have $1.5\lesssim z_{\rm1/2}\lesssim3$, $\tau=4\pm2$~Gyr and $f_{\rm in-situ}\gtrsim0.9$.
[22]  oai:arXiv.org:1905.11415  [pdf] - 1953534
A study of stellar orbit fractions: simulated IllustrisTNG galaxies compared to CALIFA observations
Comments: 12 pages, 6 figures, 3 tables, submitted to MNRAS. Comments are welcome!
Submitted: 2019-05-27
Motivated by the recently discovered kinematic "Hubble sequence" shown by the stellar orbit-circularity distribution of 260 CALIFA galaxies, we make use of a comparable galaxy sample at z = 0 with a stellar mass range from 5E9 to 5E11 solar masses, selected from the IllustrisTNG simulation and study their stellar orbit compositions in relation to a number of other fundamental galaxy properties.We find that the TNG100 simulation broadly reproduces the observed fractions of different orbital components and their stellar mass dependencies. In particular, the mean mass dependencies of the luminosity fractions for the kinematically warm and hot orbits are well reproduced within model uncertainties of the observed galaxies. The simulation also largely reproduces the observed peak and trough features at a stellar mass of 1-2E10 solar masses, in the mean distributions of the cold- and hot-orbit fractions, respectively, indicating fewer cooler orbits and more hotter orbits in both more- and less-massive galaxies beyond such a mass range. Several marginal disagreements are seen between the simulation and observations: the average cold-orbit (counter-rotating) fractions of the simulated galaxies below (above) a stellar mass of 6E10 solar masses, are systematically higher than the observational data by < 10% (absolute orbital fraction); the simulation also seems to produce more scatter for the cold-orbit fraction and less so for the non-cold orbits at any given galaxy mass. Possible causes that stem from the adopted heating mechanisms are discussed.
[23]  oai:arXiv.org:1905.10322  [pdf] - 1912875
Carbon-Chain Molecules in Molecular Outflows and Lupus I Region--New Producing Region and New Forming Mechanism
Comments: 7 figures, 8 tables, accepted by MNRAS
Submitted: 2019-05-24
Using the new equipment of the Shanghai Tian Ma Radio Telescope, we have searched for carbon-chain molecules (CCMs) towards five outflow sources and six Lupus I starless dust cores, including one region known to be characterized by warm carbon-chain chemistry (WCCC), Lupus I-1 (IRAS 15398-3359), and one TMC-1 like cloud, Lupus I-6 (Lupus-1A). Lines of HC3N J=2-1, HC5N J=6-5, HC7N J=14-13, 15-14, 16-15 and C3S J=3-2 were detected in all the targets except in the outflow source L1660 and the starless dust core Lupus I-3/4. The column densities of nitrogen-bearing species range from 10$^{12}$ to 10$^{14}$ cm$^{-2}$ and those of C$_3$S are about 10$^{12}$ cm$^{-2}$. Two outflow sources, I20582+7724 and L1221, could be identified as new carbon-chain--producing regions. Four of the Lupus I dust cores are newly identified as early quiescent and dark carbon-chain--producing regions similar to Lup I-6, which together with the WCCC source, Lup I-1, indicate that carbon-chain-producing regions are popular in Lupus I which can be regard as a Taurus like molecular cloud complex in our Galaxy. The column densities of C3S are larger than those of HC7N in the three outflow sources I20582, L1221 and L1251A. Shocked carbon-chain chemistry (SCCC) is proposed to explain the abnormal high abundances of C3S compared with those of nitrogen-bearing CCMs. Gas-grain chemical models support the idea that shocks can fuel the environment of those sources with enough $S^+$ thus driving the generation of S-bearing CCMs.
[24]  oai:arXiv.org:1904.11605  [pdf] - 1874588
Combining Stellar Populations with Orbit-Superposition Dynamical Modelling - the Formation History of the Lenticular Galaxy NGC 3115
Comments: 21 pages. 14 figures + appendices. Accepted for publication in MNRAS
Submitted: 2019-04-25, last modified: 2019-04-29
We present a combination of the Schwarzschild orbit-superposition dynamical modelling technique with the spatially-resolved mean stellar age and metallicity maps to uncover the formation history of galaxies. We apply this new approach to a remarkable 5-pointing mosaic of VLT/MUSE observations obtained by Gu\'erou et al. (2016) extending to a maximum galactocentric distance of $120''$ ($\sim 5.6\ {\rm kpc}$) along the major axis, corresponding to $\sim 2.5\ {\rm R}_e$. Our method first identifies `families' of orbits from the dynamical model that represent dynamically-distinct structures of the galaxy. Individual ages and metallicities of these components are then fit for using the stellar-population information. Our results highlight components of the galaxy that are distinct in the combined stellar dynamics/populations space, which implies distinct formation paths. We find evidence for a dynamically-cold, metal-rich disk, consistent with a gradual in-situ formation. This disk is embedded in a generally-old population of stars, with kinematics ranging from dispersion-dominated in the centre to an old, diffuse, metal-poor stellar halo at the extremities. We find also a direct correlation between the dominant dynamical support of these components, and their associated age, akin to the relation observed in the Milky Way. This approach not only provides a powerful model for inferring the formation history of external galaxies, but also paves the way to a complete population-dynamical model.
[25]  oai:arXiv.org:1904.12942  [pdf] - 1875506
Evaluating the ability of triaxial Schwarzschild modelling to estimate properties of galaxies from the Illustris simulation
Comments: 21 pages, 15 figures
Submitted: 2019-04-29
We evaluate the capabilities of Schwarzschild's orbit-superposition method by applying it to galaxies from the large scale, high resolution Illustris simulation. Nine early-type galaxies with a range of triaxiality are selected, and we create mock integral field unit data for five line-of-sight projections of each galaxy. Each of the 45 mock data sets is taken as an independent observed galaxy. Using van den Bosch's 2008 triaxial Schwarzschild implementation, we assess model estimates of various galaxy properties, covering mass profiles, intrinsic shapes, stellar orbit distributions and velocity anisotropies. Total mass within $\overline{R_{\rm e}}$ is recovered well with average deviations within $\pm15$ percent. Stellar mass is underestimated by $\sim24$ percent and dark matter overestimated by $\sim38$ percent (assuming an NFW dark matter profile and allowing for degeneracy between stellar mass and dark matter mass). Using a gNFW profile, these values improve to $\sim13$ percent for stellar mass and $\sim18$ percent for dark matter. Axis ratio estimates show a moderate bias of $\Delta (b/a)=0.07$ and $\Delta (c/a)=0.14$ ($a\ge b\ge c$). Distributions of the orbit circularities $\lambda_z$ and $\lambda_x$, representing rotation about the minor and major axes, are well reconstructed. Separating orbits into thermal categories, our models match the average fractions of these categories to within $10$ percent. Velocity anisotropy is well estimated with values matching in the inner regions but becoming slightly radially biased in the outer regions. Overall, the galaxy property estimates we obtained using Schwarzschild modelling are not implausible and are representative of the simulated galaxies we modelled.
[26]  oai:arXiv.org:1904.07221  [pdf] - 1892246
The JCMT BISTRO Survey: The Magnetic Field of the Barnard 1 Star-Forming Region
Coudé, Simon; Bastien, Pierre; Houde, Martin; Sadavoy, Sarah; Friesen, Rachel; Di Francesco, James; Johnstone, Doug; Mairs, Steve; Hasegawa, Tetsuo; Kwon, Woojin; Lai, Shih-Ping; Qiu, Keping; Ward-Thompson, Derek; Berry, David; Chen, Michael Chun-Yuan; Fiege, Jason; Franzmann, Erica; Hatchell, Jennifer; Lacaille, Kevin; Matthews, Brenda C.; Moriarty-Schieven, Gerald H.; Pon, Andy; André, Philippe; Arzoumanian, Doris; Aso, Yusuke; Byun, Do-Young; Chakali, Eswaraiah; Chen, Huei-Ru; Chen, Wen Ping; Ching, Tao-Chung; Cho, Jungyeon; Choi, Minho; Chrysostomou, Antonio; Chung, Eun Jung; Doi, Yasuo; Drabek-Maunder, Emily; Dowell, C. Darren; Eyres, Stewart P. S.; Falle, Sam; Friberg, Per; Fuller, Gary; Furuya, Ray S.; Gledhill, Tim; Graves, Sarah F.; Greaves, Jane S.; Griffin, Matt J.; Gu, Qilao; Hayashi, Saeko S.; Hoang, Thiem; Holland, Wayne; Inoue, Tsuyoshi; Inutsuka, Shu-ichiro; Iwasaki, Kazunari; Jeong, Il-Gyo; Kanamori, Yoshihiro; Kataoka, Akimasa; Kang, Ji-hyun; Kang, Miju; Kang, Sung-ju; Kawabata, Koji S.; Kemper, Francisca; Kim, Gwanjeong; Kim, Jongsoo; Kim, Kee-Tae; Kim, Kyoung Hee; Kim, Mi-Ryang; Kim, Shinyoung; Kirk, Jason M.; Kobayashi, Masato I. N.; Koch, Patrick M.; Kwon, Jungmi; Lee, Jeong-Eun; Lee, Chang Won; Lee, Sang-Sung; Li, Dalei; Li, Di; Li, Hua-bai; Liu, Hong-Li; Liu, Junhao; Liu, Sheng-Yuan; Liu, Tie; van Loo, Sven; Lyo, A-Ran; Matsumura, Masafumi; Nagata, Tetsuya; Nakamura, Fumitaka; Nakanishi, Hiroyuki; Ohashi, Nagayoshi; Onaka, Takashi; Parsons, Harriet; Pattle, Kate; Peretto, Nicolas; Pyo, Tae-Soo; Qian, Lei; Rao, Ramprasad; Rawlings, Mark G.; Retter, Brendan; Richer, John; Rigby, Andrew; Robitaille, Jean-François; Saito, Hiro; Savini, Giorgio; Scaife, Anna M. M.; Seta, Masumichi; Shinnaga, Hiroko; Soam, Archana; Tamura, Motohide; Tang, Ya-Wen; Tomisaka, Kohji; Tsukamoto, Yusuke; Wang, Hongchi; Wang, Jia-Wei; Whitworth, Anthony P.; Yen, Hsi-Wei; Yoo, Hyunju; Yuan, Jinghua; Zenko, Tetsuya; Zhang, Chuan-Peng; Zhang, Guoyin; Zhou, Jianjun; Zhu, Lei
Comments: 22 pages, 7 figures, accepted in the Astrophysical Journal
Submitted: 2019-04-15, last modified: 2019-04-23
We present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.
[27]  oai:arXiv.org:1904.05482  [pdf] - 1894293
PSR J1926-0652: A Pulsar with Interesting Emission Properties Discovered at FAST
Comments: 13pages with 12 figures
Submitted: 2019-04-10
We describe PSR J1926-0652, a pulsar recently discovered with the Five-hundred-meter Aperture Spherical radio Telescope (FAST). Using sensitive single-pulse detections from FAST and long-term timing observations from the Parkes 64-m radio telescope, we probed phenomena on both long and short time scales. The FAST observations covered a wide frequency range from 270 to 800 MHz, enabling individual pulses to be studied in detail. The pulsar exhibits at least four profile components, short-term nulling lasting from 4 to 450 pulses, complex subpulse drifting behaviours and intermittency on scales of tens of minutes. While the average band spacing P3 is relatively constant across different bursts and components, significant variations in the separation of adjacent bands are seen, especially near the beginning and end of a burst. Band shapes and slopes are quite variable, especially for the trailing components and for the shorter bursts. We show that for each burst the last detectable pulse prior to emission ceasing has different properties compared to other pulses. These complexities pose challenges for the classic carousel-type models.
[28]  oai:arXiv.org:1903.10514  [pdf] - 1905719
The Fornax 3D project: a two-dimensional view of the stellar initial mass function in the massive lenticular galaxy FCC 167
Comments: 12 pages, 7 figures. Accepted for publication in Astronomy & Astrophysics
Submitted: 2019-03-25
The stellar initial mass function (IMF) regulates the baryonic cycle within galaxies, and is a key ingredient to translate observations into physical quantities. Although for decades it was assumed to be universal, there is now growing observational evidence showing that the center of massive early-type galaxies host an enhanced population of low-mass stars compared to the expectations from the Milky Way. Moreover, these variations in the IMF have been found to be related to the radial metallicity variations in massive galaxies. We present here a two-dimensional stellar population analysis of the massive lenticular galaxy FCC 167 (NGC 1380) as part of the Fornax3D project. Using a newly developed stellar population fitting scheme, we derive a full two-dimensional IMF map of an early-type galaxy. This two-dimensional analysis allows us go further than a radial analysis, showing how the metallicity changes along a disc-like structure while the IMF follows a distinct, less disky distribution. Thus, our findings indicate that metallicity cannot be the sole driver of the observed radial IMF variations. In addition, a comparison with the orbital decomposition shows suggestive evidence of a coupling between stellar population properties and the internal dynamical structure of FCC 167, where metallicity and IMF maps seem to track the distribution of cold and warm orbits, respectively.
[29]  oai:arXiv.org:1902.07734  [pdf] - 1890373
The JCMT BISTRO Survey: The Magnetic Field In The Starless Core $\rho$ Ophiuchus C
Liu, Junhao; Qiu, Keping; Berry, David; Di Francesco, James; Bastien, Pierre; Koch, Patrick M.; Furuya, Ray S.; Kim, Kee-Tae; Coudé, Simon; Lee, Chang Won; Soam, Archana; Eswaraiah, Chakali; Li, Di; Hwang, Jihye; Lyo, A-Ran; Pattle, Kate; Hasegawa, Tetsuo; Kwon, Woojin; Lai, Shih-Ping; Ward-Thompson, Derek; Ching, Tao-Chung; Chen, Zhiwei; Gu, Qilao; Li, Dalei; Li, Hua-bai; Liu, Hong-Li; Qian, Lei; Wang, Hongchi; Yuan, Jinghua; Zhang, Chuan-Peng; Zhang, Guoyin; Zhang, Ya-Peng; Zhou, Jianjun; Zhu, Lei; André, Philippe; Arzoumanian, Doris; Aso, Yusuke; Byun, Do-Young; Chen, Michael Chun-Yuan; Chen, Huei-Ru Vivien; Chen, Wen Ping; Cho, Jungyeon; Choi, Minho; Chrysostomou, Antonio; Chung, Eun Jung; Doi, Yasuo; Drabek-Maunder, Emily; Dowell, C. Darren; Eyres, Stewart P. S.; Falle, Sam; Fanciullo, Lapo; Fiege, Jason; Franzmann, Erica; Friberg, Per; Friesen, Rachel K.; Fuller, Gary; Gledhill, Tim; Graves, Sarah F.; Greaves, Jane S.; Griffin, Matt J.; Han, Ilseung; Hatchell, Jennifer; Hayashi, Saeko S.; Hoang, Thiem; Holland, Wayne; Houde, Martin; Inoue, Tsuyoshi; Inutsuka, Shu-ichiro; Iwasaki, Kazunari; Jeong, Il-Gyo; Johnstone, Doug; Kanamori, Yoshihiro; Kang, Ji-hyun; Kang, Miju; Kang, Sung-ju; Kataoka, Akimasa; Kawabata, Koji S.; Kemper, Francisca; Kim, Gwanjeong; Kim, Jongsoo; Kim, Kyoung Hee; Kim, Mi-Ryang; Kim, Shinyoung; Kirk, Jason M.; Kobayashi, Masato I. N.; Kusune, Takayoshi; Kwon, Jungmi; Lacaille, Kevin M.; Lee, Chin-Fei; Lee, Jeong-Eun; Lee, Hyeseung; Lee, Sang-Sung; Liu, Sheng-Yuan; Liu, Tie; van Loo, Sven; Mairs, Steve; Matsumura, Masafumi; Matthews, Brenda C.; Moriarty-Schieven, Gerald H.; Nagata, Tetsuya; Nakamura, Fumitaka; Nakanishi, Hiroyuki; Ohashi, Nagayoshi; Onaka, Takashi; Parker, Josh; Parsons, Harriet; Pascale, Enzo; Peretto, Nicolas; Pon, Andy; Pyo, Tae-Soo; Rao, Ramprasad; Rawlings, Mark G.; Retter, Brendan; Richer, John; Rigby, Andrew; Robitaille, Jean-François; Sadavoy, Sarah; Saito, Hiro; Savini, Giorgio; Scaife, Anna M. M.; Seta, Masumichi; Shinnaga, Hiroko; Tamura, Motohide; Tang, Ya-Wen; Tomisaka, Kohji; Tsukamoto, Yusuke; Wang, Jia-Wei; Whitworth, Anthony P.; Yen, Hsi-Wei; Yoo, Hyunju; Zenko, Tetsuya
Comments: 22 pages, 12 figures
Submitted: 2019-02-20
We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.
[30]  oai:arXiv.org:1812.09856  [pdf] - 1806354
Photometric Investigation of the Contact Binary GU Ori with High Metallicity
Comments: 17 pages, 6 figures
Submitted: 2018-12-24
GU Ori was observed with the 1m telescope at Yunnan Observatories in 2005. To determine its physical properties, the Wilson-Devinney program is used. The results reveal that GU Ori is a W-subtype shallow contact binary with a more massive but cooler star 2. The mass of its two component stars are estimated to be $M_1 = 0.45M_\odot$, $M_2 = 1.05M_\odot$. The O'Connell effect was reported to be negative on the light curves observed in 2005. However, it changed to a positive one on the light curves observed from 2011 to 2012. The mean surface temperatures of star 2 ($T_2$) determined by the two sets of light curves were different, which may result from stellar activities. The O - C diagram shows that the period of GU Ori is decreasing at a rate of $dP/dt=-6.24\times{10^{-8}}day\cdot year^{-1}$, which may be caused by mass transfer from star 2 to star 1 with a rate of $\frac{dM_{2}}{dt}= - 2.98\times{10^{-8}}M_\odot/year$. GU Ori is a contact binary with quite high metallicity.
[31]  oai:arXiv.org:1811.08625  [pdf] - 1789691
DE CVn: an eclipsing post-common envelope binary with a circumbinary disk and a giant planet
Comments: 19 pages, 5 figures
Submitted: 2018-11-21
We present a timing analysis of the eclipsing post-common envelope binary (PCEB) DE CVn. Based on new CCD photometric observations and the published data, we found that the orbital period in DE CVn has a cyclic period oscillation with an amplitude of $28.08$ s and a period of $11.22$ years plus a rapid period decrease at a rate of $\dot{P}=-3.35\times10^{-11}ss^{-1}$. According to the evolutionary theory, secular period decreases in PCEBs arise from angular momentum losses (AMLs) driven by gravitational radiation (GR) and magnetic braking (MB). However, the observed orbital decay is too fast to be produced by AMLs via GR and MB, indicating that there could be other AML mechanism. We suggest that a circumbinary disk around DE CVn may be responsible for the additional AML. The disk mass was derived as a few$\times$$10^{-4}$-$10^{-3}$$M_{\odot}$ , which is in agreement with that inferred from previous studies in the order of magnitude. The cyclic change is most likely result of the gravitational perturbation by a circumbinary object due to the Applegate's mechanism fails to explain such a large period oscillation. The mass of the potential third body is calculated as $M_{3}\sin{i'}=0.011(\pm0.003)M_{\odot}$. Supposing the circumbinary companion and the eclipsing binary is coplanar, its mass would correspond to a giant planet. This hypothetical giant planet is moving in a circular orbit of radius $\sim5.75(\pm2.02)$ AU around its host star.
[32]  oai:arXiv.org:1810.12223  [pdf] - 1791273
A dynamical view on stellar metallicity gradient diversity across the Hubble sequence with CALIFA
Comments: 20 pages, 13 figures
Submitted: 2018-10-29
We analyze radial stellar metallicity and kinematic profiles out to 1Re in 244 CALIFA galaxies ranging from morphological type E to Sd, to study the evolutionary mechanisms of stellar population gradients. We find that linear metallicity gradients exhibit a clear correlation with galaxy morphological type - with early type galaxies showing the steepest gradients. We show that the metallicity gradients simply reflect the local mass-metallicity relation within a galaxy. This suggests that the radial stellar population distribution within a galaxys effective radius is primarily a result of the \emph{in-situ} local star formation history. In this simple picture, the dynamically derived stellar surface mass density gradient directly predicts the metallicity gradient of a galaxy. We show that this correlation and its scatter can be reproduced entirely by using independent empirical galaxy structural and chemical scaling relations. Using Schwarzschild dynamical models, we also explore the link between a galaxys local stellar populations and their orbital structures. We find that galaxies angular momentum and metallicity gradients show no obvious causal link. This suggests that metallicity gradients in the inner disk are not strongly shaped by radial migration, which is confirmed by the lack of correlation between the metallicity gradients and observable probes of radial migration in the galaxies, such as bars and spiral arms. Finally, we find that galaxies with positive metallicity gradients become increasingly common towards low mass and late morphological types - consistent with stellar feedback more efficiently modifying the baryon cycle in the central regions of these galaxies.
[33]  oai:arXiv.org:1806.09013  [pdf] - 1709561
Kinematic scaling relations of CALIFA galaxies: A dynamical mass proxy for galaxies across the Hubble sequence
Comments: accepted by MNRAS
Submitted: 2018-06-23
We used ionized gas and stellar kinematics for 667 spatially resolved galaxies publicly available from the Calar Alto Legacy Integral Field Area survey (CALIFA) 3rd Data Release with the aim of studying kinematic scaling relations as the Tully $\&$ Fisher (TF) relation using rotation velocity, $V_{rot}$, the Faber $\&$ Jackson (FJ) relation using velocity dispersion, $\sigma$, and also a combination of $V_{rot}$ and $\sigma$ through the $S_{K}$ parameter defined as $S_{K}^2 = KV_{rot}^2 + \sigma^2$ with constant $K$. Late-type and early-type galaxies reproduce the TF and FJ relations. Some early-type galaxies also follow the TF relation and some late-type galaxies the FJ relation, but always with larger scatter. On the contrary, when we use the $S_{K}$ parameter, all galaxies, regardless of the morphological type, lie on the same scaling relation, showing a tight correlation with the total stellar mass, $M_\star$. Indeed, we find that the scatter in this relation is smaller or equal to that of the TF and FJ relations. We explore different values of the $K$ parameter without significant differences (slope and scatter) in our final results with respect the case $K=0.5$ besides than a small change in the zero point. We calibrate the kinematic $S_{K}^2$ dynamical mass proxy in order to make it consistent with sophisticated published dynamical models within $0.15\ dex$. We show that the $S_{K}$ proxy is able to reproduce the relation between the dynamical mass and the stellar mass in the inner regions of galaxies. Our result may be useful in order to produce fast estimations of the central dynamical mass in galaxies and to study correlations in large galaxy surveys.
[34]  oai:arXiv.org:1805.06131  [pdf] - 1717135
Magnetic fields towards Ophiuchus-B derived from SCUBA-2 polarization measurements
Soam, Archana; Pattle, Kate; Ward-Thompson, Derek; Lee, Chang Won; Sadavoy, Sarah; Koch, Patrick M.; Kim, Gwanjeong; Kwon, Jungmi; Kwon, Woojin; Arzoumanian, Doris; Berry, David; Hoang, Thiem; Tamura, Motohide; Lee, Sang-Sung; Liu, Tie; Kim, Kee-Tae; Johnstone, Doug; Nakamura, Fumitaka; Lyo, A-Ran; Onaka, Takashi; Kim, Jongsoo; Furuya, Ray S.; Hasegawa, Tetsuo; Lai, Shih-Ping; Bastien, Pierre; Chung, Eun Jung; Kim, Shinyoung; Parsons, Harriet; Rawlings, Mark; Mairs, Steve; Graves, Sarah F.; Robitaille, J. -F.; Liu, Hong-Li; Whitworth, Anthony P.; Eswaraiah, Chakali; Rao, Ramprasad; Yoo, Hyunju; Houde, Martin; Kang, Ji-hyun; Doi, Yasuo; Choi, Minho; Kang, Miju; Coude, Simon; Li, Hua-bai; Matsumura, Masafumi; Matthews, Brenda C.; Pon, Andy; Di Francesco, James; Hayashi, Saeko S.; Kawabata, Koji S.; Inutsuka, Shu-ichiro; Qiu, Keping; Franzmann, Erica; Friberg, Per; Greaves, Jane S.; Kirk, Jason M.; Li, Di; Shinnaga, Hiroko; van Loo, Sven; Aso, Yusuke; Byun, Do-Young; Chen, Huei-Ru; Chen, Mike C. -Y.; Chen, Wen Ping; Ching, Tao-Chung; Cho, Jungyeon; Chrysostomou, Antonio; Drabek-Maunder, Emily; Eyres, Stewart P. S.; Fiege, Jason; Friesen, Rachel K.; Fuller, Gary; Gledhill, Tim; Griffin, Matt J.; Gu, Qilao; Hatchell, Jennifer; Holland, Wayne; Inoue, Tsuyoshi; Iwasaki, Kazunari; Jeong, Il-Gyo; Kang, Sung-ju; Kemper, Francisca; Kim, Kyoung Hee; Kim, Mi-Ryang; Lacaille, Kevin M.; Lee, Jeong-Eun; Li, Dalei; Liu, Junhao; Liu, Sheng-Yuan; Moriarty-Schieven, Gerald H.; Nakanishi, Hiroyuki; Ohashi, Nagayoshi; Peretto, Nicolas; Pyo, Tae-Soo; Qian, Lei; Retter, Brendan; Richer, John; Rigby, Andrew; Savini, Giorgio; Scaife, Anna M. M.; Tang, Ya-Wen; Tomisaka, Kohji; Wang, Hongchi; Wang, Jia-Wei; Yen, Hsi-Wei; Yuan, Jinghua; Zhang, Chuan-Peng; Zhang, Guoyin; Zhou, Jianjun; Zhu, Lei; Andre, Philippe; Dowell, C. Darren; Falle, Sam; Tsukamoto, Yusuke; Kanamori, Yoshihiro; Kataoka, Akimasa; Kobayashi, Masato I. N.; Nagata, Tetsuya; Saito, Hiro; Seta, Masumichi; Hwang, Jihye; Han, Ilseung; Lee, Hyeseung; Zenko, Tetsuya
Comments: 17 pages, 10 figures, Accepted for publication in ApJ
Submitted: 2018-05-16, last modified: 2018-06-07
We present the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.
[35]  oai:arXiv.org:1806.02886  [pdf] - 1705352
Morphology and kinematics of orbital components in CALIFA galaxies across the Hubble sequence
Comments: accepted by MNRAS
Submitted: 2018-06-07
Based on the stellar orbit distribution derived from orbit-superposition Schwarzschild models, we decompose each of 250 representative present-day galaxies into four orbital components: cold with strong rotation, warm with weak rotation, hot with dominant random motion and counter-rotating (CR). We rebuild the surface brightness ($\Sigma$) of each orbital component and we present in figures and tables a quantification of their morphologies using the Sersic index \textit{n}, concentration $C = \log{(\Sigma_{0.1R_e}/\Sigma_{R_e})}$ and intrinsic flattening $q_{\mathrm{Re}}$ and $q_{\mathrm{Rmax}}$, with $R_e$ the half-light-radius and $R_{\mathrm{max}}$ the CALIFA data coverage. We find that: (1) kinematic hotter components are generally more concentrated and rounder than colder components, and (2) all components become more concentrated and thicker/rounder in more massive galaxies; they change from disk-like in low mass late-type galaxies to bulge-like in high-mass early type galaxies. Our findings suggest that Sersic \textit{n} is not a good discriminator between rotating bulges and non-rotating bulges. The luminosity fraction of cold orbits $f_{\rm cold}$ is well correlated with the photometrically-decomposed disk fraction $f_{\rm disk}$ as $f_{\mathrm{cold}} = 0.14 + 0.23f_{\mathrm{\mathrm{disk}}}$. Similarly, the hot orbit fraction $f_{\rm hot}$ is correlated with the bulge fraction $f_{\rm bulge}$ as $f_{\mathrm{hot}} = 0.19 + 0.31f_{\mathrm{\mathrm{bulge}}}$. The warm orbits mainly contribute to disks in low-mass late-type galaxies, and to bulges in high-mass early-type galaxies. The cold, warm, and hot components generally follow the same morphology ($\epsilon = 1-q_{\rm Rmax}$) versus kinematics ($\sigma_z^2/\overline{V_{\mathrm{tot}}^2}$) relation as the thin disk, thick disk/pseudo bulge, and classical bulge identified from cosmological simulations.
[36]  oai:arXiv.org:1805.04978  [pdf] - 1751819
NSVS 01286630: A Detached Binary with a Close-in Companion
Comments: 12 pages, 5 figures, accepted by RAA
Submitted: 2018-05-13
New photometric observations of NSVS 01286630 were performed and two sets of four-color (B, V , Rc, Ic) light curves (LCs) were obtained. Using the 2013 version of the Wilson-Devinney (W-D) code, we analyzed these data. The photometric solutions reveal that NSVS 01286630 is an active detached eclipsing binary (EB) with a high orbit inclination. Remarkably, the temperature of the primary component (the hotter star) is higher than the secondary one, but, the value of mass ratio q(M2/M1) for NSVS 0128663 is more than 1, which can be explained that the surface of the secondary component of NSVS 01286630 is covered with big cool-star-spots. Based on our new CCD mid-eclipse times and the data published until now, the variations of the mid-eclipse times were reanalyzed in detail using a weighted least-squares method. It is discovered that the (O-C) diagram of the system shows a cyclic oscillation with a period of 3.61 yr and an amplitude of 0.001 days. The cyclic variation may be caused by the light travel time effect (LTTE) due to the presence of a third companion, whose mass we calculated is about 0.11 solar mass. The third body may affect the orbital evolution of the central binary system by transferring the angular momentum.
[37]  oai:arXiv.org:1804.09313  [pdf] - 1686875
A First Look at BISTRO Observations of The $\rho$ Oph-A core
Kwon, Jungmi; Doi, Yasuo; Tamura, Motohide; Matsumura, Masafumi; Pattle, Kate; Berry, David; Sadavoy, Sarah; Matthews, Brenda C.; Ward-Thompson, Derek; Hasegawa, Tetsuo; Furuya, Ray S.; Pon, Andy; Di Francesco, James; Arzoumanian, Doris; Hayashi, Saeko S.; Kawabata, Koji S.; Onaka, Takashi; Choi, Minho; Kang, Miju; Hoang, Thiem; Lee, Chang Won; Lee, Sang-Sung; Liu, Hong-Li; Liu, Tie; Inutsuka, Shu-Ichiro; Eswaraiah, Chakali; Bastien, Pierre; Kwon, Woojin; Lai, Shih-Ping; Qiu, Keping; Coude, Simon; Franzmann, Erica; Friberg, Per; Graves, Sarah F.; Greaves, Jane S.; Houde, Martin; Johnstone, Doug; Kirk, Jason M.; Koch, Patrick M.; Li, Di; Parsons, Harriet; Rao, Ramprasad; Rawlings, Mark; Shinnaga, Hiroko; Van Loo, Sven; Aso, Yusuke; Byun, Do-Young; Chen, Huei-Ru; Chen, Mike C. -Y.; Chen, Wen Ping; Ching, Tao-Chung; Cho, Jungyeon; Chrysostomou, Antonio; Chung, Eun Jung; Drabek-Maunder, Emily; Eyres, Stewart P. S.; Fiege, Jason; Friesen, Rachel K.; Fuller, Gary; Gledhill, Tim; Griffin, Matt J.; Gu, Qilao; Hatchell, Jennifer; Holland, Wayne; Inoue, Tsuyoshi; Iwasaki, Kazunari; Jeong, Il-Gyo; Kang, Ji-Hyun; Kang, Sung-Ju; Kemper, Francisca; Kim, Gwanjeong; Kim, Jongsoo; Kim, Kee-Tae; Kim, Kyoung Hee; Kim, Mi-Ryang; Kim, Shinyoung; Lacaille, Kevin M.; Lee, Jeong-Eun; Li, Dalei; Li, Hua-Bai; Liu, Junhao; Liu, Sheng-Yuan; Lyo, A-Ran; Mairs, Steve; Moriarty-Schieven, Gerald H.; Nakamura, Fumitaka; Nakanishi, Hiroyuki; Ohashi, Nagayoshi; Peretto, Nicolas; Pyo, Tae-Soo; Qian, Lei; Retter, Brendan; Richer, John; Rigby, Andrew; Robitaille, Jean-Franois; Savini, Giorgio; Scaife, Anna M. M.; Soam, Archana; Tang, Ya-Wen; Tomisaka, Kohji; Wang, Hongchi; Wang, Jia-Wei; Whitworth, Anthony P.; Yen, Hsi-Wei; Yoo, Hyunju; Yuan, Jinghua; Zhang, Chuan-Peng; Zhang, Guoyin; Zhou, Jianjun; Zhu, Lei; Andre, Philippe; Dowell, C. Darren; Falle, Sam; Tsukamoto, Yusuke; Nakagawa, Takao; Kanamori, Yoshihiro; Kataoka, Akimasa; Kobayashi, Masato I. N.; Nagata, Tetsuya; Saito, Hiro; Seta, Masumichi; Zenko, Tetsuya
Comments: 22 pages, 11 figures, accepted for publication in ApJ
Submitted: 2018-04-24
We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.
[38]  oai:arXiv.org:1804.06795  [pdf] - 1743695
The Fornax3D project: overall goals, galaxy sample, MUSE data analysis and initial results
Comments: Accepted for publication in Astronomy & Astrophysics. An high resolution version is available at the following link: https://www.dropbox.com/s/bext14c5udhln8n/F3D_paper.pdf?dl=0
Submitted: 2018-04-18
The Fornax cluster provides a uniquely compact laboratory to study the detailed history of early-type galaxies and the role played by environment in driving their evolution and their transformation from late-type galaxies. Using the superb capabilities of the Multi Unit Spectroscopic Explorer on the Very Large Telescope, high-quality integral-field spectroscopic data were obtained for the inner regions of all the bright ($m_B\leq15$) galaxies within the virial radius of Fornax. The stellar haloes of early-type galaxies are also covered out to about four effective radii. State-of-the-art stellar dynamical and population modelling allows to aim in particular at better characterising the disc components of fast-rotating early-type galaxies, constraining radial variations in the stellar initial-mass functions and measuring the stellar age, metallicity, and $\alpha$-element abundance of stellar haloes in cluster galaxies. This paper describes the sample selection, observations, and overall goals of the survey, and provides initial results based on the spectroscopic data, including the detailed characterisation of stellar kinematics and populations to large radii; decomposition of galaxy components directly via their orbital structure; the ability to identify globular clusters and planetary nebulae, and derivation of high-quality emission-line diagnostics in the presence of complex ionised gas.
[39]  oai:arXiv.org:1803.11083  [pdf] - 1721058
Photometric Analysis of the eclipsing Polar MN Hya
Comments:
Submitted: 2018-03-29
As an eclipsing polar with 3.39 hrs orbital period, MN Hya was going through state change when we observed it during 2009-2016. 10 new mid-eclipse times, along with others obtained from literature, allow us to give a new ephemeris. The residuals of linear fit show that period decreased during the phase of state change. It means angular momentum was lost during this phase. The X-ray observation indicates the mass accretion rate as about $3.6\times10^{-9}M_{\odot}yr^{-1}$. The period decrease gives that at least 60 percent of mass being transfered from secondary was lost, maybe in form of the spherically symmetric stellar wind. In high state, the data shows the intensity of the flickering reduced when system had higher accretion rate, and that flickering sticks out with primary timescale about 2 minutes, which implies the position of the threading point as about 30 radius of the white dwarf above the surface of it. The trend of light curves of the system in high state follows that of low state for a large fraction of phase interval from phase 0 to phase 0.4 since which the cyclotron feature is visible, and the primary intensity hump of light curves near phase 0.7 when the system is in high state did not appear on the curve when it is in low state. Those facts contradict the predictions of the two-pole model.
[40]  oai:arXiv.org:1803.02259  [pdf] - 1648953
The EDGE-CALIFA survey: validating stellar dynamical mass models with CO kinematics
Comments: 22 pages, 18 figures, Accepted for publication in MNRAS
Submitted: 2018-03-06
Deriving circular velocities of galaxies from stellar kinematics can provide an estimate of their total dynamical mass, provided a contribution from the velocity dispersion of the stars is taken into account. Molecular gas (e.g., CO) on the other hand, is a dynamically cold tracer and hence acts as an independent circular velocity estimate without needing such a correction. In this paper we test the underlying assumptions of three commonly used dynamical models, deriving circular velocities from stellar kinematics of 54 galaxies (S0-Sd) that have observations of both stellar kinematics from the CALIFA survey, and CO kinematics from the EDGE survey. We test the Asymmetric Drift Correction (ADC) method, as well as Jeans, and Schwarzschild models. The three methods each reproduce the CO circular velocity at 1Re to within 10%. All three methods show larger scatter (up to 20%) in the inner regions (R < 0.4Re) which may be due to an increasingly spherical mass distribution (which is not captured by the thin disk assumption in ADC), or non-constant stellar M/L ratios (for both the JAM and Schwarzschild models). This homogeneous analysis of stellar and gaseous kinematics validates that all three models can recover Mdyn at 1Re to better than 20%, but users should be mindful of scatter in the inner regions where some assumptions may break down.
[41]  oai:arXiv.org:1711.06728  [pdf] - 1613949
The stellar orbit distribution in present-day galaxies inferred from the CALIFA survey
Comments: Published in Nature Astronomy, 1 January 2018; doi:10.1038/s41550-017-0348-1; 22pages, 8 figures, 2 tables
Submitted: 2017-11-17, last modified: 2018-01-07
Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.
[42]  oai:arXiv.org:1707.04006  [pdf] - 1626267
LAMOST Views $\delta$ Scuti Pulsating Stars
Comments: 15 pages, 16 figures, accepted for publication in MNRAS
Submitted: 2017-07-13, last modified: 2017-12-27
About 766 $\delta$ Scuti stars were observed by LAMOST by June 16, 2017. Stellar atmospheric parameters of 525 variables were determined, while spectral types were obtained for all samples. In the paper, those spectroscopic data are catalogued. We detect a group of 131 unusual and cool variable stars (UCVs) that are distinguished from the normal $\delta$ Scuti stars (NDSTs). On the H-R diagram and the $\log g-T$ diagram, the UCVs are far beyond the red edge of pulsational instability trip. Their effective temperatures are lower than 6700\,K with periods in the range from 0.09 to 0.22\,days. NDSTs have metallicity close to that of the Sun as expected, while UCVs are slightly metal poor than NDSTs. The two peaks on the distributions of the period and stellar atmospheric parameters are all caused by the existence of UCVs. When those UCVs are excluded, it is discovered that the effective temperature, the gravitational acceleration and the metallicity all are correlated with the pulsating period for NDSTs and their physical properties and evolutionary states are discussed. The detecting of those UCVs indicates that about 25\% of the known $\delta$ Scuti stars may be misclassified. However, if some of them are confirmed to be pulsating stars, they will be a new-type pulsator and their investigations will shed light on theoretical instability domains and on the theories of interacting between the pulsation and the convection of solar-type stars. Meanwhile, 88 $\delta$ Scuti stars are detected to be the candidates of binary or multiple systems.
[43]  oai:arXiv.org:1712.09026  [pdf] - 1648771
Physical Properties and Evolutionary States of EA-type Eclipsing Binaries Observed by LAMOST
Comments: 32 pages, 12 figures, accepted for publication in ApJS
Submitted: 2017-12-25
About 3196 EA-type binaries (EAs) were observed by LAMOST by June 16, 2017 and their spectral types were derived. Meanwhile stellar atmospheric parameters of 2020 EAs were determined. In the paper, those EAs are catalogued and their physical properties and evolutionary states are investigated. The period distribution of EAs suggests that the period limit of tidal locking for the close binaries is about 6 days. It is found that the metallicity of EAs is higher than that of EWs indicating that EAs are generally younger than EWs and they are the progenitors of EWs. The metallicities of long-period EWs (0.4 < P < 1 days) are the same as those of EAs with the same periods, while their values of Log (g) are usually smaller than those of EAs. These support the evolutionary process that EAs evolve into long-period EWs through the combination of angular momentum loss (AML) via magnetic braking and case A mass transfer. For short-period EWs, their metallicities are lower than those of EAs, while their gravitational accelerations are higher. These reveal that they may be formed from cool short-period EAs through AML via magnetic braking with little mass transfer. For some EWs with high metallicities, they may be contaminated by material from the evolution of unseen neutron stars and black holes or they have third bodies that may help them to form rapidly through a short timescale of pre-contact evolution. The present investigation suggests that the modern EW populations may be formed through the combination of aforementioned mechanisms.
[44]  oai:arXiv.org:1712.09154  [pdf] - 1659643
First photometric study of ultrashort-period contact binary 1SWASP J140533.33+114639.1
Comments:
Submitted: 2017-12-25
In this paper, CCD photometric light curves for the short-period eclipsing binary 1SWASP J140533.33+114639.1 (hereafter J1405) in the $BVR$ bands are presented and analyzed using the 2013 version of the Wilson-Devinney (W-D) code. It is discovered that the J1405 is a W-subtype shallow contact binary with a contact degree of $f$ = 7.9$\pm$ 0.5\% and a mass ratio of $q$ = 1.55 $\pm$ 0.02. In order to explain the asymmetric light curves of the system, a cool star-spot on the more massive component was employed. This shallow contact eclipsing binary may be formed from a short-period detached system through the orbital shrinkage due to angular momentum loss. Based on $(O-C)$ method, the variation of the orbital period was studied using all the available times of the minimum light. The $(O-C)$ diagram reveals that the period is increasing continuously at a rate of $dP/dt=+2.09\times{10^{-7}}$days yr$^{-1}$, which can be explained by mass transfer from the less massive component to the more massive one.
[45]  oai:arXiv.org:1712.01072  [pdf] - 1648724
Multi-color light curves and orbital period research of eclipsing binary V1073 Cyg
Comments: 15 pages, 5 figures
Submitted: 2017-12-04
New Multi-color $B$ $V$ $R_c$ $I_c$ photometric observation are presented for W UMa type eclipsing binary V1073 Cyg. The multi-color light curves analysis with the Wilson-Devinney(W-D) procedure acquired the absolute parameters of this system, showing that V1073Cyg is a shallow contact binary system with fill-out factor $f=0.124(\pm0.011)$. We collected all available times of light minima spanning 119 years including CCD data to construct the O-C curve and made detailed O-C analysis. The O-C diagram shows that the period change is complex. There exist a long-term continuous decrease and a cyclic variation. The period is decreasing at a rate of $\dot P=-1.04(\pm0.18)\times 10^{-10}$ $days\cdot{cycle}^{-1}$, and with the period decrease, V1073 Cyg will evolve to deep contact stage. The cyclic variation with a period of $P_3=82.7(\pm3.6) years$ and an amplitude of $A=0.028(\pm0.002) day$ may be explained by the magnetic activity of one or both components or the light travel time effect (LTTE) caused by a distant third companion with $M_3({i'}=90^{\circ})=0.511M_\odot$.
[46]  oai:arXiv.org:1710.09349  [pdf] - 1790571
Morpho-kinematic properties of field S0 bulges in the CALIFA survey
Comments: 31 pages, 19 figures. Accepted for publication in MNRAS
Submitted: 2017-10-25
We study a sample of 28 S0 galaxies extracted from the integral-field spectroscopic (IFS) survey CALIFA. We combine an accurate two-dimensional (2D) multi-component photometric decomposition with the IFS kinematic properties of their bulges to understand their formation scenario. Our final sample is representative of S0s with high stellar masses ($M_{star}/M_{sun} > 10^{10}$). They lay mainly on the red sequence and live in relatively isolated environments similar to that of the field and loose groups. We use our 2D photometric decomposition to define the size and photometric properties of the bulges, as well as their location within the galaxies. We perform mock spectroscopic simulations mimicking our observed galaxies to quantify the impact of the underlying disc on our bulge kinematic measurements ($\lambda$ and $v/\sigma$). We compare our bulge corrected kinematic measurements with the results from Schwarzschild dynamical modelling. The good agreement confirms the robustness of our results and allows us to use bulge reprojected values of $\lambda$ and $v/\sigma$. We find that the photometric ($n$ and $B/T$) and kinematic ($v/\sigma$ and $\lambda$) properties of our field S0 bulges are not correlated. We demonstrate that this morpho-kinematic decoupling is intrinsic to the bulges and it is not due to projection effects. We conclude that photometric diagnostics to separate different types of bulges (disc-like vs classical) might not be useful for S0 galaxies. The morpho-kinematics properties of S0 bulges derived in this paper suggest that they are mainly formed by dissipation processes happening at high redshift, but dedicated high-resolution simulations are necessary to better identify their origin.
[47]  oai:arXiv.org:1709.06649  [pdf] - 1593628
Orbital decomposition of CALIFA spiral galaxies
Comments: accepted by MNRAS
Submitted: 2017-09-19
Schwarzschild orbit-based dynamical models are widely used to uncover the internal dynamics of early-type galaxies and globular clusters. Here we present for the first time the Schwarzschild models of late-type galaxies: an SBb galaxy NGC 4210 and an S0 galaxy NGC 6278 from the CALIFA survey. The mass profiles within $2\,R_e$ are constrained well with $1\sigma$ statistical error of $\sim 10\%$. The luminous and dark mass can be disentangled with uncertainties of $\sim 20\%$ and $\sim 50\%$ respectively. From $R_e$ to $2\,R_e$, the dark matter fraction increases from $14\pm10\%$ to $18\pm10\%$ for NGC 4210 and from $15\pm10\%$ to $30\pm20\%$ for NGC 6278. The velocity anisotropy profiles of both $\sigma_r/\sigma_t$ and $\sigma_z/\sigma_R$ are well constrained. The inferred internal orbital distributions reveal clear substructures. The orbits are naturally separated into three components: a cold component with near circular orbits, a hot component with near radial orbits, and a warm component in between. The photometrically-identified exponential disks are predominantly made up of cold orbits only beyond $\sim 1\,R_e$, while they are constructed mainly with the warm orbits inside. Our dynamical hot components are concentrated in the inner regions, similar to the photometrically-identified bulges. The reliability of the results, especially the orbit distribution, are verified by applying the model to mock data.
[48]  oai:arXiv.org:1709.02596  [pdf] - 1588117
A New Stellar Outburst Associated with the Magnetic Activities of the K-type Dwarf in a White-dwarf Binary
Comments: 16 pages, 6 figures, accepted for publication in ApJ
Submitted: 2017-09-08
1SWASP\,J162117.36$+$441254.2 was originally classified as an EW-type binary with a period of 0.20785\,days. However, it was detected to have undergone a stellar outburst on June 3, 2016. Although the system was latter classified as a cataclysmic variable (CV) and the event was attributed as a dwarf-nova outburst, the physical reason is still unknown. This binary has been monitored photometrically since April 19, 2016 and many light curves were obtained before, during and after the outburst. Those light and color curves observed before the outburst indicate that the system is a special CV. The white dwarf is not accreting material from the secondary and there are no accretion disks surrounding the white dwarf. By comparing the light curves obtained from April 19 to September 14, 2016, it was found that magnetic activity of the secondary is associated with the outburst. We show strong evidence that the $L_1$ region on the secondary was heavily spotted before and after the outburst and thus quench the mass transfer, while the outburst is produced by a sudden mass accretion of the white dwarf. These results suggest that J162117 is a good astrophysical laboratory to study stellar magnetic activity and its influences on CV mass transfer and mass accretion.
[49]  oai:arXiv.org:1705.03996  [pdf] - 1583324
Physical Properties and Catalogue of EW-type Eclipsing Binaries Observed by LAMOST
Comments: 20 pages, 16 figures; Accepted for publication in Research in Astronomy and Astrophysics (RAA)
Submitted: 2017-05-10, last modified: 2017-06-13
Numerous EWs were discovered by several deep photometric survey and there are about 40785 EW-type binary systems listed in the international variable star index (VSX) by March 13, 2017. 7938 of them were observed by LAMOST by November 30, 2016 and their spectral types were given. Stellar atmospheric parameters of 5363 EW-type binary stars were determined based on good spectroscopic observations. In the paper, those EWs were catalogued and their properties are analyzed. The distributions of the orbital period (P), the effect temperature (T), the gravitational acceleration (Log(g)), the metallicity ([Fe/H]) and the radial velocity (RV) are presented for those observed EW-type systems. It is shown that about 80.6\% sample stars have metallicity below zero indicating that EW-type systems are old stellar population. This is in agreement with the the conclusion that the EW binaries are formed from moderately close binaries through angular momentum loss via magnetic braking that takes a few hundred million to a few billion years. The unusual high metallicities of a few percent of EWs may be caused by contaminating of material from the evolution of unseen neutron stars and black holes in the systems. The correlations between the orbital period and the effect temperature, the gravitational acceleration and the metallicity are presented and their scatters are mainly caused by (i) the presence of the third bodies and (ii) the wrong determined periods sometimes. It is shown that some EW contain evolved component stars and the physical properties of EWs are mainly depending on their orbital periods. It is found that the extremely short-period EWs may be older than their long-period cousins because they have lower metallicities. This reveals that they have a longer timescale of pre-contact evolution and their formation and evolution are mainly driven by angular momentum loss via magnetic braking.
[50]  oai:arXiv.org:1706.03018  [pdf] - 1584478
ALMA [NII] 205 micron Imaging Spectroscopy of the Interacting Galaxy System BRI 1202-0725 at Redshift 4.7
Comments: 4 pages, 1 table, 4 figures; accepted for publication in the ApJ Letters
Submitted: 2017-06-09
We present the results from Atacama Large Millimeter/submillimeter Array (ALMA) imaging in the [NII] 205 micron fine-structure line (hereafter [NII]) and the underlying continuum of BRI 1202-0725, an interacting galaxy system at $z =$ 4.7, consisting of an optical QSO, a sub-millimeter galaxy (SMG) and two Lyman-$\alpha$ emitters (LAEs), all within $\sim$25 kpc of the QSO. We detect the QSO and SMG in both [NII] and continuum. At the $\sim$$1"$ (or 6.6 kpc) resolution, both QSO and SMG are resolved in [NII], with the de-convolved major axes of $\sim$9 and $\sim$14 kpc, respectively. In contrast, their continuum emissions are much more compact and unresolved even at an enhanced resolution of $\sim$$0.7"$. The ratio of the [NII] flux to the existing CO (7$-$6) flux is used to constrain the dust temperature ($T_{\rm dust}$) for a more accurate determination of the FIR luminosity $L_{\rm FIR}$. Our best estimated $T_{\rm dust}$ equals $43 (\pm 2)$ K for both galaxies (assuming an emissivity index $\beta = 1.8$). The resulting $L_{\rm CO(7-6)}/L_{\rm FIR}$ ratios are statistically consistent with that of local luminous infrared galaxies, confirming that $L_{\rm CO(7-6)}$ traces the star formation (SF) rate (SFR) in these galaxies. We estimate that the on-going SF of the QSO (SMG) has a SFR of 5.1 $(6.9) \times 10^3 M_{\odot}$ yr$^{-1}$ ($\pm$ 30%) assuming Chabrier initial mass function, takes place within a diameter (at half maximum) of 1.3 (1.5) kpc, and shall consume the existing 5 $(5) \times 10^{11} M_{\odot}$ of molecular gas in 10 $(7) \times 10^7$ years.
[51]  oai:arXiv.org:1705.03163  [pdf] - 1583243
Long-term photometric behaviour of the eclipsing Z Cam-type dwarf nova AY Psc
Comments: 12 pages, 10 figures
Submitted: 2017-05-08
New eclipse timings of the Z Cam-type dwarf nova AY Psc were measured and the orbital ephemeris was revised. Based on the long-term AAVSO data, moreover, the outburst behaviors were also explored. Our analysis suggests that the normal outbursts are quasi-periodic, with an amplitude of $\sim2.5(\pm0.1)$ mag and a period of $\sim18.3(\pm0.7)$ days. The amplitude vs. recurrence-time relation of AY Psc is discussed, and we concluded that this relation may represents general properties of dwarf nova (DN) outbursts. The observed standstill ends with an outburst, which is inconsistent with the general picture of Z Cam-type stars. This unusual behavior was considered to be related to the mass-transfer outbursts. Moreover, the average luminosity is brighter during standstills than during outburst cycles. The changes in brightness marks the variations in $\dot{M}_{2}$ due to the disc of AY Psc is nearly steady state. $\dot{M}_{2}$ value was limited to the range from $6.35\times10^{-9}$ to $1.18\times10^{-8}M_{\odot}yr^{-1}$. More detailed examination shows that there are a few small outbursts presence during standstills. These events with amplitudes of $\sim0.5-0.9$ mag are very similar to the stunted outbursts reported in some NLs. We discussed several possible mechanisms and suggested that the most reasonable mechanism for these stunted outbursts is a changing mass-transfer rate.
[52]  oai:arXiv.org:1705.03160  [pdf] - 1574798
Physical parameters and orbital period variation of a newly discovered cataclysmic variable GSC 4560-02157
Comments: 10 pages, 7 figures
Submitted: 2017-05-08
GSC 4560-02157 is a new eclipsing cataclysmic variable with an orbital period of $0.265359$ days. By using the published $V-$ and $R-$band data together with our observations, we discovered that the $O-C$ curve of GSC 4560-02157 may shows a cyclic variation with the period of $3.51$ years and an amplitude of $1.40$ min. If this variation is caused by a light travel-time effect via the existence of a third body, its mass can be derived as $M_{3}sini'\approx91.08M_{Jup}$, it should be a low-mass star. In addition, several physical parameters were measured. The colour of the secondary star was determined as $V-R=0.77(\pm0.03)$ which corresponds to a spectral type of K2-3. The secondary star's mass was estimated as $M_{2}=0.73(\pm0.02)M_{\odot}$ by combing the derived $V-R$ value around phase 0 with the assumption that it obeys the mass-luminosity relation of the main sequence stars. This mass is consistent with the mass$-$period relation of CV donor stars. For the white dwarf, the eclipse durations and contacts of the white dwarf yield an upper limit of the white dwarf's radius corresponding to a lower limit mass of $M_{1}\approx0.501M_{\odot}$. The overestimated radius and previously published spectral data indicate that the boundary layer may has a very high temperature.
[53]  oai:arXiv.org:1705.03155  [pdf] - 1578111
WZ Sge: an eclipsing cataclysmic variable evolving towards the period minimum
Comments: 8 pages, 5 figures
Submitted: 2017-05-08
We present the photometric results of the eclipsing cataclysmic variable (CV) WZ Sge near the period minimum ($P_{min}$). Eight new mid-eclipse times were determined and the orbital ephemeris was updated. Our result shows that the orbital period of WZ Sge is decreasing at a rate of $\dot{P}=-2.72(\pm0.23)\times{10^{-13}}\,s s^{-1}$. This secular decrease, coupled with previous detection of its donor, suggest that WZ Sge is a pre-bounce system. Further analysis indicates that the observed period decrease rate is about $1.53$ times higher than pure gravitational radiation (GR) driving. We constructed the evolutionary track of WZ Sge, which predicts that $P_{min}$ of WZ Sge is $\sim77.98 (\pm0.90)$ min. If the orbital period decreases at the current rate, WZ Sge will evolve past its $P_{min}$ after $\sim25.3$ Myr. Based on the period evolution equation we find $\dot{M}_{2}\simeq4.04(\pm0.10)\times10^{-11}M_{\odot}yr^{-1}$, which is compatible with the current concept of CV evolution at ultrashort orbital periods.
[54]  oai:arXiv.org:1705.03145  [pdf] - 1578110
Cyclic period oscillation of the eclipsing dwarf nova DV UMa
Comments: 16 pages, 6 figures
Submitted: 2017-05-08
DV UMa is an eclipsing dwarf nova with an orbital period of $\sim2.06$ h, which lies just at the bottom edge of the period gap. To detect its orbital period changes we present 12 new mid-eclipse times by using our CCD photometric data and archival data. Combining with the published mid-eclipse times in quiescence, spanning $\sim30$ yr, the latest version of the $O-C$ diagram was obtained and analyzed. The best fit to those available eclipse timings shows that the orbital period of DV UMa is undergoing a cyclic oscillation with a period of $17.58(\pm0.52)$ yr and an amplitude of $71.1(\pm6.7)$ s. The periodic variation most likely arises from the light-travel-time effect via the presence of a circumbinary object because the required energy to drive the Applegate mechanism is too high in this system. The mass of the unseen companion was derived as $M_{3}\sin{i'}=0.025(\pm0.004)M_{\odot}$. If the third body is in the orbital plane (i.e. $i'=i=82.9^{\circ}$) of the eclipsing pair, it would match to a brown dwarf. This hypothetical brown dwarf is orbiting its host star at a separation of $\sim8.6$ AU in an eccentric orbit ($e=0.44$).
[55]  oai:arXiv.org:1705.03166  [pdf] - 1583244
Double cyclic variations in orbital period of the eclipsing cataclysmic variable EX Dra
Comments: 8 pages, 5 figures
Submitted: 2017-05-08
EX Dra is a long-period eclipsing dwarf nova with $\sim2-3$ mag amplitude outbursts. This star has been monitored photometrically from November, 2009 to March, 2016 and 29 new mid-eclipse times were obtained. By using new data together with the published data, the best fit to the $O-C$ curve indicate that the orbital period of EX Dra have an upward parabolic change while undergoing double-cyclic variations with the periods of 21.4 and 3.99 years, respectively. The upward parabolic change reveals a long-term increase at a rate of $\dot{P}={+7.46}\times10^{-11}{s} {s^{-1}}$. The evolutionary theory of cataclysmic variables (CVs) predicts that, as a CV evolves, the orbital period should be decreasing rather than increasing. Secular increase can be explained as the mass transfer between the secondary and primary or may be just an observed part of a longer cyclic change. Most plausible explanation for the double-cyclic variations is a pair of light travel-time effect via the presence of two companions. Their masses are determined to be $M_{A}sini'_{A}=29.3(\pm0.6) M_{Jup}$ and $M_{B}sini'_{B}=50.8(\pm0.2) M_{Jup}$. When the two companions are coplanar to the orbital plane of the central eclipsing pair, their masses would match to brown dwarfs.
[56]  oai:arXiv.org:1704.08552  [pdf] - 1582919
First results from BISTRO -- a SCUBA-2 polarimeter survey of the Gould Belt
Ward-Thompson, Derek; Pattle, Kate; Bastien, Pierre; Furuya, Ray S.; Kwon, Woojin; Lai, Shih-Ping; Qiu, Keping; Berry, David; Choi, Minho; Coudé, Simon; Di Francesco, James; Hoang, Thiem; Franzmann, Erica; Friberg, Per; Graves, Sarah F.; Greaves, Jane S.; Houde, Martin; Johnstone, Doug; Kirk, Jason M.; Koch, Patrick M.; Kwon, Jungmi; Lee, Chang Won; Li, Di; Matthews, Brenda C.; Mottram, Joseph C.; Parsons, Harriet; Pon, Andy; Rao, Ramprasad; Rawlings, Mark; Shinnaga, Hiroko; Sadavoy, Sarah; van Loo, Sven; Aso, Yusuke; Byun, Do-Young; Chakali, Eswariah; Chen, Huei-Ru; Chen, Mike C. -Y.; Chen, Wen Ping; Ching, Tao-Chung; Cho, Jungyeon; Chrysostomou, Antonio; Chung, Eun Jung; Doi, Yasuo; Drabek-Maunder, Emily; Eyres, Stewart P. S.; Fiege, Jason; Friesen, Rachel K.; Fuller, Gary; Gledhill, Tim; Griffin, Matt J.; Gu, Qilao; Hasegawa, Tetsuo; Hatchell, Jennifer; Hayashi, Saeko S.; Holland, Wayne; Inoue, Tsuyoshi; Inutsuka, Shu-ichiro; Iwasaki, Kazunari; Jeong, Il-Gyo; Kang, Ji-hyun; Kang, Miju; Kang, Sung-ju; Kawabata, Koji S.; Kemper, Francisca; Kim, Gwanjeong; Kim, Jongsoo; Kim, Kee-Tae; Kim, Kyoung Hee; Kim, Mi-Ryang; Kim, Shinyoung; Lacaille, Kevin M.; Lee, Jeong-Eun; Lee, Sang-Sung; Li, Dalei; Li, Hua-bai; Liu, Hong-Li; Liu, Junhao; Liu, Sheng-Yuan; Liu, Tie; Lyo, A-Ran; Mairs, Steve; Matsumura, Masafumi; Moriarty-Schieven, Gerald H.; Nakamura, Fumitaka; Nakanishi, Hiroyuki; Ohashi, Nagayoshi; Onaka, Takashi; Peretto, Nicolas; Pyo, Tae-Soo; Qian, Lei; Retter, Brendan; Richer, John; Rigby, Andrew; Robitaille, Jean-François; Savini, Giorgio; Scaife, Anna M. M.; Soam, Archana; Tamura, Motohide; Tang, Ya-Wen; Tomisaka, Kohji; Wang, Hongchi; Wang, Jia-Wei; Whitworth, Anthony P.; Yen, Hsi-Wei; Yoo, Hyunju; Yuan, Jinghua; Zhang, Chuan-Peng; Zhang, Guoyin; Zhou, Jianjun; Zhu, Lei; André, Philippe; Dowell, C. Darren; Falle, Sam; Tsukamoto, Yusuke
Comments: 12 pages, 8 figures, accepted for publication in ApJ
Submitted: 2017-04-27
We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.
[57]  oai:arXiv.org:1701.01583  [pdf] - 1534216
Triaxial orbit-based modelling of the Milky Way Nuclear Star Cluster
Comments: 14 pages, 10 figures. Accepted for publication in MNRAS
Submitted: 2017-01-06
We construct triaxial dynamical models for the Milky Way nuclear star cluster using Schwarzschild's orbit superposition technique. We fit the stellar kinematic maps presented in Feldmeier et al. (2014). The models are used to constrain the supermassive black hole mass M_BH, dynamical mass-to-light ratio M/L, and the intrinsic shape of the cluster. Our best-fitting model has M_BH = (3.0 +1.1 -1.3)x10^6 M_sun, M/L = (0.90 +0.76 -0.08) M_sun/L_{sun,4.5micron}, and a compression of the cluster along the line-of-sight. Our results are in agreement with the direct measurement of the supermassive black hole mass using the motion of stars on Keplerian orbits. The mass-to-light ratio is consistent with stellar population studies of other galaxies in the mid-infrared. It is possible that we underestimate M_BH and overestimate the cluster's triaxiality due to observational effects. The spatially semi-resolved kinematic data and extinction within the nuclear star cluster bias the observations to the near side of the cluster, and may appear as a compression of the nuclear star cluster along the line-of-sight. We derive a total dynamical mass for the Milky Way nuclear star cluster of M_MWNSC = (2.1 +-0.7)x10^7 M_sun within a sphere with radius r = 2 x r_eff = 8.4 pc. The best-fitting model is tangentially anisotropic in the central r = 0.5-2 pc of the nuclear star cluster, but close to isotropic at larger radii. Our triaxial models are able to recover complex kinematic substructures in the velocity map.
[58]  oai:arXiv.org:1611.04699  [pdf] - 1521448
A low-mass-ratio and deep contact binary as the progenitor of the merger V1309 Sco
Comments: 7 pages, 4 figures
Submitted: 2016-11-14
Nova Sco 2008 (=V1309 Sco) is an example of a V838 Mon type eruption rather than a typical classical nova. This enigmatic object was recently shown to have resulted from the merger of two stars in a contact binary. It is the first stellar merger that was identified to be undergoing a common envelope transient. To understand the properties of its binary progenitor, the pre-outburst light curves were analyzed by using the W-D method. The photometric solution of the 2002 light curve shows that it is a deep contact binary (f = 89.5(+-40.5)%) with a mass ratio of 0.094. The asymmetry of the light curve is explained by the presence of a dark spot on the more massive component. The extremely high fill-out factor suggests that the merging of the contact binary is driven by dynamical mass loss from the outer Lagrange point. However, the analysis of the 2004 light curve indicates that no solutions were obtained even at an extremely low mass ratio of q = 0.03. This suggests that the common convective envelope of the binary system disappeared and the secondary component spiraled into the envelope of the primary in 2004. Finally, the ejection of the envelope of the primary produced the outburst.
[59]  oai:arXiv.org:1609.06446  [pdf] - 1531678
Stellar kinematics across the Hubble sequence in the CALIFA survey: General properties and aperture corrections
Comments: 52 pages, 44 figures, A&A accepted. FITS tables with the stellar kinematic results are available at: http://www.caha.es/CALIFA/public_html/?q=content/science-dataproducts A full-resolution version of the manuscript is also available on the same webpage
Submitted: 2016-09-21, last modified: 2016-09-22
We present the stellar kinematic maps of a large sample of galaxies from the integral-field spectroscopic survey CALIFA. The sample comprises 300 galaxies displaying a wide range of morphologies across the Hubble sequence, from ellipticals to late-type spirals. This dataset allows us to homogeneously extract stellar kinematics up to several effective radii. In this paper, we describe the level of completeness of this subset of galaxies with respect to the full CALIFA sample, as well as the virtues and limitations of the kinematic extraction compared to other well-known integral-field surveys. In addition, we provide averaged integrated velocity dispersion radial profiles for different galaxy types, which are particularly useful to apply aperture corrections for single aperture measurements or poorly resolved stellar kinematics of high-redshift sources. The work presented in this paper sets the basis for the study of more general properties of galaxies that will be explored in subsequent papers of the survey.
[60]  oai:arXiv.org:1608.08239  [pdf] - 1475522
A discrete chemo-dynamical model of the dwarf spheroidal galaxy Sculptor: mass profile, velocity anisotropy and internal rotation
Comments: Accepted by MNRAS
Submitted: 2016-08-29
We present a new discrete chemo-dynamical axisymmetric modeling technique, which we apply to the dwarf spheroidal galaxy Sculptor. The major improvement over previous Jeans models is that realistic chemical distributions are included directly in the dynamical modelling of the discrete data. This avoids loss of information due to spatial binning and eliminates the need for hard cuts to remove contaminants and to separate stars based on their chemical properties. Using a combined likelihood in position, metallicity and kinematics, we find that our models naturally separate Sculptor stars into a metal-rich and a metal-poor population. Allowing for non-spherical symmetry, our approach provides a central slope of the dark matter density of $\gamma = 0.5 \pm 0.3$. The metal-rich population is nearly isotropic (with $\beta_r^{red} = 0.0\pm0.1$) while the metal-poor population is tangentially anisotropic (with $\beta_r^{blue} = -0.2\pm0.1$) around the half light radius of $0.26$ kpc. A weak internal rotation of the metal-rich population is revealed with $v_{max}/\sigma_0 = 0.15 \pm 0.15$. We run tests using mock data to show that a discrete dataset with $\sim 6000$ stars is required to distinguish between a core ($\gamma = 0$) and cusp ($\gamma = 1$), and to constrain the possible internal rotation to better than $1\,\sigma$ confidence with our model. We conclude that our discrete chemo-dynamical modelling technique provides a flexible and powerful tool to robustly constrain the internal dynamics of multiple populations, and the total mass distribution in a stellar system.
[61]  oai:arXiv.org:1608.08238  [pdf] - 1475521
A discrete chemo-dynamical model of the giant elliptical galaxy NGC 5846: dark matter fraction, internal rotation and velocity anisotropy out to six effective radii
Comments: Accepted by MNRAS
Submitted: 2016-08-29
We construct a suite of discrete chemo-dynamical models of the giant elliptical galaxy NGC 5846. These models are a powerful tool to constrain both the mass distribution and internal dynamics of multiple tracer populations. We use Jeans models to simultaneously fit stellar kinematics within the effective radius $R_{\rm e}$, planetary nebula (PN) radial velocities out to $3\, R_{\rm e}$, and globular cluster (GC) radial velocities and colours out to $6\,R_{\rm e}$. The best-fitting model is a cored DM halo which contributes $\sim 10\%$ of the total mass within $1\,R_{\rm e}$, and $67\% \pm 10\%$ within $6\,R_{\rm e}$, although a cusped DM halo is also acceptable. The red GCs exhibit mild rotation with $v_{\rm max}/\sigma_0 \sim 0.3$ in the region $R > \,R_{\rm e}$, aligned with but counter-rotating to the stars in the inner parts, while the blue GCs and PNe kinematics are consistent with no rotation. The red GCs are tangentially anisotropic, the blue GCs are mildly radially anisotropic, and the PNe vary from radially to tangentially anisotropic from the inner to the outer region. This is confirmed by general made-to-measure models. The tangential anisotropy of the red GCs in the inner regions could stem from the preferential destruction of red GCs on more radial orbits, while their outer tangential anisotropy -- similar to the PNe in this region -- has no good explanation. The mild radial anisotropy of the blue GCs is consistent with an accretion scenario.
[62]  oai:arXiv.org:1608.06483  [pdf] - 1531339
The Sloan Digital Sky Survey Quasar Catalog: twelfth data release
Comments: Accepted for publication in A&A. The catalog is publicly available here: http://www.sdss.org/dr12/algorithms/boss-dr12-quasar-catalog
Submitted: 2016-08-23
We present the Data Release 12 Quasar catalog (DR12Q) from the Baryon Oscillation Spectroscopic Survey (BOSS) of the SDSS-III. This catalog includes all SDSS-III/BOSS objects that were spectroscopically targeted as quasar candidates during the full survey and that are confirmed as quasars via visual inspection of the spectra, have luminosities Mi[z=2]<-20.5 (in a LCDM cosmology with H_0 = 70 km/s/Mpc, Omega_M = 0.3, and Omega _L=0.7), and either display at least one emission line with a full width at half maximum (FWHM)larger than 500 km/s or, if not, have interesting/complex absorption features. The catalog also includes previously known quasars (mostly from SDSS-I and II) that were reobserved by BOSS. The catalog contains 297,301 quasars detected over 9,376 square degrees with robust identification and redshift measured by a combination of principal component eigenspectra. The number of quasars with z>2.15 is about an order of magnitude greater than the number of z>2.15 quasars known prior to BOSS. Redshifts and FWHMs are provided for the strongest emission lines (CIV, CIII], MgII). The catalog identifies 29,580 broad absorption line quasars and lists their characteristics. For each object, the catalog presents five-band (u, g, r, i, z) CCD-based photometry together with some information on the optical morphology and the selection criteria. When available, the catalog also provides information on the optical variability of quasars using SDSS and PTF multi-epoch photometry. The catalog also contains X-ray, ultraviolet, near-infrared, and radio emission properties of the quasars, when available, from other large-area surveys. The calibrated digital spectra, covering the wavelength region 3,600-10,500A at a spectral resolution in the range 1,300<R<2,500, can be retrieved from the SDSS Catalog Archive Server.
[63]  oai:arXiv.org:1608.00390  [pdf] - 1450710
Sparse aperture differential piston measurements using the pyramid wave-front sensor
Comments: 10 pages, 9 figures, Proceeding 9909-253 of the SPIE Astronomical Telescopes + Instrumentation 2016, 26 June-1 July 2016 Edinburgh, Scotland, UK
Submitted: 2016-08-01
In this paper we report on the laboratory experiment we settled in the Shanghai Astronomical Observatory (SHAO) to investigate the pyramid wavefront sensor (WFS) ability to measure the differential piston on a sparse aperture. The ultimate goal is to verify the ability of the pyramid WFS work in closed loop to perform the phasing of the primary mirrors of a sparse Fizeau imaging telescope. In the experiment we installed on the optical bench we performed various test checking the ability to flat the wave-front using a deformable mirror and to measure the signal of the differential piston on a two pupils setup. These steps represent the background from which we start to perform full closed loop operation on multiple apertures. These steps were also useful to characterize the achromatic double pyramids (double prisms) manufactured in the SHAO optical workshop.
[64]  oai:arXiv.org:1606.07448  [pdf] - 1494735
IMF shape constraints from stellar populations and dynamics from CALIFA
Comments: 6 pages, 3 figures, accepted for publication in MNRAS Letters
Submitted: 2016-06-23
In this letter we describe how we use stellar dynamics information to constrain the shape of the stellar IMF in a sample of 27 early-type galaxies from the CALIFA survey. We obtain dynamical and stellar mass-to-light ratios, $\Upsilon_\mathrm{dyn}$ and $\Upsilon_{\ast}$, over a homogenous aperture of 0.5~$R_{e}$. We use the constraint $\Upsilon_\mathrm{dyn} \ge \Upsilon_{\ast}$ to test two IMF shapes within the framework of the extended MILES stellar population models. We rule out a single power law IMF shape for 75% of the galaxies in our sample. Conversely, we find that a double power law IMF shape with a varying high-mass end slope is compatible (within 1$\sigma$) with 95% of the galaxies. We also show that dynamical and stellar IMF mismatch factors give consistent results for the systematic variation of the IMF in these galaxies.
[65]  oai:arXiv.org:1605.04795  [pdf] - 1443970
The low dark matter content of the lenticular galaxy NGC 3998
Comments: 17 pages, 21 figures, accepted by MNRAS
Submitted: 2016-05-16
We observed the lenticular galaxy NGC 3998 with the Mitchell Integral-Field Spectrograph and extracted line-of-sight velocity distributions out to 3 half-light radii. We constructed collisionless orbit models in order to constrain NGC 3998's dark and visible structure, using kinematics from both the Mitchell and SAURON instruments. We find NGC 3998 to be almost axisymmetric, seen nearly face on with a flattened intrinsic shape - i.e., a face-on fast-rotator. We find an I-band mass-to-light ratio of $4.7_{-0.45}^{+0.32}$ in good agreement with previous spectral fitting results for this galaxy. Our best-fit orbit model shows a both a bulge and a disc component, with a non-negligible counter-rotating component also evident. We find that relatively little dark matter is needed to model this galaxy, with an inferred dark mass fraction of just $(7.1^{+8.1}_{-7.1})\%$ within one half-light radius.
[66]  oai:arXiv.org:1502.07477  [pdf] - 940835
The Concept of Few-Parameter Modelling of Eclipsing Binary and Exoplanet Transit Light Curves
Comments: 5 pages, 4 figures, Astronomical Society of the Pacific Conference Series - Litomysl 2014
Submitted: 2015-02-26
We present a new few-parameter phenomenological model of light curves of eclipsing binaries and stars with transiting planets that is able to fit the observed light curves with the accuracy better than 1\% of their amplitudes. The model can be used namely for appropriate descriptions of light curve shapes, classification, mid-eclipse time determination, and fine period analyses.
[67]  oai:arXiv.org:1407.2263  [pdf] - 1215512
The Next Generation Virgo cluster Survey. V. modelling the dynamics of M87 with the Made-to-Measure method
Comments: 20 pages, 16 figures, Accepted by ApJ
Submitted: 2014-07-08, last modified: 2014-07-13
We study the dynamics of the giant elliptical galaxy M87 from the central to the outermost regions with the made-to-measure (M2M) method. We use a new catalogue of 922 globular cluster line-of- sight velocities extending to a projected radius of 180 kpc (equivalent to 25 M87 effective radii), and SAURON integral field unit data within the central 2.4 kpc. 263 globular clusters, mainly located beyond 40 kpc, are newly observed by the Next Generation Virgo Survey (NGVS). For the M2M modelling, the gravitational potential is taken as a combination of a luminous matter potential with a constant stellar mass-to-light ratio and a dark matter potential modelled as a logarithmic potential. Our best dynamical model returns a stellar mass-to-light ratio in the I band of M/LI = 6.0(+ -0.3) M_sun/L_sun with a dark matter potential scale velocity of 591(+ -50) km/s and scale radius of 42(+ -10) kpc. We determine the total mass of M87 within 180 kpc to be (1.5 + - 0.2) 10^13 M_sun. The mass within 40 kpc is smaller than previous estimates determined using globular cluster kinematics that did not extend beyond 45 kpc. With our new globular cluster velocities at much larger radii, we see that globular clusters around 40 kpc show an anomalously large velocity dispersion which affected previous results. The mass we derived is in good agreement with that inferred from ROSAT X-ray observation out to 180 kpc. Within 30 kpc our mass is also consistent with that inferred from Chandra and XMM-Newton X-ray observations, while within 120 kpc it is about 20% smaller. The model velocity dispersion anisotropy beta parameter for the globular clusters in M87 is small, varying from -0.2 at the centre to 0.2 at 40 kpc, and gradually decreasing to zero at 120 kpc.
[68]  oai:arXiv.org:1310.7861  [pdf] - 791761
Hi-fi phenomenological description of eclipsing binary light variations as the basis for their period analysis
Comments: 4 pages, 3 figures, EAS - Proceedings of the conference: Setting a new standard in the analysis of binary stars, 16 to 19 September 2013, Leuven, Belgium
Submitted: 2013-10-29
In-depth analysis of eclipsing binary (EB) observational data collected for several decades can inform us about a lot of astrophysically interesting processes taking place in the systems. We have developed a wide-ranging method for the phenomenological modelling of eclipsing binary phase curves that enables us to combine even very disparate sources of phase information. This approach is appropriate for the processing of both standard photometric series of eclipses and data from photometric surveys of all kind. We conclude that mid-eclipse times, determined using the latest version of our 'hi-fi' phenomenological light curve models, as well as their accuracy, are nearly the same as the values obtained using much more complex standard physical EB models.
[69]  oai:arXiv.org:1309.6112  [pdf] - 1944154
Sub-arcsec Observations of NGC 7538 IRS 1: Continuum Distribution and Dynamics of Molecular Gas
Comments: 19 pages, 14 figures
Submitted: 2013-09-24
We report new results based on the analysis of the SMA and CARMA observations of NGC 7538\,IRS\,1 at 1.3 and 3.4 mm with sub-arcsec resolutions. With angular resolutions $\sim$ 0\farcs7, the SMA and CARMA observations show that the continuum emission at 1.3 and 3.4 mm from the hyper-compact \ion{H}{2} region IRS\,1 is dominated by a compact source with a tail-like extended structure to the southwest of IRS\,1. With a CARMA B-array image at 1.3 mm convolved to 0\farcs1, we resolve the hyper-compact \ion{H}{2} region into two components: an unresolved hyper-compact core, and a north-south extension with linear sizes of $<270$ AU and $\sim$2000 AU, respectively. The fine structure observed with CARMA is in good agreement with the previous VLA results at centimeter wavelengths, suggesting that the hyper-compact \ion{H}{2} region at the center of IRS\,1 is associated with an ionized bipolar outflow. We image the molecular lines OCS(19-18) and CH$_3$CN(12-11) as well as $^{13}$CO(2-1) surrounding IRS\,1, showing a velocity gradient along the southwest-northeast direction. The spectral line profiles in $^{13}$CO(2-1), CO(2-1), and HCN(1-0) observed toward IRS\,1 show broad redshifted absorption, providing evidence for gas infall with rates in the range of $3-10\times10^{-3}$ M$_\odot$ yr$^{-1}$ inferred from our observations.
[70]  oai:arXiv.org:1210.7570  [pdf] - 1152530
Results from Long-Term Optical Monitoring of the Soft X-Ray Transient SAX J1810.8-2609
Comments: 21 pages, 7 figures, ApJ accepted
Submitted: 2012-10-29, last modified: 2013-08-09
In this paper we report the long-term optical observation of the faint soft X-ray transient SAX J1810.8-2609 from OGLE and MOA. We have focused on the 2007 outburst, and also did the cross-correlate between its optical light curves and the quasi-simultaneous X-ray observations from swift. Both the optical and X-ray light curves of 2007 outburst show multi-peak features. Quasi-simultaneous optical/X-ray luminosity shows that both the X-ray reprocessing and viscously thermal emission can explain the observed optical flux. There is a slightly X-ray delay of 0.6+-0.3 days during the first peak, while the X-ray emission lags the optical emission by ~2 days during the rebrightening stage, which suggests that X-ray reprocessing emission contributes significantly to the optical flux in the first peak, but the viscously-heated-disk origin dominates the optical flux during rebrightening. It implies variation of the physical environment of the outer disk, even the source stayed in low-hard state during the whole outburst. The ~2 day X-ray lag indicates a small accretion disk of the system, and the optical counterpart was not detected by OGLE and MOA during quiescence, which constrained it to be fainter than MI = 7.5 mag. There is a suspected short-time optical flare detected at MJD = 52583.5 without X-ray counterpart detected, this single flux increase may imply a magnetic loop reconnection in the outer disk as proposed by Zurita et al (2003). The observations cover all stages of the outburst, however, due to the low sensitivity of RXTE/ASM, we cannot conclude whether it is an optical precursor at the initial rise of the outburst.
[71]  oai:arXiv.org:1212.5527  [pdf] - 606454
Basic relations for the period variation models of variable stars
Comments: 4 pages, 0 figures, written for Annual Proceendings of Science and Technology at V\v{S}B--TUO
Submitted: 2012-12-21
Models of period variations are basic tools for period analyzes of variable stars. We introduce phase function and instant period and formulate basic relations and equations among them. Some simple period models are also presented.
[72]  oai:arXiv.org:1212.5519  [pdf] - 606453
Visual minima timings of eclipsing binaries: To use or not to use?
Comments: 10 pages, 3 figures, written for Central European Astrophysical Bulletin
Submitted: 2012-12-21
Plenty of mid-eclipse timings of short-periodic eclipsing binaries derived from series of visual observations appear to be an acceptable source of information for orbital period analyses, namely if they were done in time-intervals not covered by other types of observations. However, our thorough period analysis of the nearly contact eclipsing binary BS Vulpeculae proves that visually determined timings done in 1979--2003 were undoubtedly biased to accommodate the existing linear ephemeris. The heavily subjective character of visual observations disqualifies them as a source of true phase information apt for fine eclipsing binary period analyses. Consequently we warn against the use of visual timings without a preceding careful verification.
[73]  oai:arXiv.org:1204.6104  [pdf] - 505879
Optical flares and flaring oscillations on the M-type eclipsing binary CU Cnc
Comments: It will be published in MNRAS
Submitted: 2012-04-26
We report here the discovery of an optical flare observed in R band from the red-dwarf eclipsing binary CU Cnc whose component stars are at the upper boundary of full convection (M1=0.43 and M2=0.4M0, M0 is the solar mass). The amplitude of the flare is the largest among those detected in R band (~0.52mag) and the duration time is about 73 minutes. As those observed on the Sun, quasi-periodic oscillations were seen during and after the flare. Three more R-band flares were found by follow up monitoring. In total, this binary was monitored photometrically by using R filter for 79.9 hours, which reveals a R-band flare rate about 0.05 flares per hour. These detections together with other strong chromospheric and coronal activities, i.e., very strong H_alpha and H_beta emission features and an EUV and X-ray source, indicate that it has very strong magnetic activity. Therefore, the apparent faintness (~1.4 magnitude in V) of CU Cnc compared with other single red dwarfs of the same mass can be plausibly explained by the high coverage of the dark spots.
[74]  oai:arXiv.org:1201.4205  [pdf] - 1093127
A Circumbinary Planet in Orbit Around the Short-Period White-Dwarf Eclipsing Binary RR Cae
Comments: It will be published in the MNRAS
Submitted: 2012-01-19
By using six new determined mid-eclipse times together with those collected from the literature, we found that the Observed-Calculated (O-C) curve of RR Cae shows a cyclic change with a period of 11.9 years and an amplitude of 14.3s, while it undergoes an upward parabolic variation (revealing a long-term period increase at a rate of dP/dt =+4.18(+-0.20)x10^(-12). The cyclic change was analyzed for the light-travel time effect that arises from the gravitational influence of a third companion. The mass of the third body was determined to be M_3*sin i' = 4.2(+-0.4) M_{Jup} suggesting that it is a circumbinary giant planet when its orbital inclination is larger than 17.6 degree. The orbital separation of the circumbinary planet from the central eclipsing binary is about 5.3(+-0.6)AU. The period increase is opposite to the changes caused by angular momentum loss via magnetic braking or/and gravitational radiation, nor can it be explained by the mass transfer between both components because of its detached configuration. These indicate that the observed upward parabolic change is only a part of a long-period (longer than 26.3 years) cyclic variation, which may reveal the presence of another giant circumbinary planet in a wide orbit.
[75]  oai:arXiv.org:1201.0944  [pdf] - 1092717
Observations of a high-mass protostar in NGC 7538S
Comments: ApJ Accepted 2012, 13 pages, 9 figures
Submitted: 2012-01-04
We present high angular resolution continuum observations of the high-mass protostar NGC 7538S with BIMA and CARMA at 3 and 1.4 mm, VLA observations at 1.3, 2, 3.5 and 6 cm, and archive IRAC observations from the Spitzer Space Observatory, which detect the star at 4.5, 5.8, and 8 $\mu$m. The star looks rather unremarkable in the mid-IR. The excellent positional agreement of the IRAC source with the VLA free-free emission, the OH, CH$_3$OH, H$_2$O masers, and the dust continuum confirms that this is the most luminous object in the NGC 7538S core. The continuum emission at millimeter wavelengths is dominated by dust emission from the dense cold cloud core surrounding the protostar. Including all array configurations, the emission is dominated by an elliptical source with a size of ~ 8" x 3". If we filter out the extended emission we find three compact mm-sources inside the elliptical core. The strongest one, $S_A$, coincides with the VLA/IRAC source and resolves into a double source at 1.4 mm, where we have sub-arcsecond resolution. The measured spectral index, $\alpha$, between 3 and 1.4 mm is ~ 2.3, and steeper at longer wavelengths, suggesting a low dust emissivity or that the dust is optically thick. We argue that the dust in these accretion disks is optically thick and estimate a mass of an accretion disk or infalling envelope surrounding S$_A$ to be ~ 60 solar masses.
[76]  oai:arXiv.org:1112.4269  [pdf] - 1092472
Circumbinary Planets Orbiting the Rapidly Pulsating Subdwarf B-type binary NY Vir
Comments: 12 pages, 3 figures, accepted by ApJ Letters
Submitted: 2011-12-19
We report here the tentative discovery of a Jovian planet in orbit around the rapidly pulsating subdwarf B-type (sdB-type) eclipsing binary NY Vir. By using new determined eclipse times together with those collected from the literature, we detect that the observed-calculated (O-C) curve of NY Vir shows a small-amplitude cyclic variation with a period of 7.9\,years and a semiamplitude of 6.1\,s, while it undergoes a downward parabolic change (revealing a period decrease at a rate of $\dot{P}=-9.2\times{10^{-12}}$). The periodic variation was analyzed for the light-travel time effect via the presence of a third body. The mass of the tertiary companion was determined to be $M_3\sin{i^{\prime}}=2.3(\pm0.3)$\,$M_{Jupiter}$ when a total mass of 0.60\,$M_{\odot}$ for NY Vir is adopted. This suggests that it is most probably a giant circumbinary planet orbiting NY Vir at a distance of about 3.3 astronomical units (AU). Since the rate of period decrease can not be explained by true angular momentum loss caused by gravitational radiation or/and magnetic braking, the observed downward parabolic change in the O-C diagram may be only a part of a long-period (longer than 15 years) cyclic variation, which may reveal the presence of another Jovian planet ($\sim2.5$$M_{Jupiter}$) in the system.
[77]  oai:arXiv.org:1108.0713  [pdf] - 1083171
Outflow, Infall and Protostars in the Star Forming Core W3-SE
Comments: 16 pages, 11 figures
Submitted: 2011-08-02
We report new results on outflow and infall in the star forming cores W3-SE SMA-1 and SMA-2 based on analysis of $\sim2.5\arcsec$ resolution observations of the molecular lines HCN(3-2), HCO$^+$(3-2), N$_2$H$^+$(3-2) and CH$_3$OH(5$_{2,3}-4_{1,3}$) with the Submillimeter Array. A high-velocity bipolar outflow originating from the proto-stellar core SMA-1 was observed in the HCN(3-2) line, with a projected outflow axis in position angle 48$\arcdeg$. The detection of the outflow is confirmed from other molecular lines. An inverse P-Cygni profile in the HCN(3-2) line toward SMA-1 suggests that at least one of the double cores accretes matters from the molecular core. A filamentary structure in the molecular gas surrounds SMA-1 and SMA-2. Based on the SMA observations, our analysis suggests that the double pre-stellar cores SMA-1 and SMA-2 result from fragmentation in the collapsing massive molecular core W3-SE, and it is likely that they are forming intermediate to high-mass stars which will be new members of a star cluster in the W3-SE region.
[78]  oai:arXiv.org:1105.3794  [pdf] - 366688
The Five-Hundred-Meter Aperture Spherical Radio Telescope (FAST) Project
Comments: 36 pages, 28 figures, accepted for publication in International Journal of Modern Physics D
Submitted: 2011-05-19
Five-hundred-meter Aperture Spherical radio Telescope (FAST) is a Chinese mega-science project to build the largest single dish radio telescope in the world. Its innovative engineering concept and design pave a new road to realize a huge single dish in the most effective way. FAST also represents Chinese contribution in the international efforts to build the square kilometer array (SKA). Being the most sensitive single dish radio telescope, FAST will enable astronomers to jump-start many science goals, for example, surveying the neutral hydrogen in the Milky Way and other galaxies, detecting faint pulsars, looking for the first shining stars, hearing the possible signals from other civilizations, etc. The idea of sitting a large spherical dish in a karst depression is rooted in Arecibo telescope. FAST is an Arecibo-type antenna with three outstanding aspects: the karst depression used as the site, which is large to host the 500-meter telescope and deep to allow a zenith angle of 40 degrees; the active main reflector correcting for spherical aberration on the ground to achieve a full polarization and a wide band without involving complex feed systems; and the light-weight feed cabin driven by cables and servomechanism plus a parallel robot as a secondary adjustable system to move with high precision. The feasibility studies for FAST have been carried out for 14 years, supported by Chinese and world astronomical communities. The project time is 5.5 years from the commencement of work in March of 2011 and the first light is expected to be in 2016. This review intends to introduce FAST project with emphasis on the recent progress since 2006. In this paper, the subsystems of FAST are described in modest details followed by discussions of the fundamental science goals and examples of early science projects.
[79]  oai:arXiv.org:1103.2006  [pdf] - 1052632
High Fill-Out, Extreme Mass Ratio Overcontact Binary Systems. X. The new discovered binary XY Leonis Minoris
Comments: 19 pages, 4 figures, 9 tables
Submitted: 2011-03-10
The new discovered short-period close binary star, XY LMi, was monitored photometrically since 2006. It is shown that the light curves are typical EW-type and show complete eclipses with an eclipse duration of about 80 minutes. By analyzing the complete B, V, R, and I light curves with the 2003 version of the W-D code, photometric solutions were determined. It is discovered that XY LMi is a high fill-out, extreme mass ratio overcontact binary system with a mass ratio of q=0.148 and a fill-out factor of f=74.1%, suggesting that it is on the late evolutionary stage of late-type tidal-locked binary stars. As observed in other overcontact binary stars, evidence for the presence of two dark spots on both components are given. Based on our 19 epoches of eclipse times, it is found that the orbital period of the overcontact binary is decreasing continuously at a rate of dP/dt=-1.67\times10^{-7}\,days/year, which may be caused by the mass transfer from the primary to the secondary or/and angular momentum loss via magnetic stellar wind. The decrease of the orbital period may result in the increase of the fill-out, and finally, it will evolve into a single rapid-rotation star when the fluid surface reaching the outer critical Roche Lobe.
[80]  oai:arXiv.org:1103.2005  [pdf] - 1052631
Detection of a planetary system orbiting the eclipsing polar HU Aqr
Comments: 13 pages, 4 figures, 2 tables
Submitted: 2011-03-10
Using the precise times of mid-egress of the eclipsing polar HU Aqr, we discovered that this polar is orbited by two or more giant planets. The two planets detected so far have masses of at least 5.9 and 4.5\,M_{Jup}. Their respective distances from the polar are 3.6 AU and 5.4 AU with periods of 6.54 and 11.96 years, respectively. The observed rate of period decrease derived from the downward parabolic change in O-C curve is a factor 15 larger than the value expected for gravitational radiation. This indicates that it may be only a part of a long-period cyclic variation, revealing the presence of one more planet. It is interesting to note that the two detected circumbinary planets follow the Titus-Bode law of solar planets with n=5 and 6. We estimate that another 10 years of observations will reveal the presence of the predicted third planet.
[81]  oai:arXiv.org:0912.2636  [pdf] - 902739
Calibrating the correlation between black hole mass and X-ray variability amplitude: X-ray only black hole mass estimates for active galactic nuclei and ultra-luminous X-ray sources
Comments: ApJ 710, 16-23, 2010, color figures online only
Submitted: 2009-12-14, last modified: 2010-01-18
A calibration is made for the correlation between the X-ray Variability Amplitude (XVA) and Black Hole (BH) mass. The correlation for 21 reverberation-mapped Active Galactic Nuclei (AGN) appears very tight, with an intrinsic dispersion of 0.20 dex. The intrinsic dispersion of 0.27 dex can be obtained if BH masses are estimated from the stellar velocity dispersions. We further test the uncertainties of mass estimates from XVAs for objects which have been observed multiple times with good enough data quality. The results show that the XVAs derived from multiple observations change by a factor of 3. This means that BH mass uncertainty from a single observation is slightly worse than either reverberation-mapping or stellar velocity dispersion measurements; however BH mass estimates with X-ray data only can be more accurate if the mean XVA value from more observations is used. Applying this relation, the BH mass of RE J1034+396 is found to be $4^{+3}_{-2} \times 10^6$ $M_{\odot}$. The high end of the mass range follows the relationship between the 2$f_0$ frequencies of high-frequency QPO and the BH masses derived from the Galactic X-ray binaries. We also calculate the high-frequency constant $C= 2.37 M_\odot$ Hz$^{-1}$ from 21 reverberation-mapped AGN. As suggested by Gierli\'nski et al., $M_{\rm BH}=C/C_{\rm M}$, where $C_{\rm M}$ is the high-frequency variability derived from XVA. Given the similar shape of power-law dominated X-ray spectra in ULXs and AGN, this can be applied to BH mass estimates of ULXs. We discuss the observed QPO frequencies and BH mass estimates in the Ultra-Luminous X-ray source M82 X-1 and NGC 5408 X-1 and favor ULXs as intermediate mass BH systems (abridged).
[82]  oai:arXiv.org:0912.2579  [pdf] - 902733
Dust and HCO+ Gas in the Star Forming Core W3-SE
Comments: 14 pages, 8 figures; accepted by ApJ
Submitted: 2009-12-14
We report new results from CARMA observations of both continuum and HCO+(1-0) emission at 3.4 mm from W3-SE, a molecular core of intermediate mass, together with the continuum observations at 1.1 and 0.85/0.45 mm with the SMA and JCMT. A continuum emission core elongated from SE to NW (~10"), has been observed at the and further resolved into a double source with the SMA at 1.1 mm, with a separation of ~4". Together with the measurements from the Spitzer and MSX at mid-IR, we determined the SED of W3-SE and fit it with a thermal dust emission model, suggesting the presence of two dust components with different temperatures. The emission at mm/submm wavelengths is dominated by a major cold (~41 K) with a mass of ~65 Msun. In addition, there is a weaker hot component (~400 K) which accounts for emission in the mid-IR, suggesting that a small fraction of dust has been heated by newly formed stars. We also imaged the molecular core in the HCO+(1-0) line using CARMA at an angular resolution ~6". With the CARMA observations, we have verified the presence of a blue-dominated double peak profile toward this core. The line profile cannot be explained by infall alone. The broad velocity wings of the line profile suggest that other kinematics such as outflows within the central 6" of the core likely dominate the resulting spectrum. The kinematics of the sub-structures of this core suggest that the molecular gas outside the main component appears to be dominated by the bipolar outflow originated from the dust core with a dynamical age of >30000 yr. Our analysis, based on the observations at wavelengths from mm/submm to mid-IR suggests that the molecular core W3-SE hosts a group of newly formed young stars and protostars.
[83]  oai:arXiv.org:0807.3992  [pdf] - 562899
Evidence for an intermediate line region in AGN's inner torus region and its evolution from narrow to broad line Seyfert I galaxies
Comments: 58 pages, 19 figures, 1 table. Matches the published version
Submitted: 2008-07-25, last modified: 2009-08-12
A two-components model for Broad Line Region (BLR) of Active Galactic Nuclei (AGN) has been suggested for many years but not widely accepted (e.g., Hu et al. 2008; Sulentic et al. 2000; Brotherton et al. 1996; Mason et al. 1996). This model indicates that the broad line can be described with superposition of two Gaussian components (Very Broad Gaussian Component (VBGC) and InterMediate Gaussian Component (IMGC)) which are from two physically distinct regions; i.e., Very Broad Line Region (VBLR) and InterMediate Line Region (IMLR). We select a SDSS sample to further confirm this model and give detailed analysis to the geometry, density and evolution of these two regions. Micro-lensing result of BLR in J1131-1231 and some unexplained phenomena in Reverberation Mapping (RM) experiment provide supportive evidence for this model. Our results indicate that the radius obtained from the emission line RM normally corresponds to the radius of the VBLR, and the existence of the IMGC may affect the measurement of the black hole masses in AGNs. The deviation of NLS1s from the M-sigma relation and the Type II AGN fraction as a function of luminosity can be explained in this model in a coherent way. The evolution of the two emission regions may be related to the evolutionary stages of the broad line regions of AGNs from NLS1s to BLS1s. Based on the results presented here, a unified picture of hierarchical evolution of black hole, dust torus and galaxy is proposed.
[84]  oai:arXiv.org:0907.1942  [pdf] - 901504
Diagnostics of the structure of AGN's broad line regions with reverberation mapping data: confirmation of the two-component broad line region model
Comments: 7 pages, 5 figures. Accepted for publication in the Special Issue of Science in China (G) "Astrophysics of Black holes and Related Compact Objects"
Submitted: 2009-07-11
We re-examine the ten Reverberation Mapping (RM) sources with public data based on the two-component model of the Broad Line Region (BLR). In fitting their broad H-beta lines, six of them only need one Gaussian component, one of them has a double-peak profile, one has an irregular profile, and only two of them need two components, i.e., a Very Broad Gaussian Component (VBGC) and an Inter-Mediate Gaussian Component (IMGC). The Gaussian components are assumed to come from two distinct regions in the two-component model; they are Very Broad Line Region (VBLR) and Inter-Mediate Line region (IMLR). The two sources with a two-component profile are Mrk 509 and NGC 4051. The time lags of the two components of both sources satisfy $t_{IMLR}/t_{VBLR}=V^2_{VBLR}/V^2_{IMLR}$, where $t_{IMLR}$ and $t_{VBLR}$ are the lags of the two components while $V_{IMLR}$ and $V_{VBLR}$ represent the mean gas velocities of the two regions, supporting the two-component model of the BLR of Active Galactic Nuclei (AGN). The fact that most of these ten sources only have the VBGC confirms the assumption that RM mainly measures the radius of the VBLR; consequently, the radius obtained from the R-L relationship mainly represent the radius of VBLR. Moreover, NGC 4051, with a lag of about 5 days in the one component model, is an outlier on the R-L relationship as shown in Kaspi et al. (2005); however this problem disappears in our two-component model with lags of about 2 and 6 days for the VBGC and IMGC, respectively.
[85]  oai:arXiv.org:0903.1357  [pdf] - 22166
Detection of a tertiary brown dwarf companion in the sdB-type eclipsing binary HS 0705+6700
Comments: 10 pages, 3 figures To be published in Astrophysical Journal
Submitted: 2009-03-07, last modified: 2009-03-16
HS 0705+6700 is a short-period (P=2.3 hours), close binary containing a hot sdB-type primary and a fully convective secondary. We have monitored this eclipsing binary for more than 2 years and as a result, 32 times of light minimum were obtained. Based on our new eclipse times together with these compiled from the literature, it is discovered that the O-C curve of HS 0705+6700 shows a cyclic variation with a period of 7.15 years and a semiamplitude of 92.4 s. The periodic change was analyzed for the light-travel time effect that may be due to the presence of a tertiary companion. The mass of the third body is determined to be M3 sin i = 0.0377 (+/-0.0043) Msun when a total mass of 0.617 Msun for HS 0705+6700 is adopted. For orbital inclinations i >= 32.8, the mass of the tertiary component would be below the stable hydrogen-burning limit of M3~0.072 Msun, and thus it would be a brown dwarf. The third body is orbiting the sdB-type binary at a distance shorter than 3.6 astronomical units (AU). HS 0705+6700 was formed through the evolution of a common envelope after the primary becomes a red giant. The detection of a sub-stellar companion in HS 0705+6700 system at this distance from the binary could give some constraints on stellar evolution in such systems and the interactions between red giants and their companions.
[86]  oai:arXiv.org:0901.2167  [pdf] - 20354
Evidence for an InterMediate Line Region in AGN's Inner Torus Region and Its Evolution from Narrow to Broad Line Seyfert I Galaxies
Comments: 4 pages, 6 figures, 10th-APRIM-2008-Kunming-Session-3
Submitted: 2009-01-14
We have decomposed the broad H-alpha, H-beta and H-gamma lines of 90 Active Galactic Nuclei (AGNs) into a superposition of a very broad and an intermediate Gaussian components (VBGC and IMGC) and discovered that the two Gaussian components evolve with FWHM of the whole emission lines. We suggest that the VBGC and the IMGC are produced in different emission regions, namely, Very Broad Line Region (VBLR) and Intermediate Line Region (IMLR). The details of the two components of H-alpha, H-beta and H-gamma lines indicate that the radius obtained from the emission line reverberation mapping normally corresponds to the radius of the VBLR, but the radius obtained from the infrared reverberation mapping corresponding to IMLR, i.e., the inner boundary of the dusty torus. The existence of the IMGC may affect the measurement of the black hole mass in AGNs. Therefore, the deviation of NLS1s from the M-sigma relation may be explained naturally in this way. The evolution of the two emission line regions may be related to the evolutionary stages of the broad line regions of AGNs from NLS1s to BLS1s. Other evidences for the existence of the IMLR are also presented.