Source author record

Takahiro Ueda

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

25works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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.

Published work

25 published item(s)

preprint2022arXiv

A global two-layer radiative transfer model for axisymmetric, shadowed protoplanetary disks

Understanding the thermal structure of protoplanetary disks is crucial for modeling planet formation and interpreting disk observations. We present a new two-layer radiative transfer model for computing the thermal structure of axisymmetric irradiated disks. Unlike the standard two-layer model, our model accounts for the radial as well as vertical transfer of the starlight reprocessed at the disk surface. The model thus allows us to compute the temperature below "shadowed" surfaces receiving no direct starlight. Thanks to the assumed axisymmetry, the reprocessed starlight flux is given in one-dimensional integral form that can be computed at a low cost. Furthermore, our model evolves the midplane temperature using a time-dependent energy equation and can therefore treat thermal instabilities. We apply our global two-layer model to disks with a planetary induced gap and confirm that the model reproduces the disks' temperature profiles obtained from more computationally expensive Monte Carlo radiative transfer calculations to an accuracy of less than 20%. We also apply the model to study the long-term behavior of the thermal wave instability in irradiated disks. Being simple and computationally efficient, the global two-layer model will be suitable for studying the interplay between disks' thermal evolution and dust evolution.

preprint2022arXiv

Massive compact dust disk with a gap around CW Tau revealed by ALMA multi-band observations

Compact protoplanetary disks with a radius of $\lesssim$ 50 au are common around young low-mass stars. We report high resolution ALMA dust continuum observations toward a compact disk around CW Tau at Band 4 ($λ=2.2$ mm), 6 (1.3 mm), 7 (0.89 mm) and 8 (0.75 mm). The SED shows the spectral slope of $2.0\pm0.24$ between 0.75 and 1.3 mm, while it is $3.7\pm0.29$ between 2.17 and 3.56 mm. The steep slope between 2.17 and 3.56 mm is consistent with that of optically thin emission from small grains ($\lesssim$ 350 ${\rm μm}$). We perform parametric fitting of the ALMA data to characterize the dust disk. Interestingly, if the dust-to-gas mass ratio is 0.01, the Toomre's Q parameter reaches $\sim$ 1-3, suggesting that the CW Tau disk might be marginally gravitationally unstable. The total dust mass is estimated as $\sim250M_{\oplus}$ for the maximum dust size of 140 ${\rm μm}$ that is inferred from the previous Band 7 polarimetric observation and at least $80M_{\oplus}$ even for larger grain sizes. This result shows that the CW Tau disk is quite massive in spite of its smallness. Furthermore, we clearly identify a gap structure located at $\sim20$ au, which might be induced by a giant planet. In spite of these interesting characteristics, the CW Tau disk has normal disk luminosity, size and spectral index at ALMA Band 6, which could be a clue to the mass budget problem in Class II disks.

preprint2022arXiv

Probing Inner-Edge of Dead Zones in Protoplanetary Disks with ALMA and Next Generation Very Large Array

The discovery of substructures in protoplanetary disks with ALMA has provided us key insights on the formation of planets. However, observational constraints on the formation of rocky planets have been still sparse, especially because of the limited spatial resolution. The inner edge of so-called dead zone is one of the preferential sites of rocky planet formation. We investigate the capabilities of ALMA and ngVLA for observing a dust concentration expected at the inner edge of the dead-zone around a Herbig star. Herbig Ae/Be stars are useful laboratories for exploring the evolution of rocky grains in protoplanetary disks because of their high luminosity which pushes the dead-zone inner edge outward. We find that, thanks to its unprecedented angular resolution and sensitivity, ngVLA can detect the dust concentration at the dead-zone inner edge, with a reasonable integration time of 10 hrs at $λ=3,7$ mm and 1 cm. The dust concentration is expected to be optically thick at the ALMA wavelengths and cannot be spatially resolved due to its limited resolution. On the other hand, the flux density from the inner disk regions ($\sim$3--4 au) observed with current VLA is higher for disks with a dust ring, and hence would be a useful indicator that help us choose potential candidates of disks having a dust concentration at the inner most region. With these observations we can characterize the process of dust concentration in the innermost disk regions, where rocky planets can form.

preprint2021arXiv

Higgs triplet extension of GRACE

Much theoretical effort and automatization are required to confront new physics models with experimental data for many types of particle reactions at future colliders. In this context, we extend GRACE, an automatic calculation system for invariant amplitudes, to incorporate particles and interactions in the Georgi-Machacek model. With the extended GRACE system, we study fermiophobic Higgs boson production processes at $e^+ e^-$ and $e^- e^-$ colliders in the model. The results show some advantages of $e^- e^-$ colliders over $e^+ e^-$ colliders for new physics search and thus its complementary role.

preprint2020arXiv

Non-global logarithms in hadron collisions at $N_c=3$

We calculate the rapidity gap survival probability associated with the Higgs decay and Higgs plus dijet production in proton-proton collisions by resumming the leading non-global logarithms without any approximation to the number of colors. For dijet production, depending on partonic subprocesses, the probability involves various `color multipoles', i.e., the product of 4 ($qq\to qq$) or 6 ($qg\to qg$) or 8 ($gg\to gg$) Wilson lines. We calculate all these multipoles for a fixed dijet configuration and discuss the factorization of higher multipoles into lower multipoles as well as the validity of the large-$N_c$ approximation.

preprint2020arXiv

Scattering-induced intensity reduction: large mass content with small grains in the inner region of the TW Hya disk

Dust continuum observation is one of the best methods to constrain the properties of protoplanetary disks. Recent theoretical studies have suggested that the dust scattering at the millimeter wavelength potentially reduces the observed intensity, which results in an underestimate in the dust mass. We investigate whether the dust scattering indeed reduces the observed continuum intensity by comparing the ALMA archival data of the TW Hya disk at Band 3, 4, 6, 7 and 9 to models obtained by radiative transfer simulations. We find that the model with scattering by 300 ${\rm μm}$-sized grains well reproduces the observed SED of the central part of the TW Hya disk while the model without scattering is also consistent within the errors of the absolute fluxes. To explain the intensity at Band 3, the dust surface density needs to be $\sim$ 10 ${\rm g\,cm^{-2}}$ at 10 au in the model with scattering, which is 26 times more massive than previously predicted. The model without scattering needs 2.3 times higher dust mass than the model with scattering because it needs lower temperature. At Band 7, scattering reduces the intensity by $\sim$ 35% which makes the disk looks optically thin even though it is optically thick. Our study suggests the TW Hya disk is still capable of forming cores of giant planets at where the current solar system planets exist.

preprint2019arXiv

QCD Resummation in Hard Diffractive Dijet Production at the Electron-Ion Collider

Diffractive dijet production at the electron-ion collider (EIC) has been proposed to study the gluon Wigner distribution at small-$x$. We investigate the soft gluon radiation associated with the final state jets and an all order resummation formula is derived. We show that the soft gluon resummation plays an important role to describe E791 data on $π$-induced diffractive dijet production at Fermilab. Predictions for the EIC are presented, and we emphasize that the soft gluon resummation is an important aspect to explore the nucleon/nucleus tomography through these processes.

preprint2015arXiv

Exact N$^3$LO results for $q q^\prime\to H +X$

We compute the contribution to the total cross section for the inclusive production of a Standard Model Higgs boson induced by two quarks with different flavour in the initial state. Our calculation is exact in the Higgs boson mass and the partonic center-of-mass energy. We describe the reduction to master integrals, the construction of a canonical basis, and the solution of the corresponding differential equations. Our analytic result contains both Harmonic Polylogarithms and iterated integrals with additional letters in the alphabet.

preprint2015arXiv

How to solve the cake-cutting problem in sublinear time

In this paper, we show algorithms for solving the cake-cutting problem in sublinear-time. More specifically, we preassign (simple) fair portions to o(n) players in o(n)-time, and minimize the damage to the rest of the players. All currently known algorithms require Omega(n)-time, even when assigning a portion to just one player, and it is nontrivial to revise these algorithms to run in $o(n)$-time since many of the remaining players, who have not been asked any queries, may not be satisfied with the remaining cake. To challenge this problem, we begin by providing a framework for solving the cake-cutting problem in sublinear-time. Generally speaking, solving a problem in sublinear-time requires the use of approximations. However, in our framework, we introduce the concept of "eps n-victims," which means that eps n players (victims) may not get fair portions, where 0< eps =< 1 is an arbitrary constant. In our framework, an algorithm consists of the following two parts: In the first (Preassigning) part, it distributes fair portions to r < n players in o(n)-time. In the second (Completion) part, it distributes fair portions to the remaining n-r players except for the eps n victims in poly}(n)-time. There are two variations on the r players in the first part. Specifically, whether they can or cannot be designated. We will then present algorithms in this framework. In particular, an O(r/eps)-time algorithm for r =< eps n/127 undesignated players with eps n-victims, and an O~(r^2/eps)-time algorithm for r =< eps e^{sqrt{ln{n}}/{7}} designated players and eps =< 1/e with eps n-victims are presented.

preprint2015arXiv

The diamond rule for multi-loop Feynman diagrams

An important aspect of improving perturbative predictions in high energy physics is efficiently reducing dimensionally regularised Feynman integrals through integration by parts (IBP) relations. The well-known triangle rule has been used to achieve simple reduction schemes. In this work we introduce an extensible, multi-loop version of the triangle rule, which we refer to as the diamond rule. Such a structure appears frequently in higher-loop calculations. We derive an explicit solution for the recursion, which prevents spurious poles in intermediate steps of the computations. Applications for massless propagator type diagrams at three, four, and five loops are discussed.

preprint2014arXiv

Adequate bases of phase space master integrals for $gg \to h$ at NNLO and beyond

We study master integrals needed to compute the Higgs boson production cross section via gluon fusion in the infinite top quark mass limit, using a canonical form of differential equations for master integrals, recently identified by Henn, which makes their solution possible in a straightforward algebraic way. We apply the known criteria to derive such a suitable basis for all the phase space master integrals in afore mentioned process at next-to-next-to-leading order in QCD and demonstrate that the method is applicable to next-to-next-to-next-to-leading order as well by solving a non-planar topology. Furthermore, we discuss in great detail how to find an adequate basis using practical examples. Special emphasis is devoted to master integrals which are coupled by their differential equations.

preprint2013arXiv

Resummation of non-global logarithms at finite $N_c$

In the context of inter-jet energy flow, we present the first quantitative result of the resummation of non-global logarithms at finite N_c. This is achieved by refining Weigert's approach in which the problem is reduced to the simulation of associated Langevin dynamics in the space of Wilson lines. We find that, in e+e- annihilation, the exact result is rather close to the result previously obtained in the large-N_c mean field approximation. However, we observe enormous event-by-event fluctuations in the Langevin process which may have significant consequences in hadron collisions.

preprint2012arXiv

Higgs boson production at the LHC: NNLO partonic cross sections through order $ε$ and convolutions with splitting functions to N$^3$LO

We consider Higgs boson production at hadron colliders in the gluon fusion channel and compute higher order terms in the regularization parameter $ε$. In particular, the next-to-next-to-leading order cross section is evaluated including order $ε$ terms. These results are used to compute all convolutions with the splitting functions entering the next-to-next-to-next-to-leading order cross section.

preprint2011arXiv

Eccentricity and elliptic flow in proton-proton collisions from parton evolution

It has been argued that high-multiplicity proton-proton collisions at the LHC may exhibit collective phenomena usually studied in the context of heavy-ion collisions, such as elliptic flow. We study this issue using DIPSY - a Monte Carlo event generator based on the QCD dipole model. We calculate the eccentricity of the transverse area defined by the spatial distribution of produced gluons. The resulting elliptic flow is estimated to be about 6%, comparable to the value in nucleus-nucleus collisions at RHIC and the LHC. Experimentally, elliptic flow is inferred from the azimuthal correlation between hadrons, which receives contributions from collective flow, and from various other effects referred to as "nonflow". We discuss how to identify in experiments the signal of flow in the presence of large nonflow effects.

preprint2010arXiv

Heavy quark effects on parton distribution functions in the unpolarized virtual photon up to the next-to-leading order in QCD

We investigate the heavy quark mass effects on the parton distribution functions in the unpolarized virtual photon up to the next-to-leading order in QCD. Our formalism is based on the QCD-improved parton model described by the DGLAP evolution equation as well as on the operator product expansion supplemented by the mass-independent renormalization group method. We evaluate the various components of the parton distributions inside the virtual photon with the massive quark effects, which are included through the initial condition for the heavy quark distributions, or equivalently from the matrix element of the heavy quark operators. We discuss some features of our results for the heavy quark effects and their factorization-scheme dependence.

preprint2010arXiv

Resummation of large logarithms in the heavy quark effects on the parton distributions inside the virtual photon

We discuss the resummation of the large logarithmic terms appearing in the heavy quark effects on parton distribution functions inside the virtual photon. We incorporate heavy quark mass effects by changing the initial condition of the leading-order DGLAP evolution equation. In a certain kinematical limit, we recover the logarithmic terms of the next-to-leading order heavy quark effects obtained in the previous work. This method enables us to resum the large logarithmic terms due to heavy quark mass effects on the parton distributions in the virtual photon. We numerically calculate parton distributions using the formulae derived in this work, and discuss the property of the resummed heavy quark effects.

preprint2010arXiv

Sector decomposition via computational geometry

A non-iterative method is presented for the factorization step of sector decomposition method, which separates infrared divergent part from loop integration. This method is based on a classification of asymptotic behavior of polynomials. The problem is converted to ones for convex body in Euclidean space. They are solved with algorithms developed in computational geometry. A test implementation shows that this method produces less number of decomposed sectors than usual iterative sector decompositions.

preprint2010arXiv

Soft photon anomaly and gauge/string duality

Motivated by the recent DELPHI report on the anomalous photon production in e^+e^- annihilation, we exactly calculate the inclusive cross section of soft photons in the strong coupling limit of N=4 super Yang-Mills. We find that the energy distribution is that of the Bremsstrahlung, while the angular distribution is spherical. Our result elucidates a new non-perturbative source of soft photons not associated with the final state hadronic Bremsstrahlung.

preprint2009arXiv

A geometric method of sector decomposition

We propose a new geometric method of IR factorization in sector decomposition. The problem is converted into a set of problems in convex geometry. The latter problems are solved using algorithms in combinatorial geometry. This method provides a deterministic algorithm and never falls into an infinite loop. The number of resulting sectors depends on the algorithm of triangulation. Our test implementation shows smaller number of sectors comparing with other existing methods with iterations.

preprint2009arXiv

New implementation of the sector decomposition on FORM

Nowadays the sector decomposition technique, which can isolate divergences from parametric representations of integrals, becomes a quite useful tool for numerical evaluations of the Feynman loop integrals. It is used to verify the analytical results of multi-loop integrals in the Euclidean region, or in some cases practically used in the physical region by combining with other methods handling the threshold. In an intermediate stage of the sector decomposition for the multi-loop integrals, one often has to handle enormously large expressions containing tons of terms. The symbolic manipulation system FORM is originally designed to treat such huge expressions and has a strong advantage for it. In this talk, the implementation of the sector decomposition algorithm on FORM is discussed. A number of concrete examples including cases of multi-loop diagrams are also shown.

preprint2005arXiv

Motion of the Tippe Top : Gyroscopic Balance Condition and Stability

We reexamine a very classical problem, the spinning behavior of the tippe top on a horizontal table. The analysis is made for an eccentric sphere version of the tippe top, assuming a modified Coulomb law for the sliding friction, which is a continuous function of the slip velocity $\vec v_P$ at the point of contact and vanishes at $\vec v_P=0$. We study the relevance of the gyroscopic balance condition (GBC), which was discovered to hold for a rapidly spinning hard-boiled egg by Moffatt and Shimomura, to the inversion phenomenon of the tippe top. We introduce a variable $ξ$ so that $ξ=0$ corresponds to the GBC and analyze the behavior of $ξ$. Contrary to the case of the spinning egg, the GBC for the tippe top is not fulfilled initially. But we find from simulation that for those tippe tops which will turn over, the GBC will soon be satisfied approximately. It is shown that the GBC and the geometry lead to the classification of tippe tops into three groups: The tippe tops of Group I never flip over however large a spin they are given. Those of Group II show a complete inversion and the tippe tops of Group III tend to turn over up to a certain inclination angle $θ_f$ such that $θ_f<π$, when they are spun sufficiently rapidly. There exist three steady states for the spinning motion of the tippe top. Giving a new criterion for stability, we examine the stability of these states in terms of the initial spin velocity $n_0$. And we obtain a critical value $n_c$ of the initial spin which is required for the tippe top of Group II to flip over up to the completely inverted position.