Source author record

Hidekazu Tanaka

Hidekazu Tanaka appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

cond-mat.str-el astro-ph.EP astro-ph.SR hep-ph cond-mat.soft physics.chem-ph cond-mat.mtrl-sci cond-mat.stat-mech hep-ex hep-lat nucl-th physics.flu-dyn astro-ph.GA astro-ph.HE cond-mat.mes-hall hep-th nucl-ex physics.atm-clus physics.comp-ph

Catalog footprint

What is connected

75works

19topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Contrasting Magnetic Structures in SrLaCuSbO$_{6}$ and SrLaCuNbO$_{6}$: Spin-1/2 Quasi-Square-Lattice $J_1 {-} J_2$ Heisenberg Antiferromagnets

We report the magnetic properties of the double perovskites SrLaCuSbO$_6$ (SLCSO) and SrLaCuNbO$_6$ (SLCNO). The temperature dependence of the magnetic susceptibilities of both compounds shows a broad maximum characteristic of an $S\,=\,1/2$ square lattice Heisenberg antiferromagnet. Magnetic ordering occurs at $T_{\rm N}\,{=}\,13.6$ and 15.7 K for SLCSO and SLCNO, respectively. Neutron powder diffraction measurements reveal contrasting spin structures in both compounds. The spin structures of SLCSO and SLCNO below $T_{\rm N}$ are Néel antiferromagnetic and collinear antiferromagnetic, respectively. This result demonstrates that the nearest-neighbor interaction is dominant in SLCSO, whereas the next-nearest-neighbor interaction is dominant in SLCNO. The magnitude of the ordered moment was evaluated at 3.5 K to be $m\,{=}\,0.39(3)\,μ_{\rm B}$ for SLCSO and $0.37(1)\,μ_{\rm B}$ for SLCNO, which are significantly smaller than those calculated using linear spin wave theory. We infer that the small ordered moment is caused by the effect of exchange bond randomness arising from the site disorder of Sr and La ions.

preprint2022arXiv

Impacts of viscous dissipation on collisional growth and fragmentation of dust aggregates

Understanding the collisional behavior of dust aggregates consisting of submicron-sized grains is essential to unveiling how planetesimals formed in protoplanetary disks. It is known that the collisional behavior of individual dust particles strongly depends on the strength of viscous dissipation force; however, impacts of viscous dissipation on the collisional behavior of dust aggregates have not been studied in detail, especially for the cases of oblique collisions. Here we investigated the impacts of viscous dissipation on the collisional behavior of dust aggregates. We performed numerical simulations of collisions between two equal-mass dust aggregates with various collision velocities and impact parameters. We also changed the strength of viscous dissipation force systematically. We found that the threshold collision velocity for the fragmentation of dust aggregates barely depends on the strength of viscous dissipation force when we consider oblique collisions. In contrast, the size distribution of fragments changes significantly when the viscous dissipation force is considered. We obtained the empirical fitting formulae for the size distribution of fragments for the case of strong dissipation, which would be useful to study the evolution of size and spatial distributions of dust aggregates in protoplanetary disks.

preprint2022arXiv

Structures of magnetic excitations in the spin-1/2 kagome-lattice antiferromagnets Cs$_2$Cu$_3$SnF$_{12}$ and Rb$_2$Cu$_3$SnF$_{12}$

We show the structures of magnetic excitations in spin-1/2 kagome-lattice antiferromagnets Cs$_2$Cu$_3$SnF$_{12}$ and Rb$_2$Cu$_3$SnF$_{12}$ investigated by inelastic neutron scattering in wide energy and momentum ranges. For Cs$_2$Cu$_3$SnF$_{12}$, four single-magnon excitation modes were observed. Low-energy three modes are assigned to be transverse modes and the high-energy fourth mode is suggested to be an amplitude mode. It was found that the broad excitation continuum without a marked structure spreads in a wide energy range from $0.15J$ to approximately $2.5J$ in contrast to the clearly structured excitation continuum observed in the spin-1/2 triangular-lattice Heisenberg antiferromagnet. These findings strongly suggest spinon excitations as elementary excitations in Cs$_2$Cu$_3$SnF$_{12}$. In Rb$_2$Cu$_3$SnF$_{12}$, singlet-triplet excitations from the pinwheel VBS state and their ghost modes caused by the enlargement of the chemical unit cell were clearly confirmed. It was found that the excitation continuum is structured in the low-energy region approximately below $J_{\mathrm{avg}}$ and the almost structureless high-energy excitation continuum extends to approximately $2.6J_{\mathrm{avg}}$. The characteristics of the high-energy excitation continuum are common to both Cs$_2$Cu$_3$SnF$_{12}$ and Rb$_2$Cu$_3$SnF$_{12}$, irrespective of their ground states. The experimental results strongly suggest that the spin liquid component remains in the ground state as quantum fluctuations in Cs$_2$Cu$_3$SnF$_{12}$ and Rb$_2$Cu$_3$SnF$_{12}$.

preprint2021arXiv

Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback

Stellar-mass BHs (sBHs) are predicted to be embedded in active galactic nuclei (AGN) disks due to gravitational drag and in-situ star formation. However, we find that due to a high gas density in an AGN disk environment, compact objects may rapidly grow to intermediate-mass BHs and deplete matter from the AGN disk unless accretion is suppressed by some feedback process(es). These consequences are inconsistent with AGN observations and the dynamics of the Galactic center. Here we consider mechanical feedback mechanisms for the reduction of gas accretion. Rapidly accreting sBHs launch winds and/or jets via the Blandford-Znajek mechanism, which produce high-pressure shocks and cocoons. Such a shock and cocoon can spread laterally in the plane of the disk, eject the outer regions of a circum-sBH disk (CsBD) and puncture a hole in the AGN disk with horizontal size comparable to the disk scale-height. Since the depletion timescale of the bound CsBD is much shorter than the resupply timescale of gas to the sBH, the time-averaged accretion rate onto sBHs is reduced by this process by a factor of $\sim 10$--$100$. This feedback mechanism can therefore help alleviate the sBH over-growth and AGN-disk depletion problems. On the other hand, we find that cocoons of jets can unbind a large fraction of the gas accreting in the disks of less massive SMBHs, which may help explain the dearth of high-Eddington ratio AGNs with SMBH mass $\lesssim10^5{\rm M_\odot}$.

preprint2021arXiv

Eccentric gap induced by a super-Jupiter mass planet

A giant planet embedded in a protoplanetary disk opens a gap by tidal interaction, and properties of the gap strongly depend on the planetary mass and disk parameters. Many numerical simulations of this process have been conducted, but detailed simulations and analysis of gap formation by a super-Jupiter mass planet have not been thoroughly conducted. We performed two-dimensional numerical hydrodynamic simulations of the gap formation process by a super-Jupiter mass planet and examined the eccentricity of the gap. When the planet is massive, the radial motion of gas is excited, causing the eccentricity of the gap's outer edge to increase. Our simulations showed that the critical planetary mass for the eccentric gap was $\sim3~M_{\rm J}$ in a disk with $α=4.0\times10^{-3}$ and $h/r=0.05$, a finding that was consistent with that reported in a previous work. The critical planetary mass for the eccentric gap depends on the viscosity and the disk scale height. We found that the critical mass could be described by considering a dimensionless parameter related to the gap depth. The onset of gap eccentricity enhanced the surface density inside the gap, shallowing the gap more than the empirical relation derived in previous studies for a planet heavier than the critical mass. Therefore, our results suggest that the mass accretion rate, which strongly depends on the gas surface density in the gap is also enhanced for super-Jupiter mass planets. These results may substantially impact the formation and evolution processes of super-Jupiter mass planets and population synthesis calculations.

preprint2021arXiv

Magnons and Spinons in $\mathrm{Ba}_2\mathrm{CoTeO}_6 $: A Composite System of Isolated Spin-$1/2$ Triangular Heisenberg-like and Frustrated Honeycomb Ising-like Antiferromagnets

We report the neutron scattering results on magnetic orderings and excitations in $\mathrm{Ba}_2\mathrm{CoTeO}_6$ composed of two almost isolated subsystems A and B, which are described as an $S\,{=}\,1/2$ triangular Heisenberg-like antiferromagnet and a frustrated honeycomb Ising-like antiferromagnet, respectively. Stripy ordering of subsystem B was confirmed below $T_{\rm N1}\,{=}\,12.0$ K, whereas sharp streaks were observed along $(1/3, 1/3, L)$ and $(2/3, 2/3, L)$ at 0.3 K (${\ll}\,T_{\rm N2}\,{=}\,3.0$ K). This indicates the two-dimensional nature of ordering in subsystem A. It was found that the excitation spectra of both subsystems are well separated and independent of each other. The excitation spectrum of subsystem A is composed of two single-magnon branches with roton-like minima at the M point and a clearly structured intense continuum, as similarly observed in $\mathrm{Ba}_3\mathrm{CoSb}_2\mathrm{O}_9$, which is strongly indicative of spinon excitations. Dispersion curves for subsystem B can be described by linear spin wave theory within the third-neighbor exchange interaction.

preprint2020arXiv

An Experimental and Theoretical Insights into the Dielectric Properties of (Li, Nd) Co-doped ZnO Ceramics

In this work, we report the combined effect of donor (Nd) and acceptor (Li) co-doping at the Zn-site of ZnO ceramics on structural, microstructural and dielectric properties. Combining experimental observations with DFT based theoretical study, we have shown that before experimental fabrication DFT based first principles study can be used as a good indication to have prior qualitative assessment of a dielectric medium. For implementing this objective various Li and Nd co-doped ZnO ceramics have been synthesized through the conventional solid-state reaction route. Quantitative XRD analysis reveals the formation of wurtzite hexagonal structured ZnO having space group P63mc. Meanwhile, FESEM micrographs confirm the formation of randomly aligned non-uniform grains in size and shape. We show that the average grain size distribution and density of the studied compositions are two tuning factors to control the dielectric properties of these compounds. Though the value of dielectric constant is decreased with the increase in doping content, the optimum composition Nd0.005Li0.005Zn0.99O exhibits slightly lower dielectric constant (2066 at 1 KHz) than pristine ZnO but relatively very low dielectric loss (0.20 at 1 KHz) at room temperature than pure ZnO ceramics sintered at 1623 K. For understanding the dielectric relaxation mechanism in the studied ceramics, complex impedance spectra analysis have also been performed and discussed thoroughly. This study provides a new insight for further development of colossal permittivity (CP) ceramics and extends the current understanding of CP mechanism in ceramic materials.

preprint2020arXiv

Comments on "Type II migration strikes back -- An old paradigm for planet migration in discs" by Scardoni et al

In the conventional view of type II migration, a giant planet migrates inward in the viscous velocity of the accretion disc in the so-call disc-dominate case. Recent hydrodynamic simulations, however, showed that planets migrate with velocities much faster than the viscous one in massive discs. Such fast migration cannot be explained by the conventional picture. Scardoni et al. (2020) has recently argued this new picture. By carrying out similar hydrodynamic simulations, they found that the migration velocity slows down with time and eventually reaches the prediction by the conventional theory. They interpreted the fast migration as an initial transient one and concluded that the conventional type II migration is realised after the transient phase. We show that the migration velocities obtained by Scardoni et al. (2020) are consistent with the previous simulations even in the transient phase that they proposed. We also find that the transient fast migration proposed by Scardoni et al. (2020) is well described by a new model of Kanagawa et al. (2018). The new model can appropriately describe significant inward migration during the initial transient phase that Scardoni et al. (2020) termed. Hence, we conclude that the time-variation of the transient migration velocity is due to the changes of the orbital radius of the planet and its background surface density during the migration.

preprint2020arXiv

Coupled spin-1/2 antiferromagnetic chain Cs$_2$LiRuCl$_6$ with partially disordered crystal lattice

We determined the crystal structure of Cs$_2$LiRuCl$_6$, which was synthesized first in this work, and investigated its magnetic properties. Cs$_2$LiRuCl$_6$ has a hexagonal structure composed of linear chains of face-sharing RuCl$_6$ and LiCl$_6$ octahedra. In two-thirds of the structural chains, Ru$^{3+}$ and Li$^+$ sites are almost ordered, while in the other chains their sites are disordered. This situation is analogous to the ground state of the antiferromagnetic Ising model on a triangular lattice. Using electron paramagnetic resonance, we evaluated the $g$ factors of Ru$^{3+}$ with effective spin-1/2 as $g_c\,{=}\,2.72$ and $g_{ab}\,{=}\,1.50$ for magnetic fields $H$ parallel and perpendicular to the $c$ axis, respectively. Magnetization curves for $H\,{\parallel}\,c$ and $H\,{\perp}\,c$ are highly anisotropic. However, these magnetization curves approximately coincide when normalized by the $g$ factors. It was found from the magnetization and specific heat results that Cs$_2$LiRuCl$_6$ can be described as a coupled one-dimensional $S\,{=}\,1/2$ Heisenberg-like antiferromagnet with $J/k_{\rm B}\,{\simeq}\,3.7$ K. Three-dimensional ordering occurs at $T_{\rm N}\,{=}\,0.48$ K. A magnetic phase diagram for $H\,{\parallel}\,c$ is also presented.

preprint2020arXiv

Final Masses of Giant Planets III: Effect of Photoevaporation and a New Planetary Migration Model

We herein develop a new simple model for giant planet formation, which predicts the final mass of a giant planet born in a given disk, by adding the disk mass loss due to photoevaporation and a new type II migration formula to our previous model. The proposed model provides some interesting results. First, it gives universal evolution tracks in the diagram of planetary mass and orbital radius, which clarifies how giant planets migrate at growth in the runaway gas accretion stage. Giant planets with a few Jupiter masses or less suffer only a slight radial migration in the runaway gas accretion stage. Second, the final mass of giant planets is approximately given as a function of only three parameters: the initial disk mass at the starting time of runaway gas accretion onto the planet, the mass loss rate due to photoevaporation, and the starting time. On the other hand, the final planet mass is almost independent of the disk radius, viscosity, and initial orbital radius. The obtained final planet mass is similar to or less than 10% of the initial disk mass. Third, the proposed model successfully explains properties in the mass distribution of giant exoplanets with the mass distribution of observed protoplanetary disks for a reasonable range of the mass loss rate due to photoevaporation.

preprint2020arXiv

Ring formation by coagulation of dust aggregates in early phase of disk evolution around a protostar

Ring structures are observed by (sub-)millimeter dust continuum emission in various circumstellar disks from early stages of Class 0 and I to late stage of Class II young stellar objects (YSOs). In this paper, we study one of the possible scenarios of such ring formation in early stage, which is coagulation of dust aggregates. The dust grains grow in an inside-out manner because the growth timescale is roughly proportional to the orbital period. The boundary of the dust evolution can be regarded as the growth front, where the growth time is comparable to the disk age. With radiative transfer calculations based on the dust coagulation model, we find that the growth front can be observed as a ring structure because dust surface density is sharply changed at this position. Furthermore, we confirm that the observed ring positions in the YSOs with an age of $\lesssim1$ Myr are consistent with the growth front. The growth front could be important to create the ring structure in particular for early stage of the disk evolution such as Class 0 and I sources.

preprint2019arXiv

Quarks mass function at finite density in real-time formalism

Chiral symmetry restoration of quarks is investigated at finite density in quantum chromodynamics. The effective quark mass is calculated with the Schwinger-Dyson equation in the real-time formalism without the instantaneous exchange approximation. We present some properties of the quark mass functions and the quark propagator at zero temperature.

preprint2019arXiv

Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization

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.

preprint2016arXiv

Light scattering by fractal dust aggregates: I. Angular dependence of scattering

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.

preprint2016arXiv

Magnetic Structure of the $S\,{=}\,1/2$ Quasi-Two-Dimensional Square-Lattice Heisenberg Antiferromagnet Sr$_2$CuTeO$_6$

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é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.

preprint2016arXiv

Magnetic-field- and pressure-induced quantum phase transition in CsFeCl$_3$ proved via magnetization measurement

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}^ϕ$. 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}$.

preprint2016arXiv

Mass constraint for a planet in a protoplanetary disk from the gap width

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.

preprint2016arXiv

Quasi-two-dimensional Bose-Einstein condensation of spin triplets in dimerized quantum magnet Ba$_2$CuSi$_2$O$_6$Cl$_2$

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 $Δ/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.

preprint2016arXiv

Sintering-induced Dust Ring Formation in Protoplanetary Disks: Application to the HL Tau Disk

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 μ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%.

preprint2015arXiv

Detailed structure of the outer disk around HD 169142 with polarized light in H-band

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 λ= 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.

preprint2015arXiv

Electrical Oscillation in Pt/VO2 Bilayer Strips

We report on the observation of stable electrical oscillation in Pt/VO2 bilayer strips, in which the Pt overlayer serves the dual purposes of heating up the VO2 and weakening the electric field in the VO2 layer. Systematic measurements in an ultrahigh vacuum nanoprobe system show that the oscillation frequency increases with the bias current and/or with decreasing device dimension. In contrast to most VO2-based oscillators reported to date, which are electrically triggered, current-induced Joule heating in the Pt overlayer is found to play a dominant role in the generation of oscillation in Pt/VO2 bilayers. A simple model involving thermally triggered transition of VO2 on a heat sink is able to account for the experimental observations. The results in this work provide an alternative view of the triggering mechanism in VO2-based oscillators.

preprint2015arXiv

Field-driven successive phase transitions in quasi-two-dimensional frustrated antiferromagnet Ba$_2$CoTeO$_6$ and highly degenerate classical ground states

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.

preprint2015arXiv

Final Masses of Giant Planets II: Jupiter Formation in a Gas-Depleted Disk

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 $α\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.

preprint2015arXiv

Formation of a disc gap induced by a planet: Effect of the deviation from Keplerian disc rotation

The gap formation induced by a giant planet is important in the evolution of the planet and the protoplanetary disc. We examine the gap formation by a planet with a new formulation of one-dimensional viscous discs which takes into account the deviation from Keplerian disc rotation due to the steep gradient of the surface density. This formulation enables us to naturally include the Rayleigh stable condition for the disc rotation. It is found that the derivation from Keplerian disc rotation promotes the radial angular momentum transfer and makes the gap shallower than in the Keplerian case. For deep gaps, this shallowing effect becomes significant due to the Rayleigh condition. In our model, we also take into account the propagation of the density waves excited by the planet, which widens the range of the angular momentum deposition to the disc. The effect of the wave propagation makes the gap wider and shallower than the case with instantaneous wave damping. With these shallowing effects, our one-dimensional gap model is consistent with the recent hydrodynamic simulations.

preprint2015arXiv

From Planetesimals to Planets in Turbulent Protoplanetary Disks I. Onset of Runaway Growth

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.

preprint2015arXiv

Magnetic-Field-Induced Polar Phase in Chiral Magnet CsCuCl$_3$

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 μ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.

preprint2015arXiv

Mass Estimates of a Giant Planet in a Protoplanetary Disk from the Gap Structures

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 $α$ 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} α^{-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.

preprint2015arXiv

Quantum phase transition between disordered and ordered states in the spin-1/2 kagome lattice antiferromagnet (Rb$_{1-x}$Cs$_{x}$)$_2$Cu$_3$SnF$_{12}$

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é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}$.

preprint2015arXiv

Simple improvements to classical bubble nucleation models

We revisit classical nucleation theory (CNT) for the homogeneous bubble nucleation rate and improve the classical formula using a new prefactor in the nucleation rate. Most of the previous theoretical studies have used the constant prefactor determined by the bubble growth due to the evaporation process from the bubble surface. However, the growth of bubbles is also regulated by the thermal conduction, the viscosity, and the inertia of liquid motion. These effects can decrease the prefactor significantly, especially when the liquid pressure is much smaller than the equilibrium one. The deviation in the nucleation rate between the improved formula and the CNT can be as large as several orders of magnitude. Our improved, accurate prefactor and recent advances in molecular dynamics simulations and laboratory experiments for argon bubble nucleation enable us to precisely constrain the free energy barrier for bubble nucleation. Assuming the correction to the CNT free energy is of the functional form suggested by Tolman, the precise evaluations of the free energy barriers suggest the Tolman length is $\simeq 0.3 σ$ independently of the temperature for argon bubble nucleation, where $σ$ is the unit length of the Lenard-Jones potential. With this Tolman correction and our new prefactor one gets accurate bubble nucleation rate predictions in the parameter range probed by current experiments and molecular dynamics simulations.

preprint2015arXiv

Spinon, soliton and breather in the spin-1/2 antiferromagnetic chain KCuGaF$_6$

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.

preprint2015arXiv

Successive magnetic phase transitions in $\bmα$-RuCl$_3$: XY-like frustrated magnet on the honeycomb lattice

The layered compound $α$-RuCl$_3$ is composed of a honeycomb lattice of magnetic Ru$^{3+}$ ions with the $4d^5$ electronic state. We have investigated the magnetic properties of $α$-RuCl$_3$ via magnetization and specific heat measurements using single crystals. It was observed that $α$-RuCl$_3$ undergoes a structural phase transition at $T_{t}\,{\simeq}\,150$ K accompanied by fairly large hysteresis. This structural phase transition is expected to be similar to that observed in closely related CrCl$_3$. The magnetizations and magnetic susceptibilities are strongly anisotropic, which mainly arise from the anisotropic $g$-factors, i.e., $g_{ab}\,{\simeq}\,2.5$ and $g_c\,{\simeq}\,0.4$ for magnetic fields parallel and perpendicular to the $ab$ plane, respectively. These $g$-factors and the obtained entropy indicate that the effective spin of Ru$^{3+}$ is one-half, which results from the low-spin state. Specific heat data show that magnetic ordering occurs in four steps at zero magnetic field. The successive magnetic phase transitions should be ascribed to the competition among exchange interactions. The magnetic phase diagram for $H\,{\parallel}\,ab$ is obtained. We discuss the strongly anisotropic $g$-factors in $α$-RuCl$_3$ and deduce that the exchange interaction is strongly XY-like. $α$-RuCl$_3$ is magnetically described as a three-dimensionally coupled XY-like frustrated magnet on a honeycomb lattice.

preprint2014arXiv

Almost Perfect Frustration in the Dimer Magnet Ba$_2$CoSi$_2$O$_6$Cl$_2$

We determined the crystal structure of Ba$_2$CoSi$_2$O$_6$Cl$_2$, which was synthesized in this work, and investigated its quantum magnetic properties using single crystals. This compound should be described as a two-dimensionally coupled spin-1/2 XY-like spin dimer system. Ba$_2$CoSi$_2$O$_6$Cl$_2$ exhibits a stepwise magnetization process with a plateau at half of the saturation magnetization, irrespective of the field direction, although all the Co$^{2+}$ sites are equivalent. This indicates that spin triplets are localized owing to the almost perfect frustration of interdimer exchange interactions. Thus, the spin states for the zero and 1/2 magnetization-plateau states are almost exactly given by the simple product of singlet dimers and the alternate product of singlet and triplet dimers, respectively.

preprint2014arXiv

Bubble Evolution and Properties in Homogeneous Nucleation Simulations

We analyze the properties of naturally formed nano-bubbles in Lennard-Jones molecular dynamics simulations of liquid-to-vapor nucleation in the boiling and the cavitation regimes. The large computational volumes provide a realistic environment at unchanging average temperature and liquid pressure, which allows us to accurately measure properties of bubbles from their inception as stable, critically sized bubbles, to their continued growth into the constant speed regime. Bubble gas densities are up to 50$%$ lower than the equilibrium vapor densities at the liquid temperature, yet quite close to the gas equilibrium density at the lower gas temperatures measured in the simulations: The latent heat of transformation results in bubble gas temperatures up to 25$%$ below those of the surrounding bulk liquid. In the case of rapid bubble growth - typical for the cavitation regime - compression of the liquid outside the bubble leads to local temperature increases of up to 5$%$, likely significant enough to alter the surface tension as well as the local viscosity. The liquid-vapor bubble interface is thinner than expected from planar coexistence simulations by up to $50%$. Bubbles near the critical size are extremely non-spherical, yet they quickly become spherical as they grow.

preprint2014arXiv

Direct Simulations of Homogeneous Bubble Nucleation: Agreement with CNT and no Local Hot Spots

We present results from direct, large-scale molecular dynamics (MD) simulations of homogeneous bubble (liquid-to-vapor) nucleation. The simulations contain half a billion Lennard-Jones (LJ) atoms and cover up to 56 million time-steps. The unprecedented size of the simulated volumes allows us to resolve the nucleation and growth of many bubbles per run in simple direct micro-canonical (NVE) simulations while the ambient pressure and temperature remain almost perfectly constant. We find bubble nucleation rates which are lower than in most of the previous, smaller simulations. It is widely believed that classical nucleation theory (CNT) generally underestimates bubble nucleation rates by very large factors. However, our measured rates are within two orders of magnitude of CNT predictions - only at very low temperatures does CNT underestimate the nucleation rate significantly. Introducing a small, positive Tolman length leads to very good agreement at all temperatures, as found in our recent vapor-to-liquid nucleation simulations. The critical bubbles sizes derived with the nucleation theorem agree well with the CNT predictions at all temperatures. Local hot spots reported in the literature are not seen: Regions where a bubble nucleation events will occur are not above the average temperature, and no correlation of temperature fluctuations with subsequent bubble formation is seen.

preprint2014arXiv

Free energy of cluster formation and a new scaling relation for the nucleation rate

Recent very large molecular dynamics simulations of homogeneous nucleation with $(1-8) \cdot 10^9$ Lennard-Jones atoms [Diemand et al. J. Chem. Phys. {\bf 139}, 074309 (2013)] allow us to accurately determine the formation free energy of clusters over a wide range of cluster sizes. This is now possible because such large simulations allow for very precise measurements of the cluster size distribution in the steady state nucleation regime. The peaks of the free energy curves give critical cluster sizes, which agree well with independent estimates based on the nucleation theorem. Using these results, we derive an analytical formula and a new scaling relation for nucleation rates: $\ln J' / η$ is scaled by $\ln S / η$, where the supersaturation ratio is $S$, $η$ is the dimensionless surface energy, and $J'$ is a dimensionless nucleation rate. This relation can be derived using the free energy of cluster formation at equilibrium which corresponds to the surface energy required to form the vapor-liquid interface. At low temperatures (below the triple point), we find that the surface energy divided by that of the classical nucleation theory does not depend on temperature, which leads to the scaling relation and implies a constant, positive Tolman length equal to half of the mean inter-particle separation in the liquid phase.

preprint2014arXiv

Magnetic Phase Diagram of the S=1/2 Triangular-Lattice Heisenberg Antiferromagnet Ba3CoNb2O9

We report the results of low-temperature thermal and magnetic measurements on Ba$_3$CoNb$_2$O$_9$ powder, described as a uniform triangular-lattice antiferromagnet (TLAF) with a fictitious spin-1/2. Ba$_3$CoNb$_2$O$_9$ is found to undergo two-step antiferromagnetic transitions at $T_{\rm N1}=1.39~{\rm K}$ and $T_{\rm N2}=1.13~{\rm K}$. As the magnetic field is increased, both $T_{\rm N1}$ and $T_{\rm N2}$ monotonically decrease. The magnetic field vs temperature phase diagram indicates that the exchange interactions are nearly of the Heisenberg type with weak easy-axis anisotropy and that the exchange interaction between triangular lattices is crucial, in contrast to the case of the quasi-two-dimensional TLAF Ba$_3$CoSb$_2$O$_9$ [Susuki {\it et al}., Phys. Rev. Lett. \textbf{110}, 267201 (2013)].

preprint2014arXiv

Molecular Dynamics Simulations of the Nucleation of Water: Determining the Sticking Probability and Formation Energy of a Cluster

We performed molecular dynamics (MD) simulations of the nucleation of water vapor in order to test nucleation theories. Simulations were performed for a wide range of supersaturation ratios (S = 3-25) and water temperatures (Tw=300-390K). We obtained the nucleation rates and the formation free energies of a subcritical cluster from the cluster size distribution. The classical nucleation theory (CNT) and the modified classical nucleation theory (MCNT) overestimate the nucleation rates in all cases. The semi-phenomenological (SP) model, which corrects the MCNT prediction using the second virial coefficient of a vapor, reproduces the formation free energy of a cluster with the size < 20 to within 10 % and the nucleation rate and cluster size distributions to within one order of magnitude. The sticking probability of the vapor molecules to the clusters was also determined from the growth rates of the clusters. The sticking probability rapidly increases with the supersaturation ratio S, which is similar to the Lennard-Jones system.

preprint2014arXiv

Non-$γ_{5}$hermiticity minimal doubling fermion

We formulate new two dimensional fermions breaking $γ_{5}$hermiticity, based on the minimal doubling fermion. We investigate their properties: (I) Symmetries, (II) eigenvalue distributions, and (III)the number of poles. As a simple application of the fermions, the Gross-Neveu model in two dimensions is studied using the fermion. We obtain the parity broken phase diagram called Aoki phase and the chiral broken phase diagram of the model with an imaginary chemical potential.

preprint2014arXiv

Opacity of fluffy dust aggregates

Context. Dust grains coagulate to form dust aggregates in protoplanetary disks. Their porosity can be extremely high in the disks. Although disk emission may come from fluffy dust aggregates, the emission has been modeled with compact grains. Aims. We aim to reveal the mass opacity of fluffy aggregates from infrared to millimeter wavelengths with the filling factor ranging from 1 down to $10^{-4}$. Methods. We use Mie calculations with an effective medium theory. The monomers are assumed to be 0.1 ${\rm μm}$ sized grains, which is much shorter than the wavelengths that we focus on. Results. We find that the absorption mass opacity of fluffy aggregates are characterized by the product $a\times f$, where $a$ is the dust radius and $f$ is the filling factor, except for the interference structure. The scattering mass opacity is also characterized by $af$ at short wavelengths while it is higher in more fluffy aggregates at long wavelengths. We also derive the analytic formula of the mass opacity and find that it reproduces the Mie calculations. We also calculate the expected difference of the emission between compact and fluffy aggregates in protoplanetary disks with a simple dust growth and drift model. We find that compact grains and fluffy aggregates can be distinguished by the radial distribution of the opacity index $β$. The previous observation of the radial distribution of $β$ is consistent with the fluffy case, but more observations are required to distinguish between fluffy or compact. In addition, we find that the scattered light would be another way to distinguish between compact grains and fluffy aggregates.

preprint2014arXiv

Properties of Liquid Clusters in Large-scale Molecular Dynamics Nucleation Simulations

We have performed large-scale Lennard-Jones molecular dynamics simulations of homogeneous vapor-to-liquid nucleation, with $10^9$ atoms. This large number allows us to resolve extremely low nucleation rates, and also provides excellent statistics for cluster properties over a wide range of cluster sizes. The nucleation rates, cluster growth rates, and size distributions are presented in Diemand et al. [J. Chem. Phys. {\bf 139}, 74309 (2013)], while this paper analyses the properties of the clusters. We explore the cluster temperatures, density profiles, potential energies and shapes. A thorough understanding of the properties of the clusters is crucial to the formulation of nucleation models. Significant latent heat is retained by stable clusters, by as much as $ΔkT = 0.1 ε$ for clusters with size $i = 100$. We find that the clusters deviate remarkably from spherical - with ellipsoidal axis ratios for critical cluster sizes typically within $b/c = 0.7\pm 0.05$ and $a/c = 0.5 \pm 0.05$. We examine cluster spin angular momentum, and find that it plays a negligible role in the cluster dynamics. The interfaces of large, stable clusters are thiner than planar equilibrium interfaces by $10-30\%$. At the critical cluster size, the cluster central densities are between $5-30\%$ lower than the bulk liquid expectations. These lower densities imply larger-than-expected surface areas, which increase the energy cost to form a surface, which lowers nucleation rates.

preprint2014arXiv

Strong Suppression of Magnetic Ordering in an S = 1/2 Square-Lattice Heisenberg Antiferromagnet Sr2CuTeO6

We report the magnetic susceptibility and the low-temperature specific heat of a double perovskite compound, Sr2CuTeO6, which is expected as the spin-1/2 square-lattice Heisenberg antiferromagnet with nearest-neighbor J1 and next-nearest-neighbor J2 exchange interactions. The specific heat exhibits a sharp lambda-like anomaly at TN = 4.8 K indicative of magnetic ordering, which is much lower than those in isostructural A2CuMO6 with A = Ba, Sr and M = W, Mo. This indicates the strong suppression of magnetic ordering in Sr2CuTeO6. We estimate the exchange constant J1/kB ~ 80 K and J2/J1 < 0.07 from magnetic susceptibility data using Pade approximation of high temperature expansion and quantum Monte Carlo method. The strong suppression of magnetic ordering in Sr2CuTeO6 should be ascribed to the weekness of the interlayer exchange.

preprint2014arXiv

The Fermi Function Beyond The Second Order Perturbation Theory

The Fermi function is historically derived from the Dirac equation or the Schrödinger equation. However, we claim that the Fermi function should be derived from quantum field theory. Then, we obtain the following results: (1) We give the electromagnetic correction of the beta decay to the fourth order in $α/v$, where $v$ is the electron velocity. (2) Our result is consistent with the ordinary non-relativistic Fermi function. (3) We give the iterated integral representation of the Fermi function for arbitrary order in $α/v$. (4) This representation is related to the multiple L values, and they enable us to calculate the amplitude in the algebraic approach.

preprint2014arXiv

The multifaceted planetesimal formation process

Accumulation of dust and ice particles into planetesimals is an important step in the planet formation process. Planetesimals are the seeds of both terrestrial planets and the solid cores of gas and ice giants forming by core accretion. Left-over planetesimals in the form of asteroids, trans-Neptunian objects and comets provide a unique record of the physical conditions in the solar nebula. Debris from planetesimal collisions around other stars signposts that the planetesimal formation process, and hence planet formation, is ubiquitous in the Galaxy. The planetesimal formation stage extends from micrometer-sized dust and ice to bodies which can undergo run-away accretion. The latter ranges in size from 1 km to 1000 km, dependent on the planetesimal eccentricity excited by turbulent gas density fluctuations. Particles face many barriers during this growth, arising mainly from inefficient sticking, fragmentation and radial drift. Two promising growth pathways are mass transfer, where small aggregates transfer up to 50% of their mass in high-speed collisions with much larger targets, and fluffy growth, where aggregate cross sections and sticking probabilities are enhanced by a low internal density. A wide range of particle sizes, from mm to 10 m, concentrate in the turbulent gas flow. Overdense filaments fragment gravitationally into bound particle clumps, with most mass entering planetesimals of contracted radii from 100 to 500 km, depending on local disc properties. We propose a hybrid model for planetesimal formation where particle growth starts unaided by self-gravity but later proceeds inside gravitationally collapsing pebble clumps to form planetesimals with a wide range of sizes.

preprint2014arXiv

TMD Parton Distributions based on Three-Body Decay Functions in NLL Order of QCD

Three-body decay functions in space-like parton branches are implemented to evaluate transverse-momentum-dependent (TMD) parton distribution functions in the next-to-leading logarithmic (NLL) order of quantum chromodynamics (QCD). Interference contributions due to the next-to-leading order contribution are taken into account for the evaluation of the transverse momenta in initial state parton radiations. Some properties of the decay functions are also examined. As an example, we compare our results with an algorithm proposed in Ref.1), in which a transverse momentum distributions are evaluated at the last step of parton evolutions.

preprint2013arXiv

Constraints on the charged scalar effects using the forward-backward asymmetry on $B\to D^{(*)}τ\bar{ν_τ}$

The decay modes $\bar{B}\to D^{(*)}τ\barν_τ$ are sensitive to charged scalar effects, such as the charged Higgs effects. In this paper we suggest a method to determine their effects by using the ratio of branching fractions and forward-backward asymmetries. In particular, forward-backward asymmetries on $\bar{B}\to D^{(*)}τ(\to πν_τ)\barν_τ$, $\bar{B}\to D^{(*)}τ(\to ρν_τ)\barν_τ$, and $\bar{B}\to D^{(*)}τ(\to a_1ν_τ)\barν_τ$ play an important role, which discriminate the Standard Model from other New Physics scenarios.

preprint2013arXiv

Fluffy dust forms icy planetesimals by static compression

Context: In planetesimal formation theory, several barriers have been proposed, which are bouncing, fragmentation, and radial drift problems. To understand the structure evolution of dust aggregates is a key in the planetesimal formation. Dust grains become fluffy by coagulation in protoplanetary disks. However, once they become fluffy, they are not sufficiently compressed by collisional compression to form compact planetesimals. Aims: We aim to reveal the pathway of the dust structure evolution from dust grains to compact planetesimals. Methods: Using the compressive strength formula, we analytically investigate how fluffy dust aggregates are compressed by static compression due to ram pressure of the disk gas and self gravity of the aggregates in protoplanetary disks. Results: We reveal the pathway of the porosity evolution from dust grains via fluffy aggregates to form planetesimals, circumventing the barriers in planetesimal formation. The aggregates are compressed by the disk gas to the density of 10^{-3} g/cm^3 in coagulation, which is more compact than the case with collisional compression. Then, they are compressed more by self gravity to 10^{-1} g/cm^3 when the radius is 10 km. Although the gas compression decelerate the growth, they grow enough rapidly to avoid the radial drift barrier when the orbital radius is < 6 AU in a typical disk. Conclusions: We propose fluffy dust growth scenario from grains to planetesimals. It enables the icy planetesimal formation in a wide range beyond the snowline in protoplanetary disks. This result proposes a concrete initial condition of planetesimals for the later stages of the planet formation.

preprint2013arXiv

Ground state of the spin-1/2 chain of green dioptase at high fields

The gem-stone dioptase Cu6Si6O18.6H2O has a chiral crystal structure of equilateral triangular helices consisting of Cu-3d spins. It shows an antiferromagnetic order with an easy axis along c at TN = 15.5 K under zero field, and a magnetization jump at HC = 13.5 T when the field is applied along c-axis. By 29Si-NMR measurements, we have revealed that the high-field state is essentially the two sub-lattice structure, and that the component within ab-plane is collinear. The result indicates no apparent match with the geometrical pattern of helical spin chain.

preprint2013arXiv

Large Scale Molecular Dynamics Simulations of Homogeneous Nucleation

We present results from large-scale molecular dynamics (MD) simulations of homogeneous vapor-to-liquid nucleation. The simulations contain between one and eight billion Lennard-Jones (LJ) atoms, covering up to 1.2 μs (56 million time-steps). They cover a wide range of supersaturation ratios, S=1.55 to 10^4, and temperatures from kT = 0.3 to 1.0 ε (where ε is the depth of the LJ potential, and k the Boltzmann constant). We have resolved nucleation rates as low as 10^{17} cm^{-3} s^{-1} (in the argon system), and critical cluster sizes as large as 100 atoms. Recent argon nucleation experiments probe nucleation rates in an overlapping range, making the first direct comparison between laboratory experiments and molecular dynamics simulations possible: We find very good agreement within the uncertainties, which are mainly due to the extrapolations of argon and LJ saturation curves to very low temperatures. The self-consistent, modified classical nucleation model of Girshick and Chiu [J. Chem. Phys. 93, 1273 (1990)] underestimates the nucleation rates by up to 9 orders of magnitudes at low temperatures, and at kT = 1.0 ε it overestimates them by up to 10^5. The predictions from a semi-phenomenological model by Laaksonen et al. [Phys. Rev. E 49, 5517 (1994)] are much closer to our MD results, but still differ by factors of up to 104 in some cases. At low temperatures, the classical theory predicts critical clusters sizes, which match the simulation results (using the first nucleation theorem) quite well, while the semi-phenomenological model slightly underestimates them. At kT = 1.0 ε the critical sizes from both models are clearly too small. (abridged)

preprint2013arXiv

Magnetization Process and Collective Excitations in the S=1/2 Triangular-Lattice Heisenberg Antiferromagnet Ba3CoSb2O9

We have performed high-field magnetization and ESR measurements on Ba$_3$CoSb$_2$O$_9$ single crystals, as a representative two-dimensional spin-1/2 Heisenberg antiferromagnet on a uniform triangular lattice, and have determined all the magnetic parameters. For an applied magnetic field $H$ parallel to the $ab$-plane, the entire magnetization curve including the plateau at one-third of the saturation magnetization ($M_\mathrm{s}$) is in excellent agreement with the results of theoretical calculations except an anomaly near $(3/5)M_\mathrm{s}$, indicative of an additional field-induced transition. However, for $H\,{\parallel}\,c$, the magnetization curve exhibits a cusp near $M_\mathrm{s}/3$. Paramagnetic resonance signals above the Néel temperature indicate the nearly isotropic $g$-factor. A detailed analysis of the collective ESR modes observed in the ordered state, combined with the magnetization process, provides evidence of the weak easy-plane anisotropy.

preprint2013arXiv

Minimal Doubling Fermion and Hermiticity

We analyze the lattice fermion kinetic term using PT symmetry, R-hermiticity, and $γ_{5}$-hermiticity. R-hermiticity is a condition for Hermite action and it is related to $γ_{5}$-hermiticity and PT symmetry. Assuming that a translation-invariant kinetic term with continuum and periodic function does not have PT symmetry, it can have R-hermiticity or $γ_{5}$-hermiticity. We prove that a kinetic term with continuum and periodic function that is PT symmetric does not reduce doublers. As a simple example, we analyze the two-dimensional two-flavor Gross-Neveu model with minimal doubling fermions. The minimal doubling fermions break PT symmetry and R-hermiticity, hence complex or non-Hermite coupling constants are caused by quantum correction.

preprint2013arXiv

Non-$γ_{5}$hermiticity fermions in two dimensions

We construct 2D non-$γ_{5}$hermiticity fermions based on the minimal doubling fermion. We investigate symmetries, reflection positivity, eigenvalue distribution and the number of poles for our fermions. As simple tests for application to the fermion, the Gross-Neveu model in two dimensions is studied using the non-$γ_{5}$hermiticity fermion. We draw the parity broken phase diagram, called Aoki phase and the chiral broken phase diagram for the model with an imaginary chemical potential.

preprint2013arXiv

Static compression of porous dust aggregates

Context: In protoplanetary disks, dust grains coagulate with each other and grow to form aggregates. As these aggregates grow by coagulation, their filling factor ϕdecreases down to ϕ<< 1. However, comets, the remnants of these early planetesimals, have ϕ~ 0.1. Thus, static compression of porous dust aggregates is important in planetesimal formation. However, the static compression strength has been investigated only for relatively high density aggregates (ϕ> 0.1). Aims: We investigate and find the compression strength of highly porous aggregates (ϕ<< 1). Methods: We perform three dimensional N-body simulations of aggregate compression with a particle-particle interaction model. We introduce a new method of static compression: the periodic boundary condition is adopted and the boundaries move with low speed to get closer. The dust aggregate is compressed uniformly and isotropically by themselves over the periodic boundaries. Results: We empirically derive a formula of the compression strength of highly porous aggregates (ϕ<< 1). We check the validity of the compression strength formula for wide ranges of numerical parameters, such as the size of initial aggregates, the boundary speed, the normal damping force, and material. We also compare our results to the previous studies of static compression in the relatively high density region (ϕ> 0.1) and confirm that our results consistently connect to those in the high density region. The compression strength formula is also derived analytically.

preprint2012arXiv

A Derivation of the Fermi Function in Perturbative Quantum Field Theory

We postulate that the Fermi function should be derived from the amplitude, not from the solution of the Dirac equation, in the quantum field theory. Then, we obtain the following results. 1, We give the amplitude and the width of the neutron beta decay, $n \to p + e^- + \bar ν_e $ to one loop order. It is carried out by the Feynman parameter integration. 2, As the result, we find the terms which can be interpreted as the Fermi function expanded to order $α$. 3, We also give the same result using complex analysis. 4, We check that there are no such terms in the similar process, $\bar ν_e + p \to e^+ + n$. 5, We perform the Fermi function expanded to order $α^2$ using complex analysis.

preprint2012arXiv

Ferromagnetism induced in anisotropic stacked kagome-lattice antiferromagnet Cs$_2$Cu$_3$CeF$_{12}$

The magnetic properties of Cs$_2$Cu$_3$CeF$_{12}$ were investigated through magnetization and specific heat measurements. Cs$_2$Cu$_3$CeF$_{12}$ is composed of a buckled kagome lattice of Cu$^{2+}$, which is stacked along the b axis. The exchange network in the buckled kagome lattice is strongly anisotropic. Consequently, Cs$_2$Cu$_3$CeF$_{12}$ can be divided into two subsystems: alternating Heisenberg chains with strong antiferromagnetic exchange interactions and dangling spins. The dangling spins couple with one another via effective exchange interactions, which are mediated by chain spins. The dangling spins are further divided into two subsystems, DS1 and DS2. The dangling spins in DS1 undergo three-dimensional ferromagnetic ordering at 3.14 K, while those in DS2 remain paramagnetic down to 0.35 K. The effective interaction between the DS1 spins is approximately expressed by the ferromagnetic $XXZ$ model with the $z$ direction parallel to the crystallographic c axis. A magnetic phase diagram for $H {\parallel} c$ was obtained and was analyzed within the framework of the molecular field approximation. With increasing magnetic field, the dangling spins are polarized and the magnetization curve exhibits a wide plateau at one-third of the saturation magnetization.

preprint2012arXiv

Geometrical Cross Sections of Dust Aggregates and a Compression Model for Aggregate Collisions

Geometrical cross sections of dust aggregates determine their coupling with disk gas, which governs their motions in protoplanetary disks. Collisional outcomes also depend on geometrical cross sections of initial aggregates. In the previous paper, we performed three-dimensional N-body simulations of sequential collisions of aggregates composed of a number of sub-micron-sized icy particles and examined radii of gyration (and bulk densities) of the obtained aggregates. We showed that collisional compression of aggregates is not efficient and that aggregates remain fluffy. In the present study, we examine geometrical cross sections of the aggregates. Their cross sections decreases due to the compression as well as their gyration radii. It is found that a relation between the cross section and the gyration radius proposed by Okuzumi et al. is valid for the compressed aggregates. We also refine the compression model proposed in our previous paper. The refined model enables us to calculate the evolution of both gyration radii and cross sections of growing aggregates and reproduces well our numerical results of sequential aggregate collisions. The refined model can describe non-equal-mass collisions as well as equal-mass case. Although we do not take into account oblique collisions in the present study, oblique collisions would further hinder compression of aggregates.

preprint2012arXiv

Migration rates of planets due to scattering of planetesimals

Planets migrate due to the recoil they experience from scattering solid (planetesimal) bodies. To first order, the torques exerted by the interior and exterior disks cancel, analogous to the cancellation of the torques from the gravitational interaction with the gas (type I migration). Assuming the dispersion-dominated regime and power-laws characterized by indices α and β for the surface density and eccentricity profiles, we calculate the net torque on the planet. We consider both distant encounters and close (orbit-crossing) encounters. We find that the close and distant encounter torques have opposite signs with respect to their α and β dependences; and that the torque is especially sensitive to the eccentricity gradient (β). Compared to type-I migration due to excitation of density waves, the planetesimal-driven migration rate is generally lower due to the lower surface density of solids in gas-rich disk, although this may be partially or fully offset when their eccentricity and inclination are small. Allowing for the feedback of the planet on the planetesimal disk through viscous stirring, we find that under certain conditions a self-regulated migration scenario emerges, in which the planet migrates at a steady pace that approaches the rate corresponding to the one-sided torque. If the local planetesimal disk mass to planet mass ratio is low, however, migration stalls. We quantify the boundaries separating the three migration regimes.

preprint2012arXiv

Rapid Coagulation of Porous Dust Aggregates Outside the Snow Line: A Pathway to Successful Icy Planetesimal Formation

Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles' aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.

preprint2011arXiv

Determination of the Higgs CP property in Hadron Colliders

We propose three ways to determine the CP eigenvalue of the Higgs boson at the hadron collider as follows: 1. We determine the Higgs CP eigenvalue from the production cross section which is affected by the CP eigenvalue of the Higgs boson. 2. We adopt the CP selection rules to determine the Higgs CP eigenvalue. 3. We determine the CP property by the momentum distribution of the decay products of the Higgs boson. Our methods can be applied for a wide range of the Higgs mass.

preprint2011arXiv

Electrostatic Barrier against Dust Growth in Protoplanetary Disks. I. Classifying the Evolution of Size Distribution

Collisional growth of submicron-sized dust grains into macroscopic aggregates is the first step of planet formation in protoplanetary disks. These grains are expected to carry nonzero negative charges in the weakly ionized disks, but its effect on their collisional growth has not been fully understood so far. In this paper, we investigate how the charging affects the evolution of the dust size distribution properly taking into account the charging mechanism in a weakly ionized gas as well as porosity evolution through low-energy collisions. To clarify the role of the size distribution, we divide our analysis into two steps. First, we analyze the collisional growth of charged aggregates assuming a monodisperse (i.e., narrow) size distribution. We show that the monodisperse growth stalls due to the electrostatic repulsion when a certain condition is met, as is already expected in the previous work. Second, we numerically simulate dust coagulation using Smoluchowski's method to see how the outcome changes when the size distribution is allowed to freely evolve. We find that, under certain conditions, the dust undergoes bimodal growth where only a limited number of aggregates continue to grow carrying the major part of the dust mass in the system. This occurs because remaining small aggregates efficiently sweep up free electrons to prevent the larger aggregates from being strongly charged. We obtain a set of simple criteria that allows us to predict how the size distribution evolves for a given condition. In Paper II (arXiv:1009.3101), we apply these criteria to dust growth in protoplanetary disks.

preprint2011arXiv

Electrostatic Barrier against Dust Growth in Protoplanetary Disks. II. Measuring the Size of the "Frozen" Zone

Coagulation of submicron-sized dust grains into porous aggregates is the initial step of dust evolution in protoplanetary disks. Recently, it has been pointed out that negative charging of dust in the weakly ionized disks could significantly slow down the coagulation process. In this paper, we apply the growth criteria obtained in Paper I to finding out a location ("frozen" zone) where the charging stalls dust growth at the fractal growth stage. For low-turbulence disks, we find that the frozen zone can cover the major part of the disks at a few to 100 AU from the central star. The maximum mass of the aggregates is approximately 10^{-7} g at 1 AU and as small as a few monomer masses at 100 AU. Strong turbulence can significantly reduce the size of the frozen zone, but such turbulence will cause the fragmentation of macroscopic aggregates at later stages. We examine a possibility that complete freezeout of dust evolution in low-turbulence disks could be prevented by global transport of dust in the disks. Our simple estimation shows that global dust transport can lead to the supply of macroscopic aggregates and the removal of frozen aggregates on a timescale of 10^6 yr. This overturns the usual understanding that tiny dust particles get depleted on much shorter timescales unless collisional fragmentation is effective. The frozen zone together with global dust transport might explain "slow" (\sim 10^6 yr) dust evolution suggested by infrared observation of T Tauri stars and by radioactive dating of chondrites.

preprint2011arXiv

Experimental Realization of Spin-1/2 Triangular-Lattice Heisenberg Antiferromagnet

We report the results of magnetization and specific heat measurements on Ba$_3$CoSb$_2$O$_9$, in which the magnetic Co$^{2+}$ ion has a fictitious spin-1/2, and show evidence that a spin-1/2 Heisenberg antiferromagnet on a regular triangular lattice is actually realized in Ba$_3$CoSb$_2$O$_9$. We found that the entire magnetization curve including the one-third quantum magnetization plateau is in excellent agreement with theoretical calculations at a quantitative level.

preprint2011arXiv

NLL Order Contributions for Exclusive Processes in Jet-Calculus Scheme

We investigate the next-to-leading logarithmic (NLL) order contributions of the quantum chromodynamics (QCD) for exclusive processes evaluated by Monte Carlo methods. Ambiguities of the Monte Carlo calculation based on the leading-logarithmic (LL) order approximations are pointed out. To remove these ambiguities, we take into account the NLL order terms. In a model presented in this paper, interference contributions due to the NLL order terms are included for the generation of the transverse momenta in initial-state parton radiations. Furthermore, a kinematical constraint due to parton radiation, which is also a part of the NLL order contributions, is taken into account. This method guarantees a proper phase space boundary for hard scattering cross sections as well as parton radiations. As an example, cross sections for lepton pair productions mediated by a virtual photon in hadron-hadron collisions are calculated, using the jet-calculus scheme for flavor nonsinglet quarks.

preprint2011arXiv

Planetary Core Formation with Collisional Fragmentation and Atmosphere to Form Gas Giant Planets

Massive planetary cores ($\sim 10$ Earth masses) trigger rapid gas accretion to form gas giant planets \rev{such as} Jupiter and Saturn. We investigate the core growth and the possibilities for cores to reach such a critical core mass. At the late stage, planetary cores grow through collisions with small planetesimals. Collisional fragmentation of planetesimals, which is induced by gravitational interaction with planetary cores, reduces the amount of planetesimals surrounding them, and thus the final core masses. Starting from small planetesimals that the fragmentation rapidly removes, less massive cores are formed. However, planetary cores acquire atmospheres that enlarge their collisional cross section before rapid gas accretion. Once planetary cores exceed about Mars mass, atmospheres significantly accelerate the growth of cores. We show that, taking into account the effects of fragmentation and atmosphere, initially large planetesimals enable formation of sufficiently massive cores. On the other hand, because the growth of cores is slow for large planetesimals, a massive disk is necessary for cores to grow enough within a disk lifetime. If the disk with 100\,km-sized initial planetesimals is 10 times as massive as the minimum mass solar nebula, planetary cores can exceed 10 Earth masses in the Jovian planet region ($>5\,$AU).

preprint2011arXiv

Quantum Magnetization Plateau in Spin-1 Triangular-Lattice Antiferromagnet Ba$_3$NiSb$_2$O$_9$

We report the results of magnetization and specific heat measurements on Ba$_3$NiSb$_2$O$_9$, which is a quasi-two-dimensional spin-1 triangular-lattice antiferromagnet. We observed a nonclassical magnetization plateau at one-third of the saturation magnetization that is driven by spin frustration and quantum fluctuation. Exact diagonalization for a 21-site rhombic cluster was performed to analyze the magnetization process. Experimental and calculated results agree well.

preprint2011arXiv

Reply to "Comment on `Transition from Bose glass to a condensate of triplons in Tl$_{1-x}$K$_x$CuCl$_3$'"

Showing low-temperature specific heat and other experimental data and also on the basis of established physics, we argue against the comment made by Zheludev and Hüvonnen criticizing our recent study on the magnetic-field-induced spin ordering and critical behavior in Tl$_{1-x}$K$_x$CuCl$_3$, which is described as the Bose glass-condensate transition of triplons.

preprint2011arXiv

Thermodynamic properties of quantum sine-Gordon spin chain system KCuGaF6

We investigated the thermodynamic properties of the spin-1/2 one-dimensional Heisenberg antiferromagnet KCuGaF6 by measuring the specific heat in magnetic fields. When this compound is subjected to a uniform magnetic field H a transverse staggered magnetic field h is induced in this compound owing to the staggered component of the g tensor and the Dzyaloshinskii-Moriya interaction with an alternating D vector. Consequently, the quantum sine-Gordon (SG) model is an effective model of this compound in a uniform magnetic field. In three different field directions, we observed a magnetic-field-induced gap, which increases with H. We analyzed experimental results using specific heat theory based on quantum SG theory. The thermodynamic property for H // c is very well described in terms of the elementary excitations characteristic of the quantum SG model, while for the other field directions, significant contributions from other excitation modes beyond the framework of the quantum SG model were observed. For H // b, a quantum phase transition between gapless and gapped ground states was observed.

preprint2010arXiv

A New Proposal for Neutrino Mass and |V(ud)| Measurements

We introduce a new method to detect the absolute neutrino mass scale. It uses a macroscopic mass of tritium source. We explain that the neutrino mass can be measured by scaling the mass difference of the source between initial and final state, and its heat value. This method is free from the electron energy resolution limit and the statistical error. We estimate the required accuracy to measure the neutrino mass. We also report that the {u,d} component of the CKM matrix, |V(ud)| may be determined in 10^-6 accuracy as an application of this work.

preprint2010arXiv

Numerical Calculation of Schwinger-Dyson Equation with Momentum-Dependent Gauge Parameter at Finite Temperature

Chiral symmetry at finite temperature is studied using the Schwinger-Dyson equation. We calculate numerically the critical temperature using the Schwinger-Dyson equation with the gauge parameter that depends on an external momentum. The critical temperature obtained by this method is similar to that with the Landau gauge and wave function renormalization constant 1. Moreover, the gauge invariance in the ladder approximation is examined using our method.

preprint2010arXiv

Planetary Growth with Collisional Fragmentation and Gas Drag

As planetary embryos grow, gravitational stirring of planetesimals by embryos strongly enhances random velocities of planetesimals and makes collisions between planetesimals destructive. The resulting fragments are ground down by successive collisions. Eventually the smallest fragments are removed by the inward drift due to gas drag. Therefore, the collisional disruption depletes the planetesimal disk and inhibits embryo growth. We provide analytical formulae for the final masses of planetary embryos, taking into account planetesimal depletion due to collisional disruption. Furthermore, we perform the statistical simulations for embryo growth (which excellently reproduce results of direct $N$-body simulations if disruption is neglected). These analytical formulae are consistent with the outcome of our statistical simulations. Our results indicate that the final embryo mass at several AU in the minimum-mass solar nebula can reach about $\sim 0.1$ Earth mass within $10^7$ years. This brings another difficulty in formation of gas giant planets, which requires cores with $\sim 10$ Earth masses for gas accretion. However, if the nebular disk is 10 times more massive than the minimum-mass solar nebula and the initial planetesimal size is larger than 100 km, as suggested by some models of planetesimal formation, the final embryo mass reaches about 10 Earth masses at 3-4 AU. The enhancement of embryos' collisional cross sections by their atmosphere could further increase their final mass to form gas giant planets at 5-10 AU in the solar system.

preprint2010arXiv

Pressure-induced new magnetic phase in Tl(Cu$_{0.985}$Mg$_{0.015}$)Cl$_3$ probed by muon spin rotation

We carried out zero-field muon-spin-rotation (ZF-$μ$SR) measurements in hydrostatic pressures in impurity-doped quantum spin gap system Tl(Cu$_{0.985}$Mg$_{0.015}$)Cl$_3$ to investigate microscopic magnetic properties of the pressure-induced phase. The spontaneous muon spin precession, which indicates the existence of a long-range coherent order, is observed in pressures. With decreasing temperature in 3.1 kbar, the internal static magnetic field at the muon sites {\it H}$_{\rm int}$ tends to saturate to 280 Oe around 4 K, however, decreases to 240 Oe at 2.3 K. These results suggest a rearrangement of ordered spins, and we speculate that the oblique antiferromagnetic phase, which is observed in the pressure of 14 kbar on the pure TlCuCl$_3$, appears in the Mg-doped system in lower pressures.

preprint2009arXiv

Elementary excitations of S=1/2 one-dimensional antiferromagnet KCuGaF6 in magnetic field and quantum sine-Gordon model

Elementary excitations of the S=1/2 one-dimensional Heisenberg antiferromagnet KCuGaF$_6$ with exchange constant $J/k_B$=103 K were investigated by high-frequency ESR measurements combined with a pulsed high magnetic field. When an external magnetic field H is applied in KCuGaF$_6$, a staggered magnetic field h is induced perpendicular to H owing to the staggered g tensor and the Dzyaloshinsky-Moriya (DM) interaction with an alternating D vector. Consequently, KCuGaF$_6$ in a magnetic field is represented by the quantum sine-Gordon (SG) model. We observed many resonance modes including a soliton resonance, breathers, interbreather transitions and two-breather excitation. Their resonance conditions are beautifully described by the quantum SG field theory with one adjustable parameter $c_ s=h/H$. To investigate the relationship between the Curie term due to the DM interaction and the proportional coefficient $c_s$, magnetic susceptibility measurements were also performed varying the external field direction.

preprint2009arXiv

Numerical Modeling of the Coagulation and Porosity Evolution of Dust Aggregates

Porosity evolution of dust aggregates is crucial in understanding dust evolution in protoplanetary disks. In this study, we present useful tools to study the coagulation and porosity evolution of dust aggregates. First, we present a new numerical method for simulating dust coagulation and porosity evolution as an extension of the conventional Smoluchowski equation. This method follows the evolution of the mean porosity for each aggregate mass simultaneously with the evolution of the mass distribution function. This method reproduces the results of previous Monte Carlo simulations with much less computational expense. Second, we propose a new collision model for porous dust aggregates on the basis of our N-body experiments on aggregate collisions. We first obtain empirical data on porosity changes between the classical limits of ballistic cluster-cluster and particle-cluster aggregation. Using the data, we construct a recipe for the porosity change due to general hit-and-stick collisions as well as formulae for the aerodynamical and collisional cross sections. Simple coagulation simulations using the extended Smoluchowski method show that our collision model explains the fractal dimensions of porous aggregates observed in a full N-body simulation and a laboratory experiment. Besides, we discover that aggregates at the high-mass end of the distribution can have a considerably small aerodynamical cross section per unit mass compared with aggregates of lower masses. We point out an important implication of this discovery for dust growth in protoplanetary disks.

preprint2007arXiv

Magnetic-Field Induced Bose-Einstein Condensation of Magnons and Critical Behavior in Interacting Spin Dimer System TlCuCl$_3$

Magnetization measurements were performed to investigate the critical behavior of the field-induced magnetic ordering in gapped spin system TlCuCl$_3$. The critical density of the magnons was obtained as a function of temperature and the magnon-magnon interaction constant was evaluated. The experimental phase boundary for $T < 5$ K agrees almost perfectly with the magnon BEC theory based on the Hartree-Fock approximation with realistic dispersion relations. The phase boundary can be described by the power law $[H_{N}(T)-H_{c}] propto T^{phi}$. With decreasing fitting temperature range, the critical exponent ${phi}$ decreases and converges at $phi_{BEC} =3/2$ predicted by the magnon BEC theory.

preprint2005arXiv

Structural Transitions and Magnetic Structure in NH4CuCl3 via 14N-NMR

We report results of 14N-NMR experiments on NH4CuCl3 at the magnetic field of 7 T, where the 1/4-magnetization plateau is observed at low temperatures. The quadrupole splitting parameter $ν_{z}$ splits below 70 K, indicating a structural phase transition. At 4.2 K, eight N sites with distinct values of both $ν_{z}$ and the magnetic hyperfine shift $K_{z}$ are resolved in the NMR spectrum for general field directions. We then conlude that the magnetic structure in the 1/4-plateau does not break the symmetry of the crystal. Based on the NMR and the recent neutron scattering results by Ruegg et al. [Phys. Rev. Lett. 93 (2004) 037207], we propose that triplet dimers in the 1/4-plateau is formed not between the nearest neighbor pairs but over different chains.

preprint1995arXiv

The SUSY-GRACE system

We introduce a new method to treat Majorana fermions on the GRACE system which has been developed for the computation of the matrix elements for the processes of the standard model. In the standard model, we already have such particles as Dirac fermions, gauge bosons and scalar bosons in the system. On the other hand, in the SUSY models there are Majorana fermions. In the first instance, we have constructed a system for the automatic computation of cross-sections for the processes of the SUSY QED. It is remarkable that our system is also applicable to another model including Majorana fermions (e.g. MSSM) once the definition of the model file is given.

Hidekazu Tanaka

What is connected

Connect this record

See the researcher in context

Building this map preview

75 published item(s)

Contrasting Magnetic Structures in SrLaCuSbO$_{6}$ and SrLaCuNbO$_{6}$: Spin-1/2 Quasi-Square-Lattice $J_1 {-} J_2$ Heisenberg Antiferromagnets

Impacts of viscous dissipation on collisional growth and fragmentation of dust aggregates

Structures of magnetic excitations in the spin-1/2 kagome-lattice antiferromagnets Cs$_2$Cu$_3$SnF$_{12}$ and Rb$_2$Cu$_3$SnF$_{12}$

Can stellar-mass black hole growth disrupt disks of active galactic nuclei? The role of mechanical feedback

Eccentric gap induced by a super-Jupiter mass planet

Magnons and Spinons in $\mathrm{Ba}_2\mathrm{CoTeO}_6 $: A Composite System of Isolated Spin-$1/2$ Triangular Heisenberg-like and Frustrated Honeycomb Ising-like Antiferromagnets

An Experimental and Theoretical Insights into the Dielectric Properties of (Li, Nd) Co-doped ZnO Ceramics

Comments on "Type II migration strikes back -- An old paradigm for planet migration in discs" by Scardoni et al

Coupled spin-1/2 antiferromagnetic chain Cs$_2$LiRuCl$_6$ with partially disordered crystal lattice

Final Masses of Giant Planets III: Effect of Photoevaporation and a New Planetary Migration Model

Ring formation by coagulation of dust aggregates in early phase of disk evolution around a protostar

Quarks mass function at finite density in real-time formalism

Unveiling Dust Aggregate Structure in Protoplanetary Disks by Millimeter-wave Scattering Polarization

Light scattering by fractal dust aggregates: I. Angular dependence of scattering

Magnetic Structure of the $S\,{=}\,1/2$ Quasi-Two-Dimensional Square-Lattice Heisenberg Antiferromagnet Sr$_2$CuTeO$_6$

Magnetic-field- and pressure-induced quantum phase transition in CsFeCl$_3$ proved via magnetization measurement

Mass constraint for a planet in a protoplanetary disk from the gap width

Quasi-two-dimensional Bose-Einstein condensation of spin triplets in dimerized quantum magnet Ba$_2$CuSi$_2$O$_6$Cl$_2$

Sintering-induced Dust Ring Formation in Protoplanetary Disks: Application to the HL Tau Disk

Detailed structure of the outer disk around HD 169142 with polarized light in H-band

Electrical Oscillation in Pt/VO2 Bilayer Strips

Field-driven successive phase transitions in quasi-two-dimensional frustrated antiferromagnet Ba$_2$CoTeO$_6$ and highly degenerate classical ground states

Final Masses of Giant Planets II: Jupiter Formation in a Gas-Depleted Disk

Formation of a disc gap induced by a planet: Effect of the deviation from Keplerian disc rotation

From Planetesimals to Planets in Turbulent Protoplanetary Disks I. Onset of Runaway Growth

Magnetic-Field-Induced Polar Phase in Chiral Magnet CsCuCl$_3$

Mass Estimates of a Giant Planet in a Protoplanetary Disk from the Gap Structures

Quantum phase transition between disordered and ordered states in the spin-1/2 kagome lattice antiferromagnet (Rb$_{1-x}$Cs$_{x}$)$_2$Cu$_3$SnF$_{12}$

Simple improvements to classical bubble nucleation models

Spinon, soliton and breather in the spin-1/2 antiferromagnetic chain KCuGaF$_6$

Successive magnetic phase transitions in $\bmα$-RuCl$_3$: XY-like frustrated magnet on the honeycomb lattice

Almost Perfect Frustration in the Dimer Magnet Ba$_2$CoSi$_2$O$_6$Cl$_2$

Bubble Evolution and Properties in Homogeneous Nucleation Simulations

Direct Simulations of Homogeneous Bubble Nucleation: Agreement with CNT and no Local Hot Spots

Free energy of cluster formation and a new scaling relation for the nucleation rate

Magnetic Phase Diagram of the S=1/2 Triangular-Lattice Heisenberg Antiferromagnet Ba3CoNb2O9

Molecular Dynamics Simulations of the Nucleation of Water: Determining the Sticking Probability and Formation Energy of a Cluster

Non-$γ_{5}$hermiticity minimal doubling fermion

Opacity of fluffy dust aggregates

Properties of Liquid Clusters in Large-scale Molecular Dynamics Nucleation Simulations

Strong Suppression of Magnetic Ordering in an S = 1/2 Square-Lattice Heisenberg Antiferromagnet Sr2CuTeO6

The Fermi Function Beyond The Second Order Perturbation Theory

The multifaceted planetesimal formation process

TMD Parton Distributions based on Three-Body Decay Functions in NLL Order of QCD

Constraints on the charged scalar effects using the forward-backward asymmetry on $B\to D^{(*)}τ\bar{ν_τ}$

Fluffy dust forms icy planetesimals by static compression

Ground state of the spin-1/2 chain of green dioptase at high fields

Large Scale Molecular Dynamics Simulations of Homogeneous Nucleation

Magnetization Process and Collective Excitations in the S=1/2 Triangular-Lattice Heisenberg Antiferromagnet Ba3CoSb2O9

Minimal Doubling Fermion and Hermiticity

Non-$γ_{5}$hermiticity fermions in two dimensions

Static compression of porous dust aggregates

A Derivation of the Fermi Function in Perturbative Quantum Field Theory

Ferromagnetism induced in anisotropic stacked kagome-lattice antiferromagnet Cs$_2$Cu$_3$CeF$_{12}$

Geometrical Cross Sections of Dust Aggregates and a Compression Model for Aggregate Collisions

Migration rates of planets due to scattering of planetesimals

Rapid Coagulation of Porous Dust Aggregates Outside the Snow Line: A Pathway to Successful Icy Planetesimal Formation

Determination of the Higgs CP property in Hadron Colliders

Electrostatic Barrier against Dust Growth in Protoplanetary Disks. I. Classifying the Evolution of Size Distribution

Electrostatic Barrier against Dust Growth in Protoplanetary Disks. II. Measuring the Size of the "Frozen" Zone

Experimental Realization of Spin-1/2 Triangular-Lattice Heisenberg Antiferromagnet

NLL Order Contributions for Exclusive Processes in Jet-Calculus Scheme

Planetary Core Formation with Collisional Fragmentation and Atmosphere to Form Gas Giant Planets

Quantum Magnetization Plateau in Spin-1 Triangular-Lattice Antiferromagnet Ba$_3$NiSb$_2$O$_9$

Reply to "Comment on `Transition from Bose glass to a condensate of triplons in Tl$_{1-x}$K$_x$CuCl$_3$'"

Thermodynamic properties of quantum sine-Gordon spin chain system KCuGaF6

A New Proposal for Neutrino Mass and |V(ud)| Measurements

Numerical Calculation of Schwinger-Dyson Equation with Momentum-Dependent Gauge Parameter at Finite Temperature

Planetary Growth with Collisional Fragmentation and Gas Drag

Pressure-induced new magnetic phase in Tl(Cu$_{0.985}$Mg$_{0.015}$)Cl$_3$ probed by muon spin rotation

Elementary excitations of S=1/2 one-dimensional antiferromagnet KCuGaF6 in magnetic field and quantum sine-Gordon model

Numerical Modeling of the Coagulation and Porosity Evolution of Dust Aggregates

Magnetic-Field Induced Bose-Einstein Condensation of Magnons and Critical Behavior in Interacting Spin Dimer System TlCuCl$_3$

Structural Transitions and Magnetic Structure in NH4CuCl3 via 14N-NMR