We introduce new analysis methods for studying the star cluster formation processes in Orion A, especially examining the scenario of a cloud-cloud collision. We utilize the CARMA-NRO Orion survey $^{13}$CO (1-0) data to compare molecular gas to the properties of YSOs from the SDSS III IN-SYNC survey. We show that the increase of $v_{\rm 13CO} - v_{\rm YSO}$ and $\Sigma$ scatter of older YSOs can be signals of cloud-cloud collision. SOFIA-upGREAT 158$\mu$m [CII] archival data toward the northern part of Orion A are also compared to the $^{13}$CO data to test whether the position and velocity offsets between the emission from these two transitions resemble those predicted by a cloud-cloud collision model. We find that the northern part of Orion A, including regions ONC-OMC-1, OMC-2, OMC-3 and OMC-4, shows qualitative agreements with the cloud-cloud collision scenario, while in one of the southern regions, NGC1999, there is no indication of such a process in causing the birth of new stars. On the other hand, another southern cluster, L1641N, shows slight tendencies of cloud-cloud collision. Overall, our results support the cloud-cloud collision process as being an important mechanism for star cluster formation in Orion A.

We conduct a census of the high-mass protostellar population of the $\sim70,000\:M_\odot$ Infrared Dark Cloud (IRDC) G028.37+00.07, identifying 35 sources based on their $70\:\mu$m emission, as reported in the {\it Herschel} Hi-GAL catalog of Molinari et al. (2016). We perform aperture photometry to construct spectral energy distributions (SEDs), which are then fit with the massive protostar models of Zhang & Tan (2018). We find that the sources span a range of isotropic luminosities from $\sim$20 to 4,500$\:L_\odot$. The most luminous sources are predicted to have current protostellar masses of $m_{*}\sim10\:M_\odot$ forming from cores of mass $M_{c}\sim40$ to $400\:M_\odot$. The least luminous sources in our sample are predicted to be protostars with masses as low as $\sim 0.5\:M_\odot$ forming from cores with $M_{c}\sim10\:M_\odot$, which are the minimum values explored in the protostellar model grid. The detected protostellar population has a total estimated protostellar mass of $M_{*}\sim 100\:M_\odot$. Allowing for completeness corrections, which are constrained by comparison with an ALMA study in part of the cloud, we estimate a star formation efficiency per free-fall time of $\sim3\%$ in the IRDC. Finally, analyzing the spatial distribution of the sources, we find relatively low degrees of central concentration of the protostars. The protostars, including the most massive ones, do not appear to be especially centrally concentrated in the protocluster as defined by the IRDC boundary.

We study the gas kinematics and dynamics of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array using spectral lines of $C^{18}O$(2-1), $N_2D^+$(3-2), $DCO^+$(3-2) and $DCN$(3-2). On the parsec clump scale, $C^{18}O$ emission appears highly filamentary around the systemic velocity. $N_2D^+$ and $DCO^+$ are more closely associated with the dust continuum. $DCN$ is strongly concentrated towards the protocluster center, where no or only weak detection is seen for $N_2D^+$ and $DCO^+$, possibly due to this region being at a relatively evolved evolutionary stage. Spectra of 76 continuum defined dense cores, typically a few 1000 AU in size, are analysed to measure their centroid velocities and internal velocity dispersions. There are no statistically significant velocity offsets of the cores among the different dense gas tracers. Furthermore, the majority (71\%) of the dense cores have subthermal velocity offsets with respect to their surrounding, lower density $C^{18}O$ emitting gas. Within the uncertainties, the dense cores in G286 show internal kinematics that are consistent with being in virial equilibrium. On clumps scales, the core to core velocity dispersion is also similar to that required for virial equilibrium in the protocluster potential. However, the distribution in velocity of the cores is largely composed of two spatially distinct groups, which indicates that the dense molecular gas has not yet relaxed to virial equilibrium, perhaps due to there being recent/continuous infall into the system.

We present our second set of results from our mid-infrared imaging survey of Milky Way Giant HII regions. We used the FORCAST instrument on the Stratospheric Observatory For Infrared Astronomy to obtain 20 and 37$\mu$m images of the central ~10'X10' area of M17. We investigate the small- and large-scale properties of M17 using our data in conjunction with previous multi-wavelength observations. The spectral energy distributions of individual compact sources were constructed with Spitzer-IRAC, SOFIA-FORCAST, and Herschel-PACS photometry data and fitted with massive young stellar object (MYSO) models. Seven sources were found to match the criteria for being MYSO candidates, four of which are identified here for the first time, and the stellar mass of the most massive object, UC1, is determined to be 64 solar mass. We resolve the extended mid-infrared emission from the KW Object, and suggest that the angle of this extended emission is influenced by outflow. IRS5 is shown to decrease in brightness as a function of wavelength from the mid- to far-infrared, and has several other indicators that point to it being an intermediate mass Class II object and not a MYSO. We find that the large-scale appearance of emission in M17 at 20$\mu$m is significantly affected by contamination from the [SIII] emission line from the ionized gas of the Giant HII region. Finally, a number of potential evolutionary tracers yield a consistent picture suggesting that the southern bar of M17 is likely younger than the northern bar.

Dust offers a unique probe of the interstellar medium (ISM) across multiple size, density, and temperature scales. Dust is detected in outflows of evolved stars, star-forming molecular clouds, planet-forming disks, and even in galaxies at the dawn of the Universe. These grains also have a profound effect on various astrophysical phenomena from thermal balance and extinction in galaxies to the building blocks for planets, and changes in dust grain properties will affect all of these phenomena. A full understanding of dust in all of its forms and stages requires a multi-disciplinary investigation of the dust life cycle. Such an investigation can be achieved with a statistical study of dust properties across stellar evolution, star and planet formation, and redshift. Current and future instrumentation will enable this investigation through fast and sensitive observations in dust continuum, polarization, and spectroscopy from near-infrared to millimeter wavelengths.

We discuss the first results from our mid-infrared imaging survey of Milky Way Giant HII regions with our detailed analysis of W51A, which is one of the largest GHII regions in our Galaxy. We used the FORCAST instrument on SOFIA to obtain 20 and 37$\mu$m images of the central $10' \times 20'$ area, which encompasses both of the G49.5-0.4 and G49.4-0.3 sub-regions. Based on these new data, and in conjunction with previous multi-wavelength observations, we conjecture on the physical nature of several individual sources and sub-components within W51A. We find that extinction seems to play an important role in the observed structures we see in the near- to mid-infrared, both globally and locally. We used the SOFIA photometry combined with Spitzer-IRAC and Herschel-PACS photometry data to construct spectral energy distributions (SEDs) of sub-components and point sources detected in the SOFIA images. We fit those SEDs with young stellar object models, and found 41 sources that are likely to be massive young stellar objects, many of which are identified as such in this work for the first time. Close to half of the massive young stellar objects do not have detectable radio continuum emission at cm wavelengths, implying a very young state of formation. We derived luminosity-to-mass ratio and virial parameters of the extended radio sub-regions of W51A to estimate their relative ages.

We present SOFIA-upGREAT observations of [CII] emission of Infrared Dark Cloud (IRDC) G035.39-00.33, designed to trace its atomic gas envelope and thus test models of the origins of such clouds. Several velocity components of [CII] emission are detected, tracing structures that are at a wide range of distances in the Galactic plane. We find a main component that is likely associated with the IRDC and its immediate surroundings. This strongest emission component has a velocity similar to that of the $^{13}$CO(2-1) emission of the IRDC, but offset by $\sim3\:{\rm km\:s}^{-1}$ and with a larger velocity width of $\sim9\:{\rm km\:s}^{-1}$. The spatial distribution of the [CII] emission of this component is also offset predominantly to one side of the dense filamentary structure of the IRDC. The CII column density is estimated to be of the order of $\sim10^{17}-10^{18}\,{\rm cm}^{-2}$. We compare these results to the [CII] emission from numerical simulations of magnetized, dense gas filaments formed from giant molecular cloud (GMC) collisions, finding similar spatial and kinematic offsets. These observations and modeling of [CII] add further to the evidence that IRDC G035.39-00.33 has been formed by a process of GMC-GMC collision, which may thus be an important mechanism for initiating star cluster formation.

We study the core mass function (CMF) of the massive protocluster G286.21+0.17 with the Atacama Large Millimeter/submillimeter Array via 1.3~mm continuum emission at a resolution of 1.0\arcsec\ (2500~au). We have mapped a field of 5.3\arcmin$\times$5.3\arcmin\ centered on the protocluster clump. We measure the CMF in the central region, exploring various core detection algorithms, which give source numbers ranging from 60 to 125, depending on parameter selection. We estimate completeness corrections due to imperfect flux recovery and core identification via artificial core insertion experiments. For masses $M\gtrsim1\:M_\odot$, the fiducial dendrogram-identified CMF can be fit with a power law of the form ${\rm{d}}N/{\rm{d}}{\rm{log}}M\propto{M}^{-\alpha}$ with $\alpha \simeq1.24\pm0.17$, slightly shallower than, but still consistent with, the index of the Salpeter stellar initial mass function of 1.35. Clumpfind-identified CMFs are significantly shallower with $\alpha\simeq0.64\pm0.13$. While raw CMFs show a peak near $1\:M_\odot$, completeness-corrected CMFs are consistent with a single power law extending down to $\sim 0.5\:M_\odot$, with only a tentative indication of a shallowing of the slope around $\sim1\:M_\odot$. We discuss the implications of these results for star and star cluster formation theories.

Aims. Infrared dark clouds represent the earliest stages of high-mass star formation. Detailed observations of their physical conditions on all physical scales are required to improve our understanding of their role in fueling star formation. Methods. We investigate the large-scale structure of the IRDC G035.39-00.33, probing the dense gas with the classical ammonia thermometer. This allows us to put reliable constraints on the temperature of the extended, pc-scale dense gas reservoir and to probe the magnitude of its non-thermal motions. Available far-infrared observations can be used in tandem with the observed ammonia emission to estimate the total gas mass contained in G035.39-00.33. Results. We identify a main velocity component as a prominent filament, manifested as an ammonia emission intensity ridge spanning more than 6 pc, consistent with the previous studies on the Northern part of the cloud. A number of additional line-of-sight components are found, and a large scale, linear velocity gradient of ~0.2 km s$^{-1}$ pc$^{-1}$ is found along the ridge of the IRDC. In contrast to the dust temperature map, an ammonia-derived kinetic temperature map, presented for the entirety of the cloud, reveals local temperature enhancements towards the massive protostellar cores. We show that without properly accounting for the line of sight contamination, the dust temperature is 2-3 K larger than the gas temperature measured with NH$_3$. Conclusions. While both the large scale kinematics and temperature structure are consistent with that of starless dark filaments, the kinetic gas temperature profile on smaller scales is suggestive of tracing the heating mechanism coincident with the locations of massive protostellar cores.

We study the probability distribution function (PDF) of mass surface densities, $\Sigma$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since it has almost 100,000 solar masses within a radius of 8 parsecs, making it one of the most massive, dense molecular structures known and is thus a potential site for the formation of a "super star cluster." We study $\Sigma$ in two ways. First, we use a combination of NIR and MIR extinction maps that are able to probe the bulk of the cloud structure up to $\Sigma\sim1\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq200$ mag). Second, we study the FIR and sub-mm dust continuum emission from the cloud utilizing Herschel PACS and SPIRE images and paying careful attention to the effects of foreground and background contamination. We find that the PDFs from both methods, applied over a $\sim20^\prime$(30 pc)-sized region that contains $\simeq1.5\times10^5\:M_\odot$ and encloses a minimum closed contour with $\Sigma\simeq0.013\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq3$ mag), shows a log-normal shape with the peak measured at $\Sigma\simeq0.021\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq4.7$ mag). There is tentative evidence for the presence of a high-$\Sigma$ power law tail that contains from $\sim3\%$ to 8\% of the mass of the cloud material. We discuss the implications of these results for the physical processes occurring in this cloud.

We present spectroscopic tests of MIR to FIR extinction laws in IRDC G028.36+00.07, a potential site of massive star and star cluster formation. Lim & Tan (2014) developed methods of FIR extinction mapping of this source using ${\it Spitzer}$-MIPS ${\rm 24\mu m}$ and ${\it Herschel}$-PACS ${\rm 70\mu m}$ images, and by comparing to MIR ${\it Spitzer}$-IRAC $3$--${\rm 8\mu m}$ extinction maps, found tentative evidence for grain growth in the highest mass surface density regions. Here we present results of spectroscopic infrared extinction (SIREX) mapping using ${\it Spitzer}$-IRS (14 to ${\rm 38\mu m}$) data of the same IRDC. These methods allow us to first measure the SED of the diffuse Galactic ISM that is in the foreground of the IRDC. We then carry out our primary investigation of measuring the MIR to FIR opacity law and searching for potential variations as a function of mass surface density within the IRDC. We find relatively flat, featureless MIR-FIR opacity laws that lack the $\sim{\rm 12\mu m}$ and $\sim{\rm 35\mu m}$ features associated with the thick water ice mantle models of Ossenkopf & Henning (1994). Their thin ice mantle models and the coagulating aggregate dust models of Ormel et al. (2011) are a generally better match to the observed opacity laws. We also find evidence for generally flatter MIR to FIR extinction laws as mass surface density increases, strengthening the evidence for grain and ice mantle growth in higher density regions.

We present new numerical cosmological solutions of the Einstein Field Equations. The spacetime is spherically symmetric with a source of dust and radiation approximated as a perfect fluid. The dust and radiation are necessarily non-comoving due to the inhomogeneity of the spacetime. Such a model can be used to investigate non-linear general relativistic effects present during decoupling or big-bang nucleosynthesis, as well as for investigating void models of dark energy with isocurvature degrees of freedom. We describe the full evolution of the spacetime as well as the redshift and luminosity distance for a central observer. After demonstrating accuracy of the code, we consider a few example models, and demonstrate the sensitivity of the late time model to the degree of inhomogeneity of the initial radiation contrast.

Progress in understanding star formation requires detailed observational constraints on the initial conditions, i.e. dense clumps and cores in giant molecular clouds that are on the verge of gravitational instability. Such structures have been studied by their extinction of Near-Infrared (NIR) and, more recently, Mid-Infrared (MIR) background light. It has been somewhat more of a surprise to find that there are regions that appear as dark shadows at Far-Infrared (FIR) wavelengths as long as $\sim$100$\mu m$. Here we develop analysis methods of FIR images from Spitzer-MIPS and Herschel-PACS that allow quantitative measurements of cloud mass surface density, $\Sigma$. The method builds upon that developed for MIR extinction mapping (MIREX) (Butler and Tan 2012), in particular involving a search for independent saturated, i.e. very opaque, regions that allow measurement of the foreground intensity. We focus on three massive starless core/clumps in IRDC G028.37+00.07, deriving mass surface density maps from 3.5 to 70$\mu m$. A by-product of this analysis is measurement of the spectral energy distribution of the diffuse foreground emission. The lower opacity at 70$\mu m$ allows us to probe to higher $\Sigma$ values, up to $\sim1\:{\rm{g\:cm}^{-2}}$ in the densest parts of the core/clumps. Comparison of the $\Sigma$ maps at different wavelengths constrains the shape of the MIR-FIR dust opacity law in IRDCs. We find it is most consistent with the thick ice mantle models of Ossenkopf and Henning (1994). There is tentative evidence for grain ice mantle growth as one goes from lower to higher $\Sigma$ regions.

We present an analysis of the papers published in the journals Nature and Science in the years from 2006 to 2010. During this period, a total of 7788 papers were published in the two journals. This includes 544 astronomy papers that comprise 7.0% of the papers in `all' research fields and 18.9% of those in the fields of `physical sciences'. The sub-fields of research of the astronomy papers are distributed, in descending order of number of papers, in Solar System, stellar astronomy, galaxies and the universe, the Milky Way Galaxy, and exoplanets. The observational facilities used for the studies are mainly ground-based telescopes (31.1%), spacecrafts (27.0%), and space telescopes (22.8%), while 16.0% of papers did not use any noticeable facilities and 1.7% used other facilities. Korean scientists have published 86 papers (33 in Nature and 53 in Science), which is 1.10% of all the papers (N=7788) in the two journals. The share of papers by Korean astronomers among the scientific papers by Koreans is 8.14%, slightly higher than the contribution of astronomy papers (7.0%) in both journals.

A period of slow contraction with equation of state w > 1, known as an ekpyrotic phase, has been shown to flatten and smooth the universe if it begins the phase with small perturbations. In this paper, we explore how robust and powerful the ekpyrotic smoothing mechanism is by beginning with highly inhomogeneous and anisotropic initial conditions and numerically solving for the subsequent evolution of the universe. Our studies, based on a universe with gravity plus a scalar field with a negative exponential potential, show that some regions become homogeneous and isotropic while others exhibit inhomogeneous and anisotropic behavior in which the scalar field behaves like a fluid with w=1. We find that the ekpyrotic smoothing mechanism is robust in the sense that the ratio of the proper volume of the smooth to non-smooth region grows exponentially fast along time slices of constant mean curvature.

We study tilted perfect fluid cosmological models with a constant equation of state parameter in spatially homogeneous models of Bianchi type VI$_{-1/9}$ using dynamical systems methods and numerical simulations. We study models with and without vorticity, with an emphasis on their future asymptotic evolution. We show that for models with vorticity there exists, in a small region of parameter space, a closed curve acting as the attractor.

We study tilted perfect fluid cosmological models with a constant equation of state parameter in spatially homogeneous models of Bianchi type VI_h using dynamical systems methods and numerical experimentation, with an emphasis on their future asymptotic evolution. We determine all of the equilibrium points of the type VI_h state space (which correspond to exact self-similar solutions of the Einstein equations, some of which are new), and their stability is investigated. We find that there are vacuum plane-wave solutions that act as future attractors. In the parameter space, a `loophole' is shown to exist in which there are no stable equilibrium points. We then show that a Hopf-bifurcation can occur resulting in a stable closed orbit (which we refer to as the Mussel attractor) corresponding to points both inside the loophole and points just outside the loophole; in the former case the closed curves act as late-time attractors while in the latter case these attracting curves will co-exist with attracting equilibrium points. In the special Bianchi type III case, centre manifold theory is required to determine the future attractors. Comprehensive numerical experiments are carried out to complement and confirm the analytical results presented. We note that the Bianchi type VI_h case is of particular interest in that it contains many different subcases which exhibit many of the different possible future asymptotic behaviours of Bianchi cosmological models.

We investigate the effect of dark energy on the limits on the shear anisotropy in spatially homogeneous Bianchi cosmological models obtained from measurements of the temperature anisotropies in the cosmic microwave background. We shall primarily assume that the dark energy is modelled by a cosmological constant. In general, we find that there are tighter bounds on the shear than in models with no cosmological constant, although the limits are (Bianchi) model dependent. In addition, there are special spatially homogeneous cosmological models whose rate of expansion is highly anisotropic, but whose cosmic microwave background temperature is measured to be exactly isotropic at one instant of time.

We study the properties of a future singularity encountered by a perfect fluid observer in tilting spatially homogeneous Bianchi cosmologies. We derive the boost formulae for the Weyl tensor to establish that, for two observers that are asymptotically null with respect to each other, their respective Weyl parameters generally both tend to zero, constant, or infinity together. We examine three classes of typical examples and one exceptional class. Given the behaviour of the Weyl parameter, we can predict that the singularity encountered is a Weyl singularity or a kinematic singularity. The analysis suggests that the kinematic variables are also useful in indicating a singularity in these models.

We show that in tilting perfect fluid cosmological models with an ultra-radiative equation of state, generically the tilt becomes extreme at late times and, as the tilt instability sets in, observers moving with the tilting fluid will experience singular behaviour in which infinite expansion is reached within a finite proper time, similar to that of phantom cosmology (but without the need for exotic forms of matter).

We use the dynamical systems approach to investigate the Bianchi type VIII models with a tilted $\gamma$-law perfect fluid. We introduce expansion-normalised variables and investigate the late-time asymptotic behaviour of the models and determine the late-time asymptotic states. For the Bianchi type VIII models the state space is unbounded and consequently, for all non-inflationary perfect fluids, one of the curvature variables grows without bound. Moreover, we show that for fluids stiffer than dust ($1<\gamma<2$), the fluid will in general tend towards a state of extreme tilt. For dust ($\gamma=1$), or for fluids less stiff than dust ($0<\gamma< 1$), we show that the fluid will in the future be asymptotically non-tilted. Furthermore, we show that for all $\gamma\geq 1$ the universe evolves towards a vacuum state but does so rather slowly, $\rho/H^2\propto 1/\ln t$.

We use expansion-normalised variables to investigate the Bianchi type VII$_0$ model with a tilted $\gamma$-law perfect fluid. We emphasize the late-time asymptotic dynamical behaviour of the models and determine their asymptotic states. Unlike the other Bianchi models of solvable type, the type VII$_0$ state space is unbounded. Consequently we show that, for a general non-inflationary perfect fluid, one of the curvature variables diverges at late times, which implies that the type VII$_0$ model is not asymptotically self-similar to the future. Regarding the tilt velocity, we show that for fluids with $\gamma<4/3$ (which includes the important case of dust, $\gamma=1$) the tilt velocity tends to zero at late times, while for a radiation fluid, $\gamma=4/3$, the fluid is tilted and its vorticity is dynamically significant at late times. For fluids stiffer than radiation ($\gamma>4/3$), the future asymptotic state is an extremely tilted spacetime with vorticity.

In this paper we investigate asymptotic isotropization. We derive the asymptotic dynamics of spatially inhomogeneous cosmological models with a perfect fluid matter source and a positive cosmological constant near the de Sitter equilibrium state at late times, and near the flat FL equilibrium state at early times. Our results show that there exists an open set of solutions approaching the de Sitter state at late times, consistent with the cosmic no-hair conjecture. On the other hand, solutions that approach the flat FL state at early times are special and admit a so-called isotropic initial singularity. For both classes of models the asymptotic expansion of the line element contains an arbitrary spatial metric at leading order, indicating asymptotic spatial inhomogeneity. We show, however, that in the asymptotic regimes this spatial inhomogeneity is significant only at super-horizon scales.

In this letter we will show that, contrary to what is widely believed, an almost isotropic cosmic microwave background (CMB) temperature does not imply that the universe is ``close to a Friedmann-Lemaitre universe''. There are two important manifestations of anisotropy in the geometry of the universe, (i) the anisotropy in the overall expansion, and (ii) the intrinsic anisotropy of the gravitational field, described by the Weyl curvature tensor, although the former usually receives more attention than the latter in the astrophysical literature. Here we consider a class of spatially homogeneous models for which the anisotropy of the CMB temperature is within the current observational limits but whose Weyl curvature is not negligible, i.e. these models are not close to isotropy even though the CMB temperature is almost isotropic.