Order In Expansion Research Articles

A novel method for highly linear gradient coil design based on stream function.

Magnetoencephalography (MEG) and magnetic resonance imaging (MRI) are non-invasive imaging techniques that offer effective means for disease diagnosis. A more straightforward and optimized method is presented for designing gradient coils which are pivotal parts of the above imaging systems. A novel design method based on stream function combining an optimization algorithm is proposed to obtain highly linear gradient coil. Two-dimensional Fourier expansion of the current field on the surface where the coil is located and the equipotential line of the expansion term superposition according to the number of turns of the coil are used to represent the coil shape. Particle swarm optimization is utilized to optimize the coil shape while linearity and field uniformity are used as parameters to evaluate the coil performance. Through this method, the main parameters such as input current distribution region, coil turns, desired magnetic field strength, expansion order and iteration times can be combined in a given solution space to optimize coil design. Simulation results show that the maximum linearity spatial deviation of the designed bi-planar x-gradient coil compared with that of target field method is reduced from 14% to 0.54%, and that of the bi-planar z-gradient coil is reduced from 8.98% to 0.52%. Similarly, that of the cylindrical x-gradient coil is reduced from 2% to 0.1%, and that of the cylindrical z-gradient coil is reduced from 0.87% to 0.45%. The similar results are found in the index of inhomogeneity error. Moreover, it has also been verified experimentally that the result of measured magnetic field is consist with simulated result. The proposed method provides a straightforward way that simplifies the design process and improves the linearity of designed gradient coil, which could be beneficial to realize better magnetic field in engineering applications.

Numerical simulation of wave propagation by using a hybrid method with an arbitrary order accuracy in both spatial and temporal approximations

This paper introduces an innovative numerical methodology designed to achieve high precision solution of acoustic wave propagation problem in isotropic material. During the temporal discretization process, the Krylov deferred correction (KDC) technique is employed, wherein a new variable is introduced to handle the second-order time derivative in the governing equations. An improved strategy is adopted for precisely implementing boundary conditions. Following this, the arbitrary-order generalized finite difference method (GFDM) is employed to simulate the transformed boundary value problem at each time step, enabling our framework to select the Taylor series expansion order arbitrarily. Ultimately, a hybrid numerical approach for wave problems is developed to achieve arbitrary-order accuracy in both spatial and temporal approximations. The developed scheme undergoes validation through four different numerical experiments in two- or three-dimension, wherein the numerical solutions obtained are compared with either analytical solutions or results from COMSOL software.

Estimates on the convergence of expansions at finite baryon chemical potentials

Convergence of three different expansion schemes at finite baryon chemical potentials, including the conventional Taylor expansion, the Padé approximants, and the T′ expansion proposed recently in lattice QCD simulations, have been investigated in a low energy effective theory within the functional renormalization group approach. It is found that the convergence of the T′ expansion and the Padé approximants is consistent with the conventional Taylor expansion, within the expansion orders considered in this work. Furthermore, we find that the consistent regions of the three different expansions are in agreement with the convergence radius of the Lee-Yang edge singularities. Published by the American Physical Society 2024

The NNLO soft function for N-jettiness in hadronic collisions

We compute the N-jettiness soft function in hadronic collisions to next-to-next-to-leading order (NNLO) in the strong-coupling expansion. Our calculation is based on an extension of the SoftSERVE framework to soft functions that involve an arbitrary number of lightlike Wilson lines. We present numerical results for 1-jettiness and 2-jettiness, and illustrate that our formalism carries over to a generic number of jets by calculating a few benchmark points for 3-jettiness. We also perform a detailed analytic study of the asymptotic behaviour of the soft-function coefficients at the edges of phase space, where one of the jets becomes collinear to another jet or beam direction, and comment on previous calculations of the N-jettiness soft function.

Quadrupolar radiation in de Sitter: displacement memory and Bondi metric

We obtain the closed form expression for the metric perturbation around de Sitter spacetime generated by a matter source below Hubble scale both in generalized harmonic gauge and in Bondi gauge up to quadrupolar order in the multipolar expansion, including both parities (i.e. both mass and current quadrupoles). We demonstrate that such a source causes a displacement memory effect close to future infinity that originates, in the even-parity sector, from a Λ-BMS transition between the two non-radiative regions of future infinity.

The NNLO gluon beam function for jet-veto resummation

We compute the gluon beam function for jet-veto resummation to next-to-next-to-leading order (NNLO) in the strong-coupling expansion. Our calculation is based on an automated framework that was previously used for the computation of the respective quark beam function, and which we significantly extended for the present calculation. In particular, the perturbative matching kernels are directly calculated in momentum space, without the need to perform an additional Mellin transform. We present results for both gluon and quark-initiated processes, which we cross-checked with an independent semi-analytical method that exploits the similarity of the beam functions to the more familiar case of transverse-momentum resummation. Our computation is relevant for jet-veto resummations at NNLL′ accuracy.

The one-loop bispectrum of galaxies in redshift space from the Effective Field Theory of Large-Scale Structure

We derive the Effective Field Theory of Large-Scale Structure kernels and counterterms for the one-loop bispectrum of dark matter and of biased tracers in real and redshift space. This requires the expansion of biased tracers up to fourth order in fluctuations. In the process, we encounter several subtleties related to renormalization. One is the fact that, in renormalizing the momentum, a local counterterm contributes non-locally. A second subtlety is related to the renormalization of local products of the velocity fields, which need to be expressed in terms of the renormalized velocity in order to preserve Galilean symmetry. We check that the counterterms we identify are necessary and sufficient to renormalize the one-loop bispectrum at leading and subleading order in the derivative expansion. The kernels that we originally present here have already been used for the first analyses of the one-loop bispectrum in BOSS data [1,2].

Multiscale formulation for materials composed by a saturated porous matrix and solid inclusions

Despite all the progress achieved in the characterization of heterogeneous materials by using multiscale paradigms based on the Representative Volume Element concept (RVE), there are still many aspects that demand ongoing development. We mention, for instance, in-homogeneous media with internal micro-structure comprising a mixture of components that require a dissimilar number/character of primary fields to describe their physical behavior. Micro-structures constituted by a saturated porous matrix (described through the pair displacement/pore pressure fields) endowed with impermeable non-porous solid inclusions (described only by the displacement field) are a typical example among many other practical applications. This new level of heterogeneity, between the set of primal descriptors for each micro-scale constituent, claims detailed revisions and novel adjustments in multiscale RVE-based theories.In this work, we present a first-order RVE-based formulation to address with the aforementioned phenomenology at a smaller length scale. At the macro-scale, the well-known poromechanics theory is preserved, where the constitutive response is provided by homogenization of the corresponding micro-scale problem. Two important attributes feature the model: (i) it is encompassed within the general framework posited by the “Method of Multiscale Virtual Power (MMVP)”, this leads to a variationally consistent methodology; (ii) our formulation preserves the objectivity with respect to the RVE-size (an indispensable but uncommon characteristic in RVE-based multiscale approaches for saturated porous media), due to it is based in a recent contribution of the authors (the so-called Selective Order Expansion technique “SOE”).The formulation is implemented following the squared Finite Element (FE2) numerical scheme. The potentialities of the model are demonstrated through a detailed series of numerical examples, including rigorous comparisons against the Direct Numerical Simulation (DNS) procedure.

Scaling results for charged sectors of near conformal QCD

We provide the leading near conformal corrections on a cylinder to the scaling dimension of the lowest-lying fixed isospin charge Q operators defined at the lower boundary of the quantum chromodynamics conformal window, Δ˜Q=Δ˜Q*+(mσ4πν)2QΔ3B1+(mπ(θ)4πν)4Q23(1−γ)B2+O(mσ4,mπ8,mσ2mπ4). Here, Δ˜Q/r is the classical ground state energy of the theory on R×Sr3 at fixed isospin charge while Δ˜Q* is the scaling dimension at the leading order in the large charge expansion. In the conformal limit mσ=mπ=0, the state-operator correspondence implies Δ˜Q=Δ˜Q*. The near-conformal corrections are expressed in powers of the dilaton and pion masses in units of the chiral symmetry breaking scale 4πν with the θ-angle dependence encoded directly in the pion mass. The characteristic Q-scaling is dictated by the quark mass operator anomalous dimension γ and the one characterizing the dilaton potential Δ. The coefficients Bi with i=1,2 depend on the geometry of the cylinder and properties of the nearby conformal field theory. Published by the American Physical Society 2024

A Complexity Dichotomy in Spatial Reasoning via Ramsey Theory

Constraint satisfaction problems (CSPs) for first-order reducts of finitely bounded homogeneous structures form a large class of computational problems that might exhibit a complexity dichotomy, P versus NP-complete. A powerful method to obtain polynomial-time tractability results for such CSPs is a certain reduction to polynomial-time tractable finite-domain CSPs defined over k -types, for a sufficiently large k . We give sufficient conditions when this method can be applied and apply these conditions to obtain a new complexity dichotomy for CSPs of first-order expansions of the basic relations of the well-studied spatial reasoning formalism RCC5. We also classify which of these CSPs can be expressed in Datalog. Our method relies on Ramsey theory; we prove that RCC5 has a Ramsey order expansion.

VSCF/VCI theory based on the Podolsky Hamiltonian.

While the vibrational spectra of semi-rigid molecules can be computed on approaches relying on the Watson Hamiltonian, floppy molecules or molecular clusters are better described by Hamiltonians, which are capable of dealing with any curvilinear coordinates. It is the kinetic energy operator (KEO) of these Hamiltonians, which render the correlated calculations relying on them rather costly. Novel implementation of vibrational self-consistent field theory and vibrational configuration interaction theory on the basis of the Podolsky Hamiltonian are reported, in which the inverse of the metric tensor, i.e., the G matrix, is represented by an n-mode expansion expressed in terms of polynomials. An analysis of the importance of the individual terms of the KEO with respect to the truncation orders of the n-mode expansion is provided. Benchmark calculations have been performed for the cis-HOPO and methanimine, H2CNH, molecules and are compared to experimental data and to calculations based on the Watson Hamiltonian and the internal coordinate path Hamiltonian.

Higher-derivative corrections to flavoured BPS black hole thermodynamics and holography

A Cardy-like regime of the four-dimensional superconformal index has been shown to be governed by ’t Hooft anomalies and to single out a large-N saddle carrying the Bekenstein-Hawking entropy of dual supersymmetric black holes in AdS5. For the universal index where no flavour fugacities are turned on, this correspondence has been improved by matching the first subleading corrections to the saddle-point action with the four-derivative corrections to the black hole action in minimal gauged supergravity, as well as the respective corrected entropies. Here, we extend this match by including flavour symmetries. We consider five-dimensional gauged supergravity with vector multiplet and four-derivative couplings, and provide an effective theory reproducing the ’t Hooft anomalies of the R- and flavour symmetries of generic holographic superconformal field theories at next-to-leading order in the large-N expansion. Then we focus on a specific model dual to ℂ3/ℤν quiver gauge theories, where the ’t Hooft anomaly coefficients receive simple but sufficiently generic corrections. In this model, we evaluate the four-derivative corrections to the on-shell action of the supersymmetric multi-charge black hole, showing agreement with the flavoured Cardy-like formula from the index. We give a prediction for the corrected entropy of the supersymmetric black hole and discuss the general validity of our results. Taking the limit of infinite AdS5 radius, we also obtain four-derivative corrections to the action and entropy of supersymmetric asymptotically flat black holes.

Categorical-symmetry resolved entanglement in conformal field theory

We propose a symmetry resolution of entanglement for categorical noninvertible symmetries (CaT-SREE) in (1+1)-dimensional conformal field theories. The definition parallels that of grouplike invertible symmetries, employing the concept of symmetric boundary states with respect to a categorical symmetry. Our examination extends to rational conformal field theories, where the behavior of CaT-SREE mirrors that of grouplike invertible symmetries. We find that CaT-SREE can be defined if there is no obstruction to gauging the categorical symmetry, as happens in the case of grouplike symmetries. We also provide instances of the breakdown of entanglement equipartition at the next-to-leading order in the cutoff expansion. Our findings shed light on how the interplay between conformal boundary conditions and categorical symmetries lead to specific patterns in the entanglement entropy. Published by the American Physical Society 2024

Energies and a gravitational charge for massive particles in general relativity

In this paper, we investigate relations or differences among various conserved quantities which involve the matter Energy Momentum Tensor (EMT) in general relativity. These quantities include the energy with Einstein’s pseudo EMT, the generalized Komar integral, or the ADM energy, all of which can be derived from Noether’s second theorem, as well as an extra conserved charge recently proposed in general relativity. For detailed analyses, we apply definitions of these charges to a system of free massive particles. We employ the post-Newtonian (PN) expansion to make physical interpretations. We find that the generalized Komar integral is not conserved at the first non-trivial order in the PN expansion due to non-zero contributions at spatial boundaries, while the energy with Einstein’s pseudo EMT at this order agrees with a total energy of massive particles with gravitational interactions through the Newtonian potential, and thus is conserved. In addition, this total energy is shown to be identical to the ADM energy not only at this order but also all orders in the PN expansion. We next calculate an extra conserved charge for the system of massive particles, at all orders in the PN expansion, which turns out to be a total number of particles. We call it a gravitational charge, since it is clearly different from the total energy. We finally discuss an implication from a fact that there exist two conserved quantities, energy and gravitational charge, in general relativity.

High-Energy Behavior of Scattering Amplitudes in Theories with Purely Virtual Particles

We study a class of renormalizable quantum field theories with purely virtual particles that exhibits nonrenormalizable behavior in the high-energy limit of scattering cross sections, which grow as powers of the center-of-mass energy squared and seems to violate unitarity bounds. We point out that the problem should be viewed as a violation of perturbativity, instead of unitarity, and show that the resummation of self energies fixes the issue. As an explicit example, we consider a class of O(N) theories at the leading order in the large-N expansion and show that the different quantization prescription of purely virtual particles takes care of the nonrenormalizable behavior, making the resummed cross sections to decrease at high energies and the amplitudes to satisfy the unitarity bounds. We compare the results to the case of theories with ghosts, where the resummation cannot change the behavior of cross sections due to certain cancellations in the high-energy expansion of the self energies. These results are particularly relevant for quantum gravity.

Third-order relativistic fluid dynamics at finite density in a general hydrodynamic frame

The motion of water is governed by the Navier–Stokes equations, which are complemented by the continuity equation to ensure local mass conservation. In this work, we construct the relativistic generalization of these equations through a gradient expansion for a fluid with a conserved charge in a curved d-dimensional spacetime. We adopt a general hydrodynamic frame and introduce the irreducible-structure (IS) algorithm, which is based on derivatives of the expansion scalar and the shear and vorticity tensors. By this method, we systematically generate all permissible gradients up to a specified order and derive the most comprehensive constitutive relations for a charged fluid, accurate to third-order in the gradient expansion. These constitutive relations are formulated to apply to ordinary (nonconformal) and conformally invariant charged fluids. Furthermore, we examine the frame dependence of the transport coefficients for a nonconformal charged fluid up to the third order in the gradient expansion. The frame dependence of the scalar, vector, and tensor parts of the constitutive relations is obtained in terms of the (field redefinitions of the) fundamental hydrodynamic variables. Managing the frame dependencies of the constitutive relations is challenging due to their non-linear character. However, in the linear regime, the higher-order transformations become tractable, enabling the identification of a set of frame-invariant coefficients. Subsequently, the equations obtained in the linear regime are solved in momentum space, yielding dispersion relations for shear, sound, and diffusive modes for a non-conformal charged fluid, expressed in terms of a set of frame-invariant transport coefficients.

Spin relaxation rate for baryons in a thermal pion gas

We study the relaxation dynamics of the spin polarization of baryons (nucleon and Λ baryon), in a thermal pion gas as a simple model of the hadronic phase of the QCD plasma produced in relativistic heavy-ion collisions. For this purpose, we formulate the quantum kinetic theory for the spin density matrix of baryons in the leading order of the gradient expansion. Considering the baryon-pion elastic scattering processes as the dominant interaction between baryons and thermal pions, we compute the spin relaxation rate of nucleons and Λ baryons in a pion gas up to temperature 200 MeV. In the case of nucleons, we evaluate the spin relaxation rate in the s-channel resonance approximation, based on the known experimental data on Δ resonances. We also estimate the spin relaxation rate for Λ baryons, based on experimental inputs and theoretical models for the low-energy Λπ scattering, including the chiral perturbation theory. Published by the American Physical Society 2024

Solution to a class of multistate Landau-Zener model beyond integrability conditions

We study a class of multistate Landau-Zener model which cannot be solved by integrability conditions or other standard techniques. By analyzing analytical constraints on its scattering matrix and performing fitting to results from numerical simulations of the Schrödinger equation, we find nearly exact analytical expressions of all its transition probabilities for specific parameter choices. We also determine the transition probabilities up to leading orders of series expansions in terms of the inverse sweep rate (namely, in the diabatic limit) for general parameter choices. We further show that this model can describe a Su-Schrieffer-Heeger chain with couplings changing linearly in time. Our work presents a new route, i.e., analytical constraint plus fitting, to analyze those multistate Landau-Zener models which are beyond the applicability of conventional solving methods.

Vanishing of nonlinear tidal Love numbers of Schwarzschild black holes

It is well known that asymptotically flat Schwarzschild black holes in general relativity in four spacetime dimensions have vanishing induced linear tidal response. We extend this result beyond linear order for the polar sector, by solving the static nonlinear Einstein equations for the perturbations of the Schwarzschild metric and computing the quadratic corrections to the electric-type tidal Love numbers. After explicitly performing the matching with the point-particle effective theory at leading order in the derivative expansion, we show that the Love number couplings remain zero at higher order in perturbation theory.

Revisiting the Numerical Evaluation and Visualization of the Gravity Fields of Asteroid 4769 Castalia Using Polyhedron and Harmonic Expansions Models

For the convenience of comparison with previous literature, the gravity calculations are revisited for the Asteroid 4769 Castalia, but with extensions on its surface and on intersecting planes and spheres around it, using the polyhedron and harmonic expansion methods with different order and degree for different cases, especially including the gravitational accelerations inside the asteroid, which did not appear at all before. In these evaluations, a few different facts of the these methods and results are revealed, such as the fact that gravity diverges when the position radius is less than the mean radius from harmonic-expansion method, and the maximum gravity is not at the deep valley and mountain top. For a surface that intersects the asteroid, the maximum gravity on it is at the intersection lines between the asteroid surface and the spheres or planes. This means that on the sphere and the plane, the gravities inside and outside the asteroid are smaller than the gravity on the intersection, i.e., on the surface. Some analyses of these conclusions are given with many examples with different radii of the sphere and with different order and degree harmonic expansion models for the above asteroid surface and surrounding spheres. It is interesting to note that very few researchers know that the polyhedral method can also be used to calculate the gravity inside an asteroid with just some modifications of the code. Some special gravity figures on surface and planes inside the asteroid Castalia are computed and made for the first time. The calculations also include tangential gravity, potential, and gravitational slope on surface. Specifically, we find that the overall mean gravitational slope could be one kind of indicator of the density of an asteroid. The minimum overall mean slope happens when the asteroid density is about 2.9 g/cm3, which is much larger than a usually assumed value between 1.7 and 2.5 for asteroid Castalia when its period is 4.07 h, since rotation period should be a more accurate parameter than its estimated density. These conclusions about this typical prolate-like asteroid could be a benchmark for analyzing other similar asteroids.