The ATOMS, standing for ALMA Three-millimeter Observations of Massive Star-forming regions, survey has observed 146 active star-forming regions with ALMA band 3, aiming to systematically investigate the spatial distribution of various dense gas tracers in a large sample of Galactic massive clumps, to study the roles of stellar feedback in star formation, and to characterize filamentary structures inside massive clumps. In this work, the observations, data analysis, and example science of the ATOMS survey are presented, using a case study for the G9.62+0.19 complex. Toward this source, some transitions, commonly assumed to trace dense gas, including CS J = 2-1, HCO+ J = 1-0, and HCN J = 1-0, are found to show extended gas emission in low-density regions within the clump; less than 25 per cent of their emission is from dense cores. SO, CH3OH, (HCN)-C-13, and HC3N show similar morphologies in their spatial distributions and reveal well the dense cores. Widespread narrow SiO emission is present (over similar to 1 pc), which may be caused by slow shocks from large-scale colliding flows or HII regions. Stellar feedback from an expanding HII region has greatly reshaped the natal clump, significantly changed the spatial distribution of gas, and may also account for the sequential high-mass star formation in the G9.62+0.19 complex. The ATOMS survey data can be jointly analysed with other survey data, e.g. MALT90, Orion B, EMPIRE, ALMA IMF, and ALMAGAL, to deepen our understandings of 'dense gas' star formation scaling relations and massive protocluster formation.

The cryosphere of the Earth overlaps with the atmosphere, hydrosphere and lithosphere over vast areas with temperatures below 0 degrees C and pronounced H2O phase changes. In spite of its strong variability in space and time, the cryosphere plays the role of a global thermostat, keeping the thermal regime on the Earth within rather narrow limits, affording continuation of the conditions needed for the maintenance of life. Objects and processes related to cryosphere are very diverse, due to the following basic reasons: the anomalous thermodynamic and electromagnetic properties of H2O, the intermediate intensity of hydrogen bonds and the wide spread of cryogenic systems all over the Earth. However, these features attract insufficient attention from research communities. Cryology is usually understood as a descriptive discipline within physical geography, limited to glaciology and permafrost research. We emphasise its broad interdisciplinary landscape involving physical, chemical and biological phenomena related to the H2O phase transitions and various forms of ice. This paper aims to draw the attention of readers to the crucial importance of cryogenic anomalies, which make the Earth atmosphere and the entire Earth system very special, if not unique, objects in the universe.

We study interplanetary coronal mass ejections (ICMEs) measured by probes at different heliocentric distances (0.3-1 AU) to investigate the propagation of ICMEs in the inner heliosphere and determine how the generic features of ICMEs change with heliospheric distance. Using data from the MErcury Surface, Space ENvironment, GEochemistry, and Ranging (MESSENGER), Venus Express and ACE spacecraft, we analyze with the superposed epoch technique the profiles of ICME substructures, namely, the sheath and the magnetic ejecta. We determine that the median magnetic field magnitude in the sheath correlates well with ICME speeds at 1 AU, and we use this proxy to order the ICMEs at all spacecraft. We then investigate the typical ICME profiles for three categories equivalent to slow, intermediate, and fast ICMEs. Contrary to fast ICMEs, slow ICMEs have a weaker solar wind field at the front and a more symmetric magnetic field profile. We find the asymmetry to be less pronounced at Earth than at Mercury, indicating a relaxation taking place as ICMEs propagate. We also find that the magnetic field intensities in the wake region of the ICMEs do not go back to the pre-ICME solar wind intensities, suggesting that the effects of ICMEs on the ambient solar wind last longer than the duration of the transient event. Such results provide an indication of physical processes that need to be reproduced by numerical simulations of ICME propagation. The samples studied here will be greatly improved by future missions dedicated to the exploration of the inner heliosphere, such as Parker Solar Probe and Solar Orbiter.

Planar magnetic structures (PMSs) are periods in the solar wind during which interplanetary magnetic field vectors are nearly parallel to a single plane. One of the specific regions where PMSs have been reported are coronal mass ejection (CME)-driven sheaths. We use here an automated method to identify PMSs in 95 CME sheath regions observed in situ by the Wind and ACE spacecraft between 1997 and 2015. The occurrence and location of the PMSs are related to various shock, sheath, and CME properties. We find that PMSs are ubiquitous in CME sheaths; 85% of the studied sheath regions had PMSs with the mean duration of 6 h. In about one-third of the cases the magnetic field vectors followed a single PMS plane that covered a significant part (at least 67%) of the sheath region. Our analysis gives strong support for two suggested PMS formation mechanisms: the amplification and alignment of solar wind discontinuities near the CME-driven shock and the draping of the magnetic field lines around the CME ejecta. For example, we found that the shock and PMS plane normals generally coincided for the events where the PMSs occurred near the shock (68% of the PMS plane normals near the shock were separated by less than 20 degrees from the shock normal), while deviations were clearly larger when PMSs occurred close to the ejecta leading edge. In addition, PMSs near the shock were generally associated with lower upstream plasma beta than the cases where PMSs occurred near the leading edge of the CME. We also demonstrate that the planar parts of the sheath contain a higher amount of strong southward magnetic field than the non-planar parts, suggesting that planar sheaths are more likely to drive magnetospheric activity.

Observations of the submillimetre emission from Galactic dust, in both total intensity I and polarization, have received tremendous interest thanks to the Planck full-sky maps. In this paper we make use of such full-sky maps of dust polarized emission produced from the third public release of Planck data. As the basis for expanding on astrophysical studies of the polarized thermal emission from Galactic dust, we present full-sky maps of the dust polarization fraction p, polarization angle psi, and dispersion function of polarization angles ?. The joint distribution (one-point statistics) of p and N-H confirms that the mean and maximum polarization fractions decrease with increasing N-H. The uncertainty on the maximum observed polarization fraction, (max) = 22.0(-1.4)(+3.5) p max = 22 . 0 - 1.4 + 3.5 % at 353 GHz and 80 ' resolution, is dominated by the uncertainty on the Galactic emission zero level in total intensity, in particular towards diffuse lines of sight at high Galactic latitudes. Furthermore, the inverse behaviour between p and ? found earlier is seen to be present at high latitudes. This follows the ?proportional to p(-1) relationship expected from models of the polarized sky (including numerical simulations of magnetohydrodynamical turbulence) that include effects from only the topology of the turbulent magnetic field, but otherwise have uniform alignment and dust properties. Thus, the statistical properties of p, psi, and ? for the most part reflect the structure of the Galactic magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map ?xp, looking for residual trends. While the polarization fraction p decreases by a factor of 3-4 between N-H=10(20) cm(-2) and N-H=2x10(22) cm(-2), out of the Galactic plane, this product ?xp only decreases by about 25%. Because ? is independent of the grain alignment efficiency, this demonstrates that the systematic decrease in p with N-H is determined mostly by the magnetic-field structure and not by a drop in grain alignment. This systematic trend is observed both in the diffuse interstellar medium (ISM) and in molecular clouds of the Gould Belt. Second, we look for a dependence of polarization properties on the dust temperature, as we would expect from the radiative alignment torque (RAT) theory. We find no systematic trend of ?xp with the dust temperature T-d, whether in the diffuse ISM or in the molecular clouds of the Gould Belt. In the diffuse ISM, lines of sight with high polarization fraction p and low polarization angle dispersion ? tend, on the contrary, to have colder dust than lines of sight with low p and high ?. We also compare the Planck thermal dust polarization with starlight polarization data in the visible at high Galactic latitudes. The agreement in polarization angles is remarkable, and is consistent with what we expect from the noise and the observed dispersion of polarization angles in the visible on the scale of the Planck beam. The two polarization emission-to-extinction ratios, R-P/p and R-S/V, which primarily characterize dust optical properties, have only a weak dependence on the column density, and converge towards the values previously determined for translucent lines of sight. We also determine an upper limit for the polarization fraction in extinction, p(V)/E(B-V), of 13% at high Galactic latitude, compatible with the polarization fraction p approximate to 20% observed at 353 GHz. Taken together, these results provide strong constraints for models of Galactic dust in diffuse gas.

Light scattering by particles large compared to the wavelength of incident light is traditionally solved using ray optics which considers absorption inside the particle approximately, along the ray paths. To study the effects rising from this simplification, we have updated the ray-optics code SIRIS to take into account the propagation of light as inhomogeneous plane waves inside an absorbing particle. We investigate the impact of this correction on traditional ray-optics computations in the example case of light scattering by ice crystals through the extended near-infrared (NIR) wavelength regime. In this spectral range, ice changes from nearly transparent to opaque, and therefore provides an interesting test case with direct connection and applicability to atmospheric remote-sensing measurements at NIR wavelengths. We find that the correction for inhomogeneous waves systematically increases the single-scattering albedo throughout the NIR spectrum for both randomly-oriented, column-like hexagonal crystals and ice crystals shaped like Gaussian random spheres. The largest increase in the single-scattering albedo is 0.042 for hexagonal crystals and 0.044 for Gaussian random spheres, both at λ=2.725 µm. Although the effects on the 4  ×  4 scattering-matrix elements are generally small, the largest differences are seen at 2.0 µm and 3.969 µm wavelengths where the correction for inhomogeneous waves affects mostly the backscattering hemisphere of the depolarization-connected P22/P11, P33/P11, and P44/P11. We evaluated the correction for inhomogeneous waves through comparisons against the discrete exterior calculus (DEC) method. We computed scattering by hexagonal ice crystals using the DEC, a traditional ray-optics code (SIRIS3), and a ray-optics code with inhomogeneous waves (SIRIS4). Comparisons of the scattering-matrix elements from SIRIS3 and SIRIS4 against those from the DEC suggest that consideration of the inhomogeneous waves brings the ray-optics solution generally closer to the exact result and, therefore, should be taken into account in scattering by absorbing particles large compared to the wavelength of incident light.

Aims. We study whether or not rotational excitation can make a large difference to chemical models of the abundances of the H-3(+) isotopologs, including spin states, in physical conditions corresponding to starless cores and protostellar envelopes. Methods. We developed a new rate coefficient set for the chemistry of the H-3(+) isotopologs, allowing for rotational excitation, using previously published state-to-state rate coefficients. These new so-called species-to-species rate coefficients are compared with previously-used ground-state-to-species rate coefficients by calculating chemical evolution in variable physical conditions using a pseudo-time-dependent chemical code. Results. We find that the new species-to-species model produces different results to the ground state-to-species model at high density and toward increasing temperatures (T > 10 K). The most prominent difference is that the species-to-species model predicts a lower H-3(+) deuteration degree at high density owing to an increase of the rate coefficients of endothermic reactions that tend to decrease deuteration. For example at 20 K, the ground-state-to-species model overestimates the abundance of H2D+ by a factor of about two, while the abundance of D-3(+) can differ by up to an order of magnitude between the models. The spin-state abundance ratios of the various H-3(+) isotopologs are also a ffected, and the new model better reproduces recent observations of the abundances of ortho and para H2D+ and D2H+. The main caveat is that the applicability regime of the new rate coefficients depends on the critical densities of the various rotational transitions which vary with the abundances of the species and the temperature in dense clouds. Conclusions. The difference in the abundances of the H-3(+) isotopologs predicted by the species-to-species and ground state-to-species models is negligible at 10K corresponding to physical conditions in starless cores, but inclusion of the excited states is very important in studies of deuteration at higher temperatures, for example in protostellar envelopes. The species-to-species rate coefficients provide a more realistic approach to the chemistry of the H-3(+) isotopologs than the ground-state-to-species rate coefficients do, and so the former should be adopted in chemical models describing the chemistry of the H-3(+)+H-2 reacting system.

We present a hydrodynamical simulation at sub-parsec and few-solar-mass resolution of a merger between two gas-rich dwarf galaxies. Our simulation includes a detailed model for the multi-phase interstellar medium and is able to follow the entire formation history of spatially resolved star clusters, including feedback from individual massive stars. Shortly after the merger we find a population of similar to 900 stellar clusters with masses above 10(2.5) M-circle dot and a cluster mass function (CMF), which is well fitted with a power law with a slope of alpha = -1.70 +/- 0.08. We describe here in detail the formation of the three most massive clusters (M-* greater than or similar to 10(5) M-circle dot), which populate the high-mass end of the CMF. The simulated clusters form rapidly on a timescale of 6-8 Myr in converging flows of dense gas. The embedded merger phase has extremely high star formation rate surface densities of Sigma(SFR) > 10 M-circle dot yr(-1) kpc(-2) and thermal gas pressures in excess of Pth similar to 10(7) K-B cm(-3))(-1). The formation process is terminated by rapid gas expulsion driven by the first generation of supernovae, after which the cluster centers relax and both their structure and kinematics become indistinguishable from observed local globular clusters (GCs). The simulation presented here provides a general model for the formation of metal-poor GCs in chemically unevolved starbursting environments of low-mass dwarf galaxies, which are common at high redshifts.