The Event Horizon Telescope (EHT) recently released the first horizon-scale images of the black hole in M87. Combined with other astronomical data, these images constrain the mass and spin of the hole as well as the accretion rate and magnetic flux trapped on the hole. An important question for EHT is how well key parameters such as spin and trapped magnetic flux can be extracted from present and future EHT data alone. Here we explore parameter extraction using a neural network trained on high resolution synthetic images drawn from state-of-the-art simulations. We find that the neural network is able to recover spin and flux with high accuracy. We are particularly interested in interpreting the neural network output and understanding which features are used to identify, e.g., black hole spin. Using feature maps, we find that the network keys on low surface brightness features in particular.

In the standard model of solar flares, a large-scale reconnection current sheet is postulated as the central engine for powering the flare energy release and accelerating particles. However, where and how the energy release and particle acceleration occur remain unclear due to the lack of measurements for the magnetic properties of the current sheet. Here we report the measurement of spatially-resolved magnetic field and flare-accelerated relativistic electrons along a current-sheet feature in a solar flare. The measured magnetic field profile shows a local maximum where the reconnecting field lines of opposite polarities closely approach each other, known as the reconnection X point. The measurements also reveal a local minimum near the bottom of the current sheet above the flare loop-top, referred to as a "magnetic bottle". This spatial structure agrees with theoretical predictions and numerical modeling results. A strong reconnection electric field of ~4000 V/m is inferred near the X point. This location, however, shows a local depletion of microwave-emitting relativistic electrons. These electrons concentrate instead at or near the magnetic bottle structure, where more than 99% of them reside at each instant. Our observations suggest that the loop-top magnetic bottle is likely the primary site for accelerating and/or confining the relativistic electrons.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present complete results for planet-candidate hosts from the K2-HERMES survey, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra for more than 30,000 K2 stars. We present complete host-star parameters and planet-candidate radii for 224 K2 candidate planets from C1-C13. Our results cast severe doubt on 30 K2 candidates, as we derive unphysically large radii, larger than 2 $R_{Jup}$. This work highlights the importance of obtaining accurate, precise, and self-consistent stellar parameters for ongoing large planet search programs - something that will only become more important in the coming years, as TESS begins to deliver its own harvest of exoplanets.

We study gravitational lensing of the cosmic neutrino background. This signal is undetectable for the foreseeable future, but there is a rich trove of information available. At least some of the neutrinos from the early universe will be non-relativistic today, with a closer surface of last scattering (compared to the cosmic microwave background) and with larger angles of deflection. Lensing of massive neutrinos is strongly chromatic: both the amplitude of lensing and the cosmic time at which the potential is traversed depend on neutrino momentum, in principle giving access to our entire causal volume, not restricted to the light cone. As a concrete example, we focus on the case where the cosmic neutrino background would be strongly lensed when passing through halos of galaxy clusters and galaxies. We calculate the Einstein radius for cosmic neutrinos and investigate the impact of neutrino mass.

Binary stellar systems form a large fraction of the Galaxy's stars. They are useful as laboratories for studying the physical processes taking place within stars, and must be correctly taken into account when observations of stars are used to study the structure and evolution of the Galaxy. We present a sample of 12760 well-characterised double-lined spectroscopic binaries that are appropriate for statistical studies of the binary populations. They were detected as SB2s using a t-distributed stochastic neighbour embedding (t-SNE) classification and a cross-correlation analysis of GALAH spectra. This sample consists mostly of dwarfs, with a significant fraction of evolved stars and several dozen members of the giant branch. To compute parameters of the primary and secondary star ($T_{\rm eff[1,2]}$, $\log g_{[1,2]}$, [Fe/H], $V_{r[1,2]}$, $v_{\rm mic[1,2]}$, $v_{\rm broad[1,2]}$, $R_{[1,2]}$, and $E(B-V)$), we used a Bayesian approach that includes a parallax prior from Gaia DR2, spectra from GALAH, and apparent magnitudes from APASS, Gaia DR2, 2MASS, and WISE. The derived stellar properties and their distributions show trends that are expected for a population of close binaries (a $<$ 10 AU) with mass ratios $0.5 \leq q \leq 1$. The derived metallicity of these binary stars is statistically lower than that of single dwarf stars from the same magnitude-limited sample.

The possibility that in the mass range around $10^{-12}\ M_\odot$ most of dark matter constitutes of primordial black holes (PBHs) is a very interesting topic. To produce PBHs with this mass, the primordial scalar power spectrum needs to be enhanced to the order of 0.01 at the scale $k\sim 10^{12}\ \text{Mpc}^{-1}$. The enhanced power spectrum also produces large secondary gravitational waves at the mHz band. A phenomenological delta function power spectrum is usually used to discuss the production of PBHs and secondary gravitational waves. Based on G and k inflations, we propose a new mechanism to enhance the power spectrum at small scales by introducing a non-canonical kinetic term $[1-2G(\phi)]X$ with the function $G(\phi)$ having a peak. Away from the peak, $G(\phi)$ is negligible and we recover the usual slow-roll inflation which is constrained by the cosmic microwave background anisotrpy observations. Around the peak, the slow-roll inflation transiently turns to ultra slow-roll inflation. The enhancement of the power spectrum can be obtained with generic potentials, and there is no need to fine tune the parameters in $G(\phi)$. The energy spectrum $\Omega_{GW}(f)$ of secondary gravitational waves have the characteristic power law behaviour $\Omega_{GW}(f)\sim f^{n}$ and is testable by pulsar timing array and space based gravitational wave detectors.

We explore the fundamental relations governing the radial and vertical velocity dispersions of stars in the Milky Way, from combined studies of complementary surveys including GALAH, LAMOST, APOGEE, the NASA $Kepler$ and K2 missions, and $Gaia$ DR2. We find that different stellar samples, even though they target different tracer populations and employ a variety of age estimation techniques, follow the same set of fundamental relations. We provide the clearest evidence to date that, in addition to the well-known dependence on stellar age, the velocity dispersions of stars depend on orbital angular momentum $L_z$, metallicity and height above the plane $|z|$, and are well described by a multiplicatively separable functional form. The dispersions have a power-law dependence on age with exponents of 0.441$\pm 0.007$ and 0.251$\pm 0.006$ for $\sigma_z$ and $\sigma_R$ respectively, and the power law is valid even for the oldest stars. For the solar neighborhood stars, the apparent break in the power law for older stars, as seen in previous studies, is due to the anti-correlation of $L_z$ with age. The dispersions decrease with increasing $L_z$ until we reach the Sun's orbital angular momentum, after which $\sigma_z$ increases (implying flaring in the outer disc) while $\sigma_R$ flattens. The dispersions increase with decreasing metallicity, suggesting that the dispersions increase with birth radius. The dispersions also increase linearly with $|z|$. The same set of relations that work in the solar neighborhood also work for stars between $3<R/{\rm kpc}<20$. Finally, the high-[$\alpha$/Fe] stars follow the same relations as the low-[$\alpha$/Fe] stars.

The neutron star low-mass X-ray binary SAX J1712.6-3739 is known for its very long and hard thermonuclear bursts in past X-ray observations. Its thermonuclear bursts are so energetic that they can even last for tens of minutes in the hard X-ray observations of Swift/BAT (E > 15 keV). For exploring the natures of SAX J1712.6-3739 and its thermonuclear bursts, we analyzed its four bursts which were all captured by Swift/BAT. The peak flux and fluence for these bursts are derived by applying the archive observation data of Swift/BAT and Swift/XRT. Based on the peak fluxes observed by Swift, the distance of SAX J1712.6-3739 is estimated as 4.6 kpc, which is much closer than the estimate by the BeppoSAX observations. The corresponding absolute magnitude and average accretion rate agree its ultra-compact X-ray binary nature. By the measurements of effective duration, we determined that the 2010 burst corresponds to a normal X-ray burst, the 2011 burst is consistent with an intermediate-duration burst while the 2014 and the 2018 bursts are more energetic than common intermediate-duration bursts but less energetic than superbursts. Since the average mass accretion rate of SAX J1712.6-3739 was only 0.6%M_Edd, and current theory predict no carbon production in the bursters under very low accretion rate. Therefore, the 2014 and the 2018 bursts are very likely to be deep helium bursts. In the Swift observations, SAX J1712.6-3739 has shown at least three different levels of duration in thermonuclear bursts. They were possibly induced by variable accretion rates, because their ignition column depths roughly follow the prediction of ignition models for pure helium bursts at corresponding accretion rates.

Ultra-long-duration gamma-ray burst GRB 111209A was found to be associated with a very luminous supernovae (SNe) SN 2011kl. The physics of GRB 111209A/SN 2011kl has been extensively studied in the literures, but does not settle down yet. By investigating in detail the characteristics of the X-ray light curve of GRB 111209A, coupled with the temporal and spectral features observed in SN 2011kl, we argue that a short-living supramassive magnetar can be responsible for the initial shallow X-ray emission. Then the electromagnetic extraction of spin energy from a black hole results in the steeply declining X-ray flux when the magnetar collapses into a black hole (BH). A fraction of the envelope materials falls back and activates the accretion onto the newborn BH, which produces the X-ray rebrightening bump at late times. During this process, a centrifugally driven baryon-rich quasi-isotropic Blandford \& Payne outflow from the revived accretion disk deposits its kinetic energy on the SN ejecta, which powers luminous SN 2011kl. Finally, we place a limitation on the magnetar's physical parameters based on the observations.

Cryogenic phonon-scintillation detectors (CPSD) for rare event search experiments require reliable, efficient and robust photon detectors that can resolve individual photons in a scintillation event.

The Large Underground Xenon (LUX) dark matter search was a 250-kg active mass dual-phase time projection chamber that operated by detecting light and ionization signals from particles incident on a xenon target. In December 2015, LUX reported a minimum 90% upper C.L. of 6e-46 cm^2 on the spin-independent WIMP-nucleon elastic scattering cross section based on a 1.4e4 kg*day exposure in its first science run. Tension between experiments and the absence of a definitive positive detection suggest it would be prudent to search for WIMPs outside the standard spin-independent/spin-dependent paradigm. Recent theoretical work has identified a complete basis of 14 independent effective field theory (EFT) operators to describe WIMP-nucleon interactions. In addition to spin-independent and spin-dependent nuclear responses, these operators can produce novel responses such as angular-momentum-dependent and spin-orbit couplings. Here we report on a search for all 14 of these EFT couplings with data from LUX's first science run. Limits are placed on each coupling as a function of WIMP mass.

N-body simulations are essential tools in physical cosmology to understand the large-scale structure (LSS) formation of the Universe. Large-scale simulations with high resolution are important for exploring the substructure of universe and for determining fundamental physical parameters like neutrino mass. However, traditional particle-mesh (PM) based algorithms use considerable amounts of memory, which limits the scalability of simulations. Therefore, we designed a two-level PM algorithm CUBE towards optimal performance in memory consumption reduction. By using the fixed-point compression technique, CUBE reduces the memory consumption per N-body particle toward 6 bytes, an order of magnitude lower than the traditional PM-based algorithms. We scaled CUBE to 512 nodes (20,480 cores) on an Intel Cascade Lake based supercomputer with $\simeq$95\% weak-scaling efficiency. This scaling test was performed in "Cosmo-$\pi$" -- a cosmological LSS simulation using $\simeq$4.4 trillion particles, tracing the evolution of the universe over $\simeq$13.7 billion years. To our best knowledge, Cosmo-$\pi$ is the largest completed cosmological N-body simulation. We believe CUBE has a huge potential to scale on exascale supercomputers for larger simulations.

We present a new investigation of the dependence of H2O maser detection rates and properties on the mid-IR AGN luminosity, L_AGN, and the obscuring column density, N_H, based on mid-IR and hard X-ray photometry. Based on spectral energy distribution fitting that allows for decomposition of the black hole accretion and star-formation components in the mid-infrared, we show that the megamaser (disk maser) detection rate increases sharply for galaxies with 12 micron AGN luminosity L^{AGN}_{12 micron} greater than 10^42 erg/s, from ~<3%(~<2%) to ~12%(~5%). By using the ratio of the observed X-ray to mid-IR AGN luminosity as an indicator of N_H, we also find that the megamaser (disk maser) detection rates are boosted to 15%(7%) and 20%(9%) for galaxies with N_H >= 10^23 cm^{-2} and N_H >= 10^{24} cm^{-2}, respectively. Combining these column density cuts with a constraint for high L^{AGN}_{12 micron} (>=10^42 erg/s) predicts further increases in the megamaser (disk maser) detection rates to 19%(8%) and 27%(14%), revealing unprecedented potential boosts of the megamaser and disk maser detection rates by a factor of 7-15 relative to the current rates, depending on the chosen sample selection criteria. A noteworthy aspect of these new predictions is that the completeness rates are only compromised mildly, with the rates remaining at the level of ~95%(~50%) for sources with N_H >= 10^{23} cm^{-2} (N_H >= 10^24 cm^-2). Applying these selection methods to the current X-ray AGN surveys predicts the detection of >~15 new megamaser disks.

For filter-based magnetographs, the linear calibration method under the weak-field assumption is usually adopted; this leads to magnetic saturation effect in the regions with strong magnetic field. This article explores a new method to overcome the above disadvantage using a multilayer perceptron network, which we call MagMLP, based on a back-propagation algorithm with one input layer, five hidden layers, and one output layer. We use the data from the \textit{Spectropolarimeter} (SP) on board \textit{Hinode} to simulate single-wavelength observations for the model training, and take into account the influence of the Doppler velocity field and the filling factor. The training results show that the linear fitting coefficient (LFC) of the transverse field reaches above 0.91, and that of the longitudinal field is above 0.98. The generalization of the models is good because the corresponding LFCs are above 0.9 for the test subsets. Compared with the linear calibration method, the MagMLP is much more effective on dealing with the magnetic saturation effect. Analyzing an active region, the results of the linear calibration present an evident magnetic saturation effect in the umbra regions; the corresponding systematic error reaches values greater than 1000 G in most areas, or even exceeds 2000 G at some pixels. However, the results of MagMLP at these locations are very close to the inversion results, and the systematic errors are basically within 300 G. In addition, we find that there are many "bright spots" and "dark spots" on the inclination angle images from the inversion results of \textit{Hinode}/SP with values of 180 and 0 degrees, respectively, where the inversion is not reliable and does not produce a good result; the MagMLP handles these points well.

We present a novel analysis technique for liquid xenon time projection chambers that allows for a lower threshold by relying on events with a prompt scintillation signal consisting of single detected photons. The energy threshold of the LUX dark matter experiment is primarily determined by the smallest scintillation response detectable, which previously required a 2-fold coincidence signal in its photomultiplier arrays, enforced in data analysis. The technique presented here exploits the double photoelectron emission effect observed in some photomultiplier models at vacuum ultraviolet wavelengths. We demonstrate this analysis using an electron recoil calibration dataset and place new constraints on the spin-independent scattering cross section of weakly interacting massive particles (WIMPs) down to 2.5 GeV/c$^2$ WIMP mass using the 2013 LUX dataset. This new technique is promising to enhance light WIMP and astrophysical neutrino searches in next-generation liquid xenon experiments.

LUX-ZEPLIN (LZ) is a next generation dark matter direct detection experiment that will operate 4850 feet underground at the Sanford Underground Research Facility (SURF) in Lead, South Dakota, USA. Using a two-phase xenon detector with an active mass of 7 tonnes, LZ will search primarily for low-energy interactions with Weakly Interacting Massive Particles (WIMPs), which are hypothesized to make up the dark matter in our galactic halo. In this paper, the projected WIMP sensitivity of LZ is presented based on the latest background estimates and simulations of the detector. For a 1000 live day run using a 5.6 tonne fiducial mass, LZ is projected to exclude at 90% confidence level spin-independent WIMP-nucleon cross sections above $1.6 \times 10^{-48}$ cm$^{2}$ for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP. Additionally, a $5\sigma$ discovery potential is projected reaching cross sections below the existing and projected exclusion limits of similar experiments that are currently operating. For spin-dependent WIMP-neutron(-proton) scattering, a sensitivity of $2.7 \times 10^{-43}$ cm$^{2}$ ($8.1 \times 10^{-42}$ cm$^{2}$) for a 40 $\mathrm{GeV}/c^{2}$ mass WIMP is expected. With underground installation well underway, LZ is on track for commissioning at SURF in 2020.

We present isochrone ages and initial bulk metallicities ($\rm [Fe/H]_{bulk}$, by accounting for diffusion) of 163,722 stars from the GALAH Data Release 2, mainly composed of main sequence turn-off stars and subgiants ($\rm 7000 K>T_{eff}>4000 K$ and $\rm log g>3$ dex). The local age-metallicity relationship (AMR) is nearly flat but with significant scatter at all ages; the scatter is even higher when considering the observed surface abundances. After correcting for selection effects, the AMR appear to have intrinsic structures indicative of two star formation events, which we speculate are connected to the thin and thick disks in the solar neighborhood. We also present abundance ratio trends for 16 elements as a function of age, across different $\rm [Fe/H]_{bulk}$ bins. In general, we find the trends in terms of [X/Fe] vs age from our far larger sample to be compatible with studies based on small ($\sim$ 100 stars) samples of solar twins but we now extend it to both sub- and super-solar metallicities. The $\alpha$-elements show differing behaviour: the hydrostatic $\alpha$-elements O and Mg show a steady decline with time for all metallicities while the explosive $\alpha$-elements Si, Ca and Ti are nearly constant during the thin disk epoch (ages $\lessapprox $ 12 Gyr). The s-process elements Y and Ba show increasing [X/Fe] with time while the r-process element Eu have the opposite trend, thus favouring a primary production from sources with a short time-delay such as core-collapse supernovae over long-delay events such as neutron star mergers.

We describe the design and assembly of the LUX-ZEPLIN experiment, a direct detection search for cosmic WIMP dark matter particles. The centerpiece of the experiment is a large liquid xenon time projection chamber sensitive to low energy nuclear recoils. Rejection of backgrounds is enhanced by a Xe skin veto detector and by a liquid scintillator Outer Detector loaded with gadolinium for efficient neutron capture and tagging. LZ is located in the Davis Cavern at the 4850' level of the Sanford Underground Research Facility in Lead, South Dakota, USA. We describe the major subsystems of the experiment and its key design features and requirements.

The scattering of dark matter (DM) particles with sub-GeV masses off nuclei is difficult to detect using liquid xenon-based DM search instruments because the energy transfer during nuclear recoils is smaller than the typical detector threshold. However, the tree-level DM-nucleus scattering diagram can be accompanied by simultaneous emission of a Bremsstrahlung photon or a so-called "Migdal" electron. These provide an electron recoil component to the experimental signature at higher energies than the corresponding nuclear recoil. The presence of this signature allows liquid xenon detectors to use both the scintillation and the ionization signals in the analysis where the nuclear recoil signal would not be otherwise visible. We report constraints on spin-independent DM-nucleon scattering for DM particles with masses of 0.4-5 GeV/c$^2$ using 1.4$\times10^4$ kg$\cdot$day of search exposure from the 2013 data from the Large Underground Xenon (LUX) experiment for four different classes of mediators. This analysis extends the reach of liquid xenon-based DM search instruments to lower DM masses than has been achieved previously.

We apply a new deep learning technique to detect, classify, and deblend sources in multi-band astronomical images. We train and evaluate the performance of an artificial neural network built on the Mask R-CNN image processing framework, a general code for efficient object detection, classification, and instance segmentation. After evaluating the performance of our network against simulated ground truth images for star and galaxy classes, we find a precision of 92% at 80% recall for stars and a precision of 98% at 80% recall for galaxies in a typical field with $\sim30$ galaxies/arcmin$^2$. We investigate the deblending capability of our code, and find that clean deblends are handled robustly during object masking, even for significantly blended sources. This technique, or extensions using similar network architectures, may be applied to current and future deep imaging surveys such as LSST and WFIRST. Our code, Astro R-CNN, is publicly available at https://github.com/burke86/astro_rcnn.

The nature of fast radio bursts (FRBs), which occurs on millisecond time scales in the radio band, has not been well-understood. Among their unknown observational properties are their broadband spectra and persistent and transient multi-wavelength counterparts. Well-localized FRBs provide us the opportunity to address these issues in archival observations. We have performed searches for 15-150 keV hard X-ray bursts on time scales as short as millisecond in the direction of the repeating FRB 121102 (with a spacial resolution of a few arcminutes) in the archival Swift/BAT data during the period between October 2016 and September 2017. We have found no significant (5 $\sigma$) hard X-ray bursts in the direction of the repeating FRB. We have derived an upper limit of the hard X-ray (15--150 keV) flux of any X-ray bursts on 1 ms time scale of around $1.01 \times 10^{-7}$erg cm$^{-2}$s$^{-1}$, if assuming a photo-index of 2 for potential X-ray flares in X-ray band. A plausible scenario for the repeating FRB as being associated with \emph{magnetar giant flare} is still far below the upper limit.

The Orion complex is arguably the most studied star-forming region in the Galaxy. While stars are still being born in the Orion nebula, the oldest part was believed to be no more than 13 Myr old. In order to study the full hierarchy of star formation across the Orion complex, we perform a clustering analysis of the Ori OB1a region using new stellar surveys and derive robust ages for each identified stellar aggregate. We use Gaia DR2 parameters supplemented with radial velocities from the GALAH and APOGEE surveys to perform clustering of the Ori OB1a association. Five overdensities are resolved in a six-dimensional parameter space (positions, distance, proper motions, and radial velocity). Most correspond to previously known structures (ASCC 16, 25 Orionis, ASCC 20, ASCC 21). We use Gaia DR2, Pan-STARRS1 and 2MASS photometry to fit isochrones to the colour-magnitude diagrams of the identified clusters. The ages of the clusters can thus be measured with ~10% precision. While four of the clusters have ages between 11 and 13 Myr, the ASCC 20 cluster stands out at an age of 21 $\pm$ 3 Myr. This is significantly greater than the age of any previously known component of the Orion complex. To some degree, all clusters overlap in at least one of the six phase-space dimensions. We argue that the formation history of the Orion complex, and its relation to the Gould belt, must be reconsidered. A significant challenge in reconstructing the history of the Ori OB1a association is to understand the impact of the newly discovered 21 Myr old population on the younger parts of the complex, including their formation.

Gaia DR2 has revealed new small-scale and large-scale patterns in the phase-space distribution of stars in the Milky Way. In cylindrical Galactic coordinates $(R,\phi,z)$, ridge-like structures can be seen in the \vphiR{} plane and asymmetric arch-like structures in the \vphivR{} plane. We show that the ridges are also clearly present when the third dimension of the \vphiR{} plane is represented by $\langle z \rangle$, $\langle V_z \rangle$, $\langle V_R \rangle$, $\langle$[Fe/H]$\rangle$ and $\langle[\alpha/{\rm Fe}]\rangle$. The maps suggest that stars along the ridges lie preferentially close to the Galactic midplane ($|z|<0.2$ kpc), and have metallicity and $\alpha$ elemental abundance similar to that of the Sun. We show that phase mixing of disrupting spiral arms can generate both the ridges and the arches. It also generates discrete groupings in orbital energy $-$ the ridges and arches are simply surfaces of constant energy. We identify 8 distinct ridges in the \gaia{} data: six of them have constant energy while two have constant angular momentum. Given that the signature is strongest for stars close to the plane, the presence of ridges in $\langle z \rangle$ and $\langle V_z \rangle$ suggests a coupling between planar and vertical directions. We demonstrate, using N-body simulations that such coupling can be generated both in isolated discs and in discs perturbed by an orbiting satellite like the Sagittarius dwarf galaxy.

Many debris disks seen in scattered light have shapes that imply their dust grains trace highly eccentric, apsidally aligned orbits. Apsidal alignment is surprising, especially for dust. Even when born from an apse-aligned ring of parent bodies, dust grains have their periastra dispersed in all directions by stellar radiation pressure. The periastra cannot be re-oriented by planets within the short dust lifetimes at the bottom of the collisional cascade. We propose that what re-aligns dust orbits is drag exerted by second-generation gas. Gas is largely immune to radiation pressure, and when released by photodesorption or collisions within an eccentric ring of parent bodies should occupy a similarly eccentric, apse-aligned ring. Dust grains launched onto misaligned orbits cross the eccentric gas ring supersonically and can become dragged into alignment within collisional lifetimes. The resultant dust configurations, viewed nearly but not exactly edge-on, with periastra pointing away from the observer, appear moth-like, with kinked wings and even doubled pairs of wings, explaining otherwise mysterious features in HD 61005 ("The Moth") and HD 32297, including their central bulbs when we account for strong forward scattering from irregularly shaped particles. Around these systems we predict gas at Kuiper-belt-like distances to move on highly elliptical streamlines that owe their elongation, ultimately, to highly eccentric planets. Unresolved issues and an alternative explanation for apsidal alignment are outlined.

We report results from an extensive set of measurements of the \b{eta}-decay response in liquid xenon.These measurements are derived from high-statistics calibration data from injected sources of both $^{3}$H and $^{14}$C in the LUX detector. The mean light-to-charge ratio is reported for 13 electric field values ranging from 43 to 491 V/cm, and for energies ranging from 1.5 to 145 keV.

The latest Gaia data release enables us to accurately identify stars that are more luminous than would be expected on the basis of their spectral type and distance. During an investigation of the 329 best Solar twin candidates uncovered among the spectra acquired by the GALAH survey, we identified 64 such over-luminous stars. In order to investigate their exact composition, we developed a data-driven methodology that can generate a synthetic photometric signature and spectrum of a single star. By combining multiple such synthetic stars into an unresolved binary or triple system and comparing the results to the actual photometric and spectroscopic observations, we uncovered 6 definitive triple stellar system candidates and an additional 14 potential candidates whose combined spectrum mimics the Solar spectrum. Considering the volume correction factor for a magnitude limited survey, the fraction of probable unresolved triple stars with long orbital periods is ~2 %. Possible orbital configurations of the candidates were investigated using the selection and observational limits. To validate the discovered multiplicity fraction, the same procedure was used to evaluate the multiplicity fraction of other stellar types.

Asteroseismology is a promising tool to study Galactic structure and evolution because it can probe the ages of stars. Earlier attempts comparing seismic data from the {\it Kepler} satellite with predictions from Galaxy models found that the models predicted more low-mass stars compared to the observed distribution of masses. It was unclear if the mismatch was due to inaccuracies in the Galactic models, or the unknown aspects of the selection function of the stars. Using new data from the K2 mission, which has a well-defined selection function, we find that an old metal-poor thick disc, as used in previous Galactic models, is incompatible with the asteroseismic information. We show that spectroscopic measurements of [Fe/H] and [$\alpha$/Fe] elemental abundances from the GALAH survey indicate a mean metallicity of $\log (Z/Z_{\odot})=-0.16$ for the thick disc. Here $Z$ is the effective solar-scaled metallicity, which is a function of [Fe/H] and [$\alpha$/Fe]. With the revised disc metallicities, for the first time, the theoretically predicted distribution of seismic masses show excellent agreement with the observed distribution of masses. This provides an indirect verification of the asteroseismic mass scaling relation is good to within five percent. Using an importance-sampling framework that takes the selection function into account, we fit a population synthesis model of the Galaxy to the observed seismic and spectroscopic data. Assuming the asteroseismic scaling relations are correct, we estimate the mean age of the thick disc to be about 10 Gyr, in agreement with the traditional idea of an old $\alpha$-enhanced thick disc.

Based on the Chinese historical sunspots drawings, a data set consisting of the scanned images and all their digitized parameters from 1925 to 2015 have been constructed. In this paper, we briefly describe the developmental history of sunspots drawings in China. This paper describes the preliminary processing processes that strat from the initial data (inputing to the scanning equipment) to the parameters extraction, and finally summarizes the general features of this dataset. It is the first systematic project in Chinese solar-physics community that the historical observation of sunspots drawings were digitized. Our data set fills in an almost ninety years historical gap, which span 60 degrees from east to west and 50 degrees from north to south and have no continuous and detailed digital sunspot observation information. As a complementary to other sunspots observation in the world, our dataset provided abundant information to the long term solar cycles solar activity research.

We present a statistical analysis of the first 300 stars observed by the Gemini Planet Imager Exoplanet Survey (GPIES). This subsample includes six detected planets and three brown dwarfs; from these detections and our contrast curves we infer the underlying distributions of substellar companions with respect to their mass, semi-major axis, and host stellar mass. We uncover a strong correlation between planet occurrence rate and host star mass, with stars M $>$ 1.5 $M_\odot$ more likely to host planets with masses between 2-13 M$_{\rm Jup}$ and semi-major axes of 3-100 au at 99.92% confidence. We fit a double power-law model in planet mass (m) and semi-major axis (a) for planet populations around high-mass stars (M $>$ 1.5M$_\odot$) of the form $\frac{d^2 N}{dm da} \propto m^\alpha a^\beta$, finding $\alpha$ = -2.4 $\pm$ 0.8 and $\beta$ = -2.0 $\pm$ 0.5, and an integrated occurrence rate of $9^{+5}_{-4}$% between 5-13 M$_{\rm Jup}$ and 10-100 au. A significantly lower occurrence rate is obtained for brown dwarfs around all stars, with 0.8$^{+0.8}_{-0.5}$% of stars hosting a brown dwarf companion between 13-80 M$_{\rm Jup}$ and 10-100 au. Brown dwarfs also appear to be distributed differently in mass and semi-major axis compared to giant planets; whereas giant planets follow a bottom-heavy mass distribution and favor smaller semi-major axes, brown dwarfs exhibit just the opposite behaviors. Comparing to studies of short-period giant planets from the RV method, our results are consistent with a peak in occurrence of giant planets between ~1-10 au. We discuss how these trends, including the preference of giant planets for high-mass host stars, point to formation of giant planets by core/pebble accretion, and formation of brown dwarfs by gravitational instability.

We present the chemodynamic structure of the solar neighbourhood using 62 814 stars within a 500 pc sphere of the Sun observed by GALAH and with astrometric parameters from Gaia DR2. We measure the velocity dispersion for all three components (vertical, radial, and tangential) and find that it varies smoothly with [Fe/H] and [$\alpha$/Fe] for each component. The vertical component is especially clean, with $\sigma_{v_z}$ increasing from a low of $8$ km s$^{-1}$ at solar-[$\alpha$/Fe] and [Fe/H] to a high of more than 50 km s$^{-1}$ for more metal-poor and [$\alpha$/Fe] enhanced populations. We find no evidence of a large decrease in the velocity dispersion of the higher-[$\alpha$/Fe] populations as claimed in analysis prior to Gaia DR2, although the trend of increasing velocity dispersion with [$\alpha$/Fe] for the same metallicity does significantly flatten at high-[$\alpha$/Fe]. The eccentricity distribution for local stars varies most strongly as a function of [$\alpha$/Fe], where stars with [$\alpha$/Fe]$<0.1$ dex having generally circular orbits ($e<0.15$), while the median eccentricity increases rapidly for more [$\alpha$/Fe]-enhanced stellar populations up to $e\sim0.35$. These [$\alpha$/Fe]-enhanced populations have guiding radii consistent with origins in the inner Galaxy. Of the stars with metallicities much higher than the local ISM ([Fe/H]>0.1 dex), we find that more than 70\% have $e<0.2$ and are likely observed in the solar neighbourhood through churning/migration rather than blurring effects, as the epicyclic motion for these stars is not large enough to reach the radii at which they were likely born based on their metallicity.

We use the second data releases of the ESA Gaia astrometric survey and the high-resolution GALAH spectroscopic survey to analyse the structure of our Galaxy's disc components. With GALAH, we separate the alpha-rich and alpha-poor discs (with respect to Fe), which are superposed in both position and velocity space, and examine their distributions in action space. We study the distribution of stars in the zV_z phase plane, for both V_phi and V_R, and recover the remarkable "phase spiral" discovered by Gaia. We identify the anticipated quadrupole signature in zV_z of a tilted velocity ellipsoid for stars above and below the Galactic plane. By connecting our work with earlier studies, we show that the phase spiral is likely to extend well beyond the narrow solar neighbourhood cylinder in which it was found. The phase spiral is a signature of corrugated waves that propagate through the disc, and the associated non-equilibrium phase mixing. The radially asymmetric distribution of stars involved in the phase spiral reveals that the corrugation, which is mostly confined to the alpha-poor disc, grows in z-amplitude with increasing radius. We present new simulations of tidal disturbance of the Galactic disc by the Sagittarius (Sgr) dwarf. The effect on the zV_z phase plane lasts >2 Gyr but a subsequent disc crossing wipes out the coherent structure. We find that the phase spiral was excited < 0.5 Gyr ago by an object like Sgr with total mass 3 x 10^10 Msun (stripped down from 5 x 10^10 Msun when it first entered the halo) passing through the plane.

We have investigated the outburst properties of low-mass X-ray binary transients (LMXBTs) based on a comprehensive study of the outbursts observed in the past few decades. The outburst rates were estimated based on the X-ray monitoring data from Swift/BAT, RXTE/ASM and MAXI, and previous reports in the literature. We found that almost all LMXBTs with the orbital period below $\sim$12 hr showed only one outburst in these observations. There are systematic difference in the outburst rate between long-period ($P_{\rm orb} \gtrsim$ 12 hr) and short-period ($P_{\rm orb} \lesssim$ 12 hr) systems. We infer that mass transfer rate is responsible for the systematic difference, since the disk instability model (DIM) suggested that the mass transfer rate is a key factor affecting the quiescence time. The difference in outburst rate between long-period and short-period LMXBTs is probably due to the different mass transfer mechanism at different evolutionary stages of the donors. Based on the evolutionary tracks of single stars, we derived the critical orbital period for X-ray binaries that harbor a subgiant donor in various metallicity. The critical orbital period ($P_{\rm orb,crit}=$12.4 hr) is consistent with the above orbital period boundary obtained from the statistics of outburst rates. Furthermore, we found a negative correlation between the outburst rate and the orbital period in the samples for which the luminosity class of the donor star is III/IV. The best-fitting power-law index for the black hole subsamples is roughly consistent with the theoretical prediction for those systems with a donor star evolved off the main sequence.

The overlap between the spectroscopic Galactic Archaeology with HERMES (GALAH) survey & $Gaia$ provides a high-dimensional chemodynamical space of unprecedented size. We present a first analysis of a subset of this overlap, of 7066 dwarf, turn-off, & sub-giant stars. [...] We investigate correlations between chemical compositions, ages, & kinematics for this sample. Stellar parameters & elemental abundances are derived from the GALAH spectra with the spectral synthesis code SME. [...] We report Li, C, O, Na, Mg, Al, Si, K, Ca, Sc, Ti, V, Cr, Mn, Co, Ni, Cu, Zn, Y, as well as Ba & we note that we employ non-LTE calculations for Li, O, Al, & Fe. We show that the use of astrometric & photometric data improves the accuracy of the derived spectroscopic parameters, especially $\log g$. [...] we recover the result that stars of the high-$\alpha$ sequence are typically older than stars in the low-$\alpha$ sequence, the latter spanning $-0.7<$[Fe/H]$<+0.5$. While these two sequences become indistinguishable in [$\alpha$/Fe] vs. [Fe/H] at the metal-rich regime, we find that age can be used to separate stars from the extended high-$\alpha$ & the low-$\alpha$ sequence even in this regime. [...] we find that the old stars ($>8$ Gyr have lower angular momenta $L_z$ than the Sun, which implies that they are on eccentric orbits & originate from the inner disk. Contrary to some previous smaller scale studies we find a continuous evolution in the high-$\alpha$-sequence up to super-solar [Fe/H] rather than a gap, which has been interpreted as a separate "high-$\alpha$ metal-rich" population. Stars in our sample that are younger than 10 Gyr, are mainly found on the low $\alpha$-sequence & show a gradient in $L_z$ from low [Fe/H] ($L_z>L_{z,\odot}$) towards higher [Fe/H] ($L_z<L_{z,\odot}$), which implies that the stars at the ends of this sequence are likely not originating from the close solar vicinity.

If the Galaxy is axisymmetric and in dynamical equilibrium, we expect negligible fluctuations in the residual line-of-sight velocity field. Recent results using the \apg{} survey find significant fluctuations in velocity for stars in the midplane ($|z|<$0.25 kpc) out to 5 kpc, suggesting that the dynamical influence of non-axisymmetric features i.e., the Milky Way's bar, spiral arms and merger events extends out to the Solar neighborhood. Their measured power spectrum has a characteristic amplitude of 11 \kms{} on a scale of 2.5 kpc. The existence of such large-scale streaming motions has important implications for determining the Sun's motion about the Galactic Centre. Using Red Clump stars from \glh{} and \apg{}, we map the line-of-sight velocities around the Sun (d$<$5 kpc), and $|z|<$1.25 kpc from the midplane. By subtracting a smooth axisymmetric model for the velocity field, we study the residual fluctuations and compare our findings with mock survey generated by \glx{}. We find negligible large-scale fluctuations away from the plane. In the mid-plane, we reproduce the earlier \apg{} power spectrum but with 20\% smaller amplitude (9.3 \kms{}) after taking into account a few systematics (e.g., volume completeness). Using a flexible axisymmetric model the power-amplitude is further reduced to 6.3 \kms{}. Additionally, our simulations show that, in the plane, distances are underestimated for high-mass Red Clump stars which can lead to spurious power-amplitude of about 5.2 \kms{}. Taking this into account, we estimate the amplitude of real fluctuations to be $<$4.6 \kms{}, about a factor of three less than the \apg{} result.

The authors of this report met on 28-30 March 2018 at the New Jersey Institute of Technology, Newark, New Jersey, for a 3-day workshop that brought together a group of data providers, expert modelers, and computer and data scientists, in the solar discipline. Their objective was to identify challenges in the path towards building an effective framework to achieve transformative advances in the understanding and forecasting of the Sun-Earth system from the upper convection zone of the Sun to the Earth's magnetosphere. The workshop aimed to develop a research roadmap that targets the scientific challenge of coupling observations and modeling with emerging data-science research to extract knowledge from the large volumes of data (observed and simulated) while stimulating computer science with new research applications. The desire among the attendees was to promote future trans-disciplinary collaborations and identify areas of convergence across disciplines. The workshop combined a set of plenary sessions featuring invited introductory talks and workshop progress reports, interleaved with a set of breakout sessions focused on specific topics of interest. Each breakout group generated short documents, listing the challenges identified during their discussions in addition to possible ways of attacking them collectively. These documents were combined into this report-wherein a list of prioritized activities have been collated, shared and endorsed.

Magnetic reconnection is the most likely mechanism responsible for the high temperature events that are observed in strongly magnetized locations around the temperature minimum in the low solar chromosphere. This work improves upon our previous work ["Magnetic Reconnection in Strongly Magnetized Regions of the Low Solar Chromosphere", The Astrophysical Journal 852, 95 (2018) ] by using a more realistic radiative cooling model computed from the OPACITY project and the CHIANTI database. We find that the rate of ionization of the neutral component of the plasma is still faster than recombination within the current sheet region. For low $\beta$ plasmas, the ionized and neutral fluid flows are well-coupled throughout the reconnection region resembling the single-fluid Sweet-Parker model dynamics. Decoupling of the ion and neutral inflows appears in the higher $\beta$ case with $\beta_0=1.46$, which leads to a reconnection rate about three times faster than the rate predicted by the Sweet-Parker model. In all cases, the plasma temperature increases with time inside the current sheet, and the maximum value is above $2\times10^4$~K when the reconnection magnetic field strength is greater than $500$~G. While the more realistic radiative cooling model does not result in qualitative changes of the characteristics of magnetic reconnection, it is necessary for studying the variations of the plasma temperature and ionization fraction inside current sheets in strongly magnetized regions of the low solar atmosphere. It is also important for studying energy conversion during the magnetic reconnection process when the hydrogen-dominated plasma approaches full ionization.

Identifying and characterizing young stars in the Solar neighbourhood is essential to find and describe planets in the early stages of their evolution. This work seeks to identify nearby young stars showing a Lithium 6707.78$\,\AA$ absorption line in the GALAH survey. A robust, data-driven approach is used to search for corresponding templates in the pool of 434,215 measured dwarf spectra in the survey. It enables a model-free search for best-matching spectral templates for all stars, including M dwarfs with strong molecular absorption bands. 3147 stars have been found to have measurable Lithium: 1408 G and 892 K0-K5 dwarfs (EW(Li)$>$0.1$\,$\AA), 335 K5-K9 ($>$0.07$\,$\AA) and 512 M0$\sim$M4 dwarfs ($>$0.05$\,$\AA). Stars with such Lithium features are used to investigate the possibility of searching for young stars above the main sequence based merely on their parallaxes and broad-band photometry. Selection of young stars above the main sequence is highly effective for M dwarfs, moderately effective for K dwarfs and ineffective for G dwarfs. Using a combination of the Lithium information and the complete 6D kinematics from Gaia and GALAH, 305 new candidate moving group members have been found, 123 of which belong to the Scorpius-Centaurus association, 36 to the Pleiades and 25 to the Hyades clusters.

Various dark matter models predict annual and diurnal modulations of dark matter interaction rates in Earth-based experiments as a result of the Earth's motion in the halo. Observation of such features can provide generic evidence for detection of dark matter interactions. This paper reports a search for both annual and diurnal rate modulations in the LUX dark matter experiment using over 20 calendar months of data acquired between 2013 and 2016. This search focuses on electron recoil events at low energies, where leptophilic dark matter interactions are expected to occur and where the DAMA experiment has observed a strong rate modulation for over two decades. By using the innermost volume of the LUX detector and developing robust cuts and corrections, we obtained a stable event rate of 2.3$\pm$0.2~cpd/keV$_{\text{ee}}$/tonne, which is among the lowest in all dark matter experiments. No statistically significant annual modulation was observed in energy windows up to 26~keV$_{\text{ee}}$. Between 2 and 6~keV$_{\text{ee}}$, this analysis demonstrates the most sensitive annual modulation search up to date, with 9.2$\sigma$ tension with the DAMA/LIBRA result. We also report no observation of diurnal modulations above 0.2~cpd/keV$_{\text{ee}}$/tonne amplitude between 2 and 6~keV$_{\text{ee}}$.

GALAH is a large-scale magnitude-limited southern stellar spectroscopic survey. Its second data release (GALAH DR2) provides values of stellar parameters and abundances of 23 elements for 342,682 stars (Buder et al.). Here we add a description of the public release of radial velocities with a typical accuracy of 0.1 km/s for 336,215 of these stars, achievable due to the large wavelength coverage, high resolving power and good signal to noise ratio of the observed spectra, but also because convective motions in stellar atmosphere and gravitational redshift from the star to the observer are taken into account. In the process we derive medians of observed spectra which are nearly noiseless, as they are obtained from between 100 and 1116 observed spectra belonging to the same bin with a width of 50 K in temperature, 0.2 dex in gravity, and 0.1 dex in metallicity. Publicly released 1181 median spectra have a resolving power of 28,000 and trace the well-populated stellar types with metallicities between -0.6 and +0.3. Note that radial velocities from GALAH are an excellent match to the accuracy of velocity components along the sky plane derived by Gaia for the same stars. The level of accuracy achieved here is adequate for studies of dynamics within stellar clusters, associations and streams in the Galaxy. So it may be relevant for studies of the distribution of dark matter.

Open cluster members are coeval and share the same initial bulk chemical compositions. Consequently, differences in surface abundances between members of a cluster that are at different evolutionary stages can be used to study the effects of mixing and internal chemical processing. We carry out an abundance analysis of seven elements (Li, O, Na, Mg, Al, Si, Fe) in 66 stars belonging to the open cluster M67, based on high resolution GALAH spectra, 1D MARCS model atmospheres, and, for the first time for a large sample of M67 stars, non-local thermodynamic equilibrium (non-LTE) radiative transfer. From the non-LTE analysis, we find a typical star-to-star scatter in the abundance ratios of around 0.05 dex; this scatter is slightly but systematically larger when LTE is assumed instead. We find trends in the abundance ratios with effective temperature, indicating systematic differences in the surface abundances between turn-off and giant stars; these trends are more pronounced when LTE is assumed. However, in the non-LTE analysis, most of the element trends have been flattened. Two species are exceptions to this behaviour, namely Al and Si, which both clearly display remaining trends in the non-LTE analysis. We comment on the possible origin of these trends, by comparing them with recent stellar models that include atomic diffusion.

Pebble accretion refers to the assembly of rocky planet cores from particles whose velocity dispersions are damped by drag from circumstellar disc gas. Accretion cross-sections can approach maximal Hill-sphere scales for particles whose Stokes numbers approach unity. While fast, pebble accretion is also lossy. Gas drag brings pebbles to protocores but also sweeps them past; those particles with the largest accretion cross-sections also have the fastest radial drift speeds and are the most easily drained out of discs. We present a global model of planet formation by pebble accretion that keeps track of the disc's mass budget. Cores, each initialized with a lunar mass, grow from discs whose finite stores of mm-cm sized pebbles drift inward across all radii in viscously accreting gas. For every 1 $M_\oplus$ netted by a core, at least 10 $M_\oplus$ and possibly much more are lost to radial drift. Core growth rates are typically exponentially sensitive to particle Stokes number, turbulent Mach number, and solid surface density. This exponential sensitivity, when combined with disc migration, tends to generate binary outcomes from 0.1-30 AU: either sub-Earth cores remain sub-Earth, or explode into Jupiters, with the latter migrating inward to varying degrees. When Jupiter-breeding cores assemble from mm-cm sized pebbles, they do so in discs where such particles drain out in $\sim$10$^5$ yr or less; such fast-draining discs do not fit mm-wave observations.

Two-episode emission components separated by quiescent gaps in the prompt emission of gamma-ray bursts (GRBs) have been observed in the {\em Swift} era, but there is a lack of spectral information due to the narrow energy band of the {\em Swift}/Burst Alert Telescope. In this paper, a systematic analysis of the spectral and temporal properties of the prompt emission of 101 {\em Fermi}/Gamma-ray Burst Monitor detected long GRBs show the existence of two-episode emission components in the light curves, with quiescent times of up to hundreds of seconds. We focus on investigating the differences of those two emission episodes. We find that the light curves of the two emission components exhibit different behavior, e.g., a soft emission component that either precedes or follows the main prompt emission or that the intensity of the two emission episodes are comparable with each other. No statistically significant correlation in the duration of the two emission episodes can be claimed. We define a new parameter $\varepsilon$ as the ratio of the peak flux of the first and second emission episodes and find that a higher $\varepsilon$ corresponds to a larger fluence. The preferred spectral model in our analysis is a cutoff power-law model for most GRBs. The distribution of $E_p$ for episodes I and II range from tens of keV to 1000 keV with a lognormal fit and there are no significant differences between them. Moreover, we do not find significant relationships between $\varepsilon$ and $E_p$ for the two emission episodes. Those results suggest that these two-episode emission components likely share the same physical origin.

Sparse open clusters can be found at high galactic latitudes where loosely populated clusters are more easily detected against the lower stellar background. Because most star formation takes place in the thin disk, the observed population of clusters far from the Galactic plane is hard to explain. We combined spectral parameters from the GALAH survey with the Gaia DR2 catalogue to study the dynamics and chemistry of five old sparse high-latitude clusters in more detail. We find that four of them (NGC 1252, NGC 6994, NGC 7772, NGC 7826) - originally classified in 1888 - are not clusters but are instead chance projections on the sky. Member stars quoted in the literature for these four clusters are unrelated in our multi-dimensional physical parameter space; the published cluster properties are therefore irrelevant. We confirm the existence of morphologically similar NGC 1901 for which we provide a probabilistic membership analysis. An overdensity in three spatial dimensions proves to be enough to reliably detect sparse clusters, but the whole 6-dimensional space must be used to identify members with high confidence, as demonstrated in the case of NGC 1901.

We present a study using the second data release of the GALAH survey of stellar parameters and elemental abundances of 15 pairs of stars identified by Oh et al 2017. They identified these pairs as potentially co-moving pairs using proper motions and parallaxes from Gaia DR1. We find that 11 very wide (>1.7 pc) pairs of stars do in fact have similar Galactic orbits, while a further four claimed co-moving pairs are not truly co-orbiting. Eight of the 11 co-orbiting pairs have reliable stellar parameters and abundances, and we find that three of those are quite similar in their abundance patterns, while five have significant [Fe/H] differences. For the latter, this indicates that they could be co-orbiting because of the general dynamical coldness of the thin disc, or perhaps resonances induced by the Galaxy, rather than a shared formation site. Stars such as these, wide binaries, debris of past star formation episodes, and coincidental co-orbiters, are crucial for exploring the limits of chemical tagging in the Milky Way.

We report the serendipitous observations of 571 luminous supergiants in the Magellanic Clouds by the spectroscopic GALAH and TESS-HERMES surveys: 434 stars in the Large Magellanic Cloud and 137 in the Small Magellanic Cloud. We also find one star that appears associated with structured star formation in the Magellanic Bridge. Both of these surveys are aimed at the local volume of the Galaxy but have simple, magnitude-limited selection functions that mean they include some observations of luminous extra-Galactic stars. The surveys determine stellar parameter and abundances using The Cannon, a data-driven generative modelling approach. In this work, we explore the results from The Cannon when it is fed the spectra of these intrinsically luminous supergiants in the Magellanic Clouds, which are well outside the normal bounds of The Cannon's training set. We find that, although the parameters are astrophysically incorrect, the $v\sin i$ and the abundances of lithium, barium, and magnesium are excellent discriminants of these stars. It shows that in the future, with an expanded training set, it should be possible to determine accurate values for these types of stars.

The Galactic Archaeology with HERMES (GALAH) survey is a large-scale stellar spectroscopic survey of the Milky Way and designed to deliver chemical information complementary to a large number of stars covered by the $Gaia$ mission. We present the GALAH second public data release (GALAH DR2) containing 342,682 stars. For these stars, the GALAH collaboration provides stellar parameters and abundances for up to 23 elements to the community. Here we present the target selection, observation, data reduction and detailed explanation of how the spectra were analysed to estimate stellar parameters and element abundances. For the stellar analysis, we have used a multi-step approach. We use the physics-driven spectrum synthesis of Spectroscopy Made Easy (SME) to derive stellar labels ($T_\mathrm{eff}$, $\log g$, $\mathrm{[Fe/H]}$, $\mathrm{[X/Fe]}$, $v_\mathrm{mic}$, $v \sin i$, $A_{K_S}$) for a representative training set of stars. This information is then propagated to the whole survey with the data-driven method of $The~Cannon$. Special care has been exercised in the spectral synthesis to only consider spectral lines that have reliable atomic input data and are little affected by blending lines. Departures from local thermodynamic equilibrium (LTE) are considered for several key elements, including Li, O, Na, Mg, Al, Si, and Fe, using 1D MARCS stellar atmosphere models. Validation tests including repeat observations, Gaia benchmark stars, open and globular clusters, and K2 asteroseismic targets lend confidence in our methods and results. Combining the GALAH DR2 catalogue with the kinematic information from $Gaia$ will enable a wide range of Galactic Archaeology studies, with unprecedented detail, dimensionality, and scope.

The primary goal of Galactic archaeology is to learn about the origin of the Milky Way from the detailed chemistry and kinematics of millions of stars. Wide-field multi-fibre spectrographs are increasingly used to obtain spectral information for huge samples of stars. Some surveys (e.g. GALAH) are attempting to measure up to 30 separate elements per star. Stellar abundance spectroscopy is a subtle art that requires a very high degree of spectral uniformity across each of the fibres. However wide-field spectrographs are notoriously non-uniform due to the fast output optics necessary to image many fibre outputs onto the detector. We show that precise spectroscopy is possible with such instruments across all fibres by employing a photonic comb -- a device that produces uniformly spaced spots of light on the CCD to precisely map complex aberrations. Aberrations are parametrized by a set of orthogonal moments with $\sim100$ independent parameters. We then reproduce the observed image by convolving high resolution spectral templates with measured aberrations as opposed to extracting the spectra from the observed image. Such a forward modeling approach also trivializes some spectroscopic reduction problems like fibre cross-talk, and reliably extracts spectra with a resolution $\sim2.3$ times above the nominal resolution of the instrument. Our rigorous treatment of optical aberrations also encourages a less conservative spectrograph design in the future.

Using GALAH survey data of nearby stars, we look at how structure in the planar (u,v) velocity distribution depends on metallicity and on viewing direction within the Galaxy. In nearby stars, with distance d < 1 kpc, the Hercules stream is most strongly seen in higher metallicity stars [Fe/H] > 0.2. The Hercules stream peak v value depends on viewed galactic longitude, which we interpret as due to the gap between the stellar stream and more circular orbits being associated with a specific angular momentum value of about 1640 km/s kpc. The association of the gap with a particular angular momentum value supports a bar resonant model for the Hercules stream. Moving groups previously identified in Hipparcos observations are easiest to see in stars nearer than 250 pc, and their visibility and peak velocities in the velocity distributions depends on both viewing direction (galactic longitude and hemisphere) and metallicity. We infer that there is fine structure in local velocity distributions that varies over distances of a few hundred pc in the Galaxy.

We present a Bayesian implementation of isochrone fitting in deriving stellar ages and masses, incorporating absolute K magnitude ($\rm M_K$) derived from 2MASS photometry and Gaia DR1 parallax and differentiation between initial bulk metallicity and present day surface metallicity, with allowance for incorporating further constraints (e.g., asteroseismology) when available. As a test, we re-computed stellar ages and masses of $\sim 4000$ stars in the solar neighbourhood from six well-studied literature samples using both Hipparcos and TGAS parallaxes. Our ages are found to be compatible with literature values but with reduced uncertainties in general. The inclusion of parallax-based $\rm M_K$ serves as an additional constraint on the derived quantities, especially when systematic errors in stellar parameters are underestimated. We reconstructed the age-metallicity relationship in the solar neighbourhood by re-analysing the Geneva-Copenhagen Survey with the inclusion of TGAS-parallaxes and initial bulk metallicity sampling. We found a flat trend for disk stars with ages $<$11\,Gyr but with smaller scatter at all ages compared to literature.

We study the first electromagnetic pulse after the gravitational wave chirp signal, focusing on the profile and duration. It is found that the light curve, especially the steep decay (SD) phase, can be very different by adopting different viewing angle $\theta_{\rm view}$ on the jet shell. For an on-axis jet with a power-law radiation spectrum, the observed flux in the SD is proportional to $t_{\rm{obs}}^{-2-\beta}$ with $\beta$ being the spectral index and $t_{\rm{obs}}$ being the observer time. Here, $t_{\rm{obs}}=0$ is set at the observed time of the jet ejected from the central engine. The SD may become steep by increasing $\theta_{\rm view}$. We also study the bolometric luminosity $L$ from a jet shell with a non-power-law radiation spectrum. For an on-axis jet, $L{\propto}t_{\rm{obs}}^{-3}$ is found in the SD. However, the SD is steeper than $L{\propto}t_{\rm{obs}}^{-3}$ for the radiation from an off-axis jet. The higher value of $\theta_{\rm view}$ is, the steeper of SD would be. Then, we suggest that the SD phase can be used to discriminate an off-axis jet from an on-axis jet. The reason for above behaviors is discussed. In addition, we find that the duration of first electromagnetic pulse is close to its peak time, especially for $\theta_{\rm{view}}\sim20^\circ$. This result is consistent with that found in GW~170817/GRB~170817A. Thus, the jet corresponding to the prompt emission of GRB~170817A should be immediately ejected after the merger. Our results also reveal that the duration of the first electromagnetic pulse can provide the information of the time to search gravitational waves.

Differences in the elemental abundances of planet hosting stars in binary systems can give important clues and constraints about planet formation and evolution. In this study we performed a high-precision, differential elemental abundance analysis of a wide binary system, HD 80606/80607, based on high-resolution, high signal-to-noise ratio Keck/HIRES spectra. HD 80606 is known to host a four Jupiter mass giant planet while no planet has yet been detected around HD 80607. We determined stellar parameters as well as abundances for 23 elements for these two stars with extremely high precision. Our main results are: (i) we confirmed that the two components share very similar chemical compositions, but HD 80606 is marginally more metal-rich than HD 80607 with an average difference of +0.013 $\pm$ 0.002 dex ($\sigma$ = 0.009 dex) and (ii) there is no obvious trend between abundance differences and condensation temperature. Assuming this binary formed from material with the same chemical composition, it is difficult to understand how giant planet formation could produce the present-day photospheric abundances of the elements we measure. We can not exclude the possibility that HD 80606 might have accreted about 2.5 to 5 $M_{\rm Earth}$ material onto its surface, possibly from a planet destabilised by the known highly-eccentric giant.

Using data from the GALAH pilot survey, we determine properties of the Galactic thin and thick disks near the solar neighbourhood. The data cover a small range of Galactocentric radius ($7.9 \leq R_\mathrm{GC} \leq 9.5$ kpc), but extend up to 4 kpc in height from the Galactic plane, and several kpc in the direction of Galactic anti-rotation (at longitude $260 ^\circ \leq \ell \leq 280^\circ$). This allows us to reliably measure the vertical density and abundance profiles of the chemically and kinematically defined `thick' and `thin' disks of the Galaxy. The thin disk (low-$\alpha$ population) exhibits a steep negative vertical metallicity gradient, at d[M/H]/d$z=-0.18 \pm 0.01$ dex kpc$^{-1}$, which is broadly consistent with previous studies. In contrast, its vertical $\alpha$-abundance profile is almost flat, with a gradient of d[$\alpha$/M]/d$z$ = $0.008 \pm 0.002$ dex kpc$^{-1}$. The steep vertical metallicity gradient of the low-$\alpha$ population is in agreement with models where radial migration has a major role in the evolution of the thin disk. The thick disk (high-$\alpha$ population) has a weaker vertical metallicity gradient d[M/H]/d$z = -0.058 \pm 0.003$ dex kpc$^{-1}$. The $\alpha$-abundance of the thick disk is nearly constant with height, d[$\alpha$/M]/d$z$ = $0.007 \pm 0.002$ dex kpc$^{-1}$. The negative gradient in metallicity and the small gradient in [$\alpha$/M] indicate that the high-$\alpha$ population experienced a settling phase, but also formed prior to the onset of major SNIa enrichment. We explore the implications of the distinct $\alpha$-enrichments and narrow [$\alpha$/M] range of the sub-populations in the context of thick disk formation.

The Kelvin-Hemholtz (K-H) instability in the corona EUV jet is studied via 2.5 MHD numerical simulations. The jet results from magnetic reconnection due to the interation of the new emerging magnetic field and the pre-existing magnetic field in the corona. Our results show that the Alfven Mach number along the jet is about 5-14 just before the instability occurs, and it is even higher than 14 at some local areas. During the K-H instability process, several vortex-like plasma blobs of high temperature and high density appear along the jet, and magnetic fields have also been rolled up and the magnetic configuration including anti-parallel magnetic fields forms, which leads to magnetic reconnection at many X-points and current sheet fragments inside the vortex-like blob. After magnetic islands appear inside the main current sheet, the total kinetic energy of the reconnection outflows decreases, and cannot support the formation of the vortex-like blob along the jet any longer, then the K-H instability eventually disappears. We also present the results about how the guide field and the flux emerging speed affect the K-H instability. We find that the strong guide field inhibits the shock formation in the reconnection upward outflow regions but helps secondary magnetic islands appear earlier in the main current sheet, and then apparently suppresses the K-H instability. As the speed of the emerging magnetic field decreases, the K-H instability appears later, the highest temperature inside the vortex blob gets lower and the vortex structure gets smaller.

We investigate the vertical metallicity gradients of five mono-age stellar populations between 0 and 11 Gyr for a sample of 18 435 dwarf stars selected from the cross-matched Tycho-Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiment (RAVE) Data Release 5. We find a correlation between the vertical metallicity gradients and age, with no vertical metallicity gradient in the youngest population and an increasingly steeper negative vertical metallicity gradient for the older stellar populations. The metallicity at disc plane remains almost constant between 2 and 8 Gyr, and it becomes significantly lower for the $8 < \tau \leqslant 11$ Gyr population. The current analysis also reveals that the intrinsic dispersion in metallicity increases steadily with age. We discuss that our results are consistent with a scenario that (thin) disc stars formed from a flaring (thin) star-forming disc.

Magnetic reconnection in strongly magnetized regions around the temperature minimum region of the low solar atmosphere is studied by employing MHD-based simulations of a partially ionized plasma within a reactive 2.5D multi-fluid model. It is shown that in the absence of magnetic nulls in a low $\beta$ plasma the ionized and neutral fluid flows are well-coupled throughout the reconnection region. However, non-equilibrium ionization-recombination dynamics play a critical role in determining the structure of the reconnection region, lead to much lower temperature increases and a faster magnetic reconnection rate as compared to simulations that assume plasma to be in ionization-recombination equilibrium. The rate of ionization of the neutral component of the plasma is always faster than recombination within the current sheet region even when the initial plasma $\beta$ is as high as $\beta_0=1.46$. When the reconnecting magnetic field is in excess of a kilogauss and the plasma $\beta$ is lower than 0.0145, the initially weakly ionized plasmas can become fully ionized within the reconnection region and the current sheet can be strongly heated to above $2.5\times10^4$~K, even as most of the collisionally dissipated magnetic energy is radiated away. The Hall effect increases the reconnection rate slightly, but in the absence of magnetic nulls it does not result in significant asymmetries or change the characteristics of the reconnection current sheet down to meter scales.

Accurate and precise radius estimates of transiting exoplanets are critical for understanding their compositions and formation mechanisms. To know the planet, we must know the host star in as much detail as possible. We present first results from the K2-HERMES project, which uses the HERMES multi-object spectrograph on the Anglo-Australian Telescope to obtain R$\sim$28,000 spectra of up to 360 stars in one exposure. This ongoing project aims to derive self-consistent spectroscopic parameters for about half of K2 target stars. We present complete stellar parameters and isochrone-derived masses and radii for 46 stars hosting 57 K2 candidate planets in Campaigns 1-3. Our revised host-star radii cast severe doubt on three candidate planets: EPIC\,201407812.01, EPIC\,203070421.01, and EPIC\,202843107.01, all of which now have inferred radii well in excess of the largest known inflated Jovian planets.

We present our study on the spatially resolved H_alpha and M_star relation for 536 star-forming and 424 quiescent galaxies taken from the MaNGA survey. We show that the star formation rate surface density (Sigma_SFR), derived based on the H_alpha emissions, is strongly correlated with the M_star surface density (Sigma_star) on kpc scales for star- forming galaxies and can be directly connected to the global star-forming sequence. This suggests that the global main sequence may be a consequence of a more fundamental relation on small scales. On the other hand, our result suggests that about 20% of quiescent galaxies in our sample still have star formation activities in the outer region with lower SSFR than typical star-forming galaxies. Meanwhile, we also find a tight correlation between Sigma_H_alpha and Sigma_star for LI(N)ER regions, named the resolved "LI(N)ER" sequence, in quiescent galaxies, which is consistent with the scenario that LI(N)ER emissions are primarily powered by the hot, evolved stars as suggested in the literature.

Luminosity of X-ray spectral state transitions in black hole and neutron star X-ray binaries can put constraint on the critical mass accretion rate between accretion regimes. Previous studies indicate that the hard-to-soft spectral state transitions in some ultracompact neutron star LMXBs have the lowest luminosity. With X-ray monitoring observations in the past decade, we were able to identify state transitions towards the lowest luminosity limit in 4U 0614+091, 2S 0918-549 and 4U 1246-588. By analysing corresponding X-ray pointed observations with the Swift/XRT and the RXTE/PCA, we found no hysteresis of state transitions in these sources, and determined the critical mass accretion rate in the range of 0.002 - 0.04 $\dot{\rm M}_{\rm Edd}$ and 0.003 - 0.05 $\dot{\rm M}_{\rm Edd}$ for the hard-to-soft and the soft-to-hard transition, respectively, by assuming a neutron star mass of 1.4 solar masses. This range is comparable to the lowest transition luminosity measured in black hole X-ray binaries, indicating the critical mass accretion rate is not affected by the nature of the surface of the compact stars. Our result does not support the Advection-Dominated Accretion Flow (ADAF) model which predicts that the critical mass accretion rate in neutron star systems is an order of magnitude lower if same viscosity parameters are taken. The low transition luminosity and insignificant hysteresis in these ultracompact X-ray binaries provide further evidence that the transition luminosity is likely related to the mass in the disc.

We study the role of cold gas in quenching star formation in the green valley by analysing ALMA $^{12}$CO (1-0) observations of three galaxies with resolved optical spectroscopy from the MaNGA survey. We present resolution-matched maps of the star formation rate and molecular gas mass. These data are used to calculate the star formation efficiency (SFE) and gas fraction ($f_{\rm~gas}$) for these galaxies separately in the central `bulge' regions and outer disks. We find that, for the two galaxies whose global specific star formation rate (sSFR) deviates most from the star formation main sequence, the gas fraction in the bulges is significantly lower than that in their disks, supporting an `inside-out' model of galaxy quenching. For the two galaxies where SFE can be reliably determined in the central regions, the bulges and disks share similar SFEs. This suggests that a decline in $f_{\rm~gas}$ is the main driver of lowered sSFR in bulges compared to disks in green valley galaxies. Within the disks, there exist common correlations between the sSFR and SFE and between sSFR and $f_{\rm~gas}$ on kpc scales -- the local SFE or $f_{\rm~gas}$ in the disks declines with local sSFR. Our results support a picture in which the sSFR in bulges is primarily controlled by $f_{\rm~gas}$, whereas both SFE and $f_{\rm~gas}$ play a role in lowering the sSFR in disks. A larger sample is required to confirm if the trend established in this work is representative of green valley as a whole.

The LAMOST survey has acquired low-resolution spectra (R=1,800) for 5 million stars across the Milky Way, far more than any current stellar survey at a corresponding or higher spectral resolution. It is often assumed that only very few elemental abundances can be measured from such low-resolution spectra, limiting their utility for Galactic archaeology studies. However, Ting et al. (2017) used ab initio models to argue that low-resolution spectra should enable precision measurements of many elemental abundances, at least in theory. Here we verify this claim in practice by measuring the relative abundances of 14 elements from LAMOST spectra with a precision of $\lesssim$ 0.1 dex for objects with ${\rm S/N}_{\rm LAMOST}$ > 30 (per pixel). We employ a spectral modeling method in which a data-driven model is combined with priors that the model gradient spectra should resemble ab initio spectral models. This approach assures that the data-driven abundance determinations draw on physically sensible features in the spectrum in their predictions and do not just exploit astrophysical correlations among abundances. Our analysis is constrained to the number of elemental abundances measured in the APOGEE survey, which is the source of the training labels. Obtaining high quality/resolution spectra for a subset of LAMOST stars to measure more elemental abundances as training labels and then applying this method to the full LAMOST catalog will provide a sample with more than 20 elemental abundances that is an order of magnitude larger than current high-resolution surveys, substantially increasing the sample size for Galactic archaeology.

The Transiting Exoplanet Survey Satellite (TESS) will provide high precision time-series photometry for millions of stars with at least a half-hour cadence. Of particular interest are the circular regions of 12-degree radius centered around the ecliptic poles that will be observed continuously for a full year. Spectroscopic stellar parameters are desirable to characterize and select suitable targets for TESS, whether they are focused on exploring exoplanets, stellar astrophysics, or Galactic archaeology. Here, we present spectroscopic stellar parameters ($T_{\rm eff}$, $\log g$, [Fe/H], $v \sin i$, $v_{\rm micro}$) for about 16,000 dwarf and subgiant stars in TESS' southern continuous viewing zone. For almost all the stars, we also present Bayesian estimates of stellar properties including distance, extinction, mass, radius, and age using theoretical isochrones. Stellar surface gravity and radius are made available for an additional set of roughly 8,500 red giants. All our target stars are in the range $10<V<13.1$. Among them, we identify and list 227 stars belonging to the Large Magellanic Cloud. The data were taken using the the High Efficiency and Resolution Multi-Element Spectrograph (HERMES, R $\sim 28,000$) at the Anglo-Australian Telescope as part of the TESS-HERMES survey. Comparing our results with the TESS Input Catalog (TIC) shows that the TIC is generally efficient in separating dwarfs and giants, but it has flagged more than hundred cool dwarfs ($T_{\rm eff}< 4800$ K) as giants, which ought to be high-priority targets for the exoplanet search. The catalog can be accessed via http://www.physics.usyd.edu.au/tess-hermes/ , or at MAST via https://archive.stsci.edu/prepds/tess-hermes/ .

The technique of chemical tagging uses the elemental abundances of stellar atmospheres to `reconstruct' chemically homogeneous star clusters that have long since dispersed. The GALAH spectroscopic survey --which aims to observe one million stars using the Anglo-Australian Telescope -- allows us to measure up to 30 elements or dimensions in the stellar chemical abundance space, many of which are not independent. How to find clustering reliably in a noisy high-dimensional space is a difficult problem that remains largely unsolved. Here we explore t-distributed stochastic neighbour embedding (t-SNE) -- which identifies an optimal mapping of a high-dimensional space into fewer dimensions -- whilst conserving the original clustering information. Typically, the projection is made to a 2D space to aid recognition of clusters by eye. We show that this method is a reliable tool for chemical tagging because it can: (i) resolve clustering in chemical space alone, (ii) recover known open and globular clusters with high efficiency and low contamination, and (iii) relate field stars to known clusters. t-SNE also provides a useful visualization of a high-dimensional space. We demonstrate the method on a dataset of 13 abundances measured in the spectra of 187,000 stars by the GALAH survey. We recover 7 of the 9 observed clusters (6 globular and 3 open clusters) in chemical space with minimal contamination from field stars and low numbers of outliers. With chemical tagging, we also identify two Pleiades supercluster members (which we confirm kinematically), one as far as 6$^\circ$ -- one tidal radius away from the cluster centre.

We derive a lower bound on the field excursion for the tachyon inflation, which is determined by the amplitude of the scalar perturbation and the number of $e$-folds before the end of inflation. Using the relation between the observables like $n_s$ and $r$ with the slow-roll parameters, we reconstruct three classes of tachyon potentials. The model parameters are determined from the observations before the potentials are reconstructed, and the observations prefer the concave potential. We also discuss the constraints from the reheating phase preceding the radiation domination for the three classes of models by assuming the equation of state parameter $w_{re}$ during reheating is a constant. Depending on the model parameters and the value of $w_{re}$, the constraints on $N_{re}$ and $T_{re}$ are different. As $n_s$ increases, the allowed reheating epoch becomes longer for $w_{re}=-1/3$, 0 and $1/6$ while the allowed reheating epoch becomes shorter for $w_{re}=2/3$.

We present the results of measurements demonstrating the efficiency of the EDELWEISS-III array of cryogenic germanium detectors for direct dark matter searches. The experimental setup and the FID (Fully Inter-Digitized) detector array is described, as well as the efficiency of the double measurement of heat and ionization signals in background rejection. For the whole set of 24 FID detectors used for coincidence studies, the baseline resolutions for the fiducial ionization energy are mainly below 0.7 keV$_{ee}$ (FHWM) whereas the baseline resolutions for heat energies are mainly below 1.5 keV$_{ee}$ (FWHM). The response to nuclear recoils as well as the very good discrimination capability of the FID design has been measured with an AmBe source. The surface $\beta$- and $\alpha$-decay rejection power of $R_{\rm surf} < 4 \times 10^{-5}$ per $\alpha$ at 90% C.L. has been determined with a $^{210}$Pb source, the rejection of bulk $\gamma$-ray events has been demonstrated using $\gamma$-calibrations with $^{133}$Ba sources leading to a value of $R_{\gamma{\rm -mis-fid}} < 2.5 \times 10^{-6}$ at 90% C.L.. The current levels of natural radioactivity measured in the detector array are shown as the rate of single $\gamma$ background. The fiducial volume fraction of the FID detectors has been measured to a weighted average value of $(74.6 \pm 0.4)\%$ using the cosmogenic activation of the $^{65}$Zn and $^{68,71}$Ge isotopes. The stability and uniformity of the detector response is also discussed. The achieved resolutions, thresholds and background levels of the upgraded EDELWEISS-III detectors in their setup are thus well suited to the direct search of WIMP dark matter over a large mass range.

We perform two-dimensional resistive magnetohydrodynamic simulations of coronal jets driven by flux emergence along the lower boundary. The reconnection layers are susceptible the formation of blobs that are ejected in the jet. Our simulation with low plasma $\beta$ (Case I) shows that magnetic islands form easily and propagate upwards in the jet. These islands are multithermal and thus are predicted to show up in hot channels (335 \AA\, and 211 \AA) and the cool channel (304 \AA) in observations by the Atmospheric Imaging Assembly (AIA) on the \emph{Solar Dynamics Observatory}. The islands have maximum temperatures of 8 MK, lifetimes of 120 s, diameters of 6 Mm, and velocities of 200 km s$^{-1}$. These parameters are similar to the properties of blobs observed in EUV jets by AIA\@. The Kelvin-Helmholtz instability develops in our simulation with moderately high plasma $\beta$ (Case II), and leads to the formation of bright vortex-like blobs above the multiple high magneto-sonic Mach number regions that appear along the jet. These vortex-like blobs can also be identified in the AIA channels. However, they eventually move downward and disappear after the high magneto-sonic Mach number regions disappear. In the lower plasma $\beta$ case, the lifetime for the jet is shorter, the jet and magnetic islands are formed with higher velocities and temperatures, the current sheet fragments are more chaotic, and more magnetic islands are generated. Our results show that the plasmoid instability and Kelvin-Helmholtz instability along the jet are both possible causes of the formation of blobs observed at extreme ultraviolet (EUV) wavelengths.

In this Technical Design Report (TDR) we describe the LZ detector to be built at the Sanford Underground Research Facility (SURF). The LZ dark matter experiment is designed to achieve sensitivity to a WIMP-nucleon spin-independent cross section of three times ten to the negative forty-eighth square centimeters.

We report the discovery of a mysterious giant $H_{\alpha}$ blob that is $\sim 8$ kpc away from the main MaNGA target 1-24145, one component of a dry galaxy merger, identified in the first-year SDSS-IV MaNGA data. The size of the $H_{\alpha}$ blob is $\sim$ 3-4 kpc in radius, and the $H_{\alpha}$ distribution is centrally concentrated. However, there is no optical continuum counterpart in deep broadband images reaching $\sim$26.9 mag arcsec$^{-2}$ in surface brightness. We estimate that the masses of ionized and cold gases are $3.3 \times 10^{5}$ $\rm M_{\odot}$ and $< 1.3 \times 10^{9}$ $\rm M_{\odot}$, respectively. The emission-line ratios indicate that the $H_{\alpha}$ blob is photoionized by a combination of massive young stars and AGN. Furthermore, the ionization line ratio decreases from MaNGA 1-24145 to the $H_{\alpha}$ blob, suggesting that the primary ionizing source may come from MaNGA 1-24145, likely a low-activity AGN. Possible explanations of this $H_{\alpha}$ blob include AGN outflow, the gas remnant being tidally or ram-pressure stripped from MaNGA 1-24145, or an extremely low surface brightness (LSB) galaxy. However, the stripping scenario is less favoured according to galaxy merger simulations and the morphology of the $H_{\alpha}$ blob. With the current data, we can not distinguish whether this $H_{\alpha}$ blob is ejected gas due to a past AGN outburst, or a special category of `ultra-diffuse galaxy' (UDG) interacting with MaNGA 1-24145 that further induces the gas inflow to fuel the AGN in MaNGA 1-24145.

Galah is an ongoing high-resolution spectroscopic survey with the goal of disentangling the formation history of the Milky Way, using the fossil remnants of disrupted star formation sites which are now dispersed around the Galaxy. It is targeting a randomly selected, magnitude limited ($V \leq 14$) sample of stars, with the goal of observing one million objects. To date, 300,000 spectra have been obtained. Not all of them are correctly processed by parameter estimation pipelines and we need to know about them. We present a semi-automated classification scheme which identifies different types of peculiar spectral morphologies, in an effort to discover and flag potentially problematic spectra and thus help to preserve the integrity of the survey's results. To this end we employ a recently developed dimensionality reduction technique t-SNE (t-distributed Stochastic Neighbour Embedding), which enables us to represent the complex spectral morphology in a two-dimensional projection map while still preserving the properties of the local neighbourhoods of spectra. We find that the majority (178,483) of the 209,533 Galah spectra considered in this study represents normal single stars, whereas 31,050 peculiar and problematic spectra with very diverse spectral features pertaining to 28,579 stars are distributed into 10 classification categories: Hot stars, Cool metal-poor giants, Molecular absorption bands, Binary stars, H$\alpha$/H$\beta$ emission, H$\alpha$/H$\beta$ emission superimposed on absorption, H$\alpha$/H$\beta$ P-Cygni, H$\alpha$/H$\beta$ inverted P-Cygni, Lithium absorption, and Problematic. Classified spectra with supplementary information are presented in the catalogue, indicating candidates for follow-up observations and population studies of the short-lived phases of stellar evolution.

We investigated a peculiar metric type II solar radio burst with a broken lane structure, which was observed on November 13, 2012. In addition to the radio data, we also studied the data in the other wavelengths. The bursts were associated with two CMEs and two flares that originated from active region AR 11613. A long current sheet was developed in the first CME, and the second CME collided with the current sheet first and then merged with the first one. Combing information revealed by the multi-wavelength data indicated that a coronal shock accounting for the type II radio burst, and that the collision of this shock with the current sheet resulted in the broken lane of the type II radio burst. The type II burst lane resumed after the shock passed through the current sheet. We further estimated the thickness of the current sheet according to the gap on the lane of the type II burst, and found that the result is consistent with previous ones obtained for various events observed in different wavelengths by different instruments. In addition, the regular type II burst associated with the first CME/flare was also studied, and the magnetic field in each source region of the two type II bursts was further deduced in different way.

We simulate several magnetic reconnection processes in the low solar chromosphere/photosphere, the radiation cooling, heat conduction and ambipolar diffusion are all included. Our numerical results indicate that both the high temperature($ \gtrsim 8\times10^4$~K) and low temperature($\sim 10^4$~K) magnetic reconnection events can happen in the low solar atmosphere ($100\sim600$~km above the solar surface). The plasma $\beta$ controlled by plasma density and magnetic fields is one important factor to decide how much the plasma can be heated up. The low temperature event is formed in a high $\beta$ magnetic reconnection process, Joule heating is the main mechanism to heat plasma and the maximum temperature increase is only several thousand Kelvin. The high temperature explosions can be generated in a low $\beta$ magnetic reconnection process, slow and fast-mode shocks attached at the edges of the well developed plasmoids are the main physical mechanisms to heat the plasma from several thousand Kelvin to over $8\times10^4$~K. Gravity in the low chromosphere can strongly hinder the plasmoind instability and the formation of slow-mode shocks in a vertical current sheet. Only small secondary islands are formed; these islands, however, are not well developed as those in the horizontal current sheets. This work can be applied for understanding the heating mechanism in the low solar atmosphere and could possibly be extended to explain the formation of common low temperature EBs ($\sim10^4$~K) and the high tenperature IRIS bombs ($\gtrsim 8\times10^4$) in the future.

The Galactic Archaeology with HERMES (GALAH) Survey is a massive observational project to trace the Milky Way's history of star formation, chemical enrichment, stellar migration and minor mergers. Using high-resolution (R$\simeq$28,000) spectra taken with the High Efficiency and Resolution Multi-Element Spectrograph (HERMES) instrument at the Anglo-Australian Telescope (AAT), GALAH will determine stellar parameters and abundances of up to 29 elements for up to one million stars. Selecting targets from a colour-unbiased catalogue built from 2MASS, APASS and UCAC4 data, we expect to observe dwarfs at 0.3 to 3 kpc and giants at 1 to 10 kpc. This enables a thorough local chemical inventory of the Galactic thin and thick disks, and also captures smaller samples of the bulge and halo. In this paper we present the plan, process and progress as of early 2016 for GALAH survey observations. In our first two years of survey observing we have accumulated the largest high-quality spectroscopic data set at this resolution, over 200,000 stars. We also present the first public GALAH data catalogue: stellar parameters (Teff, log(g), [Fe/H], [alpha/Fe]), radial velocity, distance modulus and reddening for 10680 observations of 9860 Tycho-2 stars that may be included in the first Gaia data release.

We present the data reduction procedures being used by the GALAH survey, carried out with the HERMES fibre-fed, multi-object spectrograph on the 3.9~m Anglo-Australian Telescope. GALAH is a unique survey, targeting 1 million stars brighter than magnitude V=14 at a resolution of 28,000 with a goal to measure the abundances of 29 elements. Such a large number of high resolution spectra necessitates the development of a reduction pipeline optimized for speed, accuracy, and consistency. We outline the design and structure of the Iraf-based reduction pipeline that we developed, specifically for GALAH, to produce fully calibrated spectra aimed for subsequent stellar atmospheric parameter estimation. The pipeline takes advantage of existing Iraf routines and other readily available software so as to be simple to maintain, testable and reliable. A radial velocity and stellar atmospheric parameter estimator code is also presented, which is used for further data analysis and yields a useful verification of the reduction quality. We have used this estimator to quantify the data quality of GALAH for fibre cross-talk level ($\lesssim0.5$%) and scattered light ($\sim5$ counts in a typical 20 minutes exposure), resolution across the field, sky spectrum properties, wavelength solution reliability (better than $1$ $\mathrm{km\ s^{-1}}$ accuracy) and radial velocity precision.

We present the results of a search for elastic scattering from galactic dark matter in the form of Weakly Interacting Massive Particles (WIMPs) in the 4-30 GeV/$c^2$ mass range. We make use of a 582 kg-day fiducial exposure from an array of 800 g Germanium bolometers equipped with a set of interleaved electrodes with full surface coverage. We searched specifically for $\sim 2.5-20$ keV nuclear recoils inside the detector fiducial volume. As an illustration the number of observed events in the search for 5 (resp. 20) GeV/$c^2$ WIMPs are 9 (resp. 4), compared to an expected background of 6.1 (resp. 1.4). A 90% CL limit of $4.3\times 10^{-40}$ cm$^2$ (resp. $9.4\times 10^{-44}$ cm$^2$) is set on the spin-independent WIMP-nucleon scattering cross-section for 5 (resp. 20) GeV/$c^2$ WIMPs. This result represents a 41-fold improvement with respect to the previous EDELWEISS-II low-mass WIMP search for 7 GeV/$c^2$ WIMPs. The derived constraint is in tension with hints of WIMP signals from some recent experiments, thus confirming results obtained with different detection techniques.

The observational data on the anisotropy of the cosmic microwave background constraints the scalar spectral tilt $n_s$ and the tensor to scalar ratio $r$ which depend on the first and second derivatives of the inflaton potential. The information can be used to reconstruct the inflaton potential in the polynomial form up to some orders. However, for some classes of potentials, $n_s$ and $r$ behave as $n_s(N)$ and $r(N)$ universally in terms of the number of e-folds $N$. The universal behaviour of $n_s(N)$ can be used to reconstruct a class of inflaton potentials. By parametrizing one of the parameters $n_s(N)$, $\epsilon(N)$ and $\phi(N)$, and fitting the parameters in the models to the observational data, we obtain the constraints on the parameters and reconstruct the classes of the inflationary models which include the chaotic inflation, T-model, hilltop inflation, s-dual inflation, natural inflation and $R^2$ inflation.

G-band bright points (GBPs) are regarded as good manifestations of magnetic flux concentrations. We aim to investigate the relationship between the dynamic properties of GBPs and their longitudinal magnetic field strengths. High spatial and temporal resolution observations were recorded simultaneously with G-band filtergrams and Narrow-band Filter Imager (NFI) Stokes I and V images with Hinode /Solar Optical Telescope. The GBPs are identified and tracked in the G-band images automatically, and the corresponding longitudinal magnetic field strength of each GBP is extracted from the calibrated NFI magnetograms by a point-to-point method. After categorizing the GBPs into five groups by their longitudinal magnetic field strengths, we analyze the dynamics of GBPs of each group. The results suggest that with increasing longitudinal magnetic field strengths of GBPs correspond to a decrease in their horizontal velocities and motion ranges as well as by showing more complicated motion paths. This suggests that magnetic elements showing weaker magnetic field strengths prefer to move faster and farther along straighter paths, while stronger ones move more slowly in more erratic paths within a smaller region. The dynamic behaviors of GBPs with different longitudinal magnetic field strengths can be explained by that the stronger flux concentrations withstand the convective flows much better than weaker ones.

In this paper, we study the energy conversion and spectra in a corona current sheet by 2.5-dimensional MHD numerical simulations. Numerical results show that many Petschek-like fine structures with slow-mode shocks mediated by plasmoid instabilities develop during the magnetic reconnection process. The termination shocks can also be formed above the primary magnetic island and at the head of secondary islands. These shocks play important roles in generating thermal energy in a corona current sheet. For a numerical simulation with initial conditions close to the solar corona environment, the ratio of the generated thermal energy to the total dissipated magnetic energy is around $1/5$ before secondary islands appear. After secondary islands appear, the generated thermal energy starts to increase sharply and this ratio can reach a value about $3/5$. In an environment with a relatively lower plasma density and plasma $\beta$, the plasma can be heated to a much higher temperature. After secondary islands appear, the one dimensional energy spectra along the current sheet do not behave as a simple power law and the spectrum index increases with the wave number. The average spectrum index for the magnetic energy spectrum along the current sheet is about $1.8$. The two dimensional spectra intuitively show that part of the high energy is cascaded to large $kx$ and $ky$ space after secondary islands appear. The plasmoid distribution function calculated from numerical simulations behaves as a power law closer to $f(\psi) \sim \psi^{-1}$ in the intermediate $\psi$ regime. By using $\eta_{eff} = v_{inflow}\cdot L$, the effective magnetic diffusivity is estimated about $10^{11}\sim10^{12}$~m$^2$\,s$^{-1}$.

The design and performance of the LUX-ZEPLIN (LZ) detector is described as of March 2015 in this Conceptual Design Report. LZ is a second-generation dark-matter detector with the potential for unprecedented sensitivity to weakly interacting massive particles (WIMPs) of masses from a few GeV/c2 to hundreds of TeV/c2. With total liquid xenon mass of about 10 tonnes, LZ will be the most sensitive experiment for WIMPs in this mass region by the end of the decade. This report describes in detail the design of the LZ technical systems. Expected backgrounds are quantified and the performance of the experiment is presented. The LZ detector will be located at the Sanford Underground Research Facility in South Dakota. The organization of the LZ Project and a summary of the expected cost and current schedule are given.

Magnetic reconnection in the partially ionized solar chromosphere is studied in 2.5-dimensional magnetohydrodynamic simulations including radiative cooling and ambipolar diffusion. A Harris current sheet with and without a guide field is considered. Characteristic values of the parameters in the middle chromosphere imply a high magnetic Reynolds number of $\sim10^{6}\mbox{--}10^7$ in the present simulations. Fast magnetic reconnection then develops as a consequence of the plasmoid instability without the need to invoke anomalous resistivity enhancements. Multiple levels of the instability are followed as it cascades to smaller scales, which approach the ion inertial length. The reconnection rate, normalized to the asymptotic values of magnetic field and Alfv\'en velocity in the inflow region, reaches values in the range $\sim0.01\mbox{--}0.03$ throughout the cascading plasmoid formation and for zero as well as for strong guide field. The out-flow velocity reaches $\approx40$~km\,s$^{-1}$. Slow-mode shocks extend from the $X$-points, heating the plasmoids up to $\sim 8\times10^4$~K. In the case of zero guide field, the inclusion of ambipolar diffusion and radiative cooling both cause a rapid thinning of the current sheet (down to $\sim30$~m) and early formation of secondary islands. Both of these processes have very little effect on the plasmoid instability for a strong guide field. The reconnection rates, temperature enhancements, and upward out-flow velocities from the vertical current sheet correspond well to their characteristic values in chromospheric jets.

We numerically study the detailed evolutionary features of the wave-like disturbance and its propagation in the eruption. This work is a follow-up to Wang et al., using significantly upgraded new simulations. We focus on the contribution of the velocity vortices and the fast shock reflection and refraction in the solar corona to the formation of the EUV waves. Following the loss of equilibrium in the coronal magnetic structure, the flux rope exhibits rapid motions and invokes the fast-mode shock forward of the rope, which then produces the type II radio burst. The expansion of the fast shock, which is associated with outward motion, takes place in various directions, and the downward expansion shows the reflection and the refraction as a result of the non-uniform background plasma. The reflected component of the fast shock propagates upward and the refracted component propagates downward. As the refracted component reaches the boundary surface, a weak echo is excited. The Moreton wave is invoked as the fast shock touches the bottom boundary, so the Moreton wave lags the type II burst. A secondary echo occurs in the area where reflection of the fast shock encounters the slow-mode shock, and the nearby magnetic field lines are further distorted because of the interaction between the secondary echo and the velocity vortices. Our results indicate that the EUV wave may arise from various processes that are revealed in the new simulations.

The GALAH survey is a large high-resolution spectroscopic survey using the newly commissioned HERMES spectrograph on the Anglo-Australian Telescope. The HERMES spectrograph provides high-resolution (R ~28,000) spectra in four passbands for 392 stars simultaneously over a 2 degree field of view. The goal of the survey is to unravel the formation and evolutionary history of the Milky Way, using fossil remnants of ancient star formation events which have been disrupted and are now dispersed throughout the Galaxy. Chemical tagging seeks to identify such dispersed remnants solely from their common and unique chemical signatures; these groups are unidentifiable from their spatial, photometric or kinematic properties. To carry out chemical tagging, the GALAH survey will acquire spectra for a million stars down to V~14. The HERMES spectra of FGK stars contain absorption lines from 29 elements including light proton-capture elements, alpha-elements, odd-Z elements, iron-peak elements and n-capture elements from the light and heavy s-process and the r-process. This paper describes the motivation and planned execution of the GALAH survey, and presents some results on the first-light performance of HERMES.

We present a study of the largest available sample of near-infrared selected (i.e., stellar mass selected) dynamically close pairs of galaxies at low redshifts ($z<0.3$). We combine this sample with new estimates of the major-merger pair fraction for stellar mass selected galaxies at $z<0.8$, from the Red Sequence Cluster Survey (RCS1). We construct our low-redshift $K-$band selected sample using photometry from the UKIRT Infrared Deep Sky Survey (UKIDSS) and the Two Micron All Sky Survey (2MASS) in the $K-$band ($\sim 2.2~\mu$m). Combined with all available spectroscopy, our $K-$band selected sample contains $\sim 250,000$ galaxies and is $> 90\%$ spectroscopically complete. The depth and large volume of this sample allow us to investigate the low-redshift pair fraction and merger rate of galaxies over a wide range in $K-$band luminosity. We find the major-merger pair fraction to be flat at $\sim 2\%$ as a function of $K-$band luminosity for galaxies in the range $10^8 - 10^{12} L_{\odot}$, in contrast to recent results from studies in the local group that find a substantially higher low-mass pair fraction. This low-redshift major-merger pair fraction is $\sim 40-50\%$ higher than previous estimates drawn from $K-$band samples, which were based on 2MASS photometry alone. Combining with the RCS1 sample we find a much flatter evolution ($m = 0.7 \pm 0.1$), in the relation $f_{\rm{pair}} \propto (1+z)^m$, than indicated in many previous studies. These results indicate that a typical $L\sim L^*$ galaxy has undergone $\sim 0.2-0.8$ major mergers since $z=1$ (depending on the assumptions of merger timescale and percentage of pairs that actually merge).

We study the motions of G band bright points (GBPs) in the quiet Sun to obtain the characteristics of different motion types. A high resolution image sequence taken with the Hinode/Solar Optical Telescope (SOT) is used, and GBPs are automatically tracked by segmenting 3D evolutional structures in a space time cube. After putting the GBPs that do not move during their lifetimes aside, the non stationary GBPs are categorized into three types based on an index of motion type. Most GBPs that move in straight or nearly straight lines are categorized into a straight motion type, a few moving in rotary paths into a rotary motion, and the others fall into a motion type we called erratic. The mean horizontal velocity is 2.18 km/s, 1.63 km/s and 1.33 km/s for straight, erratic and rotary motion type, respectively. We find that a GBP drifts at a higher and constant velocity during its whole life if it moves in a straight line. However, it has a lower and variational velocity if it moves in a rotary path. The diffusive process is ballistic, super and sub diffusion for straight, erratic and rotary motion type, respectively. The corresponding diffusion index and coefficients are 2.13 and 850 km2/s, 1.82 and 331 km2/s, 0.73 and 13 km2/s. In terms of direction of motion, it is homogeneous and isotropical, and usually persists between neighbouring frames, no matter what motion type a GBP belongs to.

The onset of a solar eruption is formulated here as either a magnetic catastrophe or as an instability. Both start with the same equation of force balance governing the underlying equilibria. Using a toroidal flux rope in an external bipolar or quadrupolar field as a model for the current-carrying flux, we demonstrate the occurrence of a fold catastrophe by loss of equilibrium for several representative evolutionary sequences in the stable domain of parameter space. We verify that this catastrophe and the torus instability occur at the same point; they are thus equivalent descriptions for the onset condition of solar eruptions.

We investigate equilibrium height of the flux rope, and its internal equilibrium in a realistic plasma environment by carrying out numerical simulations of the evolution of systems including a current-carrying flux rope. We find that the equilibrium height of the flux rope is approximately a power-law function of the relative strength of the background field. Our simulations indicate that the flux rope can escape more easily from a weaker background field. This further confirms the catastrophe in the magnetic configuration of interest can be triggered by decrease of strength of the background field. Our results show that it takes some time to reach internal equilibrium depending on the initial state of the flux rope. The plasma flow inside the flux rope due to the adjustment for the internal equilibrium of the flux rope remains small and does not last very long when the initial state of the flux rope commences from the stable branch of the theoretical equilibrium curve. This work also confirms the influence of the initial radius of flux rope on its evolution, the results indicate that the flux rope with larger initial radius erupts more easily. In addition, by using the realistic plasma environment and much higher resolution in our simulations, we notice some different characters compared to previous studies in Forbes (1990).

Giant flares on soft gamma-ray repeaters that are thought to take place on magnetars release enormous energy in a short time interval. Their power can be explained by catastrophic instabilities occurring in the magnetic field configuration and the subsequent magnetic reconnection. By analogy with the coronal mass ejection (CME) events on the Sun, we develop a theoretical model via an analytic approach for magnetar giant flares. In this model, the rotation and/or displacement of the crust causes the field to twist and deform, leading to flux rope formation in the magnetosphere and energy accumulation in the related configuration. When the energy and helicity stored in the configuration reach a threshold, the system loses its equilibrium, the flux rope is ejected outward in a catastrophic way, and magnetic reconnection helps the catastrophe develop to a plausible eruption. By taking SGR 1806 - 20 as an example, we calculate the free magnetic energy released in such an eruptive process and find that it is more than $10^{47}$ ergs, which is enough to power a giant flare. The released free magnetic energy is converted into radiative energy, kinetic energy and gravitational energy of the flux rope. We calculated the light curves of the eruptive processes for the giant flares of SGR 1806 - 20, SGR 0526-66 and SGR 1900+14, and compared them with the observational data. The calculated light curves are in good agreement with the observed light curves of giant flares.

In this paper we investigate, by means of two-dimensional magnetohydrodynamic simulations, the impact of temperature-dependent resistivity and thermal conduction on the development of plasmoid instabilities in reconnecting current sheets in the solar corona. We find that the plasma temperature in the current sheet region increases with time and it becomes greater than that in the inflow region. As secondary magnetic islands appear, the highest temperature is not always found at the reconnection $X$-points, but also inside the secondary islands. One of the effects of anisotropic thermal conduction is to decrease the temperature of the reconnecting $X-$points and transfer the heat into the $O-$points, the plasmoids, where it gets trapped. In the cases with temperature-dependent magnetic diffusivity, $\eta \sim T^{-3/2}$, the decrease in plasma temperature at the $X-$points leads to: (i) increase in the magnetic diffusivity until the characteristic time for magnetic diffusion becomes comparable to that of thermal conduction; (ii) increase in the reconnection rate; and, (iii) more efficient conversion of magnetic energy into thermal energy and kinetic energy of bulk motions. These results provide further explanation of the rapid release of magnetic energy into heat and kinetic energy seen during flares and coronal mass ejections. In this work, we demonstrate that the consideration of anisotropic thermal conduction and Spitzer-type, temperature-dependent magnetic diffusivity, as in the real solar corona, are crucially important for explaining the occurrence of fast reconnection during solar eruptions.

We study the effects of non-trivial initial quantum states for inflationary fluctuations within the context of the effective field theory for inflation constructed by Cheung et al. which allows us to discriminate between different initial states in a model-independent way. We develop a Green's function/path integral based formulation that incorporates initial state effects and use it to address questions such as how state-dependent is the consistency relation for the bispectrum, how many e-folds beyond the minimum required to solve the cosmological fine tunings of the big bang are we allowed so that some information from the initial state survives until late times, among others. We find that the so-called consistency condition relating the local limit of the bispectrum and the slow-roll parameter is a state-dependent statement that can be avoided for physically consistent initial states either with or without initial non-Gaussianities.

By incorporating the curvature $\Omega_k$ as a free parameter, it has been found that the tension between the high redshift CMB shift parameter $R(z^{\ast})$ data and the low redshift SNIa and BAO data from the combination of SDSS and 2dFGRS can be ameliorated, and both SNIa+BAO and SNIa+BAO+CMB favor that the decelerating parameter $q(z)$ shows a rapid variation in sign at small redshift. In this paper, with the MCMC method, we re-examine the evolutionary behavior of $q(z)$ using the latest observational data including the Union2 SNIa, BAO, and CMB data ($R(z^{\ast})$, $l_{A}(z^{\ast})$, $z^{\ast}$) from WMAP7. For BAO data, four different data sets obtained from the 6dFGS, the combination of SDSS and 2dFGRS, the WiggleZ dark energy survey and the BOSS, are used. Except for the spatially flat case constrained by SNIa+ the WiggleZ BAO data, both SNIa and other BAO+SNIa favor that the present cosmic acceleration is slowing down irrespective of whether the spatial curvature is included or not. However, once the WMAP7 CMB data is included, observations favor strongly the $\Lambda$CDM model, a spatially flat universe, and a speeding-up of the cosmic acceleration. Therefore, the inclusion of spatial curvature seems to has no effect on alleviating the tension between SNIa+BAO and CMB in clear contrast to the previous work in the literature.

The Swift/XRT light curve of the ultraluminous X-ray (ULX) source NGC 5408 X-1 was re-analyzed with two new numerical approaches, Weighted Wavelet $Z$-transform (WWZ) and CLEANest, that are different from previous studies. Both techniques detected a prominent periodicity with a time scale of $115.5\pm1.5$ days, in excellent agreement with the detection of the same periodicity first reported by Strohmayer (2009). Monte Carlo simulation was employed to test the statisiticak confidence of the 115.5-day periodicity, yielding a statistical significance of $> 99.98%$ (or $>3.8\sigma$). The robust detection of the 115.5-day quasi-periodic oscillations (QPOs), if it is due to the orbital motion of the binary, would infer a mass of a few thousand $M_\odot$ for the central black hole, implying an intermediate-mass black hole in NGC 5408 X-1.

We present two-dimensional resistive magnetohydrodynamic simulations of line-tied asymmetric magnetic reconnection in the context of solar flare and coronal mass ejection current sheets. The reconnection process is made asymmetric along the inflow direction by allowing the initial upstream magnetic field strengths and densities to differ, and along the outflow direction by placing the initial perturbation near a conducting wall boundary that represents the photosphere. When the upstream magnetic fields are asymmetric, the post-flare loop structure is distorted into a characteristic skewed candle flame shape. The simulations can thus be used to provide constraints on the reconnection asymmetry in post-flare loops. More hard X-ray emission is expected to occur at the footpoint on the weak magnetic field side because energetic particles are more likely to escape the magnetic mirror there than at the strong magnetic field footpoint. The footpoint on the weak magnetic field side is predicted to move more quickly because of the requirement in two dimensions that equal amounts of flux must be reconnected from each upstream region. The X-line drifts away from the conducting wall in all simulations with asymmetric outflow and into the strong magnetic field region during most of the simulations with asymmetric inflow. There is net plasma flow across the X-line for both the inflow and outflow directions. The reconnection exhaust directed away from the obstructing wall is significantly faster than the exhaust directed towards it. The asymmetric inflow condition allows net vorticity in the rising outflow plasmoid which would appear as rolling motions about the flux rope axis.

In this paper, a periodicity analysis of the radio light curves of the blazar NRAO 530 at 14.5, 8.0, and 4.8 GHz is presented employing an improved Phase Dispersion Minimization (PDM) technique. The result, which shows two persistent periodic components of $ \sim 6$ and $ \sim 10$ years at all three frequencies, is consistent with the results obtained with the Lomb-Scargle periodogram and weighted wavelet Z-transform algorithms. The reliability of the derived periodicities is confirmed by the Monte Carlo numerical simulations which show a high statistical confidence. (Quasi-)Periodic fluctuations of the radio luminosity of NRAO 530 might be associated with the oscillations of the accretion disk triggered by hydrodynamic instabilities of the accreted flow. \keywords{methods: statistical -- galaxies: active -- galaxies: quasar: individual: NRAO 530}

Taking advantage of both the high temporal and spatial resolution of the Atmospheric Imaging Assembly (AIA) on board the Solar Dynamics Observatory (SDO), we studied a limb coronal shock wave and its associated extreme ultraviolet (EUV) wave that occurred on 2010 June 13. Our main findings are (1) the shock wave appeared clearly only in the channels centered at 193 \AA and 211 \AA as a dome-like enhancement propagating ahead of its associated semi-spherical CME bubble; (2) the density compression of the shock is 1.56 according to radio data and the temperature of the shockis around 2.8 MK; (3) the shock wave first appeared at 05:38 UT, 2 minutes after the associated flare has started and 1 minute after its associated CME bubble appeared;(4) the top of the dome-like shock wave set out from about 1.23 R\odot and the thickness of the shocked layer is ~ 2\times10^4 km; (5) the speed of the shock wave is consistent with a slight decrease from about 600 km/s to 550 km/s; (6) the lateral expansion of the shock wave suggests a constant speed around 400 km/s, which varies at different heights and directions. Our findings support the view that the coronal shock wave is driven by the CME bubble, and the on-limb EUV wave is consistent with a fast wave or at least includes the fast wave component.

We have computed an extensive grid of binary evolution tracks to represent low- and intermediate mass X-ray binaries (LMXBs and IMXBs). The grid includes 42,000 models which covers 60 initial donor masses over the range of 1-4 solar masses and, for each of these, 700 initial orbital periods over the range of 10-250 hours. These results can be applied to understanding LMXBs and IMXBs: those that evolve analogously to CVs; that form ultracompact binaries with orbital periods in the range of 6-50 minutes; and that lead to wide orbits with giant donors. We also investigate the relic binary recycled radio pulsars into which these systems evolve. To evolve the donor stars in this study, we utilized a newly developed stellar evolution code called "MESA" that was designed, among other things, to be able to handle very low-mass and degenerate donors. This first application of the results is aimed at an understanding of the newly discovered pulsar PSR J1614-2230 which has a 1.97 solar masses neutron star, orbital period = 8.7 days, and a companion star of 0.5 solar mass. We show that (i) this system is a cousin to the LMXB Cyg X-2; (ii) for neutron stars of canonical birth mass 1.4 solar masses, the initial donor stars which produce the closest relatives to PSR J1614-2230 have a mass between 3.4-3.8 solar masses; (iii) neutron stars as massive as 1.97 solar masses are not easy to produce in spite of the initially high mass of the donor star, unless they were already born as relatively massive neutron stars; (iv) to successfully produce a system like PSR J1614-2230 requires a minimum initial neutron star mass of at least 1.6+-0.1 solar masses, as well as initial donor masses and orbital period of ~ 4.25+-0.10 solar masses and ~49+-2 hrs, respectively; and (v) the current companion star is largely composed of CO, but should have a surface H abundance of ~10-15%.

X-ray observations of EXO 0748-676 during thermonuclear bursts revealed a set of narrow absorption lines that potentially originate from the stellar photosphere. The identification of these lines with particular atomic transitions led to the measurement of the surface gravitational redshift of the neutron star and to constraints on its mass and radius. However, the recent detection of 552 Hz oscillations at 15% rms amplitude revealed the spin frequency of the neutron star and brought into question the consistency of such a rapid spin with the narrow width of the absorption lines. Here, we calculate the amplitudes of burst oscillations and the width of absorption lines emerging from the surface of a rapidly rotating neutron star for a wide range of model parameters. We show that no combination of neutron-star and geometric parameters can simultaneously reproduce the narrowness of the absorption lines, the high amplitude of the oscillations, and the observed flux at the time the oscillations were detected. We, therefore, conclude that the observed absorption lines are unlikely to originate from the surface of this neutron star.

We report the detection of pulsations at 552 Hz in the rising phase of two type-I (thermonuclear) X-ray bursts observed from the accreting neutron star EXO 0748-676 in 2007 January and December, by the Rossi X-ray Timing Explorer. The fractional amplitude was 15% (rms). The dynamic power density spectrum for each burst revealed an increase in frequency of approx. 1-2 Hz while the oscillation was present. The frequency drift, the high significance of the detections and the almost identical signal frequencies measured in two bursts separated by 11 months, confirms this signal as a burst oscillation similar to those found in 13 other sources to date. We thus conclude that the spin frequency in EXO 0748-676 is within a few Hz of 552 Hz, rather than 45 Hz as was suggested from an earlier signal detection by Villarreal & Strohmayer (2004). Consequently, Doppler broadening must significantly affect spectral features arising from the neutron star surface, so that the narrow absorption features previously reported from an XMM-Newton spectrum could not have arisen there. The origin of both the previously reported 45 Hz oscillation and the X-ray absorption lines is now uncertain.

The low mass X-ray binary 4U 2129+47 was discovered during a previous X-ray outburst phase and was classified as an accretion disk corona source. A 1% delay between two mid-eclipse epochs measured ~22 days apart was reported from two XMM-Newton observations taken in 2005, providing support to the previous suggestion that 4U 2129+47 might be in a hierarchical triple system. In this work we present timing and spectral analysis of three recent XMM-Newton observations of 4U 2129+47, carried out between November 2007 and January 2008. We found that absent the two 2005 XMM-Newton observations, all other observations are consistent with a linear ephemeris with a constant period of 18857.63s; however, we confirm the time delay reported for the two 2005 XMM-Newton observations. Compared to a Chandra observation taken in 2000, these new observations also confirm the disappearance of the sinusoidal modulation of the lightcurve as reported from two 2005 XMM-Newton observations. We further show that, compared to the Chandra observation, all of the XMM-Newton observations have 40% lower 0.5--2 keV absorbed fluxes, and the most recent XMM-Newton observations have a combined 2--6 keV flux that is nearly 80% lower. Taken as a whole, the timing results support the hypothesis that the system is in a hierarchical triple system (with a third body period of at least 175 days). The spectral results raise the question of whether the drop in soft X-ray flux is solely attributable to the loss of the hard X-ray tail (which might be related to the loss of sinusoidal orbital modulation), or is indicative of further cooling of the quiescent neutron star after cessation of residual, low-level accretion.

The origin of the Moreton wave observed in the chromosphere and the EIT wave observed in the corona during the eruption remains being an active research subject for a while. We investigate numerically in this work the evolutionary features of the magnetic configuration that includes a current-carrying flux rope, which is used to model the filament, after the loss of equilibrium in the system takes place in a catastrophic fashion. Rapid motions of the flux rope following the catastrophe invokes the velocity vortices behind the rope, and may invoke as well slow and fast mode shocks in front of the rope. The velocity vortices at each side of the flux rope propagate roughly horizontally away from the area where it is produced, and both shocks expand toward the flank of the flux rope. The fast one may eventually reach the bottom boundary and produces two echoes moving back into the corona, but the slow one and the vortices totally decay somewhere in the lower corona before arriving the bottom boundary. The interaction of the fast shock with the boundary leads to disturbance that accounts for the Moreton wave observed in Halpha, and the disturbance in the corona caused by the slow shock and the velocity vortices should account for the EIT wave whose speed is about 40% that of the Moreton wave. Implication of these results to the observed correlation of the type II radio burst to the fast and the slow mode shocks, and that of EIT waves to CMEs and flares have also been discussed.

Episodic ejection of plasma blobs have been observed in many black hole systems. While steady, continuous jets are believed to be associated with large-scale open magnetic fields, what causes the episodic ejection of blobs remains unclear. Here by analogy with the coronal mass ejection on the Sun, we propose a magnetohydrodynamical model for episodic ejections from black holes associated with the closed magnetic fields in an accretion flow. Shear and turbulence of the accretion flow deform the field and result in the formation of a flux rope in the disk corona. Energy and helicity are accumulated and stored until a threshold is reached. The system then loses its equilibrium and the flux rope is thrust outward by the magnetic compression force in a catastrophic way. Our calculations show that for parameters appropriate for the black hole in our Galactic center, the plasmoid can attain relativistic speeds in about 35 minutes.

On 2007 December 7, a small filament located in a small active region AR 10977 erupted and led to a B1.4 flare. An EIT wave associated with this eruption was observed both by SOHO/EIT and by EUVI on board STEREO. According to the observations from SOHO/LASCO and STEREO/COR A, we found that there was no CME associated with the EIT wave. This seems to challenge the argument that the cause of EIT waves is CME. However the data from STEREO/COR B indicated that there was a narrow CME associated with the EIT wave. This suggests that studying CMEs by investigating observations made in one direction alone may not be able to guarantee the reliability of the results.

Magnetic reconnection plays a crucial role in violent energy conversion occurring in the environments of high electrical conductivity, such as the solar atmosphere, magnetosphere, and fusion devices. We focus on the morphological features of the process in two different environments, the solar atmosphere and the geomagnetic tail. In addition to indirect evidence that indicates reconnection in progress or having just taken place, such as auroral manifestations in the magnetosphere and the flare loop system in the solar atmosphere, more direct evidence of reconnection in the solar and terrestrial environments is being collected. Such evidence includes the reconnection inflow near the reconnecting current sheet, and the outflow along the sheet characterized by a sequence of plasmoids. Both turbulent and unsteady Petschek-type reconnection processes could account for the observations. We also discuss other relevant observational consequences of both mechanisms in these two settings. While on face value, these are two completely different physical environments, there emerge many commonalities, for example, an Alfven speed of the same order of magnitude, a key parameter determining the reconnection rate. This comparative study is meant as a contribution to current efforts aimed at isolating similarities in processes occurring in very different contexts in the heliosphere, and even in the universe.

We summarize the unique capabilities of UV coronagraphic spectroscopy for determining the detailed plasma properties (e.g., density, temperature, outflow speed, composition) of the source regions of both transient phenomena such as CMEs, flares, and solar energetic particles (SEPs) and more time-steady solar wind streams. UVCS/SOHO observations have provided the first detailed diagnostics of the physical conditions of CME plasma in the extended corona. It provided new insights into the roles of shock waves, reconnection, and magnetic helicity in CME eruptions. We summarize past observations and discuss the diagnostic potential of UV coronagraphic spectroscopy for characterizing two possible sites of SEP production: CME shocks and reconnection current sheets. UVCS/SOHO has also led to fundamentally new views of the acceleration region of the solar wind. Understanding the physical processes in this region, which ranges from the low corona (1.1 to 1.5 Rsun) past the sonic points (greater than 5 Rsun), is key to linking the results of solar imaging to in situ particle and field detection. Despite the advances that have resulted from UVCS/SOHO, more advanced instrumentation could determine properties of additional ions with a wider sampling of charge/mass combinations. This would provide much better constraints on the specific kinds of waves that are present as well as the specific collisionless damping modes. Electron temperatures and departures from Maxwellian velocity distributions could also be measured. The instrumentation capable of making the above observations will be described.

Glitches are common phenomena in pulsars. After each glitch, there is usually a permanent increase in the pulsar's spin-down rate. Therefore a pulsar's present spin-down rate may be much higher than its initial value. Thus the characteristic age of a pulsar based on its present spin-down rate and period may be shorter than a pulsar's age. At the same time, the permanent increase of its spin-down rate implies that the pulsar's surface magnetic field is increased after each glitch. Consequently after many glitches some radio pulsars may evolve into AXPs and SGRs, i.e., strongly magnetized and slowly rotating neutron stars.

It is known that some observed gamma-ray bursts (GRBs) are produced at cosmological distances and that the GRB production rate may follow the star formation rate. We model the BATSE-detected intensity distribution of long GRBs in order to determine their space density distribution and opening angle distribution. Our main results are: the lower and upper distance limits to the GRB production are z 0.24 and >10, respectively; the GRB opening angle follows an exponential distribution and the mean opening angle is about 0.03 radians; and the peak luminosity appears to be a better standard candle than the total energy of a GRB.