Next-to-leading power two-loop soft functions for the Drell-Yan process at threshold

We calculate the generalized soft functions at $$ \mathcal{O} $$ O ($$ {\alpha}_s^2 $$ α s 2 ) at next-to-leading power accuracy for the Drell-Yan process at threshold. The operator definitions of these objects contain explicit insertions of soft gauge and matter fields, giving rise to a dependence on additional convolution variables with respect to the leading power result. These soft functions constitute the last missing ingredient for the validation of the bare factorization theorem to NNLO accuracy. We carry out the calculations by reducing the soft squared amplitudes into a set of canonical master integrals and we employ the method of differential equations to evaluate them. We retain the exact d-dimensional dependence of the convolution variables at the integration boundaries in order to regulate the fixed-order convolution integrals. After combining the soft functions with the relevant collinear functions, we perform checks of the results at the cross-section level against the literature and expansion-by-regions calculations, at NNLO and partly at N3LO, finding agreement.

Convolution Integral: How a Graphical Type of Solution Can Help Minimize Misconceptions

A physical system, Low Pass Filter (LPF) RC Circuit, which serves as an impulse response and a square wave input signal are utilized to derive the continuous time convolution (convolution integrals). How to set up the limits of integration correctly and how the excitation source convolves with the impulse response are explained using a graphical type of solution. This in turn, help minimize the students’ misconceptions about the convolution integral. Further, the effect of varying the circuit elements on the shape of the convolution output plot is presented allowing students to see the connection between a convolution integral and a physical system. PSpice simulation and experiment results are incorporated and are compared with those of the analytical solution associated with the convolution integral.

Application of the Efros Theorem to the Function Represented by the Inverse Laplace Transform of s−μ exp(−sν)

Using a special case of the Efros theorem which was derived by Wlodarski, and operational calculus, it was possible to derive many infinite integrals, finite integrals and integral identities for the function represented by the inverse Laplace transform. The integral identities are mainly in terms of convolution integrals with the Mittag–Leffler and Volterra functions. The integrands of determined integrals include elementary functions (power, exponential, logarithmic, trigonometric and hyperbolic functions) and the error functions, the Mittag–Leffler functions and the Volterra functions. Some properties of the inverse Laplace transform of s−μexp(−sν) with μ≥0 and 0<ν<1 are presented.

Factorization at subleading power and endpoint divergences in h → γγ decay. Part II. Renormalization and scale evolution

Building on the recent derivation of a bare factorization theorem for the b-quark induced contribution to the h → γγ decay amplitude based on soft-collinear effective theory, we derive the first renormalized factorization theorem for a process described at subleading power in scale ratios, where λ = mb/Mh « 1 in our case. We prove two refactorization conditions for a matching coefficient and an operator matrix element in the endpoint region, where they exhibit singularities giving rise to divergent convolution integrals. The refactorization conditions ensure that the dependence of the decay amplitude on the rapidity regulator, which regularizes the endpoint singularities, cancels out to all orders of perturbation theory. We establish the renormalized form of the factorization formula, proving that extra contributions arising from the fact that “endpoint regularization” does not commute with renormalization can be absorbed, to all orders, by a redefinition of one of the matching coefficients. We derive the renormalization-group evolution equation satisfied by all quantities in the factorization formula and use them to predict the large logarithms of order $$ {\alpha \alpha}_s^2{L}^k $$ αα s 2 L k in the three-loop decay amplitude, where $$ L=\ln \left(-{M}_h^2/{m}_b^2\right) $$ L = ln − M h 2 / m b 2 and k = 6, 5, 4, 3. We find perfect agreement with existing numerical results for the amplitude and analytical results for the three-loop contributions involving a massless quark loop. On the other hand, we disagree with the results of previous attempts to predict the series of subleading logarithms $$ \sim {\alpha \alpha}_s^n{L}^{2n+1} $$ ∼ αα s n L 2 n + 1 .

IMPROVED THERMAL SCATTERING DATA PREPARATION

Convenient access to accurate nuclear data, particularly data describing low-energy neutrons, is crucial for trustworthy simulations of thermal nuclear systems. Obtaining the scattering kernel for thermal neutrons (i.e., neutrons with energy ~1 eV or less) can be a difficult problem, since the neutron energy is not sufficient to break molecular bonds, and thus the neutrons must often interact with a much larger structure. The “scattering law” S(α; β), which is a function of unitless momentum α and energy β transfer, is used to relate the material’s phonon frequency distribution to the scattering kernel. LEAPR (a module of NJOY) and GASKET are two nuclear data processing codes that can be used to prepare the scattering law and use different approaches to approximate the same equations. LEAPR uses the “phonon expansion method” which involves iterative convolution. Iteratively solving convolution integrals is an expensive calculation to perform (to ease this calculation, LEAPR uses trapezoidal integration for the convolution). GASKET uses a more direct approach that, while avoiding the iterative convolutions, can become numerically unstable for some α; β combinations. When both methods are properly converged, they tend to agree quite well. The agreement and departure from agreement is presented here.

Fresnel diffraction of multiple disks on axis

Aims. We seek to study the Fresnel diffraction of external occulters that differ from a single mask in a plane. Such occulters have been used in previous space missions and are planned for the future ESA Proba 3 ASPIICS coronagraph. Methods. We studied the shading efficiency of double on-axis disks and generalized results to a 3D occulter. We used standard Fourier optics in an analytical approach. We show that the Fresnel diffraction of two and three disks on axis can be expressed using a Babinet-like approach. Results are obtained in the form of convolution integrals that can be written as Bessel-Hankel integrals; these are difficult to compute numerically for large Fresnel numbers found in solar coronagraphy. Results. We show that the shading efficiency of two disks is well characterized by the intensity of the residual Arago spot, a quantity that is easier to compute and therefore allows an interesting parametric study. Very simple conditions are derived for optimal sizes and positions of two disks to produce the darkest structure around the Arago spot. These conditions are inspired from empirical experiments performed in the sixties. A differential equation is established to give the optimal envelope for a multiple-disk occulter. The solution takes the form of a simple law, the approximation of which is a conical occulter, a shape already used in the SOHO Mission. Conclusions. In addition to quantifying expected results, the present study highlights unfortunate configurations of disks and spurious diffractions that may increase the stray light. Particular attention is paid to the possible issues of the future occulter spacecraft of ASPIICS.

Differentiation of the Mittag-Leffler Functions with Respect to Parameters in the Laplace Transform Approach

In this work, properties of one- or two-parameter Mittag-Leffler functions are derived using the Laplace transform approach. It is demonstrated that manipulations with the pair direct–inverse transform makes it far more easy than previous methods to derive known and new properties of the Mittag-Leffler functions. Moreover, it is shown that sums of infinite series of the Mittag-Leffler functions can be expressed as convolution integrals, while the derivatives of the Mittag-Leffler functions with respect to their parameters are expressible as double convolution integrals. The derivatives can also be obtained from integral representations of the Mittag-Leffler functions. On the other hand, direct differentiation of the Mittag-Leffler functions with respect to parameters produces an infinite power series, whose coefficients are quotients of the digamma and gamma functions. Closed forms of these series can be derived when the parameters are set to be integers.

Turbulence as a Network of Fourier Modes

Turbulence is the duality of chaotic dynamics and hierarchical organization of a field over a large range of scales due to advective nonlinearities. Quadratic nonlinearities (e.g., advection) in real space, translates into triadic interactions in Fourier space. Those interactions can be computed using fast Fourier transforms, or other methods of computing convolution integrals. However, more generally, they can be interpreted as a network of interacting nodes, where each interaction is between a node and a pair. In this formulation, each node interacts with a list of pairs that satisfy the triadic interaction condition with that node, and the convolution becomes a sum over this list. A regular wavenumber space mesh can be written in the form of such a network. Reducing the resolution of a regular mesh and combining the nearby nodes in order to obtain the reduced network corresponding to the low resolution mesh, we can deduce the reduction rules for such a network. This perspective allows us to develop network models as approximations of various types of turbulent dynamics. Various examples, such as shell models, nested polyhedra models, or predator–prey models, are briefly discussed. A prescription for setting up a small world variants of these models are given.

Large Amplitude Time Domain Seakeeping Simulations of KVLCC2 in Head Seas Taking Into Account Forward Speed Effect

Time-domain analysis is important for the development of performance-based criteria for the intact stability of ship (level 3 criteria in the 2nd Generation Intact Stability Criteria). It can be implemented to assess directly the vulnerability of ships against various modes of intact stability failure and ensure a sufficient level of safety. In this context, ELIGMOS, a novel time domain simulation code combining a maneuvering model and a seakeeping model, is under development. As the maneuvering model has been the subject of previous paper, the seakeeping part is presented herein and validated in terms of linear and nonlinear vertical motions in head seas. Heave and pitch motions are augmented near the resonant periods, testing the ability of the numerical tool to capture accurately the ship’s behaviour, as it plays a dominant role for her direct stability assessment afterwards. Surge motion is considered uncoupled and it is excluded from the system of equations, keeping the forward speed constant. Comparison regarding the seakeeping performance of KVLCC2 is presented against experimental results and results obtained by a 3D-potential flow, frequency domain software as well. Ship’s geometry is modelled by means of quadrilateral panels whilst radiation forces were incorporated by means of memory functions by adopting the well-known concept of convolution integrals. The ability of ELIGMOS to capture the effect of geometrical nonlinearities on the vertical motions is demonstrated through preliminary simulation of large amplitude motions of a C11-class containership.