A consensus view on the formation of planetesimals is now exposed to a threat, since recent numerical studies on the mechanical properties of dust aggregates tend to dispute the conceptual picture that submicrometer-sized grains conglomerate into planetesimals in protoplanetary disks. With the advent of precise laboratory experiments and extensive computer simulations on the interaction between elastic spheres comprising dust aggregates, we revisit a model for the tensile strength of dust aggregates consisting of small elastic grains. In the framework of contact mechanics and fracture mechanics, we examine outcomes of computer simulations and laboratory experiments on the tensile strength of dust aggregates. We provide a novel analytical formula that explicitly incorporates the volume effect on the tensile strength, namely, the dependence of tensile strength on the volume of dust aggregates. We find that our model for the tensile strength of dust aggregates well reproduces results of computer simulations and laboratory experiments, if appropriate values are adopted for the elastic parameters used in the model. Moreover, the model with dust aggregates of submicrometer-sized grains is in good harmony with the tensile strength of cometary dust and meteoroids derived from astronomical observations. Therefore, we reaffirm the commonly believed idea that the formation of planetesimals begins with conglomeration of submicrometer-sized grains condensed in protoplanetary disks.

We present optical observations of type Ia supernova (SN) 2019ein, starting at 2 days after the estimated explosion date. The spectra and the light curves show that SN 2019ein belongs to the High-Velocity (HV) and Bload Line groups with relatively rapid decline in the light curves (Delta m15(B) = 1.36 +- 0.02 mag) and the short rise time (15.37 +- 0.55 days). The Si II 6355 velocity, associated with a photospheric component but not with a detached high-velocity feature, reached ~ 20,000 km s-1 at 12 days before the B-band maximum. The line velocity however decreased very rapidly and smoothly toward the maximum light, where it was ~ 13,000 km s-1 as relatively low among HV SNe. This indicates that the speed of the spectral evolution of HV SNe Ia is correlated not only to the velocity at the maximum light, but also to the light curve decline rate like the case for Normal-Velocity (NV) SNe Ia. Spectral synthesis modeling shows that the outermost layer at > 17,000 km s-1 is well described by the O-Ne-C burning layer extending to at least 25,000 km s-1, and there is no unburnt carbon below 30,000 km s-1; these properties are largely consistent with the delayed detonation scenario, and are shared with the prototypical HV SN 2002bo despite the large difference in Delta m15(B). This structure is strikingly different from that derived for the well-studied NV SN 2011fe. We suggest that the relation between the mass of 56Ni (or Delta m15) and the extent of the O-Ne-C burning layer provides an important constraint on the explosion mechanism(s) of HV and NV SNe.

The Cometary Secondary Ion Mass Analyzer (COSIMA) onboard ESA's Rosetta orbiter has revealed that dust particles in the coma of Comet 67P/Churyumov-Gerasimenko are aggregates of small grains. We study the morphological, elastic, and electric properties of dust aggregates in the coma of Comet 67P/Churyumov-Gerasimenko using optical microscopic images taken by the COSIMA instrument. Dust aggregates in COSIMA images are well represented as fractals in harmony with morphological data from MIDAS (Micro-Imaging Dust Analysis System) and GIADA (Grain Impact Analyzer and Dust Accumulator) onboard Rosetta. COSIMA's images, together with the data from the other Rosetta's instruments such as MIDAS and GIADA do not contradict the so-called rainout growth of $10~\mu\mathrm{m}$-sized particles in the solar nebula. The elastic and electric properties of dust aggregates measured by COSIMA suggest that the surface chemistry of cometary dust is well represented as carbonaceous matter rather than silicates or ices, consistent with the mass spectra, and that organic matter is to some extent carbonized by solar radiation, as inferred from optical and infrared observations of various comets. Electrostatic lofting of cometary dust by intense electric fields at the terminator of its parent comet is unlikely, unless the surface chemistry of the dust changes from a dielectric to a conductor. Our findings are not in conflict with our current understanding of comet formation and evolution, which begin with the accumulation of condensates in the solar nebula and follow with the formation of a dust mantle in the inner solar system.

The JAXA/ISAS spacecraft DESTINY$^+$ will be launched to the active asteroid (3200) Phaethon in 2022. Among the proposed core payload is the DESTINY+ Dust Analyzer (DDA) which is an upgrade of the Cosmic Dust Analyzer flown on the Cassini spacecraft to Saturn (Srama et al. 2011). We use two up-to-date computer models, the ESA Interplanetary Meteoroid Engineering Model (IMEM, Dikarev et al. 2005), and the interstellar dust module of the Interplanetary Meteoroid environment for EXploration model (IMEX; Sterken2013 et al., Strub et al. 2019) to study the detection conditions and fluences of interplanetary and interstellar dust with DDA. Our results show that a statistically significant number of interplanetary and interstellar dust particles will be detectable with DDA during the 4-years interplanetary cruise of DESTINY+. The particle impact direction and speed can be used to descriminate between interstellar and interplanetary particles and likely also to distinguish between cometary and asteroidal particles.

The recent discovery of hot dust grains in the vicinity of main-sequence stars has become a hot issue among the scientific community of debris disks. Hot grains must have been enormously accumulated near their sublimation zones, but it is a mystery how such a high concentration of hot grains is sustained. The most difficult conundrum is that the size of hot dust grains is estimated to lie in the submicrometer range, while submicrometer-sized grains are instantly swept away from near-stellar environments by stellar radiation pressure. One and only mechanism proposed for prolonging the residence time of hot grains in the near-stellar environments is trapping of charged nanoparticles by stellar magnetic fields. We revisit the model of magnetic grain trapping around main-sequence stars of various spectral classes by taking into account sublimation and electric charging of the grains. The model of magnetic grain trapping predicts that hot dust grains are present in the vicinity of main-sequence stars with high rotation velocities and intermediate magnetic-field strengths. On the contrary, we find that the detection of hot dust grains has no correlation with the rotation velocities of central stars nor the magnetic field strengths of the stars. Our numerical evaluation of electric grain charging indicates that the surface potential of submicrometer-sized grains in the vicinity of main-sequence stars is typically 4-5 V, which is one order of magnitude smaller than the value assumed by the model of magnetic grain trapping. On the basis of our numerical simulation on sublimation of dust grains in the vicinity of a star, it turns out that their lives end due to sublimation in a timescale much shorter than the period of one revolution at the gyroradius. It is, therefore, infeasible to dynamically extend the dwell time of hot grains inside the sublimation zone by magnetic trapping, ...

This is an editorial to the special issue on Cosmic Dust IX.

Recent space missions have provided information on the physical and chemical properties of interstellar grains such as the ratio $\beta$ of radiation pressure to gravity acting on the grains in addition to the composition, structure, and size distribution of the grains. Numerical simulation on the trajectories of interstellar grains captured by Stardust and returned to Earth constrained the $\beta$ ratio for the Stardust samples of interstellar origin. However, recent accurate calculations of radiation pressure cross sections for model dust grains have given conflicting stories in the $\beta$ ratio of interstellar grains. The $\beta$ ratio for model dust grains of so-called "astronomical silicate" in the femto-kilogram range lies below unity, in conflict with $\beta \sim 1$ for the Stardust interstellar grains. Here, I tackle this conundrum by re-evaluating the $\beta$ ratio of interstellar grains on the assumption that the grains are aggregated particles grown by coagulation and composed of amorphous MgSiO$_{3}$ with the inclusion of metallic iron. My model is entirely consistent with the depletion and the correlation of major rock-forming elements in the Local Interstellar Cloud surrounding the Sun and the mineralogical identification of interstellar grains in the Stardust and Cassini missions. I find that my model dust particles fulfill the constraints on the $\beta$ ratio derived from not only the Stardust mission but also the Ulysses and Cassini missions. My results suggest that iron is not incorporated into silicates but exists as metal, contrary to the majority of interstellar dust models available to date.

This is an editorial to the special issue on Cosmic Dust VIII.

Photoelectron emission is crucial to electric charging of dust particles around main-sequence stars and gas heating in various dusty environments. An estimate of the photoelectric processes contains an ill-defined parameter called the photoelectric quantum yield, which is the total number of electrons ejected from a dust particle per absorbed photon. Here we revisit the so-called small particle effect of photoelectron emission and provide an analytical model to estimate photoelectric quantum yields of small dust particles in sizes down to nanometers. We show that the small particle effect elevates the photoelectric quantum yields of nanoparticles up to by a factor of $10^3$ for carbon, water ice, and organics, and a factor of $10^2$ for silicate, silicon carbide, and iron. We conclude the surface curvature of the particles is a quantity of great importance to the small particle effect, unless the particles are submicrometers in radius or larger.

This is an editorial to the special issue on Cosmic Dust VII.

This review focuses on numerical approaches to deducing the light-scattering and thermal-emission properties of primitive dust particles in planetary systems from astronomical observations. The particles are agglomerates of small grains with sizes comparable to visible wavelength and compositions being mainly magnesium-rich silicates, iron-bearing metals, and organic refractory materials in pristine phases. These unique characteristics of primitive dust particles reflect their formation and evolution around main-sequence stars of essentially solar composition. The development of light-scattering theories has been offering powerful tools to make a thorough investigation of light scattering and thermal emission by primitive dust agglomerates in such a circumstellar environment. In particular, the discrete dipole approximation, the T-matrix method, and effective medium approximations are the most popular techniques for practical use in astronomy. Numerical simulations of light scattering and thermal emission by dust agglomerates of submicrometer-sized constituent grains have a great potential to provide new state-of-the-art knowledge of primitive dust particles in planetary systems. What is essential to this end is to combine the simulations with comprehensive collections of relevant results from not only astronomical observations, but also in-situ data analyses, laboratory sample analyses, laboratory analogue experiments, and theoretical studies on the origin and evolution of the particles.

Adhesion forces between submicrometer-sized silicate grains play a crucial role in the formation of silicate dust agglomerates, rocky planetesimals, and terrestrial planets. The surface energy of silicate dust particles is the key to their adhesion and rolling forces in a theoretical model based on the contact mechanics. Here we revisit the cohesion of amorphous silica spheres by compiling available data on the surface energy for hydrophilic amorphous silica in various circumstances. It turned out that the surface energy for hydrophilic amorphous silica in a vacuum is a factor of 10 higher than previously assumed. Therefore, the previous theoretical models underestimated the critical velocity for the sticking of amorphous silica spheres, as well as the rolling friction forces between them. With the most plausible value of the surface energy for amorphous silica spheres, theoretical models based on the contact mechanics are in harmony with laboratory experiments. Consequently, we conclude that silicate grains with a radius of $0.1~\mu$m could grow to planetesimals via coagulation in a protoplanetary disk. We argue that the coagulation growth of silicate grains in a molecular cloud is advanced either by organic mantles rather than icy mantles or, if there are no mantles, by nanometer-sized grain radius.

We carried out the photometric observations of the SU UMa-type dwarf nova ER UMa during 2011 and 2012, which showed the existence of persistent negative superhumps even during the superoutburst. We performed two-dimensional period analysis of its light curves by using a method called "least absolute shrinkage and selection operator" (Lasso) and "phase dispersion minimization" (PDM) analysis, and we found that the period of negative superhumps systematically changed between a superoutburst and the next superoutburst. The trend of the period change can beinterpreted as reflecting the change of the disk radius. This change of the disk radius is in good agreement with the predicted change of the disk radius by the thermal-tidal instability (TTI) model. The normal outbursts within a supercycle showed a general trend that the rising rate to maximum becomes slower as the next superoutburst approaches. The change can be interpreted as the consequence of the increased gas-stream flow onto the inner region of the disk as the result of the tilted disk. Some of the superoutbursts were found to be triggered by a precursor normal outburst when the positive superhumps appeared to develop. The positive and negative superhumps co-existed during the superoutburst. The positive superhumps were prominent only during four or five days after the supermaximum, while the signal of the negative superhumps became strong after the middle phase of the superoutburst plateau. A simple combination of the positive and negative superhumps was found to be insufficient in reproducing the complex profile variation. We were able to detect the developing phase of positive superhumps (stage A superhumps) for the first time in ER UMa-type dwarf novae. Using the period of stage A superhumps, we obtained a mass ratio of 0.100(15), which indicates that ER UMa is on the ordinary evolutional track of CVs.

The trajectories of dust particles ejected from a comet are affected by solar radiation pressure as a function of their ratios of radiation pressure cross section to mass. Therefore, a study on the orbital evolution of the particles caused by the radiation pressure reveals the physical properties of dust on the surface of the comet nucleus. In the course of NASA's Deep Impact mission, the ejecta plume evolved under the influence of the radiation pressure. From the evolution and shape of the plume, we have succeeded in obtaining $\beta \approx 0.4$, where $\beta$ is the ratio of the radiation pressure to the solar gravity. Taking into account $\beta \approx 0.4$ as well as the observational constraints of a high color temperature and a small silicate-feature strength, dust particles ejected from the surface of comet 9P/Tempel 1 are likely compact dust aggregates of sizes $\approx 20\,\mu$m (mass $\sim 10^{-8}$\,g). This is comparable to the major dust on the surface of comet 1P/Halley ($\sim 10\mu$m) inferred from in-situ measurements and theoretical considerations. Since such dust aggregates with $\beta \approx 0.4$ must have survived on the surface against jets due to ice sublimation, the temperature of ice in the nucleus must be kept below 145\,K, which is much lower than equilibrium temperature determined by solar irradiation and thermal emission. These facts indicate that 9P/Tempel 1 has a dust mantle composed of $20\,\mu$m-sized dust aggregates with low thermal conductivities $\sim 1 \, {\rm erg\, cm}^{-1} \, {\rm K}^{-1}\,{\rm s}^{-1}$.

Dust particles in orbit around a star drift toward the central star by the Poynting-Robertson effect and pile up by sublimation. We analytically derive the pile-up magnitude, adopting a simple model for optical cross sections. As a result, we find that the sublimation temperature of drifting dust particles plays the most important role in the pile-up rather than their optical property does. Dust particles with high sublimation temperature form a significant dust ring, which could be found in the vicinity of the sun through in-situ spacecraft measurements. While the existence of such a ring in a debris disk could not be identified in the spectral energy distribution (SED), the size of a dust-free zone shapes the SED. Since we analytically obtain the location and temperature of sublimation, these analytical formulae are useful to find such sublimation evidences.

We construct a theoretical model for low-temperature crystallization of amorphous silicate grains induced by exothermic chemical reactions. As a first step, the model is applied to the annealing experiments, in which the samples are (1) amorphous silicate grains and (2) amorphous silicate grains covered with an amorphous carbon layer. We derive the activation energies of crystallization for amorphous silicate and amorphous carbon from the analysis of the experiments. Furthermore, we apply the model to the experiment of low-temperature crystallization of amorphous silicate core covered with an amorphous carbon layer containing reactive molecules. We clarify the conditions of low-temperature crystallization due to exothermic chemical reactions. Next, we formulate the crystallization conditions so as to be applicable to astrophysical environments. We show that the present crystallization mechanism is characterized by two quantities: the stored energy density Q in a grain and the duration of the chemical reactions \tau . The crystallization conditions are given by Q > Q_{min} and \tau < \tau _{cool} regardless of details of the reactions and grain structure, where \tau _{cool} is the cooling timescale of the grains heated by exothermic reactions, and Q_{min} is minimum stored energy density determined by the activation energy of crystallization. Our results suggest that silicate crystallization occurs in wider astrophysical conditions than hitherto considered.

The most successful model of comet dust presents comet particles as aggregates of submicron grains. It qualitatively explains the spectral and angular change in the comet brightness and polarization and is consistent with the thermal infrared data and composition of the comet dust obtained {\it in situ} for comet 1P/Halley. However, it experiences some difficulties in providing a quantitative fit to the observational data. Here we present a model that considers comet dust as a mixture of aggregates and compact particles. The model is based on the Giotto and Stardust mission findings that both aggregates (made mainly of organics, silicates, and carbon) and solid silicate particles are present in the comet dust. We simulate aggregates as {\bf Ballistic Cluster-Cluster Aggregates (BCCA)} and compact particles as polydisperse spheroids with some distribution of the aspect ratio. The particles follow a power-law size distribution with the power -3 that is close to the one obtained for comet dust {\it in situ}, at studies of the Stardust returned samples, and the results of ground-based observations of comets. The model provides a good fit to the angular polarization curve. It also reproduces the positive spectral gradient of polarization, red color of the dust, and {\bf low albedo. It also has the ratio of compact to fluffy particles close to the one found {\it in situ} for comet 1P/Halley} and the mass ratio of silicate to carbonaceous materials equal to unity that is in accordance with the elemental abundances of Halley's dust found by Giotto mission.

According to our common understandings, the original surface of a short-period comet nucleus has been lost by sublimation processes during its close approaches to the Sun. Sublimation results in the formation of a dust mantle on the retreated surface and in chemical differentiation of ices over tens or hundreds of meters below the mantle. In the course of NASA's Deep Impact mission, optical and infrared imaging observations of the ejecta plume were conducted by several researchers, but their interpretations of the data came as a big surprise: (1) The nucleus of comet 9P/Tempel 1 is free of a dust mantle, but maintains its pristine crust of submicron-sized carbonaceous grains; (2) Primordial materials are accessible already at a depth of several tens of cm with abundant silicate grains of submicrometer sizes. In this study, we demonstrate that a standard model of cometary nuclei explains well available observational data: (1) A dust mantle with a thickness of ~1-2 m builds up on the surface, where compact aggregates larger than tens of micrometers dominate; (2) Large fluffy aggregates are embedded in chemically differentiated layers as well as in the deepest part of the nucleus with primordial materials. We conclude that the Deep Impact results do not need any peculiar view of a comet nucleus.

Comets can be divided into two groups: type I, characterized by high gas/dust ratio, low polarization, and a weak or absent 10 micron silicate feature, and type II, for which a low gas/dust ratio, high polarization, and strong silicate feature are typical. We show that the low polarization is the apparent result of depolarization by gas contamination at low dust concentration, which, in turn, results from the dust in type I comets being concentrated near the nucleus. The simulations of thermal emission show that for more porous particles (BCCA), the silicate feature is more pronounced than more compact ones (BPCA), for which it even vanishes as the particles become larger. We also show that in both types of comets the main contribution to light scattering and emission comes from particles larger than 10 micron. Conclusions: .The strength of the silicate feature in the cometary infrared spectra suggests that the dust in type II comets consists of high-porosity aggregates, whereas the dust of type I comets contains low-porosity ones. This is consistent with the polarimetric features of these comets, which indicate that the dust in type I comets tends to concentrate near the nucleus. This may result from the predominance of highly processed particles in type I comets, whereas in type II comets we see pristine or slightly-processed dust. This conclusion is in accordance with the orbital characteristics of the comets. We have found that the strength of the silicate feature correlates with the semi-major axis of periodic comets and, for short-period comets, with the perihelion distance. Thus, the silicate feature weakens due to compaction of aggregate particles if a comet spends more time in the vicinity of the Sun, which allows the comet to accumulate a mantle on the surface of its nucleus.

Wide-field (~8' x 8') and deep near-infrared (JHKs bands) polarization images of the Orion nebulae (IRNe) around young stellar objects (YSOs), both massive and low-mass. We found the IRNe around both IRc2 and BN to be very extensive, suggesting that there might be two extended (>0.7 pc) bipolar/monopolar IRNe in these sources. We discovered at least 13 smaller-scale (~0.01-0.1 pc) IRNe around less-massive YSOs including the famous source theta^2 Ori C. We also suggest the presence of many unresolved (<690 AU) systems around low-mass YSOs and young brown dwarfs showing possible intrinsic polarizations. Wide-field infrared polarimetry is thus demonstrated to be a powerful technique in revealing IRNe and hence potential disk/outflow systems among high-mass to substellar YSOs.

High-resolution K band imaging polarimetry of the beta Pic dust disk has been conducted with adaptive optics and a coronagraph using the Subaru 8.2-m telescope. Polarization of ~10 % is detected out to r ~ 120 AU with a centro-symmetric vector pattern around the central star, confirming that the disk is seen as an infrared reflection nebula. We have modeled our near-infrared and previous optical polarization results in terms of dust scattering in the disk and have found that both the degrees of polarization and the radial intensity profiles are well reproduced. We argue that the observed characteristics of the disk dust are consistent with the presence of ice-filled fluffy aggregates consisting of submicron grains in the beta Pic system. There is a gap around 100 AU in both the intensity and polarization profiles, which suggests a paucity of planetesimals in this region. The radial intensity profile also shows ripple-like structures, which are indicative of the presence of multiple planetesimal belts, as in the case of the M-type Vega-like star AU Mic.

An extensive photometric monitoring of KH 15D, an enigmatic variable in the young star cluster NGC 2264, has been conducted. Simultaneous and accurate near-infrared (JHKs-bands) photometry is presented between 2003 December and 2005 March covering most of the variable phase. The infrared variability is characterized by large-amplitude and long-lasting eclipse, as observed at optical. The period of variability is 48.3 +/- 0.2 days, the maximum photometric amplitude of variability is ~4.2 mag, and the eclipse duration is \~0.5 in phase units. These are consistent with the most recent period, amplitude, and duration at optical. The blueing of the J-H color (~0.16 mag) during the eclipse, which has been suggested before, is unambiguously confirmed; a similar blueing at H-Ks is less clear but is probably present at a similar level. The overall shape of the JHKs light curves is very similar to the optical one, including a fair time-symmetry and a less stable flux during the eclipse with a slight hump near the zero phase. Most of these variability features of KH 15D observed at near-infrared wavelengths can be explained with the recent model employing an eclipse by the inclined, precessing disk and an outer scattering region around a pre-main-sequence binary.