Source author record

Hidenori Fukaya

Hidenori Fukaya 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

34works
13topics
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

34 published item(s)

preprint2026arXiv

Time-reversal invariant vortex in topological superconductors and gravitational $\mathbb{Z}_2$ topology

We study a time-reversal invariant vortex, namely a spin vortex, in helical superconductors by focusing on its emergent gravitational structure. The topology of the time-reversal invariant vortex is classified by a $\mathbb{Z}_2$ invariant: helical Majorana zero modes appear at the vortex core when the winding number is odd, while no such zero modes exist when it is even. We provide a formal mapping to the theory of gravity to describe this $\mathbb{Z}_2$ topological structure. Identifying a superconducting order parameter as a vielbein in the theory of gravity, we explicitly convert the Bogoliubov-de-Genne Hamiltonian into the Dirac Hamiltonian coupled to a nontrivial gravitational field. Then we find that a gravitational curvature is induced at the vortex core, with its total flux quantized in integer multiples of $π$, reflecting the $\mathbb{Z}_2$ topology. Although the curvature vanishes everywhere except at the vortex core, the energy spectrum remains sensitive to the total curvature flux, owing to the gravitational Aharonov-Bohm effect. We further demonstrate that our gravitational framework can be applied to the topological phase transition driven by the vortex-linking precess in three-dimensional helical superconductors such as the He-B phase.

preprint2023arXiv

Curved domain-wall fermion and its anomaly inflow

We investigate the effect of $U (1)$ gauge field on lattice fermion systems with a curved domain-wall mass term. In the same way as the conventional flat domain-wall fermion, the chiral edge modes appear localized at the wall, whose Dirac operator contains the induced gravitational potential as well as the $U(1)$ vector potential. In the case of $S^1$ domain-wall fermion on a two-dimensional flat lattice, we find a competition between the Aharonov-Bohm(AB) effect and gravitational gap in the Dirac eigenvalue spectrum, which leads to anomaly of the time-reversal ($T$) symmetry. Our numerical result shows a good consistency with the Atiyah-Patodi-Singer index theorem on a disk inside the $S^1$ domain-wall, which describes the cancellation of the $T$ anomaly between the bulk and edge. When the $U(1)$ flux is squeezed inside one plaquette, and the AB phase takes a quantized value $π$ mod $2π\mathbb{Z}$, the anomaly inflow drastically changes: the strong flux creates another domain-wall around the flux to make the two zero modes coexist. This phenomenon is also observed in the $S^2$ domain-wall fermion in the presence of a magnetic monopole. We find that the domain-wall creation around the monopole microscopically explains the Witten effect.

preprint2022arXiv

Axial U(1) symmetry at high temperatures in $N_f=2+1$ lattice QCD with chiral fermions

We study the $U(1)_A$ anomaly in the high-temperature phase of $N_f=2+1$ lattice QCD with chiral fermions. Gauge ensembles are generated with Möbius domain-wall (MDW) fermions, and in the measurements the determinant is reweighted to that of overlap fermions. We report the results for the overlap Dirac spectrum, $U(1)_A$ susceptibility, and topological susceptibility at $T=204$ and $175$ MeV.

preprint2022arXiv

Curved domain-wall fermions

We consider fermion systems on a square lattice with a mass term having a curved domain-wall. Similarly to the conventional flat domain-wall fermions, massless and chiral edge states appear on the wall. In the cases of $S^1$ and $S^2$ domain-walls embedded into flat hypercubic lattices, we find that these edge modes feel gravity through the induced Spin or Spin$^c$ connections. The gravitational effect is encoded in the Dirac eigenvalue spectrum as a gap from zero. In the standard continuum extrapolation of the square lattice, we find a good agreement with the analytic prediction in the continuum theory. We also find that the rotational symmetry of the edge modes is automatically recovered in the continuum limit.

preprint2020arXiv

A lattice formulation of the Atiyah-Patodi-Singer index

Atiyah-Singer index theorem on a lattice without boundary is well understood owing to the seminal work by Hasenfratz et al. But its extension to the system with boundary (the so-called Atiyah- Patodi-Singer index theorem), which plays a crucial role in T-anomaly cancellation between bulk- and edge-modes in 3+1 dimensional topological matters, is known only in the continuum theory and no lattice realization has been made so far. In this work, we try to non-perturbatively define an alternative index from the lattice domain-wall fermion in 3+1 dimensions. We will show that this new index in the continuum limit, converges to the Atiyah-Patodi-Singer index defined on a manifold with boundary, which coincides with the surface of the domain-wall.

preprint2020arXiv

A physicist-friendly reformulation of the Atiyah-Patodi-Singer index and its mathematical justification

The Atiyah-Patodi-Singer index theorem describes the bulk-edge correspondence of symmetry protected topological insulators. The mathematical setup for this theorem is, however, not directly related to the physical fermion system, as it imposes on the fermion fields a non-local and unnatural boundary condition known as the "APS boundary condition" by hand. In 2017, we showed that the same integer as the APS index can be obtained from the $η$ invariant of the domain-wall Dirac operator. Recently we gave a mathematical proof that the equivalence is not a coincidence but generally true. In this contribution to the proceedings of LATTICE 2019, we try to explain the whole story in a physicist-friendly way.

preprint2020arXiv

Atiyah-Patodi-Singer index on a lattice

We propose a non-perturbative formulation of the Atiyah-Patodi-Singer(APS) index in lattice gauge theory, in which the index is given by the $η$ invariant of the domain-wall Dirac operator. Our definition of the index is always an integer with a finite lattice spacing. To verify this proposal, using the eigenmode set of the free domain-wall fermion, we perturbatively show in the continuum limit that the curvature term in the APS theorem appears as the contribution from the massive bulk extended modes, while the boundary $η$ invariant comes entirely from the massless edge-localized modes.

preprint2020arXiv

Axial U(1) symmetry and mesonic correlators at high temperature in $N_f=2$ lattice QCD

We investigate the high-temperature phase of QCD using lattice QCD simulations with $N_f = 2$ dynamical Möbius domain-wall fermions. On generated configurations, we study the axial $U(1)$ symmetry, overlap-Dirac spectra, screening masses from mesonic correlators, and topological susceptibility. We find that some of the observables are quite sensitive to lattice artifacts due to a small violation of the chiral symmetry. For those observables, we reweight the Möbius domain-wall fermion determinant by that of the overlap fermion. We also check the volume dependence of observables. Our data near the chiral limit indicates a strong suppression of the axial $U(1)$ anomaly at temperatures $\geq$ 220 MeV.

preprint2020arXiv

The Anomaly Inflow of the domain-wall fermion in odd dimension

In 1985, Callan and Harvey showed a view of gauge anomaly as a missing current into an extra-dimension, and the total contribution, including the Chern-Simons current in the bulk, is conserved. However in their computation, the edge and bulk contributions are separately evaluated and their cross correlations, which should be relevant at boundary, are simply ignored. This issue has been solved in many approaches. In this work, we revisit this issue with a complete set of eigenstates of free domain-wall Hamiltonian and give the systematic evaluation, easy to take in the higher mass correction and extend to the higher dimension.

preprint2020arXiv

TKNN formula for general Hamiltonian

Topological insulators in odd dimensions are characterized by topological numbers. We prove the well-known relation between the topological number given by the Chern character of the Berry curvature and the Chern-Simons level of the low energy effective action for a general class of Hamiltonians bilinear in the fermion with general U(1) gauge interactions including non-minimal couplings by an explicit calculation. A series of Ward-Takahashi identities are crucial to relate the Chern-Simons level to a winding number, which could then be directly reduced to Chern character of Berry curvature by carrying out the integral over the temporal momenta.

preprint2020arXiv

TKNN formula for general lattice Hamiltonian in odd dimensions

Topological insulators in odd dimensions are characterized by topological numbers. We prove the well-known relation between the topological number given by the Chern character of the Berry curvature and the Chern-Simons level of the low energy effective action for a general class of Hamiltonians bilinear in the fermion with general U(1) gauge interactions including non-minimal couplings by an explicit calculation. A series of Ward-Takahashi identities are crucial to relate the Chern-Simons level to a winding number, which could then be directly reduced to Chern character of Berry curvature by carrying out the integral over the temporal momenta.

preprint2020arXiv

Towards fully non-perturbative computation of inelastic $\ell N$ scattering cross sections from lattice QCD

We propose a fully non-perturbative method to compute inelastic lepton-nucleon ($\ell N$) scattering cross sections using lattice QCD. The method is applicable even at low energies, such as the energy region relevant for the recent and future neutrino-nucleon ($νN$) scattering experiments, for which perturbative analysis is invalidated. The basic building block is the forward Compton-scattering amplitude, or the hadronic tensor, computed on a Euclidean lattice. Total cross section is constructed from the hadronic tensor by multiplying a phase space factor and integrating over the energy and momentum of final hadronic states. The energy integral that induces a sum over all possible final states is performed implicitly by promoting the phase space factor to an operator written in terms of the transfer matrix on the lattice. The formalism is imported from that of the inclusive semileptonic B meson decay [P. Gambino, S. Hashimoto, arXiv:2005.13730]. It can be generalized to compute the $\ell N$ scattering cross sections and their moments, as well as the virtual correction to the nuclear $β$-decay. Necessary quark-line contractions for two current insertions corresponding to the Compton amplitude to be computed on the lattice are summarized.

preprint2019arXiv

Axial U(1) symmetry, topology, and Dirac spectra at high temperature in $N_f=2$ lattice QCD

Using lattice QCD simulations with $N_f = 2$ dynamical fermions, we study the axial $U(1)$ symmetry, topological charge, and Dirac eigenvalue spectra in the high-temperature phase in which the chiral symmetry is restored. Our gauge ensembles are generated with Möbius domain-wall fermions, but the measurements such as susceptibilities are reweighted to those for the overlap fermions by using overlap/domain-wall reweighting technique. We find that the $U(1)_A$ and topological susceptibilities are strongly suppressed in the small quark mass region, which is related to the reduction of chiral-zero and low-nonzero modes on the Dirac spectra. We also examine their volume dependence.

preprint2018arXiv

Axial U(1) symmetry and Dirac spectra in high-temperature phase of $N_f=2$ lattice QCD

The axial $U(1)$ symmetry in the high-temperature phase is investigated with $N_f = 2$ lattice QCD simulations. The gauge ensembles are generated with Möbius domain-wall fermions, and the overlap/domain-wall reweighting is applied. We find that the $U(1)_A$ susceptibility evaluated from the spectrum of overlap-Dirac eigenvalues is strongly suppressed in the chiral limit. We also study its volume dependence.

preprint2016arXiv

Stochastic calculation of the Dirac spectrum on the lattice and a determination of chiral condensate in 2+1-flavor QCD

We compute the chiral condensate in 2+1-flavor QCD through the spectrum of low-lying eigenmodes of Dirac operator. The number of eigenvalues of the Dirac operator is evaluated using a stochastic method with an eigenvalue filtering technique on the background gauge configurations generated by lattice QCD simulations including the effects of dynamical up, down and strange quarks described by the Mobius domain-wall fermion formulation. The low-lying spectrum is related to the chiral condensate, which is one of the leading order low-energy constants in chiral effective theory, as dictated by the Banks-Casher relation. The spectrum shape and its dependence on the sea quark masses calculated in numerical simulations are consistent with the expectation from one-loop chiral perturbation theory. After taking the chiral limit as well as the continuum limit using the data at three lattice spacings ranging 0.080-0.045 fm, we obtain $Σ^{1/3}$(2 GeV) = 270.0(4.9) MeV, with the error combining those from statistical and from various sources of systematic errors. Finite volume effect is confirmed to be under control by a direct comparison of the results from two different volumes at the lightest available sea quarks corresponding to 230 MeV pions.

preprint2015arXiv

On the axial $U(1)$ symmetry at finite temperature

We study the $U(1)_A$ anomaly in two-flavor lattice QCD at finite temperature with the Möbius domain-wall Dirac operator. We generate gauge configurations in the temperature range $(0.9, 1.2) T_c$ on different physical volumes, $L=$ 2--4 fm, and lattice spacings. We measure the difference of the susceptibilities of the flavor non-singlet scalar ($χ_δ$) and pseudoscalar ($χ_π$) mesons. They are related by an axial $U(1)$ transformation and the difference vanishes if the axial $U(1)$ symmetry is respected. We identify the source of axial $U(1)$ symmetry breaking at finite temperature in the lowest eigenmodes, for the observable $χ_π- χ_δ$. We then reweight the Möbius domain-wall fermion partition function to that of the overlap-Dirac operator to fully recover chiral symmetry. Our data show a significant discrepancy in the results coming from the Möbius domain-wall valence quarks, the overlap valence quarks on our DWF configurations and the reweighted ones that have full chiral symmetry. After recovering full chiral symmetry we conclude that the difference $χ_π- χ_δ$ shows a suppression in the chiral limit that is compatible with an effective restoration of $U(1)_A$ at $T \gtrsim T_c$ in the scalar meson channels.

preprint2015arXiv

Renormalization of Extended QCD$_2$

Extended QCD (XQCD) proposed by Kaplan [1] is an interesting reformulation of QCD with additional bosonic auxiliary fields. While its partition function is kept exactly the same as that of original QCD, XQCD naturally contains properties of low energy hadronic models. We analyze the renormalization group flow of two-dimensional (X)QCD, which is solvable in the limit of large number of colors Nc, to understand what kind of roles the auxiliary degrees of freedom play and how the hadronic picture emerges in the low energy region.

preprint2015arXiv

Renormalization of two-dimensional XQCD

Recently, Kaplan proposed an interesting extension of QCD named Extended QCD or XQCD with bosonic auxiliary fields [1]. While its partition function is kept exactly the same as that of QCD, XQCD naturally contains properties of low-energy hadrons. We apply this extension to the two-dimensional QCD in the large $N_c$ limit ('t Hooft model) [2]. In this solvable model, it is possible to directly examine the hadronic picture of the 2d XQCD and analyze its renormalization group flow to understand how the auxiliary degrees of freedom behave in the low energy region. We confirm that the additional scalar fields can become dynamical acquiring the kinetic term, and its parity-odd part becomes dominant in the low energy region. This renomalization of XQCD provides an "extension" of the renormalization scheme of QCD, inserting different field variables from those in the original theory, without any changes in physical observables.

preprint2015arXiv

Violation of chirality of the Möbius domain-wall Dirac operator from the eigenmodes

We investigate the effects of the violation of the Ginsparg-Wilson (GW) relation in the Möbius domain-wall fermion formulation on the lattice with finite fifth dimension. Using a decomposion in terms of the eigenmodes of its four-dimensional effective Dirac operator, we isolate the GW-violating terms for various physical quantities including the residual mass and the meson susceptibilities relevant for the effective restoration of the axial U(1) symmetry at finite temperature. Numerical result shows that the GW-violating effect is more significant, or even overwhelming, for the quantities that are dominated by the low-lying eigenmodes.

preprint2014arXiv

Axial U(1) symmetry at finite temperature with Möbius domain-wall fermions

We investigate the axial $U(1)$ symmetry restoration at finite temperature in two flavor QCD. We employ the Möbius domain-wall formalism that is designed to achieve good chiral symmetry. We show the measurements of a difference of meson susceptibilities, sensitive to the $U(1)_A$ symmetry breaking. The signal is dominated by zero and near-zero modes. By reweighting the measure to that of overlap fermions we find a suppression of the $U(1)_A$ breaking effects above the chiral transition temperature.

preprint2014arXiv

Effects of near-zero Dirac eigenmodes on axial U(1) symmetry at finite temperature

We study the axial U(1)A symmetry of Nf = 2 QCD at finite temperature using the Dirac eigenvalue spectrum. The gauge configurations are generated employing the Mobius domain-wall fermion action on 16^3x8 and 32^3x8 lattices. The physical spatial size of these lattices is around 2 fm and 4 fm, respectively, and the simulated temperature is around 200 MeV, which is slightly above the critical temperature of the chiral phase transition. Although the Mobius domain-wall Dirac operator is expected to have a good chiral symmetry and our data actually show small values of the residual mass, we observe significant violation of the Ginsparg-Wilson relation for the low- lying eigenmodes of the Mobius domain-wall Dirac operator. Using the reweighting technique, we compute the overlap-Dirac operator spectrum on the same set of configurations and find a significant difference of the spectrum between the two Dirac operators for the low-lying eigenvalues. The overlap-Dirac spectrum shows a gap from zero, which is insensitive to the spacial volume.

preprint2014arXiv

Extracting the electro-magnetic pion form factor from QCD in a finite volume

We consider finite volume effects on the electromagnetic form factor of the pion. We compute the peudoscalar-vector-pseudoscalar correlator in the $ε$ expansion of chiral perturbation theory up to the next-to-leading order and find a way to remove the dominant part, which comes from a contribution of the pion zero-mode. Inserting non-zero momentum to relevant operators (or taking a subtraction of the correlators at different time-slices), and taking an appropriate ratio of them, one can automatically cancel the zero-mode's contribution, which becomes non-perturbatively large $\sim \mathcal{O}(100 \%)$ in the $ε$ regime. The remaining finite volume dependence, which comes from the non-zero momentum modes, is shown to be perturbatively small even in such an extremal case. Since the zero-mode's dominance is universal in any finite volume scaling, and we do not rely on any particular feature of the $ε$ expansion, our method has a wide application to many other correlators of QCD.

preprint2014arXiv

Finite volume scaling of the electro-magnetic pion form factor in the $ε$ regime

We consider finite volume effects on the electro-magnetic pion form factor near the chiral limit, in the so-called $ε$ regime. The pseudoscalar-vector-pseudoscalar three-point function is calculated in the $ε$ expansion of chiral perturbation theory to the next-to-leading order. In the $ε$ regime, finite volume effects are non-perturbatively large in general. However, we find a way to remove its dominant part, by inserting momenta to the correlators, and taking an appropriate ratio of them. The subleading contribution is, then, shown to be perturbatively small, and one can extract the form factor as in a similar way to that in the $p$ regime.

preprint2013arXiv

1st or 2nd; the order of finite temperature phase transition of Nf=2 QCD from effective theory analysis

In the previous work, we have shown that the SU(2) chiral symmetry recovered above the critical temperature gives a strong constraint on the Dirac eigenvalue spectrum and this constraint is strong enough for a set of anomalous U(1) chiral symmetry breaking operators to vanish in the thermodynamical and chiral limits. We use this condition as an input and impose a constraint on the Landau low energy effective theory of QCD. The only constraint we can set is that the mass splitting term between the pion and eta meson should vanish. All the singlet/non-singlet scalar/pseudo-scalar mesons contribute to the effective theory. We evaluate the renormalization group $β$-function for the effective theory using the $ε$-expansion at one loop level, but find no stable infra-red fixed point except for the trivial Gaussian one. The chiral phase transition seems to be of first order.

preprint2013arXiv

Chiral symmetry restoration, eigenvalue density of Dirac operator and axial U(1) anomaly at finite temperature

We reconsider constraints on the eigenvalue density of the Dirac operator in the chiral symmetric phase of 2 flavor QCD at finite temperature. To avoid possible ultra-violet(UV) divergences, we work on a lattice, employing the overlap Dirac operator, which ensures the exact "chiral" symmetry at finite lattice spacings. Studying multi-point correlation functions in various channels and taking their thermodynamical limit (and then taking the chiral limit), we obtain stronger constraints than those found in the previous studies: both the eigenvalue density at the origin and its first and second derivatives vanish in the chiral limit of 2 flavor QCD. In addition we show that the axial U(1) anomaly becomes invisible in susceptibilities of scalar and pseudo scalar mesons, suggesting that the 2nd order chiral phase transition with the O(4) scaling is not realized in 2 flavor QCD. Possible lattice artifacts when non-chiral lattice Dirac operator is employed are briefly discussed.

preprint2013arXiv

Finite temperature study of the axial U(1) symmetry on the lattice with overlap fermion formulation

We examine the axial U(1) symmetry near and above the finite temperature phase transition in two-flavor QCD using lattice QCD simulations. Although the axial U(1) symmetry is always violated by quantization, (i.e.) the chiral anomaly, the correlation functions may manifest effective restoration of the symmetry in the high temperature phase. We explicitly study this possibility by calculating the meson correlators as well as the Dirac operator spectral density near the critical point. Our numerical simulations are performed on a $16^3\times 8$ lattice with two flavors of dynamical quarks represented by the overlap fermion formalism. Chiral symmetry and its violation due to the axial anomaly is manifestly realized with this formulation, which is a prerequisite for the study of the effective restoration of the axial U(1) symmetry. In order to avoid discontinuity in the gauge configuration space, which occurs for the exactly chiral lattice fermions, the simulation is confined in a fixed topological sector. It induces finite volume effect, which is well described by a formula based on the Fourier transform from the $θ$-vacua. We confirm this formula at finite temperature by calculating the topological susceptibility in the quenched theory. Our two flavor simulations show degeneracy of the meson correlators and a gap in the Dirac operator spectral density, which implies that the axial U(1) symmetry is effectively restored in the chirally symmetric phase.

preprint2012arXiv

Chiral symmetry restoration and eigenvalue density of Dirac operator

We reinvestigate constraints on the eigenvalue density of the Dirac operator in the chiral symmetric phase of 2 flavor QCD at finite temperature, employing the overlap Dirac operator with the exact chiral symmetry at finite lattice spacings to avoid possible ultra-violet(UV) divergences. Studying multi-point correlation functions in various channels in the thermodynamical limit, we obtain stronger constraints than those found in the previous studies that not only the eigenvalue density at the origin but also its first and second derivatives vanish in the chiral limit of 2 flavor QCD. In addition we show that the axial U(1) anomaly becomes invisible in susceptibilities of scalar and pseudo scalar mesons.

preprint2012arXiv

The neutral pion decay and the chiral anomaly on the lattice

We perform a lattice QCD calculation of the π^0\toγγtransition form factor and the associated decay width. We use a Euclidean time integral of the relevant three-point function to compute the decay amplitude for two-photon final state, which is a non-QCD state. We use the all-to-all quark propagator technique to carry out this integral as well as to include the disconnected quark diagram contributions. We execute the calculation using the overlap fermion formulation, which ensures the exact chiral symmetry on the lattice and produces the chiral anomaly through the Jacobian of the chiral transformation. We examine various sources of systematic effects except for possible discretization effect. Our final results for the form factor and the decay width reproduce the ABJ anomaly in the chiral limit and agree with the experimental measurements at the physical pion mass with a precision of a few percent.

preprint2012arXiv

Two-photon decay of the neutral pion in lattice QCD

We perform non-perturbative calculation of the π^0 to γγ transition form factor and the associated decay width using lattice QCD. The amplitude for two-photon final state, which is not an eigenstate of QCD, is extracted through an Euclidean time integral of the relevant three-point function. We utilize the all-to-all quark propagator technique to carry out this integral as well as to include the disconnected quark diagram contributions. The overlap fermion formulation is employed on the lattice to ensure exact chiral symmetry on the lattice. After examining various sources of systematic effects except for possible discretization effect, we obtain Γ=7.83(31)(49) eV for the pion decay width, where the first error is statistical and the second is our estimate of the systematic error.

preprint2011arXiv

Interpolation between the epsilon and p regimes

We reconsider chiral perturbation theory in a finite volume and develop a new computational scheme which smoothly interpolates the conventional epsilon and p regimes. The counting rule is kept essentially the same as in the p expansion. The zero-momentum modes of Nambu-Goldstone bosons are, however, treated separately and partly integrated out to all orders as in the epsilon expansion. In this new scheme, the theory remains infra-red finite even in the chiral limit, while the chiral-logarithmic effects are kept present. We calculate the two-point function in the pseudoscalar channel and show that the correlator has a constant contribution in addition to the conventional hyperbolic cosine function of time t. This constant term rapidly disappears in the p regime but it is indispensable for a smooth convergence of the formula to the epsilon regime result. Our calculation is useful to precisely estimate the finite volume effects in lattice QCD simulations on the pion mass Mpi and kaon mass MK, as well as their decay constants Fpi and FK.

preprint2010arXiv

Chiral perturbation theory in a theta vacuum

We consider chiral perturbation theory (ChPT) with a non-zero theta term. Due to the CP violating term, the vacuum of chiral fields is shifted to a non-trivial element on the SU(N_f) group manifold. The CP violation also provides mixing of different CP eigenstates, between scalar and pseudoscalar, or vector and axialvector operators. We investigate upto O(theta^2) effects on the mesonic two point correlators of ChPT to the one-loop order. We also address the effects of fixing topology, by using saddle point integration in the Fourier transform with respect to theta.