Table of contents

Using detailed numerical simulations, we present the design of an experiment that will generate samples of iron under extreme conditions of density and pressure believed to exist in the interior of the Earth and interior of extrasolar Earth-like planets. In the proposed experiment design, an intense uranium beam is used to implode a multilayered cylindrical target that consists of a thin Fe cylinder enclosed in a thick massive W shell. Such intense uranium beams will be available at the heavy-ion synchrotron, SIS100, at the Facility for Antiprotons and Ion Research (FAIR), at Darmstadt, which is under construction and will become operational in the next few years. It is expected that the beam intensity will increase gradually over a couple of years to its maximum design value. Therefore, in our studies, we have considered a wide range of beam parameters, from the initial beam intensity ("Day One") to the maximum specified value. It is also worth noting that two different focal spot geometries have been used. In one case, a circular focal spot with a Gaussian transverse intensity distribution is considered, whereas in the other case, an annular focal spot is used. With these two beam geometries, one can access different parts of the Fe phase diagram. For example, heating the sample with a circular focal spot generates a hot liquid state, while an annular focal spot can produce a highly compressed liquid or a highly compressed solid phase depending on the beam intensity.

We present estimates of stellar age and mass for 0.93 million Galactic-disk main-sequence turnoff and subgiant stars from the LAMOST Galactic Spectroscopic Surveys. The ages and masses are determined by matching with stellar isochrones using a Bayesian algorithm, utilizing effective temperature ${T}_{\mathrm{eff}}$, absolute magnitude ${M}_{V}$, metallicity [Fe/H], and α-element to iron abundance ratio [α/Fe] deduced from the LAMOST spectra. Extensive examinations suggest the age and mass estimates are robust. Overall, the sample stars have a median error of 34% for the age estimates, and half of the stars older than 2 Gyr have age uncertainties of only 20%–30%. The median error for the mass estimates of the whole sample of stars is ∼8%. The huge data set demonstrates good correlations among stellar age, [Fe/H] ([α/H]), and [α/Fe]. Particularly, double-sequence features are revealed in both the age–[α/Fe] and age–[Fe/H]([α/H]) spaces. In the [Fe/H]–[α/Fe] space, stars of 8–10 Gyr exhibit both the thin and thick disk sequences, while younger (older) stars show only the thin (thick) disk sequence, indicating that the thin disk became prominent 8–10 Gyr ago, while the thick disk formed earlier and was almost quenched 8 Gyr ago. Stellar ages exhibit positive vertical and negative radial gradients across the disk, and the outer disk of R ≳ 9 kpc exhibits a strong flare in stellar age distribution.

Spectral line surveys are an indispensable tool for exploring the physical and chemical evolution of astrophysical environments due to the vast amount of data that can be obtained in a relatively short amount of time. We present deep, broadband spectral line surveys of 30 interstellar clouds using two broadband λ = 1.3 mm receivers at the Caltech Submillimeter Observatory. This information can be used to probe the influence of physical environment on molecular complexity. We observed a wide variety of sources to examine the relative abundances of organic molecules as they relate to the physical properties of the source (i.e., temperature, density, dynamics, etc.). The spectra are highly sensitive, with noise levels ≤25 mK at a velocity resolution of ∼0.35 km s⁻¹. In the initial analysis presented here, column densities and rotational temperatures have been determined for the molecular species that contribute significantly to the spectral line density in this wavelength regime. We present these results and discuss their implications for complex molecule formation in the interstellar medium.

We combine Gaia DR1, PS1, Sloan Digital Sky Survey (SDSS), and 2MASS astrometry to measure proper motions for 350 million sources across three-fourths of the sky down to a magnitude of ${m}_{r}\sim 20$. Using positions of galaxies from PS1, we build a common reference frame for the multi-epoch PS1, single-epoch SDSS and 2MASS data, and calibrate the data in small angular patches to this frame. As the Gaia DR1 excludes resolved galaxy images, we choose a different approach to calibrate its positions to this reference frame: we exploit the fact that the proper motions of stars in these patches are linear. By simultaneously fitting the positions of stars at different epochs of—Gaia DR1, PS1, SDSS, and 2MASS—we construct an extensive catalog of proper motions dubbed GPS1. GPS1 has a characteristic systematic error of less than 0.3 $\mathrm{mas}\,{\mathrm{yr}}^{-1}$ and a typical precision of 1.5–2.0 $\mathrm{mas}\,{\mathrm{yr}}^{-1}$. The proper motions have been validated using galaxies, open clusters, distant giant stars, and QSOs. In comparison with other published faint proper motion catalogs, GPS1's systematic error ($\lt 0.3\,\mathrm{mas}\,{\mathrm{yr}}^{-1}$) should be nearly an order of magnitude better than that of PPMXL and UCAC4 ($\gt 2.0\,\mathrm{mas}\,{\mathrm{yr}}^{-1}$). Similarly, its precision (∼1.5 $\mathrm{mas}\,{\mathrm{yr}}^{-1}$) is a four-fold improvement relative to PPMXL and UCAC4 (∼6.0 $\mathrm{mas}\,{\mathrm{yr}}^{-1}$). For QSOs, the precision of GPS1 is found to be worse (∼2.0–3.0 $\mathrm{mas}\,{\mathrm{yr}}^{-1}$), possibly due to their particular differential chromatic refraction. The GPS1 catalog will be released online and be available via the VizieR Service and VO Service.

There have been several attempts in the past to understand the nature of the collision of individual cases of interacting coronal mass ejections (CMEs). We selected eight cases of interacting CMEs and estimated their propagation and expansion speeds, and direction of impact and masses, by exploiting coronagraphic and heliospheric imaging observations. Using these estimates while ignoring the errors therein, we find that the nature of collisions is perfectly inelastic for two cases (i.e., 2012 March and November), inelastic for two cases (i.e., 2012 June and 2011 August), elastic for one case (i.e., 2013 October), and super-elastic for three cases (i.e., 2011 February, 2010 May, and 2012 September). Including the large uncertainties in the estimated directions, angular widths, and pre-collision speeds, the probability of a perfectly inelastic collision for the 2012 March and November cases drops from 98% to 60% and 100% to 40%, respectively, increasing the probability for other types of collision. Similarly, the probability of an inelastic collision drops from 95% to 50% for the 2012 June case, 85% to 50% for the 2011 August case, and 75% to 15% for the 2013 October case. We note that the probability of a super-elastic collision for the 2011 February, 2010 May, and 2012 September CMEs drops from 90% to 75%, 60% to 45%, and 90% to 50%, respectively. Although the sample size is small, we find good dependence of the nature of collision on the CME parameters. The crucial pre-collision parameters of the CMEs responsible for increasing the probability of a super-elastic collision are, in descending order of priority, their lower approaching speed, expansion speed of the following CME higher than the preceding one, and a longer duration of the collision phase.

The relative orientation between filamentary structures in molecular clouds and the ambient magnetic field provides insight into filament formation and stability. To calculate the relative orientation, a measurement of filament orientation is first required. We propose a new method to calculate the orientation of the one-pixel-wide filament skeleton that is output by filament identification algorithms such as filfinder. We derive the local filament orientation from the direction of the intensity gradient in the skeleton image using the Sobel filter and a few simple post-processing steps. We call this the "Sobel-gradient method." The resulting filament orientation map can be compared quantitatively on a local scale with the magnetic field orientation map to then find the relative orientation of the filament with respect to the magnetic field at each point along the filament. It can also be used for constructing radial profiles for filament width fitting. The proposed method facilitates automation in analyses of filament skeletons, which is imperative in this era of "big data."

The results of 8 yr R-band photopolarimetric data of blazar Mrk 421 collected from 2008 February to 2016 May are presented, along with extensive multiwavelength observations covering radio to TeV γ-rays around the flares observed in 2008 May, 2010 March, and 2013 April. The most important results are found in 2013, when the source displayed in the R band a very high brightness state of 11.29 ± 0.03 mag (93.60 ± 1.53 mJy) on April 10 and a polarization degree of 11.00% ± 0.44% on May 13. The analysis of the optical data shows that the polarization variability is due to the superposition of two polarized components that might be produced in two distinct emitting regions. An intranight photopolarimetric variability study carried out over seven nights after the 2013 April maximum found flux and polarization variations on the nights of April 14, 15, 16, and 19. In addition, the flux shows a minimum variability timescale of Δt = 2.34 ± 0.12 hr, and the polarization degree presents variations of ∼1%–2% on a timescale of Δ t ∼ minutes. Also, a detailed analysis of the intranight data shows a coherence length of the large-scale magnetic field of ${l}_{B}\simeq 0.3$ pc, which is the same order of magnitude as the distance traveled by the relativistic shocks. This result suggests that there is a connection between the intranight polarimetric variations and spatial changes of the magnetic field. Analysis of the complete R-band data along with the historical optical light curve found for this object shows that Mrk 421 varies with a period of 16.26 ± 1.78 yr.

We present a detailed spectral analysis of the brightest active galactic nuclei (AGNs) identified in the 7Ms Chandra Deep Field-South (CDF-S) survey over a time span of 16 years. Using a model of an intrinsically absorbed power-law plus reflection, with possible soft excess and narrow Fe Kα line, we perform a systematic X-ray spectral analysis, both on the total 7Ms exposure and in four different periods with lengths of 2–21 months. With this approach, we not only present the power-law slopes, column densities ${N}_{{\rm{H}}}$, observed fluxes, and absorption-corrected 2–10 keV luminosities L_X for our sample of AGNs, but also identify significant spectral variabilities among them on timescales of years. We find that the ${N}_{{\rm{H}}}$ variabilities can be ascribed to two different types of mechanisms, either flux-driven or flux-independent. We also find that the correlation between the narrow Fe line EW and ${N}_{{\rm{H}}}$ can be well explained by the continuum suppression with increasing ${N}_{{\rm{H}}}$. Accounting for the sample incompleteness and bias, we measure the intrinsic distribution of ${N}_{{\rm{H}}}$ for the CDF-S AGN population and present reselected subsamples that are complete with respect to ${N}_{{\rm{H}}}$. The ${N}_{{\rm{H}}}$-complete subsamples enable us to decouple the dependences of ${N}_{{\rm{H}}}$ on L_X and on redshift. Combining our data with those from C-COSMOS, we confirm the anticorrelation between the average ${N}_{{\rm{H}}}$ and L_X of AGN, and find a significant increase of the AGN-obscured fraction with redshift at any luminosity. The obscured fraction can be described as ${f}_{\mathrm{obscured}}\approx 0.42\ {(1+z)}^{0.60}$.

Spontaneous breaking of Lorentz symmetry at energies on the order of the Planck energy or lower is predicted by many quantum gravity theories, implying non-trivial dispersion relations for the photon in vacuum. Consequently, gamma-rays of different energies, emitted simultaneously from astrophysical sources, could accumulate measurable differences in their time of flight until they reach the Earth. Such tests have been carried out in the past using fast variations of gamma-ray flux from pulsars, and more recently from active galactic nuclei and gamma-ray bursts. We present new constraints studying the gamma-ray emission of the galactic Crab Pulsar, recently observed up to TeV energies by the Major Atmospheric Gamma-ray Imaging Cherenkov (MAGIC) collaboration. A profile likelihood analysis of pulsar events reconstructed for energies above 400 GeV finds no significant variation in arrival time as their energy increases. Ninety-five percent CL limits are obtained on the effective Lorentz invariance violating energy scale at the level of ${E}_{{\mathrm{QG}}_{1}}\gt 5.5\times {10}^{17}\,\mathrm{GeV}$ ($4.5\times {10}^{17}\,\mathrm{GeV}$) for a linear, and ${E}_{{\mathrm{QG}}_{2}}\gt 5.9\times {10}^{10}\,\mathrm{GeV}$ ($5.3\times {10}^{10}\,\mathrm{GeV}$) for a quadratic scenario, for the subluminal and the superluminal cases, respectively. A substantial part of this study is dedicated to calibration of the test statistic, with respect to bias and coverage properties. Moreover, the limits take into account systematic uncertainties, which are found to worsen the statistical limits by about 36%–42%. Our constraints would have been much more stringent if the intrinsic pulse shape of the pulsar between 200 GeV and 400 GeV was understood in sufficient detail and allowed inclusion of events well below 400 GeV.

We investigate the nature of the extragalactic unresolved γ-ray background (UGRB) by cross-correlating several galaxy catalogs with sky maps of the UGRB built from 78 months of Pass 8 Fermi-Large Area Telescope data. This study updates and improves similar previous analyses in several aspects. First, the use of a larger γ-ray data set allows us to investigate the energy dependence of the cross-correlation in more detail, using up to eight energy bins over a wide energy range of [0.25,500] GeV. Second, we consider larger and deeper catalogs (2MASS Photometric Redshift catalog, 2MPZ; WISE × SuperCOSMOS, WI×SC; and SDSS DR12 photometric redshift data set) in addition to the ones employed in the previous studies (NVSS and SDSS QSOs). Third, we exploit the redshift information available for the above catalogs to divide them into redshift bins and perform the cross-correlation separately in each of them. Our results confirm, with higher statistical significance, the detection of cross-correlation signals between the UGRB maps and all the catalogs considered, on angular scales smaller than 1°. Significances range from $16.3\sigma$ for NVSS, $7\sigma$ for SDSS DR12 and WI×SC, to $5\sigma$ for 2MPZ and $4\sigma$ for SDSS QSOs. Furthermore, including redshift tomography, the significance of the SDSS DR12 signal strikingly rises up to $\sim 12\sigma$ and that of WI×SC to $\sim 10.6\sigma$. We offer a simple interpretation of the signal in the framework of the halo model. The precise redshift and energy information allows us to clearly detect a change over redshift in the spectral and clustering behavior of the γ-ray sources contributing to the UGRB.

We present the second and final data release of the Siding Spring Southern Seyfert Spectroscopic Snapshot Survey (S7). Data are presented for 63 new galaxies not included in the first data release, and we provide 2D emission-line fitting products for the full S7 sample of 131 galaxies. The S7 uses the WiFeS instrument on the ANU 2.3 m telescope to obtain spectra with a spectral resolution of R = 7000 in the red (540–700 nm) and R = 3000 in the blue (350–570 nm), over an integral field of 25 × 38 arcsec² with 1 × 1 arcsec² spatial pixels. The S7 contains both the largest sample of active galaxies and the highest spectral resolution of any comparable integral field survey to date. The emission-line fitting products include line fluxes, velocities, and velocity dispersions across the WiFeS field of view, and an artificial neural network has been used to determine the optimal number of Gaussian kinematic components for emission-lines in each spaxel. Broad Balmer lines are subtracted from the spectra of nuclear spatial pixels in Seyfert 1 galaxies before fitting the narrow lines. We bin nuclear spectra and measure reddening-corrected nuclear fluxes of strong narrow lines for each galaxy. The nuclear spectra are classified on optical diagnostic diagrams, where the strength of the coronal line [Fe vii] λ6087 is shown to be correlated with [O iii]/Hβ. Maps revealing gas excitation and kinematics are included for the entire sample, and we provide notes on the newly observed objects.

Tracking the thermal evolution of plasmas, characterized by an n-distribution, using numerical simulations, requires the determination of the emission spectra and of the radiative losses due to free–free emission from the corresponding temperature-averaged and total Gaunt factors. Detailed calculations of the latter are presented and associated with n-distributed electrons with the parameter n ranging from 1 (corresponding to the Maxwell–Boltzmann distribution) to 100. The temperature-averaged and total Gaunt factors with decreasing n tend toward those obtained with the Maxwell–Boltzmann distribution. Radiative losses due to free–free emission in a plasma evolving under collisional ionization equilibrium conditions and composed by H, He, C, N, O, Ne, Mg, Si, S, and Fe ions, are presented. These losses decrease with a decrease in the parameter n, reaching a minimum when n = 1, and thus converge with the loss of thermal plasma. Tables of the thermal-averaged and total Gaunt factors calculated for n-distributions, and a wide range of electron and photon energies, are presented.

We performed simultaneous observations of SiO v = 1, 2, ²⁹SiO v = 0, J = 1–0, and H₂O ${6}_{12}\to {5}_{23}$ maser lines toward 252 OH/IR stars using the individual 21 m telescopes of the Korean VLBI Network (KVN). The observations for studying SiO and H₂O maser properties associated with the different evolutionary stages of OH/IR stars were carried out from 2011 November to 2012 July. Both H₂O and SiO masers were detected from 50 sources with a detection rate of 20% in one epoch of observation. One-sided SiO maser emissions without H₂O were detected from 108 sources and H₂O maser emission was detected from 11 sources, of which the detection rates were 43% and 4%, respectively. The overall detection rate of the SiO maser was 63%, and that of the H₂O maser was 24%. There were 65 new detections in the SiO maser lines, 22 new detections in the H₂O maser line, and 4 new detections in the ²⁹SiO maser line. For the H₂O and/or SiO maser- detected sources, mutual relations between SiO and H₂O maser properties (including peak and integrated antenna temperatures, and full widths at zero power, etc.) are investigated based on a statistical analysis. We also investigate these maser properties on an IRAS two-color diagram related to stellar evolutionary sequences. In particular, a large number of SiO v = 2-only detected sources appear among the SiO-only detected sources compared to those of both H₂O and SiO maser detected sources and also appear in the later evolutionary stages of asymptotic giant branch AGB stars in the IRAS two-color diagram. These results may be associated with the development of a hot and thick dust envelope at later stages of AGB evolution and the different excitation conditions of SiO v = 1 and v = 2 masers. Our observational results will be useful for statistical studies of circumstellar envelopes of OH/IR stars related to their late evolution and future very long baseline interferometry (VLBI) observations.

Will, in 1974, treated the perturbation of a Schwarzschild black hole due to a slowly rotating, light, concentric thin ring by solving the perturbation equations in terms of a multipole expansion of the mass-and-rotation perturbation series. In the Schwarzschild background, his approach can be generalized to perturbation by a thin disk (which is more relevant astrophysically), but, due to rather bad convergence properties, the resulting expansions are not suitable for specific (numerical) computations. However, we show that Green's functions, represented by Will's result, can be expressed in closed form (without multipole expansion), which is more useful. In particular, they can be integrated out over the source (a thin disk in our case) to yield good converging series both for the gravitational potential and for the dragging angular velocity. The procedure is demonstrated, in the first perturbation order, on the simplest case of a constant-density disk, including the physical interpretation of the results in terms of a one-component perfect fluid or a two-component dust in a circular orbit about the central black hole. Free parameters are chosen in such a way that the resulting black hole has zero angular momentum but non-zero angular velocity, as it is just carried along by the dragging effect of the disk.

We report the results of searches for strong gravitational lens systems in the Dark Energy Survey (DES) Science Verification and Year 1 observations. The Science Verification data span approximately 250 sq. deg. with a median i-band limiting magnitude for extended objects (10σ) of 23.0. The Year 1 data span approximately 2000 sq. deg. and have an i-band limiting magnitude for extended objects (10σ) of 22.9. As these data sets are both wide and deep, they are particularly useful for identifying strong gravitational lens candidates. Potential strong gravitational lens candidate systems were initially identified based on a color and magnitude selection in the DES object catalogs or because the system is at the location of a previously identified galaxy cluster. Cutout images of potential candidates were then visually scanned using an object viewer and numerically ranked according to whether or not we judged them to be likely strong gravitational lens systems. Having scanned nearly 400,000 cutouts, we present 374 candidate strong lens systems, of which 348 are identified for the first time. We provide the R.A. and decl., the magnitudes and photometric properties of the lens and source objects, and the distance (radius) of the source(s) from the lens center for each system.

Based on an extensive spectral study of a photometrically confirmed sample of Mira variables, we find a relationship between the relative Balmer emission-line strengths and spectral temperatures of O-rich Mira stars. The ${F}_{{\rm{H}}\delta }/{F}_{{\rm{H}}\gamma }$ flux ratio increases from less than unity to five as stars cool down from M0 to M10, which is likely driven by increasing TiO absorption above the deepest shock-emitting regions. We also discuss the relationship between the equivalent widths of the Balmer emission lines and the photometric luminosity phase of our Mira sample stars. Using our 291 Mira spectra as templates for reference, 191 Mira candidates are newly identified from the LAMOST DR4 catalog. We summarize the criteria adopted to select Mira candidates based on emission-line indices and molecular absorption bands. This enlarged spectral sample of Mira variables has the potential to contribute significantly to our knowledge of the optical properties of Mira stars and will facilitate further studies of these late-type, long-period variables.

We present our broadband (2–250 keV) time-averaged spectral analysis of 388 bursts from SGR J1550−5418, SGR 1900+14, and SGR 1806−20 detected with the Rossi X-ray Timing Explorer (RXTE) here and as a database in a companion web-catalog. We find that two blackbody functions (BB+BB), the sum of two modified blackbody functions (LB+LB), the sum of a blackbody function and a power-law function (BB+PO), and a power law with a high-energy exponential cutoff (COMPT) all provide acceptable fits at similar levels. We performed numerical simulations to constrain the best fitting model for each burst spectrum and found that 67.6% of burst spectra with well-constrained parameters are better described by the Comptonized model. We also found that 64.7% of these burst spectra are better described with the LB+LB model, which is employed in the spectral analysis of a soft gamma repeater (SGR) for the first time here, than with the BB+BB and BB+PO models. We found a significant positive lower bound trend on photon index, suggesting a decreasing upper bound on hardness, with respect to total flux and fluence. We compare this result with bursts observed from SGR and AXP (anomalous X-ray pulsar) sources and suggest that the relationship is a distinctive characteristic between the two. We confirm a significant anticorrelation between burst emission area and blackbody temperature, and find that it varies between the hot and cool blackbody temperatures differently than previously discussed. We expand on the interpretation of our results in the framework of a strongly magnetized neutron star.