skip to main content

This content will become publicly available on June 1, 2024

Title: Cauchy-type integral method for solving the linearized one-dimensional Vlasov-Poisson equation
We present a method for solving the linearized Vlasov-Poisson equation, based on analyticity properties of the equilibrium and initial condition through Cauchy-type integrals, that produces algebraic expressions for the distribution and field, i.e., the solution is expressed without integrals. Standard extant approaches involve deformations of the Bromwich contour that give erroneous results for certain physically reasonable configurations or eigenfunction expansions that are misleading as to the temporal structure of the solution. Our method is more transparent, lacks these defects, and predicts previously unrecognized behavior.  more » « less
Award ID(s):
Author(s) / Creator(s):
Publisher / Repository:
American Physical Society
Date Published:
Journal Name:
Physical Review E
Page Range / eLocation ID:
Subject(s) / Keyword(s):
["plasma waves","plasma kinetic theory","perturbative methods","electrostatic waves and oscillations"]
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    Abstract This paper presents the Eshelby’s tensor of a polygonal inclusion with a polynomial eigenstrain, which can provide an elastic solution to an arbitrary, convex inclusion with a continuously distributed eigenstrain by the Taylor series approximation. The Eshelby’s tensor for plane strain problem is derived from the fundamental solution of isotropic Green’s function with the Hadmard regularization, which is composed of the integrals of the derivatives of the harmonic and biharmonic potentials over the source domain. Using the Green’s theorem, they are converted to two line (contour) integrals over the polygonal cross section. This paper evaluates them by direct analytical integrals. Following Mura’s work, this paper formulates the method to derive linear, quadratic, and higher order of the Eshelby’s tensor in the polynomial form for arbitrary, convex polygonal shapes of inclusions. Numerical case studies were performed to verify the analytic results with the original Eshelby’s solution for a uniform eigenstrain in an ellipsoidal domain. It is of significance to consider higher order terms of eigenstrain for the polygon-shape inclusion problem because the eigenstrain distribution is generally non-uniform when Eshelby’s equivalent inclusion method is used. The stress disturbance due to a triangle particle in an infinite domain is demonstrated by comparison with the results of the finite element method (FEM). The present solution paves the way to accurately simulate the particle-particle, partial-boundary interactions of polygon-shape particles. 
    more » « less
  2. Abstract

    In clinical research and practice, landmark models are commonly used to predict the risk of an adverse future event, using patients' longitudinal biomarker data as predictors. However, these data are often observable only at intermittent visits, making their measurement times irregularly spaced and unsynchronized across different subjects. This poses challenges to conducting dynamic prediction at any post‐baseline time. A simple solution is the last‐value‐carry‐forward method, but this may result in bias for the risk model estimation and prediction. Another option is to jointly model the longitudinal and survival processes with a shared random effects model. However, when dealing with multiple biomarkers, this approach often results in high‐dimensional integrals without a closed‐form solution, and thus the computational burden limits its software development and practical use. In this article, we propose to process the longitudinal data by functional principal component analysis techniques, and then use the processed information as predictors in a class of flexible linear transformation models to predict the distribution of residual time‐to‐event occurrence. The measurement schemes for multiple biomarkers are allowed to be different within subject and across subjects. Dynamic prediction can be performed in a real‐time fashion. The advantages of our proposed method are demonstrated by simulation studies. We apply our approach to the African American Study of Kidney Disease and Hypertension, predicting patients' risk of kidney failure or death by using four important longitudinal biomarkers for renal functions.

    more » « less
  3. null (Ed.)
    Forward-backward stochastic differential equation (FBSDE) systems were introduced as a probabilistic description for parabolic type partial differential equations. Although the probabilistic behavior of the FBSDE system makes it a natural mathematical model in many applications, the stochastic integrals contained in the system generate uncertainties in the solutions which makes the solution estimation a challenging task. In this paper, we assume that we could receive partial noisy observations on the solutions and introduce an optimal filtering method to make a data informed solution estimation for FBSDEs. 
    more » « less
  4. Various of the single scale quantities in massless and massive QCD up to 3-loop order can be expressed by iterative integrals over certain classes of alphabets, from the harmonic polylogarithms to root-valued alphabets. Examples are the anomalous dimensions to 3-loop order, the massless Wilson coefficients and also different massive operator matrix elements. Starting at 3-loop order, however, also other letters appear in the case of massive operator matrix elements, the so called iterative non-iterative integrals, which are related to solutions based on complete elliptic integrals or any other special function with an integral representation that is definite but not a Volterra-type integral. After outlining the formalism leading to iterative non-iterative integrals,we present examples for both of these cases with the 3-loop anomalous dimension $\gamma^{(2)}_{qg}$ and the structure of the principle solution in the iterative non-interative case of the 3-loop QCD corrections to the $\rho$-parameter. 
    more » « less
  5. Abstract

    We present a calculation of the master integrals (MI’s) required for the calculation of the Electroweak corrections to$$gg\rightarrow \gamma \gamma $$ggγγproduction in which the process contains a light quark loop. The integrals can be broken down into five categories based on the flow of the heavy vector bosons throughout the loop. Three of the families are planar, and two are non-planar. We determine a canonical basis for each family which allows an efficient solution of the resulting differential equations via iterated integrals. We calculate the families in relevant physical kinematics and obtain an efficient numerical evaluation based on an implementation of Chen-iterated integrals.

    more » « less