• We report a muon spin rotation/relaxation ($\mu$SR) study of single-crystalline samples of the $\alpha$-RuCl$_3$ honeycomb magnet, which is presumed to be a model compound for the Kitaev-Heisenberg interaction. It is inferred from magnetic susceptibility and specific-heat measurements that the present samples exhibit successive magnetic transitions at different critical temperatures $T_{\rm N}$ with decreasing temperature, eventually falling into the $T_{\rm N}=7$ K antiferromagnetic (7 K) phase that has been observed in only single-crystalline specimens with the least stacking fault. Via $\mu$SR measurements conducted under a zero external field, we show that such behavior originates from a phase separation induced by the honeycomb plane stacking fault, yielding multiple domains with different $T_{\rm N}$'s. We also perform $\mu$SR measurements under a transverse field in the paramagnetic phase to identify the muon site from the muon-Ru hyperfine parameters. Based on a comparison of the experimental and calculated internal fields at the muon site for the two possible spin structures inferred from neutron diffraction data, we suggest a modulated zig-zag spin structure for the 7 K phase, with the amplitude of the ordered magnetic moment being significantly reduced from that expected for the orbital quenched spin-1/2 state.
  • Optical properties of dust aggregates become 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 becomes 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 once multiple scattering becomes important, MFT is not able to reproduce the scattering and absorption opacities. 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 discuss 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.
  • Micro- and nano-electromechanical resonators are a fundamental building block of modern technology, used in environmental monitoring, robotics, medical tools as well as fundamental science. These devices rely on dedicated electronics to generate their driving signal, resulting in an increased complexity and size. Here, we present a new paradigm to achieve high-frequency mechanical actuation based on the metal-insulator transition of VO$\mathrm{_2}$, where the steep variation of its electronic properties enables to realize high-frequency electrical oscillations. The dual nature of this phase change, which is both electronic and structural, turns the electrical oscillations into an intrinsic actuation mechanism, powered by a small DC voltage and capable to selectively excite the different mechanical modes of a microstructure. Our results pave the way towards the realization of micro- and nano-electro-mechanical systems with autonomous actuation from integrated DC power sources such as solar cells or micro-batteries.
  • CsFeBr$_3$ is an $S\,{=}\,1$ hexagonal antiferromagnet that has a singlet ground state owing to its large easy-plane single-ion anisotropy. The critical behavior of the magnetic-field-induced phase transition for a magnetic field parallel to the $c$ axis, which can be described by the Bose-Einstein condensation (BEC) of magnons under the $U(1)$ symmetry, was investigated via magnetization and specific heat measurements down to 0.1$\,$K. For the specific heat measurement, we have developed a method of effectively suppressing the torque acting on a sample with strong anisotropy that uses the spin dimer compound Ba$_2$CoSi$_2$O$_6$Cl$_2$ with large and anisotropic Van Vleck paramagnetism. The temperature dependence of the transition field $H_{\rm c}(T)$ was found to follow the power-law $H_{\rm c}(T)\,{-}\,H_{\rm c}\,{\propto}\,T^{\phi}$ with a critical exponent of ${\phi}\,{=}\,1.50\,{\pm}\,0.02$ and critical field of $H_{\rm c}\,{=}\,2.60$ T. This result verifies the universality of the three-dimensional BEC of magnons described by ${\phi}_{\rm BEC}\,{=}\,3/2$.
  • We investigate the properties of quark mass functions in quantum chromodynamics calculated by the Schwinger-Dyson equation in the strong coupling region, in which the loop integration is performed in Minkowski space. The calculated results are compared with those obtained by integration in Euclidean space.
  • We have investigated the longitudinal thermal conductivity of $\alpha$-RuCl$_{3}$, the magnetic state of which is considered to be proximate to a Kitaev honeycomb model, along with the spin susceptibility and magnetic specific heat. We found that the temperature dependence of the thermal conductivity exhibits an additional peak around 100 K, which is well above the phonon peak temperature ($\sim$ 50 K). The higher-temperature peak position is comparable to the temperature scale of the Kitaev couplings rather than the N\'eel temperatures below 15 K. The additional heat conduction was observed for all five samples used in this study, and was found to be rather immune to a structural phase transition of $\alpha$-RuCl$_{3}$, which suggests its different origin from phonons. Combined with experimental results of the magnetic specific heat, our transport measurement suggests strongly that the higher-temperature peak in the thermal conductivity is attributed to itinerant spin excitations associated with the Kitaev couplings of $\alpha$-RuCl$_{3}$. A kinetic approximation of the magnetic thermal conductivity yields a mean free path of $\sim$ 20 nm at 100 K, which is well longer than the nearest Ru-Ru distance ($\sim$ 3 \AA), suggesting the long-distance coherent propagation of magnetic excitations driven by the Kitaev couplings.
  • K$_2$CuF$_4$ is magnetically described as a spin-1/2, quasi-two-dimensional (2D), square-lattice XXZ ferromagnet with weak easy-plane anisotropy. The magnetic ordering for an applied magnetic field $H$ parallel to the $c$ axis is equivalent to the Bose\,--\,Einstein condensation (BEC) of lattice bosons, as discussed by Matsubara and Matsuda [Prog. Theor. Phys. \textbf{16}, 569 (1956)]. Magnetization and specific heat measurements were performed to obtain the temperature vs magnetic field phase diagram for $H\,{\parallel}\,c$. The phase boundary between polarized and ordered phases was found to be expressed by the power law $H_{\rm c}(T)\,{-}\,H_{\rm c}(0)\,{\propto}\,T^\phi$ with exponent $\phi\,{\approx}\,1.0$ in a wide temperature range, in agreement with the theory of quasi-2D BEC.
  • We investigate the properties of fermion mass functions in quantum electrodynamics calculated by the Schwinger-Dyson equation in the strong coupling region, in which the loop integration is performed in Minkowski space. The calculated results without the instantaneous exchange approximation are compared with those obtained by integration in Euclidean space.
  • Refractory organic compounds formed in molecular clouds are among the building blocks of the solar system objects and could be the precursors of organic matter found in primitive meteorites and cometary materials. However, little is known about the evolutionary pathways of molecular cloud organics from dense molecular clouds to planetary systems. In this study, we focus on the evolution of the morphological and viscoelastic properties of molecular cloud refractory organic matter. We found that the organic residue, experimentally synthesized at about 10 K from UV-irradiated H2O-CH3OH-NH3 ice, changed significantly in terms of its nanometer- to micrometer-scale morphology and viscoelastic properties after UV irradiation at room temperature. The dose of this irradiation was equivalent to that experienced after short residence in diffuse clouds (equal or less than 10,000 years) or irradiation in outer protoplanetary disks. The irradiated organic residues became highly porous and more rigid and formed amorphous nanospherules. These nanospherules are morphologically similar to organic nanoglobules observed in the least-altered chondrites, chondritic porous interplanetary dust particles, and cometary samples, suggesting that irradiation of refractory organics could be a possible formation pathway for such nanoglobules. The storage modulus (elasticity) of photo-irradiated organic residues is about 100 MPa irrespective of vibrational frequency, a value that is lower than the storage moduli of minerals and ice. Dust grains coated with such irradiated organics would therefore stick together efficiently, but growth to larger grains might be suppressed due to an increase in aggregate brittleness caused by the strong connections between grains.
  • A programmable micromechanical resonator based on a VO2 thin film is reported. Multiple mechanical eigenfrequency states are programmed using Joule heating as local power source, gradually driving the phase transition of VO2 around its Metal-Insulator transition temperature. Phase coexistence of domains is used to tune the stiffness of the device via local control of internal stresses and mechanical properties. This study opens perspectives for developing mechanically configurable nanostructure arrays.
  • We investigate the resistive switching behaviour of $\mathrm{VO_2}$ microbridges under current bias as a function of temperature and thermal coupling with the heat bath. Upon increasing the electrical current bias, the formation of the metallic phase can progress smoothly or through sharp jumps. The magnitude and threshold current values of these sharp resistance drops show random behaviour and are dramatically influenced by thermal dissipation conditions. Our results also evidence how the propagation of the metallic phase induced by electrical current in $\mathrm{VO_2}$, and thus the shape of the resulting high-conductivity path, are not predictable. We discuss the origin of the switching events through a simple electro-thermal model based on the domain structure of $\mathrm{VO_2}$ films that can be useful to improve the stability and controllability of future $\mathrm{VO_2}$-based devices.
  • We synthesized single crystals of composition Ba$_2$CuSi$_2$O$_6$Cl$_2$ and investigated its quantum magnetic properties. The crystal structure is closely related to that of the quasi-two-dimensional (2D) dimerized magnet BaCuSi$_2$O$_6$ also known as Han purple. Ba$_2$CuSi$_2$O$_6$Cl$_2$ has a singlet ground state with an excitation gap of ${\Delta}/k_{\rm B}\,{=}\,20.8$ K. The magnetization curves for two different field directions almost perfectly coincide when normalized by the $g$-factor except for a small jump anomaly for a magnetic field perpendicular to the $c$ axis. The magnetization curve with a nonlinear slope above the critical field is in excellent agreement with exact-diagonalization calculations based on a 2D coupled spin-dimer model. Individual exchange constants are also evaluated using density functional theory (DFT). The DFT results demonstrate a 2D exchange network and weak frustration between interdimer exchange interactions, supported by weak spin-lattice coupling implied from our magnetostriction data. The magnetic-field-induced spin ordering in Ba$_2$CuSi$_2$O$_6$Cl$_2$ is described as the quasi-2D Bose-Einstein condensation of triplets.
  • We have performed magnetization measurements of the gapped quantum magnet CsFeCl$_3$ at temperatures ($T$) down to 0.5\,K at ambient pressure and down to 1.8\,K at hydrostatic pressures ($P$) of up to 1.5\,GPa. The lower-field ($H$) phase boundary of the field-induced ordered phase at ambient pressure is found to follow the power-law behavior expressed by the formula $H_{\rm N}(T)$\,$-$\,$H_{\rm c}$\,$\propto$\,$T_{\rm N}^{\phi}$. The application of pressure extends the phase boundary to both a lower field and higher temperature. Above the critical pressure $P_{\rm c}$\,$\sim$\,0.9\,GPa, the transition field $H_{\rm N}$ associated with the excitation gap becomes zero, and a signature of the magnetic phase transition is found in the $T$-dependence of magnetization in a very low applied field. This suggests that CsFeCl$_3$ exhibits a pressure-induced magnetic phase transition at $P_{\rm c}$.
  • 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.
  • A giant planet creates a gap in a protoplanetary disk, which might explain the observed gaps in protoplanetary disks. The width and depth of the gaps depend on the planet mass and disk properties. We have performed two--dimensional hydrodynamic simulations for various planet masses, disk aspect ratios and viscosities, to obtain an empirical formula for the gap width. The gap width is proportional to the square root of the planet mass, -3/4 power of the disk aspect ratio and -1/4 power of the viscosity. This empirical formula enables us to estimate the mass of a planet embedded in the disk from the width of an observed gap. We have applied the empirical formula for the gap width to the disk around HL~Tau, assuming that each gap observed by ALMA observations is produced by planets, and discussed the planet masses within the gaps. The estimate of planet masses from the gap widths is less affected by the observational resolution and dust filtration than that from the gap depth.
  • The magnetic structure of the double perovskite compound Sr$_2$CuTeO$_6$ was determined from neutron powder diffraction data. This material is magnetically described as an $S\,{=}\,1/2$ quasi-two-dimensional square-lattice Heisenberg model with antiferromagnetic nearest-neighbor and next-nearest-neighbor interactions. Sr$_2$CuTeO$_6$ undergoes a magnetic phase transition at $T_{\rm N}\,{\simeq}\,29\,$ K. The spin structure below $T_{\rm N}$ is N\'{e}el antiferromagnetic on the square lattice, which means that the nearest-neighbor interaction ($J_1$) is stronger than the next-nearest-neighbor interaction ($J_2$), in contrast to other isostructural compounds such as Ba$_2$CuWO$_6$ and Sr$_2$CuWO$_6$, for which $|J_1|\,{<}\,|J_2|$ is realized.
  • The latest observation of HL Tau by ALMA revealed spectacular concentric dust rings in its circumstellar disk. We attempt to explain the multiple ring structure as a consequence of aggregate sintering. Sintering is known to reduce the sticking efficiency of dust aggregates and occurs at temperatures slightly below the sublimation point of their constituent material. We here present a dust growth model incorporating sintering and use it to simulate global dust evolution due to sintering, coagulation, fragmentation, and radial inward drift in a modeled HL Tau disk. We show that aggregates consisting of multiple species of volatile ices experience sintering, collisionally disrupt, and pile up at multiple locations slightly outside the snow lines of the volatiles. At wavelengths of 0.87--1.3 mm, these sintering zones appear as bright, optically thick rings with a spectral slope of $\approx 2$, whereas the non-sintering zones as darker, optically thinner rings of a spectral slope of $\approx$ 2.3--2.5. The observational features of the sintering and non-sintering zones are consistent with those of the major bright and dark rings found in the HL Tau disk, respectively. Radial pileup and vertical settling occur simultaneously if disk turbulence is weak and if monomers constituting the aggregates are $\sim 1~{\rm \mu m}$ in radius. For the radial gas temperature profile of $T = 310(r/1~{\rm AU})^{-0.57}~{\rm K}$, our model perfectly reproduces the brightness temperatures of the optically thick bright rings, and reproduces their orbital distances to an accuracy of $\lesssim$ 30%.
  • When planetesimals grow via collisions in a turbulent disk, stirring through density fluctuation caused by turbulence effectively increases the relative velocities between planetesimals, which suppresses the onset of runaway growth. We investigate the onset of runaway growth in a turbulent disk through simulations that calculate the mass and velocity evolution of planetesimals. When planetesimals are small, the average relative velocity between planetesimals, $v_{\rm r}$, is much greater than their surface escape velocity, $v_{\rm esc}$, so that runaway growth does not occur. As planetesimals become large via collisional growth, $v_{\rm r}$ approaches $v_{\rm esc}$. When $v_{\rm r} \approx 1.5 v_{\rm esc}$, runaway growth of the planetesimals occurs. During the oligarchic growth subsequent to runaway growth, a small number of planetary embryos produced via runaway growth become massive through collisions with planetesimals with radii of that at the onset of runaway growth, $r_{\rm p,run}$. We analytically derive $r_{\rm p,run}$ as a function of the turbulent strength. Growing $\sim 10\,M_\oplus$ embryos that are suitable to become the cores of Jupiter and Saturn requires $r_{\rm p,run} \sim 100$\,km, which is similar to the proposed fossil feature in the size distribution of main belt asteroids. In contrast, the formation of Mars as quickly as suggested from Hf-W isotope studies requires small planetesimals at the onset of runaway growth. Thus, the conditions required to form Mars, Jupiter, and Saturn and the size distribution of the main-belt asteroids indicate that the turbulence increased in amplitude relative to the sound speed with increasing distance from the young Sun.
  • We report the results of magnetization and specific heat measurements of Ba$_2$CoTeO$_6$ composed of two subsystems A and B, which are magnetically described as an $S\,{=}\,1/2$ triangular-lattice Heisenberg-like antiferromagnet and a $J_1-J_2$ honeycomb-lattice Ising-like antiferromagnet, respectively. These two subsystems were found to be approximately decoupled. Ba$_2$CoTeO$_6$ undergoes magnetic phase transitions at $T_{\rm N1}\,{=}\,12.0$ K and $T_{\rm N2}\,{=}\,3.0$ K, which can be interpreted as the orderings of subsystems B and A, respectively. Subsystem A exhibits a magnetization plateau at one-third of the saturation magnetization for the magnetic field $H$ perpendicular to the $c$ axis owing to the quantum order-by-disorder, whereas for $H\,{\parallel}\,c$, subsystem B shows three-step metamagnetic transitions with magnetization plateaus at zero, one-third and one-half of the saturation magnetization. The analysis of the magnetization process for subsystem B shows that the classical ground states at these plateaus are infinitely degenerate within the Ising model.
  • Elementary excitations of the S=1/2 one-dimensional antiferromagnet KCuGaF_6 were investigated by inelastic neutron scattering in zero and finite magnetic fields perpendicular to the (1, 1, 0) plane combined with specific heat measurements. KCuGaF$_6$ exhibits no long-range magnetic ordering down to 50 mK despite the large exchange interaction J/k_B=103 K. At zero magnetic field, well-defined spinon excitations were observed. The energy of the des Cloizeaux and Pearson mode of the spinon excitations is somewhat larger than that calculated with the above exchange constant. This discrepancy is mostly ascribed to the effective XY anisotropy arising from the large Dzyaloshinsky-Moriya interaction with an alternating D vector. KCuGaF_6 in a magnetic field is represented by the quantum sine-Gordon model, for which low-energy elementary excitations are composed of solitons, antisolitons and their bound states called breathers. Unlike the theoretical prediction, it was found that the energy of a soliton is smaller than that of the first breather, although the energy of the first breather coincides with that observed in a previous ESR measurement.
  • Firstly, we study the final masses of giant planets growing in protoplanetary disks through capture of disk gas, by employing an empirical formula for the gas capture rate and a shallow disk gap model, which are both based on hydrodynamical simulations. The shallow disk gaps cannot terminate growth of giant planets. For planets less massive than 10 Jupiter masses, their growth rates are mainly controlled by the gas supply through the global disk accretion, rather than their gaps. The insufficient gas supply compared with the rapid gas capture causes a depletion of the gas surface density even at the outside of the gap, which can create an inner hole in the protoplanetary disk. Our model can also predict the depleted gas surface density in the inner hole for a given planet mass. Secondly, our findings are applied to the formation of our solar system. For the formation of Jupiter, a very low-mass gas disk with a few or several Jupiter masses is required at the beginning of its gas capture because of the non-stopping capture. Such a low-mass gas disk with sufficient solid material can be formed through viscous evolution from an initially $\sim$10AU-sized compact disk with the solar composition. By the viscous evolution with a moderate viscosity of $\alpha \sim 10^{-3}$, most of disk gas accretes onto the sun and a widely spread low-mass gas disk remains when the solid core of Jupiter starts gas capture at $t \sim 10^7$ yrs. The depletion of the disk gas is suitable for explaining the high metallicity in giant planets of our solar system. A very low-mass gas disk also provides a plausible path where type I and II planetary migrations are both suppressed significantly. In particular, we also show that the type II migration of Jupiter-size planets becomes inefficient because of the additional gas depletion due to the rapid gas capture by themselves.
  • Magnetoelectric effects in the chiral magnet CsCuCl$_3$ have been investigated through the magnetization and electric polarization measurements in pulsed high magnetic fields. If the magnetic field is applied normal to the spin chiral c-axis, a plateau and a jump in magnetization are observed. Between the plateau and the jump in magnetization, a small but significant electric polarization of about 0.25 {\mu}C/m$^2$ along the a-axis is observed with the field applied almost perpendicular to the ac-plane. In addition, the paramagnetoelectric effect expected from the point group symmetry is confirmed. The emergence of the electric polarization is explained by the cooperation of the local electric polarization on the chiral solitonic spin arrangement and the paramagnetoelectric effect. Hence, we interpret this novel multiferroic phase in CsCuCl$_3$ as a polar solitonic phase.
  • We have systematically investigated the variation of the exchange parameters and the ground state in the $S = 1/2$ kagome-lattice antiferromagnet (Rb$_{1-x}$Cs$_{x}$)$_2$Cu$_3$SnF$_{12}$, via magnetic measurements using single crystals. One of the parent compounds, Rb$_2$Cu$_3$SnF$_{12}$, which has a distorted kagome lattice accompanied by four sorts of nearest-neighbor exchange interaction, has a disordered ground state described by a pinwheel valence-bond-solid state. The other parent compound, Cs$_2$Cu$_3$SnF$_{12}$, which has a uniform kagome lattice at room temperature, has an ordered ground state with the $q$ = 0 spin structure. The analysis of magnetic susceptibilities shows that with increasing cesium concentration $x$, the exchange parameters increase with the tendency to be uniform. It was found that the ground state is disordered for $x$ < 0.53 and ordered for $x$ > 0.53. The pseudogap observed for $x$ < 0.53 and the N\'{e}el temperature for $x$ > 0.53 approach zero at $x_{\rm c}$ $\simeq$ 0.53. This is indicative of the occurrence of a quantum phase transition at $x_{\rm c}$.
  • Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at \lambda = 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.
  • A giant planet embedded in a protoplanetary disk forms a gap. An analytic relationship among the gap depth, planet mass $M_{p}$, disk aspect ratio $h_p$, and viscosity $\alpha$ has been found recently, and the gap depth can be written in terms of a single parameter $K= (M_{p}/M_{\ast})^2 h_p^{-5} \alpha^{-1}$. We discuss how observed gap features can be used to constrain the disk and/or planet parameters based on the analytic formula for the gap depth. The constraint on the disk aspect ratio is critical in determining the planet mass so the combination of the observations of the temperature and the image can provide a constraint on the planet mass. We apply the formula for the gap depth to observations of HL~Tau and HD~169142. In the case of HL~Tau, we propose that a planet with $\gtrsim 0.3$ is responsible for the observed gap at $30$~AU from the central star based on the estimate that the gap depth is $\lesssim 1/3$. In the case of HD~169142, the planet mass that causes the gap structure recently found by VLA is $\gtrsim 0.4 M_J$. We also argue that the spiral structure, if observed, can be used to estimate the lower limit of the disk aspect ratio and the planet mass.