BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Thomas Luu

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

ResearcherUnclaimed source record
hep-latnucl-thhep-phcond-mat.str-elApplicationsastro-ph.COComputationcond-mat.dis-nncond-mat.quant-gascond-mat.stat-mechnucl-exphysics.comp-phq-fin.CPq-fin.RMq-fin.ST

Catalog footprint

What is connected

17works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

17 published item(s)

preprint2022arXiv

Self-consistent formulations for stochastic nonlinear neuronal dynamics

Neural dynamics is often investigated with tools from bifurcation theory. However, many neuron models are stochastic, mimicking fluctuations in the input from unknown parts of the brain or the spiking nature of signals. Noise changes the dynamics with respect to the deterministic model; in particular bifurcation theory cannot be applied. We formulate stochastic neuronal dynamics in the Martin-Siggia-Rose de Dominicis-Janssen (MSRDJ) formalism and present the fluctuation expansion of the effective action and the functional renormalization group (fRG) as two systematic ways to incorporate corrections to the mean dynamics and time-dependent statistics due to fluctuations in the presence of nonlinear neuronal gain. To formulate self-consistency equations, we derive a fundamental link between the effective action in the Onsager-Machlup(OM) formalism, which allows the study of phase transitions, and the MSRDJ effective action, which is computationally advantageous. These results in particular allow the derivation of an OM effective action for systems with non-Gaussian noise. This approach naturally leads to effective deterministic equations for the first moment of the stochastic system; they explain how nonlinearities and noise cooperate to produce memory effects. Moreover, the MSRDJ formulation yields an effective linear system that has identical power spectra and linear response. Starting from the better known loopwise approximation, we then discuss the use of the fRG as a method to obtain self-consistency beyond the mean. We present a new efficient truncation scheme for the hierarchy of flow equations for the vertex functions by adapting the Blaizot, Méndez and Wschebor approximation from the derivative expansion to the vertex expansion. The methods are presented by means of the simplest possible example of a stochastic differential equation that has generic features of neuronal dynamics.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

Gell-Mann-Low criticality in neural networks

Criticality is deeply related to optimal computational capacity. The lack of a renormalized theory of critical brain dynamics, however, so far limits insights into this form of biological information processing to mean-field results. These methods neglect a key feature of critical systems: the interaction between degrees of freedom across all length scales, which allows for complex nonlinear computation. We present a renormalized theory of a prototypical neural field theory, the stochastic Wilson-Cowan equation. We compute the flow of couplings, which parameterize interactions on increasing length scales. Despite similarities with the Kardar-Parisi-Zhang model, the theory is of a Gell-Mann-Low type, the archetypal form of a renormalizable quantum field theory. Here, nonlinear couplings vanish, flowing towards the Gaussian fixed point, but logarithmically slowly, thus remaining effective on most scales. We show this critical structure of interactions to implement a desirable trade-off between linearity, optimal for information storage, and nonlinearity, required for computation.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Avoiding Ergodicity Problems in Lattice Discretizations of the Hubbard Model

The Hubbard model arises naturally when electron-electron interactions are added to the tight-binding descriptions of many condensed matter systems. For instance, the two-dimensional Hubbard model on the honeycomb lattice is central to the ab initio description of the electronic structure of carbon nanomaterials, such as graphene. Such low-dimensional Hubbard models are advantageously studied with Markov chain Monte Carlo methods, such as Hybrid Monte Carlo (HMC). HMC is the standard algorithm of the lattice gauge theory community, as it is well suited to theories of dynamical fermions. As HMC performs continuous, global updates of the lattice degrees of freedom, it provides superior scaling with system size relative to local updating methods. A potential drawback of HMC is its susceptibility to ergodicity problems due to so-called exceptional configurations, for which the fermion operator cannot be inverted. Recently, ergodicity problems were found in some formulations of HMC simulations of the Hubbard model. Here, we address this issue directly and clarify under what conditions ergodicity is maintained or violated in HMC simulations of the Hubbard model. We study different lattice formulations of the fermion operator and provide explicit, representative calculations for small systems, often comparing to exact results. We show that a fermion operator can be found which is both computationally convenient and free of ergodicity problems.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Renormalization of a Contact Interaction on a Lattice

Contact interactions can be used to describe a system of particles at unitarity, contribute to the leading part of nuclear interactions and are numerically non-trivial because they require a proper regularization and renormalization scheme. We explain how to tune the coefficient of a contact interaction between non-relativistic particles on a discretized space in 1, 2, and 3 spatial dimensions such that we can remove all discretization artifacts. By taking advantage of a latticized Lüscher zeta function, we can achieve a momentum-independent scattering amplitude at any finite lattice spacing.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The Semimetal-Mott Insulator Quantum Phase Transition of the Hubbard Model on the Honeycomb Lattice

We take advantage of recent improvements in the grand canonical Hybrid Monte Carlo algorithm, to perform a precision study of the single-particle gap in the hexagonal Hubbard model, with on-site electron-electron interactions. After carefully controlled analyses of the Trotter error, the thermodynamic limit, and finite-size scaling with inverse temperature, we find a critical coupling of $U_c/κ=3.834(14)$ and the critical exponent $zν=1.185(43)$. Under the assumption that this corresponds to the expected anti-ferromagnetic Mott transition, we are also able to provide a preliminary estimate $β=1.095(37)$ for the critical exponent of the order parameter. We consider our findings in view of the $SU(2)$ Gross-Neveu, or chiral Heisenberg, universality class. We also discuss the computational scaling of the Hybrid Monte Carlo algorithm, and possible extensions of our work to carbon nanotubes, fullerenes, and topological insulators.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Confirming the Existence of the strong CP Problem in Lattice QCD with the Gradient Flow

We calculate the electric dipole moment of the nucleon induced by the QCD theta term. We use the gradient flow to define the topological charge and use $N_f = 2+1$ flavors of dynamical quarks corresponding to pion masses of $700$, $570$, and $410$ MeV, and perform an extrapolation to the physical point based on chiral perturbation theory. We perform calculations at $3$ different lattice spacings in the range of $0.07~{\rm fm} < a < 0.11$ fm at a single value of the pion mass, to enable control on discretization effects. We also investigate finite size effects using $2$ different volumes. A novel technique is applied to improve the signal-to-noise ratio in the form factor calculations. The very mild discretization effects observed suggest a continuum-like behavior of the nucleon EDM towards the chiral limit. Under this assumption our results read $d_{n}=-0.00152(71)\ \barθ e~\text{fm}$ and $d_{p}=0.0011(10)\ \barθ e~\text{fm}$. Assuming the theta term is the only source of CP violation, the experimental bound on the neutron electric dipole moment limits $\left|\barθ\right| < 1.98\times 10^{-10}$ ($90\%$ CL). A first attempt at calculating the nucleon Schiff moment in the continuum resulted in $S_{p} = 0.50(59)\times 10^{-4}\ \barθ e~\text{fm}^3$ and $S_{n} = -0.10(43)\times 10^{-4}\ \barθ e~\text{fm}^3$.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Quantum Monte Carlo Calculations for Carbon Nanotubes

We show how lattice Quantum Monte Carlo can be applied to the electronic properties of carbon nanotubes in the presence of strong electron-electron correlations. We employ the path-integral formalism and use methods developed within the lattice QCD community for our numerical work. Our lattice Hamiltonian is closely related to the hexagonal Hubbard model augmented by a long-range electron-electron interaction. We apply our method to the single-quasiparticle spectrum of the (3,3) armchair nanotube configuration, and consider the effects of strong electron-electron correlations. Our approach is equally applicable to other nanotubes, as well as to other carbon nanostructures. We benchmark our Monte Carlo calculations against the two- and four-site Hubbard models, where a direct numerical solution is feasible.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Ab initio alpha-alpha scattering

Processes involving alpha particles and alpha-like nuclei comprise a major part of stellar nucleosynthesis and hypothesized mechanisms for thermonuclear supernovae. In an effort towards understanding alpha processes from first principles, we describe in this letter the first ab initio calculation of alpha-alpha scattering. We use lattice effective field theory to describe the low-energy interactions of nucleons and apply a technique called the adiabatic projection method to reduce the eight-body system to an effective two-cluster system. We find good agreement between lattice results and experimental phase shifts for S-wave and D-wave scattering. The computational scaling with particle number suggests that alpha processes involving heavier nuclei are also within reach in the near future.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Applying Twisted Boundary Conditions for Few-body Nuclear Systems

We describe and implement twisted boundary conditions for the deuteron and triton systems within finite-volumes using the nuclear lattice EFT formalism. We investigate the finite-volume dependence of these systems with different twists angles. We demonstrate how various finite-volume information can be used to improve calculations of binding energies in such a framework. Our results suggests that with appropriate twisting of boundaries, infinite-volume binding energies can be reliably extracted from calculations using modest volume sizes with cubic length $L\approx8-14$ fm. Of particular importance is our derivation and numerical verification of three-body analogue of `i-periodic' twist angles that eliminate the leading order finite-volume effects to the three-body binding energy.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Nuclear Lattice Simulations using Symmetry-Sign Extrapolation

Projection Monte Carlo calculations of lattice Chiral Effective Field Theory suffer from sign oscillations to a varying degree dependent on the number of protons and neutrons. Hence, such studies have hitherto been concentrated on nuclei with equal numbers of protons and neutrons, and especially on the alpha nuclei where the sign oscillations are smallest. Here, we introduce the "symmetry-sign extrapolation" method, which allows us to use the approximate Wigner SU(4) symmetry of the nuclear interaction to systematically extend the Projection Monte Carlo calculations to nuclear systems where the sign problem is severe. We benchmark this method by calculating the ground-state energies of the $^{12}$C, $^6$He and $^6$Be nuclei, and discuss its potential for studies of neutron-rich halo nuclei and asymmetric nuclear matter.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The nucleon electric dipole moment with the gradient flow: the $θ$-term contribution

We propose a new method to calculate electric dipole moments induced by the strong QCD $θ$-term. The method is based on the gradient flow for gauge fields and is free from renormalization ambiguities. We test our method by computing the nucleon electric dipole moments in pure Yang-Mills theory at several lattice spacings, enabling a first-of-its-kind continuum extrapolation. The method is rather general and can be applied for any quantity computed in a $θ$ vacuum. This first application of the gradient flow has been successful and demonstrates proof-of-principle, thereby providing a novel method to obtain precise results for nucleon and light nuclear electric dipole moments.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Two-Nucleon Systems in a Finite Volume: (II) 3S1-3D1 Coupled Channels and the Deuteron

The energy spectra of two nucleons in a cubic volume provide access to the two phase shifts and one mixing angle that define the S-matrix in the 3S1-3D1 coupled channels containing the deuteron. With the aid of recently derived energy quantization conditions for such systems, and the known scattering parameters, these spectra are predicted for a range of volumes. It is found that extractions of the infinite-volume deuteron binding energy and leading scattering parameters, including the S-D mixing angle at the deuteron pole, are possible from Lattice QCD calculations of two-nucleon systems with boosts of |P| <= 2pi sqrt{3}/L in volumes with 10 fm <~ L <~ 14 fm. The viability of extracting the asymptotic D/S ratio of the deuteron wavefunction from Lattice QCD calculations is discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Extracting Scattering Phase-Shifts in Higher Partial-Waves from Lattice QCD Calculations

Lüscher's method is routinely used to determine meson-meson, meson-baryon and baryon-baryon s-wave scattering amplitudes below inelastic thresholds from Lattice QCD calculations - presently at unphysical light-quark masses. In this work we review the formalism and develop the requisite expressions to extract phase-shifts describing meson-meson scattering in partial-waves with angular-momentum l<=6 and l=9. The implications of the underlying cubic symmetry, and strategies for extracting the phase-shifts from Lattice QCD calculations, are presented, along with a discussion of the signal-to-noise problem that afflicts the higher partial-waves.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quantile Mechanics II: Changes of Variables in Monte Carlo methods and GPU-Optimized Normal Quantiles

This article presents differential equations and solution methods for the functions of the form $Q(x) = F^{-1}(G(x))$, where $F$ and $G$ are cumulative distribution functions. Such functions allow the direct recycling of Monte Carlo samples from one distribution into samples from another. The method may be developed analytically for certain special cases, and illuminate the idea that it is a more precise form of the traditional Cornish-Fisher expansion. In this manner the model risk of distributional risk may be assessed free of the Monte Carlo noise associated with resampling. Examples are given of equations for converting normal samples to Student t, and converting exponential to hyperbolic, variance gamma and normal. In the case of the normal distribution, the change of variables employed allows the sampling to take place to good accuracy based on a single rational approximation over a very wide range of the sample space. The avoidance of any branching statement is of use in optimal GPU computations as it avoids the effect of {\it warp divergence}, and we give examples of branch-free normal quantiles that offer performance improvements in a GPU environment, while retaining the best precision characteristics of well-known methods. We also offer models based on a low-probability of warp divergence. Comparisons of new and old forms are made on the Nvidia Quadro 4000, GTX 285 and 480, and Tesla C2050 GPUs. We argue that in single-precision mode, the change-of-variables approach offers performance competitive with the fastest existing scheme while substantially improving precision, and that in double-precision mode, this approach offers the most GPU-optimal Gaussian quantile yet, and without compromise on precision for Monte Carlo applications, working twice as fast as the CUDA 4 library function with increased precision.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Quark mass variation constraints from Big Bang nucleosynthesis

We study the impact on the primordial abundances of light elements created by a variation of the quark masses at the time of Big Bang nucleosynthesis (BBN). In order to navigate through the particle and nuclear physics required to connect quark masses to binding energies and reaction rates in a model-independent way, we use lattice QCD data and a hierarchy of effective field theories. We find that the measured Helium-4 abundances put a bound of -1 % <~ d m_q/m_q <~ 0.7 % on a possible variation of quark masses. The effect of quark mass variations on the deuterium abundances can be largely compensated by changes of the baryon-to-photon ratio eta. Including bounds on the variation of eta coming from WMAP results and adding some additional assumptions further narrows the range of allowed values of d m_q/m_q.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Nucleon-Nucleon Scattering in a Harmonic Potential

The discrete energy-eigenvalues of two nucleons interacting with a finite-range nuclear force and confined to a harmonic potential are used to numerically reconstruct the free-space scattering phase shifts. The extracted phase shifts are compared to those obtained from the exact continuum scattering solution and agree within the uncertainties of the calculations. Our results suggest that it might be possible to determine the amplitudes for the scattering of complex systems, such as n-d, n-t or n-alpha, from the energy-eigenvalues confined to finite volumes using ab-initio bound-state techniques.

0 citations0 reviewsNo code linkedOpen access