Order In Perturbative Expansion Research Articles

Pole decomposition of BFKL eigenvalue at zero conformal spin and the real part of digamma function

We consider the powers of leading order eigenvalue of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation at zero conformal spin. Using reflection identities of harmonic sums we demonstrate how involved generalized polygamma functions are introduced by pole separation of a rather simple digamma function. This generates higher weight generalized polygamma functions at any given order of perturbative expansion. As a byproduct of our analysis we develop a general technique for calculating powers of the real part of digamma function in a pole separated form.

Analytical modeling of variable thickness cylindrical shallow shells using extended Kantorovich method

The objective of the current paper is to provide novel numerical and analytical solutions for the variable cross section shallow shells under static pressure loads. First, the governing equations and the corresponded boundary conditions of the system are derived using the principle of virtual work. Then using the extended Kantorovich method (EKM), the partial differential equations of the system are transformed into two sets of variable coefficient ordinary differential equations, which are solved iteratively using both a novel numerical finite difference scheme, as well as an analytical perturbation tool. The convergence of the developed finite difference scheme is analyzed and it is observed that considering 600 nodes in the normalized domain is sufficient for achieving highly accurate results. Similarly, analytical simulations reveal that a fourth order perturbation expansion leads to an error of less than 0.2% in calculating the pick transverse response. The convergence of the EKM is also investigated and it is observed that the maximum transverse displacement can be predicted with less than 0.01% error using only two iterates. Additionally, variations of the stress and strain components in the mid-surface are studied using contour plots. Finally, parametric studies are carried out to characterize the effect of the shell's aspect ratio on the maximum transverse deflection. The analytical and numerical findings reported in this paper are closely validated with FE results as well as other available data in the literature.

Generalized master equation for charge transport in a molecular junction: Exact memory kernels and their high order expansion.

We derive a set of generalized master equations (GMEs) to study charge transport dynamics in molecular junctions using the Nakajima-Zwanzig-Mori projection operator approach. In the new GME, time derivatives of population on each quantum state of the molecule, as well as the tunneling current, are calculated as the convolution of time non-local memory kernels with populations on all system states. The non-Markovian memory kernels are obtained by combining the hierarchical equations of motion (HEOM) method and a previous derived Dyson relation for the exact kernel. A perturbative expansion of these memory kernels is then calculated using the extended HEOM developed in our previous work [M. Xu et al., J. Chem. Phys. 146, 064102 (2017)]. By using the resonant level model and the Anderson impurity model, we study properties of the exact memory kernels and analyze convergence properties of their perturbative expansions with respect to the system-bath coupling strength and the electron-electron repulsive energy. It is found that exact memory kernels calculated from HEOM exhibit short memory times and decay faster than the population and current dynamics. The high order perturbation expansion of the memory kernels can give converged results in certain parameter regimes. The Padé and Landau-Zener resummation schemes are also found to give improved results over low order perturbation theory.

Transverse momentum dependent parton densities in processes with heavy quark generations

To study the heavy quark production processes, we use the transverse momentum dependent (TMD, or unintegrated) gluon distribution function in a proton obtained recently using the Kimber-Martin-Ryskin prescription from the Bessel-inspired behavior of parton densities at small Bjorken $x$ values. We obtained a good agreement of our results with the latest HERA experimental data for reduced cross sections $\sigma^{c\overline{c}}_{\rm red}(x,Q^2)$ and $\sigma^{b\overline{b}}_{\rm red}(x,Q^2)$, and also for deep inelastic structure functions $F_2^c(x,Q^2)$ and $F_2^b(x,Q^2)$ in a wide range of $x$ and $Q^2$ values. Comparisons with the predictions based on the Ciafaloni-Catani-Fiorani-Marchesini evolution equation and with the results of conventional pQCD calculations performed at first three orders of perturbative expansion are presented.

Perturbation approaches for describing dipolar fluids and electrolyte solutions.

This work proposes perturbation approaches for describing dipolar fluids as well as model and aqueous electrolyte solutions. The electrostatic pair potentials are split into short- and long-ranged contributions, whereas a third order perturbation expansion is applied for the short-ranged potentials. This circumvents the problem of divergent correlation integrals. The dipolar perturbation terms are represented through a [2,1]-Padé approximation to resum the poorly convergent series. For the remaining charge-charge and charge-dipole contributions, we present a new approximant, which provides a (quasi)linear dependence of the Helmholtz energy. The underlying correlation integrals are adjusted to results from molecular simulations. The long-ranged contribution to the electrostatic interactions is treated through an analytic expression developed by Rodgers and Weeks [J. Chem. Phys. 131, 244108 (2010)]. Theoretical predictions of our perturbation theory are compared to results from a widely used integral equation theory, namely, the mean spherical approximation, and we find that our perturbation theory provides much more accurate results. Furthermore, the theory shows some quantities in rather good agreement with reference data, namely, Helmholtz energies, internal energies, and densities at higher densities of solutions. Limitations of the approach, however, are observed for several other partial molar quantities, such as the mean activity coefficient.

An efficient spectral method for numerical time-dependent perturbation theory.

We develop the Fourier-Laplace Inversion of the Perturbation Theory (FLIPT), a novel numerically exact "black box" method to compute perturbative expansions of the density matrix with rigorous convergence conditions. Specifically, the FLIPT method is extremely well-suited to simulate multiphoton pulsed laser experiments with complex pulse shapes. The n-dimensional frequency integrals of the nth order perturbative expansion are evaluated numerically using tensor products. The N-point discretized integrals are computed in O(N2) operations, a significant improvement over the O(Nn) scaling of standard quadrature methods.

Disentangling Electronic and Vibrational Effects in the Prediction of Band Shapes for Singlet–Triplet Transitions

The band shape of phosphorescent emission of benzophenone has been computed by using the first order perturbative expansion of singlet and triplet states with the spin–orbit coupling operator as perturbation and by evaluating Franck–Condon integrals with an efficient strategy for handling the whole set of vibrational coordinates. The computed band shape compares well with the experimental one, showing that modern computational tools yield reliable spin–orbit couplings to be used for evaluating the rates of singlet–triplet transitions in modern optoelectronic devices.

Synchronization of stochastic hybrid oscillators driven by a common switching environment.

Many systems in biology, physics, and chemistry can be modeled through ordinary differential equations (ODEs), which are piecewise smooth, but switch between different states according to a Markov jump process. In the fast switching limit, the dynamics converges to a deterministic ODE. In this paper, we suppose that this limit ODE supports a stable limit cycle. We demonstrate that a set of such oscillators can synchronize when they are uncoupled, but they share the same switching Markov jump process. The latter is taken to represent the effect of a common randomly switching environment. We determine the leading order of the Lyapunov coefficient governing the rate of decay of the phase difference in the fast switching limit. The analysis bears some similarities to the classical analysis of synchronization of stochastic oscillators subject to common white noise. However, the discrete nature of the Markov jump process raises some difficulties: in fact, we find that the Lyapunov coefficient from the quasi-steady-state approximation differs from the Lyapunov coefficient one obtains from a second order perturbation expansion in the waiting time between jumps. Finally, we demonstrate synchronization numerically in the radial isochron clock model and show that the latter Lyapunov exponent is more accurate.

Higher Order Algorithm for Solving Lambert’s Problem

This work presents a high-order perturbation expansion method for solving Lambert’s problem. The necessary condition for the problem is defined by a fourth-order Taylor expansion of the terminal error vector. The Taylor expansion partial derivative models are generated by Computational Differentiation (CD) tools. A novel derivative enhanced numerical integration algorithm is presented for computing nonlinear state transition tensors, where only the equation of motion is coded. A high-order successive approximation algorithm is presented for inverting the problems nonlinear necessary condition. Closed-form expressions are obtained for the first, second,third, and fourth order perturbation expansion coefficients. Numerical results are presented that compare the convergence rate and accuracy of first-through fourth-order expansions. The initial p-iteration starting guess is used as the Lambert’s algorithm initial condition. Numerical experiments demonstrate that accelerated convergence is achieved for the second-, third-, and fourth-order expansions, when compared to a classical first-order Newton method.

Rotations versus perturbative expansions for calculating neutrino oscillation probabilities in matter

We further develop a simple and compact technique for calculating the three flavor neutrino oscillation probabilities in uniform matter density. By performing additional rotations instead of implementing a perturbative expansion we significantly decrease the scale of the perturbing Hamiltonian and therefore improve the accuracy of zeroth order. We explore the relationship between implementing additional rotations and that of performing a perturbative expansion. Based on our analysis, independent of the size of the matter potential, we find that the first order perturbation expansion can be replaced by two additional rotations and a second order perturbative expansion can be replaced by one more rotation. Numerical tests have been applied and all the exceptional features of our analysis have been verified.

Exact generator and its high order expansions in time-convolutionless generalized master equation: Applications to spin-boson model and excitation energy transfer

The time-convolutionless (TCL) quantum master equation provides a powerful tool to simulate reduced dynamics of a quantum system coupled to a bath. The key quantity in the TCL master equation is the so-called kernel or generator, which describes effects of the bath degrees of freedom. Since the exact TCL generators are usually hard to calculate analytically, most applications of the TCL generalized master equation have relied on approximate generators using second and fourth order perturbative expansions. By using the hierarchical equation of motion (HEOM) and extended HEOM methods, we present a new approach to calculating the exact TCL generator and its high order perturbative expansions. The new approach is applied to the spin-boson model with different sets of parameters, to investigate the convergence of the high order expansions of the TCL generator. We also discuss circumstances where the exact TCL generator becomes singular for the spin-boson model, and a model of excitation energy transfer in the Fenna-Matthews-Olson complex.

All-Order Volume Conjecture for Closed 3-Manifolds from Complex Chern–Simons Theory

We propose an extension of the recently-proposed volume conjecture for closed hyperbolic 3-manifolds, to all orders in perturbative expansion. We first derive formulas for the perturbative expansion of the partition function of complex Chern-Simons theory around a hyperbolic flat connection, which produces infinitely-many perturbative invariants of the closed oriented 3-manifold. The conjecture is that this expansion coincides with the perturbative expansion of the Witten-Reshetikhin-Turaev invariants at roots of unity $q=e^{2 \pi i/r}$ with $r$ odd, in the limit $r \to \infty$. We provide numerical evidence for our conjecture.

Dark states in spin-polarized transport through triple quantum dot molecules

We study the spin-polarized transport through a triple quantum dot molecule weakly coupled to ferromagnetic leads. The analysis is performed by means of the real-time diagrammatic technique including up to the second order of perturbation expansion with respect to the tunnel coupling. The emphasis is put on the impact of dark states on spin-resolved transport characteristics. It is shown that the interplay of coherent population trapping and cotunneling processes results in a highly nontrivial behavior of the tunnel magnetoresistance, which can take negative values. Moreover, a super-Poissonian shot noise is found in transport regimes where the current is blocked by the formation of dark states, which can be additionally enhanced by spin-dependence of tunneling processes, depending on magnetic configuration of the device. The mechanisms leading to those effects are thoroughly discussed.

Analytical study on motion and shape of creeping Boger drops falling through viscoelastic media

In the present study, an analytical solution based on a double perturbation method is obtained to simulate the creeping motion of a Boger drop descending in an immiscible viscoelastic quiescent medium. The perturbation parameters are considered to be the Deborah number defined based on the relaxation times of interior and exterior media and the Capillary number. The Oldroyd-B model is implemented as a constitutive equation for both interior and exterior phases. The analytical solutions are obtained up to second order of perturbation expansion for both phases. The motion of the droplet in the surrounding medium is considered to be creeping. Effects of different parameters including, Deborah numbers, viscosity ratios, elasto-capillary number, and viscosity ratio between the drop and external flows are investigated on the shape, motion, and stream function in detail. In contrast to the generally spherical shape of creeping drops in Newtonian regimes, it is shown that the effect of elasticity in a non-Newtonian scenario can feature prolate or oblate shapes. It is shown that an increment in elastic properties of interior fluid is followed by decrement in the terminal velocity of the droplet, while an increment in elasticity of exterior fluid has a contrary effect and increases the terminal velocity. The results have shown suitable agreements with previously reported experimental results.

Investigation of Bose-Einstein Condensates in q-Deformed Potentials with First Order Perturbation Theory

The Gross-Pitaevskii equation, which is the governor equation of Bose-Einstein condensates, is solved by first order perturbation expansion under various q-deformed potentials. Stationary probability distributions reveal one and two soliton behavior depending on the type of the q-deformed potential. Additionally a spatial shift of the probability distribution is found for the dark soliton solution, when the q parameter is changed.

Hamiltonian approach to second order gauge invariant cosmological perturbations

In view of growing interest in tensor modes and their possible detection, we clarify the definition of tensor modes up to 2nd order in perturbation theory within the Hamiltonian formalism. Like in gauge theory, in cosmology the Hamiltonian is a suitable and consistent approach to reduce the gauge degrees of freedom. In this paper we employ the Faddeev-Jackiw method of Hamiltonian reduction. An appropriate set of gauge invariant variables that describe the dynamical degrees of freedom may be obtained by suitable canonical transformations in the phase space. We derive a set of gauge invariant variables up to 2nd order in perturbation expansion and for the first time we reduce the 3rd order action without adding gauge fixing terms. In particular, we are able to show the relation between the uniform-$\phi$ and Newtonian slicings, and study the difference in the definition of tensor modes in these two slicings.

Second-order singular pertubative theory for gravitational lenses

The extension of the singular perturbative approach to the second order is presented in this paper. The general expansion to the second order is derived. The second order expansion is considered as a small correction to the first order expansion. Using this approach it is demonstrated that the second order expansion is reducible to a first order expansion via a re-definition of the first order pertubative fields. Even if in practice the second order correction is small the reducibility of the second order expansion to the first order expansion indicates a degeneracy problem. In general this degeneracy is hard to break. A useful and simple second order approximation is the thin source approximation which offers a direct estimation of the correction. The practical application of the corrections derived in this paper are illustrated by using an elliptical NFW lens model. The second order pertubative expansion provides a noticeable improvement, even for the simplest case of thin source approximation. To conclude it is clear that for accurate modelisation of gravitational lenses using the perturbative method the second order perturbative expansion should be considered. In particular an evaluation of the degeneracy due to the second order term should be performed, for which the thin source approximation is particularly useful.

Dynamics of a mass–spring–beam with 0:1:1 internal resonance using the analytical and continuation method

The nonlinear harmonic response of an autoparametric system comprised of a linear oscillator with a vertically attached flexural beam is investigated and the capability of the beam as a vibration absorber is assessed. A weak torsional spring is used for constraining the rotation of the beam giving rise to an almost non-flexural rotational mode with a low frequency. The system parameters are also tuned to enforce the zero-to-one-to-one internal resonance condition. The Lagrange’s formulation accompanied by the assumed-mode method is used to derive the discretized equations based on the order three nonlinear Euler–Bernoulli beam theory. An analytical solution is developed based on the method of multiple scales where the generalized coordinate corresponding to the non-flexural rotational mode is approximated by higher order perturbation expansion than the other coordinates, due to much larger contribution of the non-flexural rotation to the response. Comprehensive response and bifurcation analysis are performed using analytical and direct numerical solutions. The results are obtained for vertically-aligned and also initially inclined beams and various interesting behaviors are recognized for different non-dimensional system parameters. Different types of bifurcations such as the Pitch-fork, Hopf, Period-doubling and symmetry breaking bifurcations are observed in the solution of slow-flow equations and some of them are found to be beneficial for vibration absorption in a limited range of excitation amplitudes and frequencies.

Convergence of high order perturbative expansions in open system quantum dynamics.

We propose a new method to directly calculate high order perturbative expansion terms in open system quantum dynamics. They are first written explicitly in path integral expressions. A set of differential equations are then derived by extending the hierarchical equation of motion (HEOM) approach. As two typical examples for the bosonic and fermionic baths, specific forms of the extended HEOM are obtained for the spin-boson model and the Anderson impurity model. Numerical results are then presented for these two models. General trends of the high order perturbation terms as well as the necessary orders for the perturbative expansions to converge are analyzed.

Design and construction of ultra-relativistic collision PV panel and its application into building sector to mitigate total energy demand

To meet total energy demand for a building, a sustainable technology has been developed to double up the photon production from a single one by ultra-relativistic collision into the solar panel, which is actually a part of building's exterior skin curtain wall. Subsequently, mathematical analysis has been performed considering the reaction between two bare ions which is activated by quark-lepton complementarily particle's interaction into the solar panel. To confirm this pair photon production, a further mathematical modeling has been implemented by assuming the creation of pair electron into a bound ionic state (e+e−) which is employed by two approaches (i) the equivalent photon approximation, and (ii) the first order perturbation and multipole expansion of wave functions. With the analysis of these approaches, a detailed calculations have been conducted considering two bound free e+e− pairs as well as of bound free e+e− and free μ+μ− pairs created by collisions of bare ions which suggest that the double photon production processes from a single photon is feasible into the ultra-relativistic collision solar panel that interestingly allow to double up the energy production simultaneously into the solar panel. The energy production efficiently, is thereafter, been estimated which revealed that if only a mere of 12.5% of building exterior curtain wall is used as ultra-relativistic collision capable solar panel it will meet the total energy demand for a building which is also 100% environmentally friendly.