We report the detection of four new hot corino sources, G211.47-19.27S, G208.68-19.20N1, G210.49-19.79W and G192.12-11.10 from a survey study of Planck Galactic Cold Clumps in the Orion Molecular Cloud Complex with the Atacama Compact Array (ACA). Three sources had been identified as low mass Class 0 protostars in the Herschel Orion Protostar Survey (HOPS). One source in the lambda Orionis region is firstly reported as a protostellar core. We have observed abundant complex organic molecules (COMs), primarily methanol but also other oxygen-bearing COMs (in G211.47-19.27S and G208.68-19.20N1) and the molecule of prebiotic interest NH2CHO (in G211.47-19.27S), signifying the presence of hot corinos. While our spatial resolution is not sufficient for resolving most of the molecular emission structure, the large linewidth and high rotational temperature of COMs suggest that they likely reside in the hotter and innermost region immediately surrounding the protostar. In G211.47-19.27S, the D/H ratio of methanol ([CH2DOH]/[CH3OH]) and the 12C/13C ratio of methanol ([CH3OH]/[13CH3OH]) are comparable to those of other hot corinos. Hydrocarbons and long carbon-chain molecules such as c-C3H2 and HCCCN are also detected in the four sources, likely tracing the outer and cooler molecular envelopes.

One of the unsolved questions in gamma-ray astronomy is whether the extragalactic gamma-ray background is of the discrete-source origin. To respond to this question, one first needs to reduce the data by differentiating charged particles from gamma rays. This procedure is usually performed on the basis of the detector responses. In this paper, we showed that the geomagnetic shielding effect at GeV energies can, to some extent, be used for this purpose for gamma-ray telescopes in a low Earth orbit. We illustrated this method by applying it to the Fermi-LAT data. To partially decompose the charge-filtered background, we examined the contribution from star-forming galaxies by implying a radio/gamma connection in consideration of next generation radio surveys.

A systematic search for cyclical sources of gamma-ray emission on the period range from days to years in the Fermi-LAT sky is performed. Looking for cyclical emission, the sky is binned into equal-area pixels and the generalised Lomb-Scargle periodogram is computed for each of these pixels. The search on the period range between 2.5 and 30 days in the Galactic plane confirms periodicities of three binaries, LSI +61 303, LS 5039, and 1FGL J1018.6-5856. The all-sky search on the period range between 30 days and 2.5 years confirms periodicities of three blazars, PG 1553+113, PKS 2155-304, and BL Lacertae. Evidence for periodic behaviours of four blazars, 4C +01.28, S5 0716+71, PKS 0805-07, and PKS 2052-47, are presented. Three of these blazars, 4C +01.28, PKS 0805-07, and PKS 2052-47, are located at high redshifts. These three sources are potential candidates to binary systems of supermassive black holes provided that major galaxy mergers are more frequent and that galaxies are more gas-rich at high redshifts.

We develop a method to search for pair halos around active galactic nuclei (AGN) through a temporal analysis of gamma-ray data. The basis of our method is an analysis of the spatial distributions of photons coming from AGN flares and from AGN quiescent states and a further comparison of these two spatial distributions. This method can also be used for a reconstruction of a point spread function (PSF). We found no evidence for a pair halo component through this method by applying it to the Fermi-LAT data in the energy bands of 4.5-6, 6-10, and >10 GeV and set upper limits on the fraction of photons attributable to a pair halo component. An illustration of how to reconstruct the PSF of Fermi-LAT is given. We demonstrate that the PSF reconstructed by using this method is in good agreement with that which was obtained by using the gamma-ray data taken by LAT in the direction of the Crab pulsar and nebula.

Recent high-resolution high-sensitivity observations of protostellar jets have shown many to possess deviations to their trajectories. HH 211 is one such example where sub-mm observations with the SMA have revealed a clear reflection-symmetric wiggle. The most likely explanation is that the HH 211 jet source could be moving as part of a protobinary system. Here we test this assumption by simulating HH 211 through 3D hydrodynamic jet propagation simulations using the PLUTO code with a molecular chemistry and cooling module, and initial conditions based on an analytical model derived from SMA observations. Our results show the reflection-symmetric wiggle can be recreated through the assumption of a jet source perturbed by binary motion at its base, and that a regular sinusoidal velocity variation in the jet beam can be close to matching the observed knot pattern. However, a more complex model with either additional heating from the protostar, or a shorter period velocity pulsation may be required to account for enhanced emission near the source, and weaker knot emission downstream. Position velocity diagrams along the pulsed jet beam show a complex structure with detectable signatures of knots and show caution must be exercised when interpreting radial velocity profiles through observations. Finally, we make predictions for future HH 211 observations with ALMA.

In order to understand the formation of the multipolar structures of the pre-planetary nebula (PPN) CRL 618, we perform 3D simulations using a multi-directional bullet model. The optical lobes of CRL 618 and fast molecular outflows at the tips of the lobes have been found to have similar expansion ages of ~ 100 yr. Additional fast molecular outflows were found near the source along the outflow axes with ages of ~ 45 yr, suggesting a second episode of bullet ejections. Thus, in our simulations, two episodes of bullet ejections are assumed. The shaping process is simulated using the ZEUS-3D hydrodynamics code that includes molecular and atomic cooling. In addition, molecular chemistry is also included to calculate the CO intensity maps. Our results show the following: (1) Multi-epoch bullets interacting with the toroidal dense core can produce the collimated multiple lobes as seen in CRL 618. The total mass of the bullets is ~ 0.034 solar mass, consistent with the observed high-velocity CO emission in fast molecular outflows. (2) The simulated CO J=3-2 intensity maps show that the low-velocity cavity wall and the high-velocity outflows along the lobes are reasonably consistent with the observations. The position-velocity diagram of the outflows along the outflow axes shows a linear increase of velocity with distance, similar to the observations. The ejections of these bullets could be due to magneto-rotational explosions or nova-like explosions around a binary companion.

We investigate the power spectra of outflow-driven turbulence through high-resolution three-dimensional isothermal numerical simulations where the turbulence is driven locally in real-space by a simple spherical outflow model. The resulting turbulent flow saturates at an average Mach number of ~2.5 and is analysed through density and velocity power spectra, including an investigation of the evolution of the solenoidal and compressional components. We obtain a shallow density power spectrum with a slope of ~-1.2 attributed to the presence of a network of localised dense filamentary structures formed by strong shock interactions. The total velocity power spectrum slope is found to be ~-2.0, representative of Burgers shock dominated turbulence model. The density weighted velocity power spectrum slope is measured as ~-1.6, slightly less than the expected Kolmogorov scaling value (slope of -5/3) found in previous works. The discrepancy may be caused by the nature of our real space driving model and we suggest there is no universal scaling law for supersonic compressible turbulence. We find that on average, solenoidal modes slightly dominate in our turbulence model as the interaction between strong curved compressible shocks generates solenoidal modes, and compressible modes decay faster.

Globular Cluster (GC) Color-Magnitude Diagrams (CMDs) are reasonably well understood in terms of standard stellar-evolution. However, there are still some open issues, such as fully accounting for the Horizontal Branch (HB) morphology in terms of chemical and dynamical parameters. Mass-loss on the Red Giant Branch (RGB) shapes the mass-distribution of the HB stars, and the color distribution in turn. The physical mechanisms driving mass-loss are still unclear, as direct observations fail to reveal a clear correlation between mass-loss rate and stellar properties. The horizontal-branch mass-distribution is further complicated by Helium-enhanced multiple stellar populations, because of differences in the evolving mass along the HB. We present a simple analytical mass-loss model, based on tidal stripping through Roche-Lobe OverFlow (RLOF) during stellar encounters. Our model naturally results in a non-gaussian mass-loss distribution, with high skewness, and contains only two free parameters. We fit it to the HB mass distribution of four Galactic GCs, as obtained from fitting the CMD with Zero Age HB (ZAHB) models. The best-fit model accurately reproduces the observed mass-distribution. If confirmed on a wider sample of GCs, our results would account for the effects of dynamics in RGB mass-loss processes and provide a physically motivated procedure for synthetic CMDs of GCs. Our physical modeling of mass-loss may result in the ability to disentangle the effects of dynamics and helium-enhanced multiple-populations on the HB morphology and is instrumental in making HB morphology a probe of the dynamical state of GCs, leading to an improved understanding of their evolution.

Context. Stellar collision rates in globular clusters (GCs) do not appear to correlate with horizontal branch (HB) morphology, sug- gesting that dynamics does not play a role in the second-parameter problem. However, core densities and collision rates derived from surface-brightness may be significantly underestimated as the surface-brightness profile of GCs is not necessarily a good indicator of the dynamical state of GC cores. Core-collapse may go unnoticed if high central densities of dark remnants are present. Aims. We test whether GC HB morphology data supports a dynamical contribution to the so-called second-parameter effect. Methods. To remove first-parameter dependence we fitted the maximum effective temperature along the HB as a function of metal- licity in a sample of 54 Milky Way GCs. We plotted the residuals to the fit as a function of second-parameter candidates, namely dynamical age and total luminosity. We considered dynamical age (i.e. the ratio between age and half-light relaxation time) among possible second-parameters. We used a set of direct N-body simulations, including ones with dark remnants to illustrate how core density peaks, due to core collapse, in a dynamical-age range similar to that in which blue HBs are overabundant with respect to the metallicity expectation, especially for low-concentration initial conditions. Results. GC total luminosity shows nonlinear behavior compatible with the self-enrichment picture. However, the data are amenable to a different interpretation based on a dynamical origin of the second-parameter effect. Enhanced mass-stripping in the late red-giant- branch phase due to stellar interactions in collapsing cores is a viable candidate mechanism. In this picture, GCs with HBs bluer than expected based on metallicity are those undergoing core-collapse.

Protostellar jets and outflows are signatures of star formation and promising mechanisms for driving supersonic turbulence in molecular clouds. We quantify outflow-driven turbulence through three-dimensional numerical simulations using an isothermal version of the robust total variation diminishing code. We drive turbulence in real-space using a simplified spherical outflow model, analyse the data through density probability distribution functions (PDF), and investigate the Core Formation Rate per free-fall time (CFR_ff). The real-space turbulence driving method produces a negatively skewed density PDF possessing an enhanced tail on the low-density side. It deviates from the log-normal distributions typically obtained from Fourier-space turbulence driving at low densities, but can provide a good fit at high-densities, particularly in terms of mass weighted rather than volume weighted density PDF. Due to this fact, we suggest that the CFR_ff determined from a Fourier-driven turbulence model could be comparable to that of our particular real-space driving model, which has a ratio of solenoidal to compressional components from the resulting turbulence velocity fields of ~0.6.

We present a study of the observational properties of Millisecond Pulsars (MSPs) by way of their magnetic fields, spin periods and masses. These measurements are derived through the scenario of Accretion Induced Collapse (AIC) of white dwarfs (WDs) in stellar binary systems, in order to provide a greater understanding of the characteristics of MSP populations. In addition, we demonstrate a strong evolutionary connection between neutron stars and WDs with binary companions from a stellar binary evolution perspective via the AIC process.

We perform relativistic hydrodynamic simulations of the formation and evolution of AGN cocoons produced by very light powerful jets. We calculate the intensity maps of the Sunyaev-Zel'dovich (SZ) effect at high frequencies for the simulated AGN cocoons using the relativistically correct Wright formalism. Our fully relativistic calculations demonstrate that the contribution from the high temperature gas (kb Te ~ 100 keV) to the SZ signal from AGN cocoons at high frequencies is stronger than that from the shocked ambient intercluster medium owing to the fact that the relativistic spectral functions peak at these temperature values. We present simulations of the SZ effect from AGN cocoons at various frequencies, and demonstrate that SZ observations at 217 GHz and at higher frequencies, such as 857 GHz, will provide us with knowledge about the dynamically-dominant component of AGN cocoons.