The magnetic field structure of a star-forming Bok globule BHR 71 was determined based on near-infrared polarimetric observations of background stars. The magnetic field in BHR 71 was mapped from 25 stars. By using a simple 2D parabolic function, the plane-of-sky magnetic axis of the core was found to be $\theta_{\rm mag} = 125^{\circ} \pm 11^{\circ}$. The plane-of-sky mean magnetic field strength of BHR 71 was found to be $B_{\rm pos} = 8.8 - 15.0$ $\mu$G, indicating that the BHR 71 core is magnetically supercritical with $\lambda = 1.44 - 2.43$. Taking into account the effect of thermal/turbulent pressure and the plane-of-sky magnetic field component, the critical mass of BHR 71 was $M_{\rm cr} = 14.5-18.7$ M$_{\odot}$, which is consistent with the observed core mass of $M_{\rm core} \approx 14.7$ M$_{\odot}$ (Yang et al. 2017). We conclude that BHR 71 is in a condition close to a kinematically critical state, and the magnetic field direction lies close to the plane of sky. Since BHR 71 is a star-forming core, a significantly subcritical condition (i.e., the magnetic field direction deviating from the plane of sky) is unlikely, and collapsed from a condition close to a kinematically critical state. There are two possible scenarios to explain the curved magnetic fields of BHR 71, one is an hourglass-like field structure due to mass accumulation and the other is the Inoue \& Fukui (2013) mechanism, which proposes the interaction of the core with a shock wave to create curved magnetic fields wrapping around the core.

Observed extinction curves of active galactic nuclei (AGNs) are significantly different from those observed in the Milky Way. The observations require preferential removal of small grains at the AGN environment; however, the physics for this remains unclear. In this paper, we propose that dust destruction by charging, or Coulomb explosion, may be responsible for AGN extinction curves. Harsh AGN radiation makes a dust grain highly charged through photoelectric emission, and grain fission via the Coulomb explosion occurs when the electrostatic tensile stress of a charge grain exceeds its tensile strength. We show that the Coulomb explosion can preferentially remove both small silicate and graphite grains and successfully reproduce both flat extinction curves and the absence of 2175\AA~bump.

Recent mid-infrared high spatial resolution observations have revealed that active galactic nuclei (AGNs) may host a polar dust region with the size of several pc, and such dust may be carried by radiation from the central engine. The polar dust emission often exhibits very weak or absence of the silicate 10-$\mu$m emission feature. A possible explanation is that the polar dust is dominated by micron-sized large grains because these grains do not show the silicate feature, while it remains unclear how large grains are preferentially supplied to the polar region. We here propose a new scenario describing the prevalence of large grains at the polar region. We show that grains are accelerated to the hypersonic drift velocity by the radiation pressure from AGN, and the hypersonic drift results in dust destruction via kinetic sputtering. Sputtering destroys small grains faster than the large ones, and thus larger grains will be preferentially blown over longer distance. Although the hypersonic drift, or kinetic sputtering, tends to be suppressed for very small grains due to the Coulomb drag, they might also be disrupted by Coulomb explosion. Removal of small grains and/or survival of large grains may explain the lack of a silicate 10-$\mu$m emission feature in polar dust emission.

In this study, the detailed magnetic field structure of the dense protostellar core Barnard 335 (B335) was revealed based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. Magnetic fields pervading B335 were mapped using 24 stars after subtracting unrelated ambient polarization components, for the first time revealing that they have an axisymmetrically distorted hourglass-shaped structure toward the protostellar core. On the basis of simple two- and three-dimensional magnetic field modeling, magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $90^{\circ} \pm 7^{\circ}$ and $50^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B335 was determined to be $30.2 \pm 17.7$ $\mu {\rm G}$. The critical mass of B335, evaluated using both magnetic and thermal/turbulent support against collapse, was determined to be $M_{\rm cr} = 3.37 \pm 0.94$ ${\rm M}_{\odot}$, which is identical to the observed core mass of $M_{\rm core}=3.67$ M$_{\odot}$. We thus concluded that B335 started its contraction from a condition near equilibrium. We found a linear relationship in the polarization versus extinction diagram, up to $A_V \sim 15$ mag toward the stars with the greatest obscuration, which verified that our observations and analysis provide an accurate depiction of the core.

Transmission spectrum surveys have suggested the ubiquity of high-altitude clouds in exoplanetary atmospheres. Theoretical studies have investigated the formation processes of the high-altitude clouds; however, cloud particles have been commonly approximated as compact spheres, which is not always true for solid mineral particles that likely constitute exoplanetary clouds. Here, we investigate how the porosity of cloud particles evolve in exoplanetary atmospheres and influence the cloud vertical profiles. We first construct a porosity evolution model that takes into account the fractal aggregation and the compression of cloud particle aggregates. Using a cloud microphysical model coupled with the porosity model, we demonstrate that the particle internal density can significantly decrease during the cloud formation. As a result, fluffy-aggregate clouds ascend to altitude much higher than that for compact-sphere clouds assumed so far. We also examine how the fluffy-aggregate clouds affect transmission spectra. We find that the clouds largely obscure the molecular features and produce a spectral slope originated by the scattering properties of aggregates. Finally, we compare the synthetic spectra with the observations of GJ1214 b and find that its flat spectrum could be explained if the atmospheric metallicity is sufficiently high ($\ge100\times$ solar) and the monomer size is sufficiently small ($r_{\rm mon}<1~{\rm {\mu}m}$). The high-metallicity atmosphere may offer the clues to explore the gas accretion processes onto past GJ1214b.

The detailed magnetic field structure of the starless dense core CB81 (L1774, Pipe 42) in the Pipe Nebula was determined based on near-infrared polarimetric observations of background stars to measure dichroically polarized light produced by magnetically aligned dust grains in the core. The magnetic fields pervading CB81 were mapped using 147 stars and axisymmetrically distorted hourglass-like fields were identified. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles in the plane-of-sky and line-of-sight directions were determined to be $4^{\circ} \pm 8^{\circ}$ and $20^{\circ} \pm 20^{\circ}$, respectively. The total magnetic field strength of CB81 was found to be $7.2 \pm 2.3$ $\mu{\rm G}$. Taking into account the effects of thermal/turbulent pressure and magnetic fields, the critical mass of CB81 was calculated to be $M_{\rm cr}=4.03 \pm 0.40$ M$_{\odot}$, which is close to the observed core mass of $M_{\rm core}=3.37 \pm 0.51$ M$_{\odot}$. We thus conclude that CB81 is in a condition close to the critical state. In addition, a spatial offset of $92''$ was found between the center of magnetic field geometry and the dust extinction distribution; this offset structure could not have been produced by self-gravity. The data also indicate a linear relationship between polarization and extinction up to $A_V \sim 30$ mag going toward the core center. This result confirms that near-infrared polarization can accurately trace the overall magnetic field structure of the core.

Observational data for the hourglass-like magnetic field toward the starless dense core FeSt 1-457 were compared with a flux freezing magnetic field model (Myers et al. 2018). Fitting of the observed plane-of-sky magnetic field using the flux freezing model gave a residual angle dispersion comparable with the results based on a simple three-dimensional parabolic model. The best-fit parameters for the flux freezing model were a line-of-sight magnetic inclination angle of $\gamma_{\rm mag} = 35^{\circ} \pm 15^{\circ}$ and a core center to ambient (background) density contrast of $\rho_{\rm c} / \rho_{\rm bkg} = 75$. The initial density for core formation ($\rho_0$) was estimated to be $\rho_{\rm c} / 75 = 4670$ cm$^{-3}$, which is about one order of magnitude higher than the expected density ($\sim 300$ cm$^{-3}$) for the inter-clump medium of the Pipe Nebula. FeSt 1-457 is likely to have been formed from the accumulation of relatively dense gas, and the relatively dense background column density of $A_V \simeq 5$ mag supports this scenario. The initial radius (core formation radius) $R_0$ and the initial magnetic field strength $B_0$ were obtained to be 0.15 pc ($1.64 R$) and $10.8-14.6$ $\mu$G, respectively. We found that the initial density $\rho_0$ is consistent with the mean density of the nearly critical magnetized filament with magnetic field strength $B_0$ and radius $R_0$. The relatively dense initial condition for core formation can be naturally understood if the origin of the core is the fragmentation of magnetized filaments.

The magnetic field structure, kinematical stability, and evolutionary status of the starless dense core Barnard 68 (B68) are revealed based on the near-infrared polarimetric observations of background stars, measuring the dichroically polarized light produced by aligned dust grains in the core. After subtracting unrelated ambient polarization components, the magnetic fields pervading B68 are mapped using 38 stars and axisymmetrically distorted hourglass-like magnetic fields are obtained, although the evidence for the hourglass field is not very strong. On the basis of simple 2D and 3D magnetic field modeling, the magnetic inclination angles on the plane-of-sky and in the line-of-sight direction are determined to be $47^{\circ} \pm 5^{\circ}$ and $20^{\circ} \pm 10^{\circ}$, respectively. The total magnetic field strength of B68 is obtained to be $26.1 \pm 8.7$ $\mu {\rm G}$. The critical mass of B68, evaluated using both magnetic and thermal/turbulent support, is $M_{\rm cr} = 2.30 \pm 0.20$ ${\rm M}_{\odot}$, which is consistent with the observed core mass of $M_{\rm core}=2.1$ M$_{\odot}$, suggesting nearly critical state. We found a relatively linear relationship between polarization and extinction up to $A_V \sim 30$ mag toward the stars with deepest obscuration. Further theoretical and observational studies are required to explain the dust alignment in cold and dense regions in the core.

Dust coagulation in a protoplanetary disk is the first step of planetesimal formation. However, a pathway from dust aggregates to planetesimals remains unclear. Both numerical simulations and laboratory experiments have suggested the importance of dust structure in planetesimal formation, but it is not well constrained by observations. We study how dust structure and porosity alters polarimetric images at millimeter wavelength by performing 3D radiative transfer simulations. Aggregates with different porosity and fractal dimension are considered. As a result, we find that dust aggregates with lower porosity and/or higher fractal dimension are favorable to explain observed millimeter-wave scattering polarization of disks. Aggregates with extremely high porosity fail to explain the observations. In addition, we also show that particles with moderate porosity show weak wavelength dependence of scattering polarization, indicating that multi-wavelength polarimetry is useful to constrain dust porosity. Finally, we discuss implications for dust evolution and planetesimal formation in disks.

Recent (sub)millimeter polarimetric observations toward the young star HL Tau have successfully detected polarization emission from its circumstellar disk. The polarization pattern observed at 0.87 mm is uniform and parallel to the disk's minor axis, consistent with the self-scattering of thermal emission by dust particles whose maximum radius is $\approx 100~\mu m$. However, this maximum size is considerably smaller than anticipated from dust evolution models that assume a high sticking efficiency for icy particles. Here, we propose that the unexpectedly small particle size can be explained if CO$_2$ ice covers the particles in the outer region of the HL Tau disk. CO$_2$ ice is one of the most major interstellar ices, and laboratory experiments show that it is poorly sticky. Based on dust evolution models accounting for CO$_2$ ice mantles as well as aggregate sintering and post-processing radiative transfer, we simulate the polarimetric observation of HL Tau at 0.87 mm. We find that the models with CO$_2$ ice mantles better match the observation. These models also predict that only particles lying between the H$_2$O and CO$_2$ snow lines can grow to millimeter to centimeter sizes, and that their rapid inward drift results in a local dust gap similar to the 10 au gap of the HL Tau disk. We also suggest that the millimeter spectral index for the outer part of the HL Tau disk is largely controlled by the optical thickness of this region and does not necessarily indicate dust growth to millimeter sizes.

We find that the polarimetric observations of protoplanetary disks are useful to search for ultra-light axion dark matter. Axion dark matter predicts the rotation of the linear polarization plane of propagating light, and protoplanetary disks are ideal targets to observe it. We show that a recent observation puts the tightest constraint on the axion-photon coupling constant for axion mass $m\lesssim10^{-21}$eV.

Polarized dust emission outside of disks reveal the magnetic field morphology of molecular clouds. Within disks, however, polarized dust emission can arise from very different mechanisms (e.g., self-scattering), and each of them are useful for constraining physical properties in the disk. For example, these mechanisms allow us to constrain the disk grain size distributions and grain/disk geometries, independent from current methods of measuring these parameters. To accurately model these features and disentangle the various polarization mechanisms, multiwavelength observations at very high resolution and sensitivity are required. With significant upgrades to current interferometric facilities, we can understand how grains evolve in disks during the planet formation process.

Dust grains emit intrinsic polarized emission if they are elongated and aligned in the same direction. The direction of the grain alignment is determined by external forces, such as magnetic fields, radiation, and gas flow against the dust grains. In this letter, we apply the concept of the grain alignment by gas flow, which is called mechanical alignment, to the situation of a protoplanetary disk. We assume that grains have a certain helicity, which results in the alignment with the minor axis parallel to the grain velocity against the ambient disk gas and discuss the morphology of polarization vectors in a protoplanetary disk. We find that the direction of the polarization vectors depends on the Stokes number, which denotes how well grains are coupled to the gas. If the Stokes number is less than unity, orientation of polarization is in the azimuthal direction since the dust velocity against the gas is in the radial direction. If the Stokes number is as large as unity, the polarization vectors show a leading spiral pattern since the radial and azimuthal components of the gas velocity against the dust grains are comparable. This suggests that if the observed polarization vectors show a leading spiral pattern, it would indicate that Stokes number of dust grains is around unity, which is presumably radially drifting.

We study scattered light properties of protoplanetary discs at near-infrared wavelengths for various dust size and structure by performing radiative transfer simulations. We show that different dust structures might be probed by measuring disk polarisation fraction as long as the dust radius is larger than the wavelength. When the radius is larger than the wavelength, disc scattered light will be highly polarised for highly porous dust aggregates, whereas more compact dust structure tends to show low polarisation fraction. Next, roles of monomer radius and fractal dimension for scattered light colours are studied. We find that, outside the Rayleigh regime, as fractal dimension or monomer radius increases, colours of the effective albedo at near-infrared wavelengths vary from blue to red. Our results imply that discs showing grey or slightly blue colours and high polarisation fraction in near-infrared wavelengths might be explained by the presence of large porous aggregates containing sub-microns sized monomers.

The relationship between submillimeter (submm) dust emission polarization and near-infrared (NIR) $H$-band polarization produced by dust dichroic extinction was studied for the cold starless dense core FeSt 1-457. The distribution of polarization angles ($90^{\circ}$-rotated for submm) and degrees were found to be very different between at submm and NIR wavelengths. The mean polarization angles for FeSt 1-457 at submm and NIR wavelengths are $132.1^{\circ} \pm 22.0^{\circ}$ and $2.7^{\circ} \pm 16.2^{\circ}$, respectively. The correlation between $P_H$ and $A_V$ was found to be linear from outermost regions to relatively dense line of sight of $A_V \approx 25$ mag, indicating that NIR polarization reflects overall polarization (magnetic field) structure of the core at least in this density range. The flat $P_H/A_V$ versus $A_V$ correlations were confirmed, and the polarization efficiency was found to be comparable to the observational upper limit (Jones 1989). On the other hand, as reported by Alves et al., submm polarization degrees show clear linearly decreasing trend against $A_V$ from $A_V \approx 20$ mag to the densest center ($A_V \approx 41$ mag), appearing as "polarization hole" structure. The power law index for the $P_{\rm submm}$ versus $A_V$ relationship was obtained to be $\approx -1$, indicating that the alignment for the submm sensitive dust is lost. These very different polarization distributions at submm and NIR wavelengths suggest that (1) there is different radiation environment at these wavelengths or (2) submm-sensitive dust is localized or the combination of them.

Optical properties of dust aggregates are important at various astrophysical environments. To find a reliable approximation method for optical properties of dust aggregates, we calculate the opacity and the asymmetry parameter of dust aggregates by using a rigorous numerical method, the T-Matrix Method (TMM), and then the results are compared to those obtained by approximate methods; the Rayleigh-Gans-Debye (RGD) theory, the effective medium theory (EMT), and the distribution of hollow spheres method (DHS). First of all, we confirm that the RGD theory breaks down when multiple scattering is important. In addition, we find that both EMT and DHS fail to reproduce the optical properties of dust aggregates with fractal dimension of 2 when the incident wavelength is shorter than the aggregate radius. In order to solve these problems, we test the mean field theory (MFT), where multiple scattering can be taken into account. We show that the extinction opacity of dust aggregates can be well reproduced by MFT. However, it is also shown that MFT is not able to reproduce the scattering and absorption opacities when multiple scattering is important. We successfully resolve this weak point of MFT, by newly developing a modified mean field theory (MMF). Hence, we conclude that MMF can be a useful tool to investigate radiative transfer properties of various astrophysical environments. We also point out an enhancement of the absorption opacity of dust aggregates in the Rayleigh domain, which would be important to explain the large millimeter-wave opacity inferred from observations of protoplanetary disks.

We estimate the dust torus cooling timescale once the active galactic nucleus (AGN) is quenched. In a clumpy torus system, once the incoming photons are suppressed, the cooling timescale of one clump from $T_{\rm dust}=1000$ K to several $10$ K is less than $10$ years, indicating that the dust torus cooling time is mainly governed by the light crossing time of the torus from the central engine. After considering the light crossing time of the torus, the AGN torus emission at $12~\mu$m becomes over two orders of magnitude fainter within $100$ years after the quenching. We also propose that those "dying" AGN could be found using the AGN indicators with different physical scale $R$ such as $12~\mu$m band luminosity tracing AGN torus ($R \sim 10$ pc) and the optical [OIII]$\lambda5007$ emission line narrow line regions ($R=10^{2-4}$ pc).

We apply the theory of radiative torque (RAT) alignment for studying protoplanetary disks around a T-Tauri star and perform 3D radiative transfer calculations to provide the expected maps of polarized radiation to be compared with observations, such as with ALMA. We revisit the issue of grain alignment for large grains expected in the protoplanetary disks and find that mm-sized grains at midplane do not align with the magnetic field as the Larmor precession timescale for such large grains becomes longer than the gaseous damping timescale. Hence, for these grains the RAT theory predicts that the alignment axis is determined by the grain precession with respect to the radiative flux. As a result, we expect that the polarization will be in the azimuthal direction for a face-on disk. It is also shown that if dust grains have superparamagnetic inclusions, magnetic field alignment is possible for (sub-)micron grains at the surface layer of disks, and this can be tested by mid-infrared polarimetric observations.

In protoplanetary disks, micron-sized dust grains coagulate to form highly porous dust aggregates. Because the optical properties of these aggregates are not completely understood, it is important to investigate how porous dust aggregates scatter light. In this study, the light scattering properties of porous dust aggregates were calculated using a rigorous method, the T-matrix method, and the results were then compared with those obtained using the Rayleigh-Gans-Debye (RGD) theory and Mie theory with the effective medium approximation (EMT). The RGD theory is applicable to moderately large aggregates made of nearly transparent monomers. This study considered two types of porous dust aggregates, ballistic cluster-cluster agglomerates (BCCAs) and ballistic particle-cluster agglomerates (BPCAs). First, the angular dependence of the scattered intensity was shown to reflect the hierarchical structure of dust aggregates; the large-scale structure of the aggregates is responsible for the intensity at small scattering angles, and their small-scale structure determines the intensity at large scattering angles. Second, it was determined that the EMT underestimates the backward scattering intensity by multiple orders of magnitude, especially in BCCAs, because the EMT averages the structure within the size of the aggregates. It was concluded that the RGD theory is a very useful method for calculating the optical properties of BCCAs.

We study the dust motion at the surface layer of protoplanetary disks. Dust grains in surface layer migrate outward due to angular momentum transport via gas-drag force induced by the stellar radiation pressure. In this study, we calculate mass flux of the outward motion of compact grains and porous dust aggregates by the radiation pressure. The radiation pressure force for porous dust aggregates is calculated using the T-Matrix Method for the Clusters of Spheres. First, we confirm that porous dust aggregates are forced by strong radiation pressure even if they grow to be larger aggregates in contrast to homogeneous and spherical compact grains to which efficiency of radiation pressure becomes lower when their sizes increase. In addition, we find that the outward mass flux of porous dust aggregates with monomer size of 0.1 $\mu$m is larger than that of compact grains by an order of magnitude at the disk radius of 1 AU, when their sizes are several microns. This implies that large compact grains like calcium-aluminum rich inclusions (CAIs) are hardly transported to outer region by stellar radiation pressure, whereas porous dust aggregates like chondritic-porous interplanetary dust particles (CP-IDPs) are efficiently transported to comet formation region. Crystalline silicates are possibly transported in porous dust aggregates by stellar radiation pressure from inner hot region to outer cold cometary region in the protosolar nebula.