Asymptotic Form Research Articles

Non-Gaussian diffusive fluctuations in Dirac fluids

Dirac fluids-interacting systems obeying particle-hole symmetry and Lorentz invariance-are among the simplest hydrodynamic systems; they have also been studied as effective descriptions of transport in strongly interacting Dirac semimetals. Direct experimental signatures of the Dirac fluid are elusive, as its charge transport is diffusive as in conventional metals. In this paper, we point out a striking consequence of fluctuating relativistic hydrodynamics: The full counting statistics (FCS) of charge transport is highly non-Gaussian. We predict the exact asymptotic form of the FCS, which generalizes a result previously derived for certain interacting integrable systems. A consequence is that, starting from quasi-one-dimensional nonequilibrium initial conditions, charge noise in the hydrodynamic regime is parametrically enhanced relative to that in conventional diffusive metals.

Non-equilibrium dynamics of symmetry-resolved entanglement and entanglement asymmetry: exact asymptotics in Rule 54*

Abstract Symmetry resolved entanglement and entanglement asymmetry are two measures of quantum correlations sensitive to symmetries of the system. Here we discuss their non-equilibrium dynamics in the Rule 54 cellular automaton, a simple, yet interacting, integrable model. Both quantities can be expressed in terms of the more analytically tractable ‘charged moments’, i.e. traces of powers of a suitably deformed density matrix, via a replica trick. We express them in terms of a tensor network, which we contract in space using a system of local algebraic relations. This gives the asymptotic form for the charged moments, valid in the regime of large but finite time that is shorter than all the relevant subsystem sizes. In this regime the charge moments decay exponentially with the rate given by the leading solution to a cubic equation.

Two-dimensional Coulomb gas on an elliptic annulus

Abstract It is well-known that two-dimensional Coulomb gases at 
a special inverse temperature $\beta = 2$ can be analyzed 
by using the orthogonal polynomial method borrowed from the 
theory of random matrices. In this paper, such Coulomb gas molecules 
are studied when they are distributed on an elliptic annulus, and the 
asymptotic forms of the molecule correlation functions in the thermodynamic 
limit are evaluated. For that purpose, two-dimensional orthogonality 
relations of the Chebyshev polynomials on an elliptic annulus are utilized. 

Analytical expression for continuum–continuum transition amplitude of hydrogen-like atoms with angular-momentum dependence

Abstract Attosecond chronoscopy typically utilises interfering two-photon transitions to access the phase information. Simulating these two-photon transitions is challenging due to the continuum–continuum transition term. The hydrogenic approximation within second-order perturbation theory has been widely used due to the existence of analytical expressions of the wave functions. So far, only (partially) asymptotic results have been derived, which fail to correctly describe the low-kinetic-energy behaviour, especially for high angular-momentum states. Here, we report an analytical expression that overcomes these limitations. It is based on the Appell’s F 1 function and uses the confluent hypergeometric function of the second kind as the intermediate state. We show that the derived formula quantitatively agrees with the numerical simulations using the time-dependent Schrödinger equation for various angular-momentum states, which improves the accuracy compared to the other analytical approaches that were previously reported. Furthermore, we give an angular-momentum-dependent asymptotic form of the outgoing wavefunction and the corresponding continuum–continuum dipole transition amplitudes.

Deflection of light by a compact object with electric charge and magnetic dipole in Einstein–Born–Infeld gravity

We consider the bending of light around a compact astrophysical object with both the electric field and the magnetic field in Einstein–Born–Infeld theory. From the null geodesic of a light ray passing a massive object with electric charge and magnetic dipole, the effective metric was obtained from the light cone condition reflecting the nonlinear electromagnetic effects. We found the asymptotic form of the effective metric up to the first order in gravitational constant G and Born–Infeld parameter 1/β2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$1/\\beta ^2$$\\end{document} on the equatorial plane. Then we compute the bending angle of light from the geodesic equation. The result includes particular cases where only one type of field is present taking the appropriate limits.

Drawing of fibres composed of shear-thinning or shear-thickening fluid with internal holes

We explore the drawing of a shear-thinning or shear-thickening thread with an axisymmetric hole that evolves due to axial drawing, inertia and surface tension effects. The stress is assumed to be proportional to the shear rate raised to the $n$ th power. The presence of non-Newtonian rheology and surface tension forces acting on the hole introduces radial pressure gradients that make the derivation of long-wavelength equations significantly more challenging than either a Newtonian thread with a hole or shear-thinning and shear-thickening threads without a hole. In the case of weak surface tension, we determine the steady-state profiles. Our results show that for negligible inertia the hole size at the exit becomes smaller as $n$ is decreased (i.e. strong shear-thinning effects) above a critical draw ratio, but surprisingly the opposite is true below this critical draw ratio. We determine an accurate estimate of the critical draw ratio and also discuss how inertia affects this process. We further show that the dynamics of hole closure is dominated by a different limit, and we determine the asymptotic forms of the hole closure process for shear-thinning and shear-thickening fluids with inertia. A linear instability analysis is conducted to predict the onset of draw resonance. We show that increased shear thinning, surface tension and inlet hole size all act to destabilise the flow. We also show that increasing shear-thinning effects reduce the critical Reynolds number required for unconditional stability. Our study provides valuable insights into the drawing process and its dependence on the physical effects.

Integrability and Inverse Scattering Transform of the Modified Benjamin-Ono Equation

Abstract This paper presents a B{"a}cklund transformation, the Lax representation, and conserved quantities for the modified Benjamin--Ono equation. The initial problem of the modified Benjamin--Ono equation on the line was studied by the inverse scattering transform method, presenting a nonlocal Riemann--Hilbert problem in the inverse problem to reconstruct the explicit potential function. Further, the exact $N$-soliton solutions and long--time asymptotic behavior are provided, respectively. We also graphically show that the propagation of soliton solutions is consistent with the result of large-time asymptotic forms. As a revised version, the mBO equation incorporates the solution features of the BO equation, and its solutions are presented in logarithmic form, which seeks to describe related physical phenomena more accurately. Our results will contribute to further exploration of physical and mathematical problems related to the mBO equation.

Boosting the HP filter for trending time series with long-range dependence

This article extends recent asymptotic theory developed for the Hodrick Prescott (HP) filter and boosted HP (bHP) filter to long-range dependent time series that have fractional Brownian motion (fBM) limit processes after suitable standardization. Under general conditions, it is shown that the asymptotic form of the HP filter is a smooth curve, analogous to the finding in Phillips and Jin for integrated time series and series with deterministic drifts. Boosting the filter using the iterative procedure suggested in Phillips and Shi leads under well-defined rate conditions to a consistent estimate of the fBM limit process or the fBM limit process with an accompanying deterministic drift when that is present. A stopping criterion is used to automate the boosting algorithm, giving a data-determined method for practical implementation. The theory is illustrated in simulations and two real data examples that highlight the differences between simple HP filtering and the use of boosting. The analysis is assisted by employing a uniformly and almost surely convergent trigonometric series representation of fBM.

One-dimensional Fermi polaron after a kick: Two-sided singularity of the momentum distribution, Bragg reflection and other exact results

A mobile impurity particle immersed in a quantum fluid forms a polaron - a quasiparticle consisting of the impurity and a local disturbance of the fluid around it. We ask what happens to a one-dimensional polaron after a kick, i.e. an abrupt application of a force that instantly delivers a finite impulse to the impurity. In the framework of an integrable model describing an impurity in a one-dimensional gas of fermions or hard-core bosons, we calculate the distribution of the polaron momentum established when the post-kick relaxation is over. A remarkable feature of this distribution is a two-sided power-law singularity. It emerges due to one of two processes. In the first process, the whole impulse is transferred to the polaron, without creating phonon-like excitations of the fluid. In the second process, the impulse is shared between the polaron and the center-of-mass motion of the fluid, again without creating any fluid excitations. The latter process is, in fact, a Bragg reflection at the edge of the emergent Brillouin zone. We carefully analyze the conditions for each of the two processes. The asymptotic form of the distribution in the vicinity of the singularity is derived.

Heavy states in 3d gravity and 2d CFT

We discuss correlators of light fields in heavy states in AdS3 gravity and holographic 2d CFTs. In the bulk, the propagator of free fields in AdS backgrounds containing a conical defect or a BTZ black hole can be obtained by solving a wave equation, as well as by the method of images. On the boundary, these geometries are sourced by heavy operator insertions, and the propagator is dual to a heavy-light (HHLL) correlator. By matching its expansion in Virasoro blocks to our bulk results, we determine the OPE coefficients of all contributing states in both the s and t channels. In the s channel, these states are excitations of the light field on top of the heavy state, and their OPE coefficients are the amplitudes to create them. The t-channel OPE is dominated by the Virasoro vacuum block, but there is also an infinite family of light two-particle states that contribute to the correlator. The OPE coefficients that couple these states to heavy operators represent their expectation values in heavy backgrounds. We determine them exactly, derive their asymptotic form at large twist, and discuss their behavior near and above the BTZ threshold, where they become thermal one-point functions.

Quantum Dynamics Studies of the Li + Na2 (V = 0, j = 0) → Na + NaLi Reaction on a New Neural Network Potential Energy Surface.

The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na2Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na2Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na2 (v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10-7 to 80 cm-1. In the high collision energy range of 8 to 80 cm-1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward-backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10-7 to 8 cm-1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10-10 cm3/s.

Viscous fluid flow and heat transfer past a permeable wall jet with convective boundary conditions

PurposeThe steady laminar wall jet flow over a stretching/shrinking surface in the presence of lateral suction or injection with a convective boundary condition is considered.Design/methodology/approachThe partial differential equations for mass, momentum and energy conservation are changed to the system of ordinary differential equations through similarity solution transformations. Solutions, both numerical and asymptotic, to these similarity equations are found in some new ranges of parameters in the governing equations.FindingsThe equations are solved both asymptotically and numerically for a range of the transpiration parameter S and the flow parameter λ given in Mahros et al. (2023), thus greatly extending the range of these previous solutions. Asymptotic solutions for both large and small values of the Prandtl number σ are derived, showing good agreement with additional numerical integrations. It should be noted that in Mahros et al. (2023), only the case when σ=1 was treated. A solution for large λ when S=1 is obtained, showing a different asymptotic form to the case when S>0 in Mahros et al. (2023). Multiple solutions were seen by them for S<0 and the nature of the lower solution branch as S→0 from below is discussed. The question as to whether the lower branch solutions join as λ>0 when S<0 is resolved through obtaining an asymptotic solution λ small.Originality/valueThe accuracy of the solutions has been checked through a detailed comparison between the solutions obtained numerically and analytically, where excellent agreement has been found. This study is important for scientists working in the area of jet flows to become familiar with the flow properties and behaviour of jets.

Deflection of light by a Reissner–Nordström black hole and Painlevé VI equation

We consider the bending angle of the trajectory of a photon incident from and deflected to infinity around a Reissner–Nordström black hole. We treat the bending angle as a function of the squared reciprocal of the impact parameter and the squared electric charge of the background normalized by the mass of the black hole. It is shown that the bending angle satisfies a system of two inhomogeneous linear partial differential equations with polynomial coefficients. This system can be understood as an isomonodromic deformation of the inhomogeneous Picard–Fuchs equation satisfied by the bending angle in the Schwarzschild spacetime, where the deformation parameter is identified as the background electric charge. Furthermore, the integrability condition for these equations is found to be a specific type of the Painlevé VI equation that allows an algebraic solution. We solve the differential equations both at the weak and strong deflection limits. In the weak deflection limit, the bending angle is expressed as a power series expansion in terms of the squared reciprocal of the impact parameter and we obtain the explicit full-order expression for the coefficients. In the strong deflection limit, we obtain the asymptotic form of the bending angle that consists of the divergent logarithmic term and the finite O(1) term supplemented by linear recurrence relations which enable us to straightforwardly derive higher order coefficients. In deriving these results, the isomonodromic property of the differential equations plays an important role. Lastly, we briefly discuss the applicability of our method to other types of spacetimes such as a spinning black hole.

Fermion-soliton scattering in a modified CP1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbb{C}\\mathbb{P}^{1}$$\\end{document} model

The scattering of fermions in the background field of a topological soliton of the modified (2+1)\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$(2 + 1)$$\\end{document}-dimensional CP1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbb{C}\\mathbb{P}^{1}$$\\end{document} model is studied here both analytically and numerically. Unlike the original CP1\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$\\mathbb{C}\\mathbb{P}^{1}$$\\end{document} model, the Lagrangian of the modified model contains a potential term. Due to this, a dilatation zero mode of the topological soliton disappears, which results in stability of the fermion-soliton system. The symmetry properties of the fermion-soliton system are established, and the asymptotic forms of fermionic radial wave functions are studied. Questions related to the bound states of the fermion-soliton system are then discussed. General formulae describing the scattering of fermions are presented. The amplitudes of the fermion-soliton scattering are obtained in an analytical form within the framework of the Born approximation, and their symmetry properties and asymptotic forms are studied. The energy levels of the fermionic bound states and the partial phase shifts of fermionic scattering are obtained by numerical methods, and the ultrarelativistic limits of the partial phase shifts are found.

Bifurcation of solutions of Ornstein–Zernike equation in the gas–liquid coexistence region of simple fluids

One objective of this paper is to investigate bifurcation of solutions of Ornstein–Zernike equation in liquid state theory for the gas–liquid coexistence region of simple fluids. The analysis uses a discrete form of the Ornstein–Zernike equation in real space, together with a closure relation which provides nearly self consistent thermodynamic properties of fluids governed by interaction potentials like the Lennard-Jones potential. Accurate asymptotic forms of correlation functions are incorporated in the real space algorithm so as to mitigate the influence of cutoff length used in the discrete approximations. The global error reduction of this algorithm is similar to that of Simpson’s rule. It is found that there is no spinodal on the vapor side on a sub-critical isotherm. However, there is a density at which compressibility vanishes, but this lies on a nonphysical branch. There are fold bifurcation points on the vapor and liquid sides together with a ‘no-real-solution-region’ in between. This is followed by a region of negative compressibility extending up to the spinodal on the liquid side. Importantly, it is found that in regions where compressibility is non-positive, there are infinite number of solutions to OZE which violate an integral constraint that is needed for applications to thermodynamics. The emerging picture is somewhat similar to that in the case of the hypernetted-chain-closure. It is also found that complex solutions connect the physical solutions that exist outside the ‘no-real-solution-region’ as well as in region of negative compressibility. The second objective of the paper, in the light of bifurcation scenario, is to evolve an interpretation of the coupling-parameter expansion of correlation functions pertaining to the region of negative compressibility. It is shown that the expansion provides an accurate description of results of the restricted ensembles, implied in mean field theories and simulation methods. Bifurcation of singlet density to an in-homogeneous distribution at a spinodal is also established using a simple equation, which uses the direct correlation function. Thus the paper re-establishes the fact that restriction to a homogeneous singlet density is the root cause of all inconsistencies found in the gas–liquid coexistence region.

Amplification of the anomalous scaling in the Kazantsev-Kraichnan model with finite-time correlations and spatial parity violation.

By using the field theoretic renormalization group technique together with the operator product expansion, simultaneous influence of the spatial parity violation and finite-time correlations of an electrically conductive turbulent environment on the inertial-range scaling behavior of correlation functions of a passively advected weak magnetic field is investigated within the corresponding generalized Kazantsev-Kraichnan model in the second order of the perturbation theory (in the two-loop approximation). The explicit dependence of the anomalous dimensions of the leading composite operators on the fixed point value of the parameter that controls the presence of finite-time correlations of the turbulent field as well as on the parameter that drives the amount of the spatial parity violation (helicity) in the system is found even in the case with the presence of the large-scale anisotropy. In accordance with the Kolmogorov's local isotropy restoration hypothesis, it is shown that, regardless of the amount of the spatial parity violation, the scaling properties of the model are always driven by the anomalous dimensions of the composite operators near the isotropic shell. The asymptotic (inertial-range) scaling form of all single-time two-point correlation functions of arbitrary order of the passively advected magnetic field is found. The explicit dependence of the corresponding scaling exponents on the helicity parameter as well as on the parameter that controls the finite-time velocity correlations is determined. It is shown that, regardless of the amount of the finite-time correlations of the given Gaussian turbulent environment, the presence of the spatial parity violation always leads to more negative values of the scaling exponents, i.e., to the more pronounced anomalous scaling of the magnetic correlation functions. At the same time, it is shown that the stronger the violation of spatial parity, the larger the anomalous behavior of magnetic correlations.

Random vector functional link networks for function approximation on manifolds

The learning speed of feed-forward neural networks is notoriously slow and has presented a bottleneck in deep learning applications for several decades. For instance, gradient-based learning algorithms, which are used extensively to train neural networks, tend to work slowly when all of the network parameters must be iteratively tuned. To counter this, both researchers and practitioners have tried introducing randomness to reduce the learning requirement. Based on the original construction of Igelnik and Pao, single layer neural-networks with random input-to-hidden layer weights and biases have seen success in practice, but the necessary theoretical justification is lacking. In this study, we begin to fill this theoretical gap. We then extend this result to the non-asymptotic setting using a concentration inequality for Monte-Carlo integral approximations. We provide a (corrected) rigorous proof that the Igelnik and Pao construction is a universal approximator for continuous functions on compact domains, with approximation error squared decaying asymptotically like O(1/n) for the number n of network nodes. We then extend this result to the non-asymptotic setting, proving that one can achieve any desired approximation error with high probability provided n is sufficiently large. We further adapt this randomized neural network architecture to approximate functions on smooth, compact submanifolds of Euclidean space, providing theoretical guarantees in both the asymptotic and non-asymptotic forms. Finally, we illustrate our results on manifolds with numerical experiments.

Variational Bayesian Neural Networks via Resolution of Singularities

In this work, we advocate for the importance of singular learning theory (SLT) as it pertains to the theory and practice of variational inference in Bayesian neural networks (BNNs). To begin, we lay to rest some of the confusion surrounding discrepancies between downstream predictive performance measured via the test log predictive density and the variational objective. Next, we use the SLT-corrected asymptotic form for singular posterior distributions to inform the design of the variational family itself. Specifically, we build upon the idealized variational family introduced in Bhattacharya, Pati, and Plummer which is theoretically appealing but practically intractable. Our proposal takes shape as a normalizing flow where the base distribution is a carefully-initialized generalized gamma. We conduct experiments comparing this to the canonical Gaussian base distribution and show improvements in terms of variational free energy and variational generalization error. Supplemental appendices and code for the article are available online.

Kinematical higher-twist corrections in γ*→M1M2γ : Charged meson production

Generalized distribution amplitudes (GDAs) of mesons can be probed by the reactions e−e+→M1M2γ, which are accessible at electron-positron colliders such as BESIII and Belle II. After discussing the neutral meson production case in our previous paper, we discuss here the complementary case of the charged meson (M+M−) production, where one can extract the complete information on GDAs from the interference of the amplitudes of the two competing processes where the photon is emitted either in the initial or in the final state. Considering the importance of the charged meson production, we present a complete expression for the interference term of the cross section, which is experimentally accessible thanks to its charge conjugation specific property. We adopt two types of models for leading-twist ππ GDAs to estimate the size of the interference term in the process e−e+→π+π−γ numerically, namely, a model extracted from previous experimental results on γ*γ→π0π0 at Belle and the asymptotic form predicted by QCD evolution equations. We include in the calculation the kinematical power suppressed (sometimes called kinematical higher-twist) corrections up to 1/Q2 for the helicity amplitudes. Both models of GDAs indicate that the kinematical corrections are not negligible for the interference term of the cross section measured at BESIII, thus it is necessary to include them if we try to extract the GDAs precisely. On the other side, the kinematical corrections are very tiny for the measurements at Belle II, and the leading twist-two formula of the interference term will be good enough to describe the charge conjugation odd part of the differential cross section. Published by the American Physical Society 2024

Natural orbitals and their occupation numbers for free anyons in the magnetic gauge

We investigate the properties of natural orbitals and their occupation numbers of the ground state of two noninteracting anyons characterized by the fractional statistics parameter α and confined in a harmonic trap. We work in the boson magnetic gauge where the anyons are modeled as composite bosons with magnetic flux quanta attached to their positions. We derive an asymptotic form of the weakly occupied natural orbitals, and show that their corresponding (ordered descendingly) occupation numbers decay according to the power law n−(4+2α), where n is the index of the natural orbital. We find remarkable numerical agreement of the theory with the natural orbitals and their occupation numbers computed from the spectral decomposition of the system's wave function. We explain that the same results apply to the fermion magnetic gauge. Published by the American Physical Society 2024