Higher Order Correction Terms Research Articles

Wrinkling of differentially growing bilayers with similar film and substrate moduli

Growth-induced surface wrinkling in constrained bilayers comprising a thin film attached to a thick substrate is a canonical model for understanding pattern formation in many biological systems. While the bilayer model has received much prior attention, the nonlinear behaviour for arrangements with similar film and substrate properties, or substrate growth that outpaces film growth, remains poorly understood. This paper therefore focuses on these cases in which the substrate’s elasticity dominates surface wrinkling. By combining analytical modelling and finite element simulations incorporating advanced path-following techniques, we characterise critical wrinkling states and explore the initial and advanced post-buckling behaviour for a wide range of film-to-substrate modulus and growth rate ratios. It is shown that the classical leading-order asymptotic expression for the critical strain is not of sufficient accuracy for film-to-substrate modulus ratios below 50, and higher order correction terms are presented. Leveraging a weakly nonlinear analytical model, we introduce a phase diagram categorising stable and unstable post-buckling regimes. Advanced path-following simulations are employed to unveil the evolution of wrinkling patterns, from initial sinusoidal instabilities to complex formations involving period doubling, quadrupling, and eventual surface creasing. These comprehensive phase diagrams, parameterised by film-to-substrate modulus and growth rate ratios, provide new insights into the rich dynamics of surface wrinkling. Finally, we demonstrate the existence of multi-stability in the advanced post-buckling regimes for bilayers experiencing substrate-dominated growth. This work contributes to the understanding of the mechanics underlying pattern formation in growing bilayers over the entire parameter regime, with potential implications for explaining biological morphogenesis and informing the development of novel diagnostic tools and artificial flexible electronic skins.

Reduced model and nonlinear analysis of localized instabilities of residually stressed cylinders under axial stretch

In this paper, we present a dimensional reduction to obtain a one-dimensional model to analyze localized necking or bulging in a residually stressed circular cylindrical solid. The nonlinear theory of elasticity is first specialized to obtain the equations governing the homogeneous deformation. Then, to analyze the nonhomogeneous part, we include higher-order correction terms of the axisymmetric displacement components leading to a three-dimensional form of the total potential energy functional. Details of the reduction to the one-dimensional form are given. We focus on a residually stressed Gent material and use numerical methods to solve the governing equations. Two loading conditions are considered. First, the residual stress is maintained constant, while the axial stretch is used as the loading parameter. Second, we keep the pre-stretch constant and monotonically increase the residual stress until bifurcation occurs. We specify initial conditions, find the critical values for localized bifurcation, and compute the change in radius during localized necking or bulging growth. Finally, we optimize material properties and use the one-dimensional model to simulate necking or bulging until the Maxwell values of stretch are reached.

The Potential Roles of Transacylation in Intracellular Lipolysis and Related Qssa Approximations

Fatty acids (FAs) are crucial energy metabolites, signalling molecules, and membrane building blocks for a wide range of organisms. Adipose triglyceride lipase (ATGL) is the first and presumingly most crucial regulator of FA release from triacylglycerols (TGs) stored within cytosolic lipid droplets. However, besides the function of releasing FAs by hydrolysing TGs into diacylglycerols (DGs), ATGL also promotes the transacylation reaction of two DG molecules into one TG and one monoacylglycerol molecule. To date, it is unknown whether DG transacylation is a coincidental byproduct of ATGL-mediated lipolysis or whether it is physiologically relevant. Experimental evidence is scarce since both, hydrolysis and transacylation, rely on the same active site of ATGL and always occur in parallel in an ensemble of molecules. This paper illustrates the potential roles of transacylation. It shows that, depending on the kinetic parameters but also on the state of the hydrolytic machinery, transacylation can increase or decrease downstream products up to 80% respectively 30%. We provide an extensive asymptotic analysis including quasi-steady-state approximations (QSSA) with higher order correction terms and provide numerical simulation. We also argue that when assessing the validity of QSSAs one should include parameter sensitivity derivatives. Our results suggest that the transacylation function of ATGL is of biological relevance by providing feedback options and altogether stability to the lipolytic machinery in adipocytes.

Impact of higher-order effects on the dynamics of soliton solutions in the (3+1)D cubic-quintic-septic complex Ginzburg–Landau equation with higher-order dispersion terms

We investigate the influence of the higher-order correction term of the nonlinear refractive index and the nonlinear-gain absorption on the dynamics of soliton solutions in a three-dimensional (3D) dissipative medium described by the higher-order (3+1)-dimensional cubic-quintic-septic complex Ginzburg–Landau [(3+1)D CQS-CGL] equation with the viscous (spectral-filtering) term, self steepening, Raman effect, dispersion terms up to six, diffraction and cubic-quintic-septic nonlinearities. By mean of direct numerical simulation of the dynamical model, we show that these two parameters can have a significant impact on the dynamics of solitons govern by the above mentioned equation, since their formation deeply depends of the value of each of the above cited parameters. Some interesting particular cases have been detected concerning the combined action of the two parameters or the isolated action of each of them. Among those solutions, we have bell-shaped, snail and creeping dissipative solitons. Also, we have found that the propagation distance z is another important parameter which has allowed us to appreciate the impact of each parameter. At a given distance z or for different values of z, the formation of soliton depends on the value of the higher-order correction term and the nonlinear-gain absorption involved. We have also shown that the above solutions can be self-trapped over a huge propagation distance even in the presence of random perturbations. The discovered nonlinear dependence of the higher-order correction term and the nonlinear-gain absorption on the soliton solution can also be used for the better understanding of the (3+1)D dissipative light bullets in a doped nonlinear kerr medium.

Identification and Modeling Method of Longitudinal Stall Aerodynamic Parameters of Civil Aircraft Based on Improved Kirchhoff Stall Aerodynamic Model

The stall aerodynamic model based on Kirchhoff’s theory of flow separation is widely used in the identification and modeling of stall aerodynamic parameters. However, it has two defects. First, its model structure is significantly different from the pre-stall model used for a small attack angle, meaning the identification results cannot be combined with the pre-stall model to form the full flight envelope model. Second, the pitching moment model, which is used in conjunction with the Kirchhoff lift model, cannot accurately describe the aircraft stall pitching moment characteristics. To ensure the compatibility of the two models, this paper proposes a method to determine some unknown parameters in the stall model. The mechanism for the pitching moment generation of the aircraft stall is analyzed, and two high-order correction terms are added to the pitching moment model to better describe the longitudinal stall aerodynamic characteristics. Based on the identification of aerodynamic parameters, a longitudinal stall aerodynamic modeling method used for the aircraft stall process is developed. The identification and simulation validation results based on the quasi-steady stall flight data of a civil aircraft show that the improved stall model can accurately describe the quasi-steady stall pitching moment. The established stall aerodynamic model can accurately characterize the longitudinal quasi-steady stall aerodynamic characteristics of aircraft under different stall degrees.

Optical soliton propagation under the influence of higher order nonlinearities

We numerically investigate the propagation of optical soliton under the influence of higher-order nonlinearities in the anomalous dispersion regime by solving the complex cubic–quintic Ginzburg–Landau equation which is a generalized nonlinear Schrödinger equation. The nonlinear gain-absorption processes, the higher-order correction term to the intensity-dependence of refractive index and intra-pulse Raman scattering are considered in the generalized nonlinear Schrödinger equation. With the help of the numerical simulation, we note that the higher-order dispersive term and nonlinearities can perturb a soliton, resulting in a rich variety of solitons such as period-doubling bifurcation of dissipative soliton, Raman-induced spectral shifts, the formation of multi-branch pulsating solitons, and the creation of pulsating soliton with different periodicity.

Response to “Comment on ‘Even-order harmonic generation from nonlinear Thomson backscatter in a tightly focused Gaussian laser pulse’” [Phys. Plasmas 29, 114703 (2022)

In a recent Comment, Chang et al. [Phys. Plasmas 29, 114703 (2022)] have studied the electron dynamics for an electron accelerating in a tightly focused Gaussian laser pulse with the higher order correction terms, and it is found that the initial position of the electron is the reason why the electron is pushed to the −z axis at the end. In this Reply, the electron dynamics and its nonlinear Thomson backscattering in a tightly focused Gaussian laser pulse with the higher order correction terms are presented, and it is found that the result is consistent with the lowest order approximation case in our paper [Hong et al., Phys. Plasmas 29, 043102 (2022)]. Meanwhile, it is also found that when the longitudinal deceleration effect of the ponderomotive force on the electron introduced by the falling part of the tightly focused laser pulse is greater than the longitudinal acceleration effect of the rising part, the electron that is initially stationary will be slightly pushed to the −z axis at the end, which can well explain the “return” phenomenon of the electron in the longitudinal direction both in our paper and in the Comment.

Stability of nonparaxial gap-soliton bullets in waveguide gratings

We investigate the formation and propagation of gap-soliton bullets in nonlinear periodic waveguides at frequencies close to the gap for Bragg reflection beyond the paraxial approximation. Using a multiple-scales analysis, we derive a two-dimensional (2D) nonlinear Schrödinger equation with higher-order correction terms that consider the nonparaxial regimes in the slowly-varying envelope approximation. In addition, a fully numerical simulation of the newly derived model equation demonstrates that the mutual balancing between Kerr, dimensionality, higher-order dispersions and nonparaxiality allows shape-preserving propagation of gap-soliton bullets in a grating waveguide.

Nonlinear characteristics and corrections of near-field underwater explosion shock waves

The shock wave characteristics within the near-field are one of the most challenging aspects of understanding an underwater explosion. The latest numerical and experimental techniques were utilized to investigate the near-field pressure distribution and decay features after a shock disturbance. The governing equations in the numerical simulation were discretized with a fifth-order weighted essentially non-oscillatory scheme in space and a third-order Runge–Kutta scheme in time, and multi-medium interactions were defined and resolved via the modified ghost fluid method. The test system consisted of a synchronized high-speed framing camera and polyvinylidene fluoride (PVDF) sensors. Three identical spherical composition B charges were examined under the same test conditions, and the raw data from the high-speed camera were processed with edge detection and circle fitting techniques. The comparison showed that the high-speed camera image data, the PVDF signals, and the numerical computation results were highly consistent with each other. Higher-order correction terms were added to the pressure peak distribution model and the pressure decay model as nonlinear corrections based on further comprehensive and insightful analysis of the verified results. The corrected models not only fit with the near-field data but had better accuracy under the far-field condition as well.

Higher-order correction terms to the nonlinear amplification or absorption, the nonlinear refractive index, and the intrapulse Raman scattering.

In this paper we have investigated through the numerical solution of the basic equation as well as through the dynamic model the influence of higher-order correction terms to the nonlinear amplification (absorption) and to the nonlinear refractive index on the self-frequency shift of Raman dissipative solitons. We have found a nonlinear dependence of the self-frequency shift of Raman dissipative solitons on the parameter describing intrapulse Raman scattering in the presence of the saturation of the nonlinear gain. With the increase of the absolute value of the saturation of the nonlinear gain, the maximum absolute value of the frequency shift decreases and its position moves to larger values of the parameter describing intrapulse Raman scattering. The increase in the value of the nonlinear gain leads to an increase in the maximum absolute value of the frequency shift, without changing its position. We have also observed the nonlinear dependence of the absolute value of the frequency shift on the parameter describing intrapulse Raman scattering in the presence of higher-order correction term to the nonlinear refractive index. The discovered nonlinear dependence of the self-frequency shift on the value of the saturation of the nonlinear gain as well as on the higher-order correction term to the nonlinear refractive index can be used for the better understanding and control of the spectral characteristics of Raman dissipative solitons. The dynamic model correctly describes all the features of the observed phenomena.

A new conservative/dissipative time integration scheme for nonlinear mechanical systems

We present a conservative/dissipative time integration scheme for nonlinear mechanical systems. Starting from a weak form, we derive algorithmic forces and velocities that guarantee the desired conservation/dissipation properties. Our approach relies on a collection of linearly constrained quadratic programs defining high order correction terms that modify, in the minimum possible way, the classical midpoint rule so as to guarantee the strict energy conservation/dissipation properties. The solution of these programs provides explicit formulas for the algorithmic forces and velocities which can be easily incorporated into existing implementations. Similarities and differences between our approach and well-established methods are discussed as well. The approach, suitable for reduced-order models, finite element models, or multibody systems, is tested and its capabilities are illustrated by means of several examples.

Pair and many-body interactions between ligated Au nanoparticles.

We report the results of molecular dynamics simulations of the properties of a pseudo-atom (united atom) model of dodecane thiol ligated 5-nm diameter gold nanoparticles (AuNPs) in a vacuum as a function of ligand coverage and particle separation in three states of aggregation, namely, the isolated AuNPs, the isolated pair of AuNPs, and a square lattice of four AuNPs. Our calculations show that the ligand density along a radius emanating from the core of an isolated AuNP has the same gross features for all values of the coverage; it oscillates around a constant value up to a distance along the chain corresponding to the position of the fourth pseudo-atom and then smoothly decays to zero, reflecting both the restricted conformations of the chain near the core surface and the larger numbers of conformations available further from the core. Interaction between two AuNPs generates changes in the ligand distributions of each. We examine the structure and general shape of the ligand envelope as a function of the coverage and demonstrate that the equilibrium structure of the envelope and the deformation of that envelope generated by interaction between the NPs are coverage-dependent so that the shape, depth, and position of the minimum of the potential of mean force display a systematic dependence on the ligand coverage. We propose an accurate analytical description of the calculated potential of mean force as a function of a set of parameters that scale linearly with the ligand coverage. Noting that the conformational freedom of the ligands implies that multiparticle induced deviations from additivity of the pair potential of mean force are likely important; we define and calculate a "bond stretching" effective pair potential of mean force for a square lattice of particles that contains, implicitly, both the three- and four-NP contributions. We find that the bond stretching effective pair potential of mean force in this cluster has a different minimum and a different well depth from the isolated pair potential of mean force. Previous work has found that the three-particle contribution to deviation from pair additivity is monotonically repulsive, whereas we find that the combined three- and four-particle contributions have an attractive well, implying that the three- and four-particle contributions are of comparable magnitude but opposite sign, thereby suggesting that even higher order correction terms likely play a significant role in the behavior of dense assemblies of many nanoparticles.

Thermodynamics of black holes with higher order corrected entropy

For analyzing the thermodynamical behavior of two well-known black holes, such as Reissner–Nordström – anti-de Sitter (RN-AdS) black hole with global monopole and f(R) black hole, we consider the higher order logarithmic corrected entropy. We develop various thermodynamical properties, such as entropy, specific heat, pressure, and Gibbs and Helmhotz free energies for both black holes in the presence of corrected entropy. A versatile study on the stability of black holes is made by using various frameworks, such as the ratio of heat capacities (γ), grand canonical and canonical ensembles, and phase transition in view of higher order logarithmic corrected entropy. It is observed that both black holes exhibit more stability (locally as well as globally) for growing values of cosmological constant and higher order correction terms.

Analysis of multichannel optical soliton collisions with asymmetric energy in strongly dispersion managed transmission networks considering higher-order effects

We study the interaction properties of optical soliton pulses propagating in a two-channel wavelength-division multiplexed strongly dispersion managed communication system. We analyze the effects of third order dispersion, Raman scattering and self-steepening using an ordinary differential equations model obtained by a variational method applied to the Generalized Nonlinear Schrodinger Equation. The validity of the model is assessed against the integration of the full nonlinear partial differential equations. The variational equations are initially solved for a single pulse in order to identify the launching parameters for each pulse in the first DM cell of the system. One pulse per wavelength is then injected in the transmission link. We then systematically study the evolution of multiplexed soliton trains and analyse the transmission system in terms of residual frequency shifts and interchannel interactions as the map strength is varied, focusing on the ratio of dissimilar peak powers in a broad range of dispersion difference values, concluding that the transmission characteristics improve by using specific values of unequal energies and considering higher-order correction terms. The work presented is an extension of previous analysis where isolated intra-channel two-pulse interactions had been addressed considering pulses with dissimilar energies.

Generalisation of Gilbert damping and magnetic inertia parameter as a series of higher-order relativistic terms

The phenomenological Landau–Lifshitz–Gilbert (LLG) equation of motion remains as the cornerstone of contemporary magnetisation dynamics studies, wherein the Gilbert damping parameter has been attributed to first-order relativistic effects. To include magnetic inertial effects the LLG equation has previously been extended with a supplemental inertia term; the arising inertial dynamics has been related to second-order relativistic effects. Here we start from the relativistic Dirac equation and, performing a Foldy–Wouthuysen transformation, derive a generalised Pauli spin Hamiltonian that contains relativistic correction terms to any higher order. Using the Heisenberg equation of spin motion we derive general relativistic expressions for the tensorial Gilbert damping and magnetic inertia parameters, and show that these tensors can be expressed as series of higher-order relativistic correction terms. We further show that, in the case of a harmonic external driving field, these series can be summed and we provide closed analytical expressions for the Gilbert and inertial parameters that are functions of the frequency of the driving field.

An Extended Guinier Analysis for Intrinsically Disordered Proteins

Guinier analysis allows model-free determination of the radius of gyration (Rg) of a biomolecule from X-ray or neutron scattering data, in the limit of very small scattering angles. Its range of validity is well understood for globular proteins, but is known to be more restricted for unfolded or intrinsically disordered proteins (IDPs). We have used ensembles of disordered structures from molecular dynamics simulations to investigate which structural properties cause deviations from the Guinier approximation at small scattering angles. We find that the deviation from the Guinier approximation is correlated with the polymer scaling exponent ν describing the unfolded ensemble. We therefore introduce an empirical, ν-dependent, higher-order correction term, to augment the standard Guinier analysis. We test the new fitting scheme using all-atom simulation data for several IDPs and experimental data for both an IDP and a destabilized mutant of a folded protein. In all cases tested, we achieve an accuracy of the inferred Rg within ∼3% of the true Rg. The method is straightforward to implement and extends the range of validity to a maximum qRg of ∼2 versus ∼1.1 for Guinier analysis. Compared with the Guinier or Debye approaches, our method allows data from wider angles with lower noise to be used to analyze scattering data accurately. In addition to Rg, our fitting scheme also yields estimates of the scaling exponent ν in excellent agreement with the reference ν determined from the underlying molecular ensemble.

Holographic conductivity of holographic superconductors with higher-order corrections

We analytically and numerically disclose the effects of the higher-order correction terms in the gravity and in the gauge field on the properties of s-wave holographic superconductors. On the gravity side, we consider the higher curvature Gauss–Bonnet corrections and on the gauge field side, we add a quadratic correction term to the Maxwell Lagrangian. We show that, for this system, one can still obtain an analytical relation between the critical temperature and the charge density. We also calculate the critical exponent and the condensation value both analytically and numerically. We use a variational method, based on the Sturm–Liouville eigenvalue problem for our analytical study, as well as a numerical shooting method in order to compare with our analytical results. For a fixed value of the Gauss–Bonnet parameter, we observe that the critical temperature decreases with increasing the nonlinearity of the gauge field. This implies that the nonlinear correction term to the Maxwell electrodynamics makes the condensation harder. We also study the holographic conductivity of the system and disclose the effects of the Gauss–Bonnet and nonlinear parameters alpha and b on the superconducting gap. We observe that, for various values of alpha and b, the real part of the conductivity is proportional to the frequency per temperature, omega /T, as the frequency is large enough. Besides, the conductivity has a minimum in the imaginary part which is shifted toward greater frequency with decreasing temperature.

Quantum interaction between two gravitationally polarizable objects in the presence of boundaries

We investigate, in the framework of the linearized quantum gravity and the leading-order perturbation theory, the quantum correction to the classical Newtonian interaction between a pair of gravitationally polarizable objects in the presence of both Neumann and Dirichlet boundaries. We obtain general results for the interaction potential and find that the presence of a boundary always strengthens in the leading-order the interaction as compared with the case in absence of boundaries. But different boundaries yield a different degree of strengthening. In the limit when one partner of the pair is placed very close to the Neumann boundary, the interaction potential is larger when the pair is parallel with the boundary than when it is perpendicular to, which is just opposite to the case when the boundary is Dirichlet where the latter is larger than the former. In addition, we find that the pair-boundary separation dependence of the higher-order correction term is determined by the orientation of the pair with respect to boundary, with the parallel case giving a quadratic behavior and the perpendicular case a linear one.

Near integrability of kink lattice with higher order interactions**Supported by Shandong Provincial Natural Science Foundation (ZR2014AQ007), National Natural Science Foundation of China (11403015, U1531105), S. He is supported by Max-Planck fellowship in Germany and National Natural Science Foundation of China (11305235)

We make use of Manton’s analytical method to investigate the force between kinks and anti-kinks at large distances in 1+1 dimensional field theory. The related potential has infinite order corrections of exponential pattern, and the coefficients for each order are determined. These coefficients can also be obtained by solving the equation of the fluctuations around the vacuum. At the lowest order, the kink lattice represents the Toda lattice. With higher order correction terms, the kink lattice can represent one kind of generic Toda lattice. With only two sites, the kink lattice is classically integrable. If the number of sites of the lattice is larger than two, the kink lattice is not integrable but is a near integrable system. We make use of Flaschka’s variables to study the Lax pair of the kink lattice. These Flaschka’s variables have interesting algebraic relations and non-integrability can be manifested. We also discuss the higher Hamiltonians for the deformed open Toda lattice, which has a similar result to the ordinary deformed Toda.

Towards a quantitative comparison between global and local stability analysis

A methodology is proposed here to estimate the stability characteristics of bluff-body wakes using local analysis under the assumption of weakly non-parallel flows. In this connection, a generalisation of the classic spatio-temporal stability analysis for fully three-dimensional flows is first described. Secondly, an additional higher-order correction term with respect to the common saddle-point global frequency estimation is included in the analysis. The proposed method is first validated for the case of the flow past a circular cylinder and then applied to two fully three-dimensional flows: the boundary layer flow over a wall-mounted hemispherical body and the wake flow past a fixed sphere. In all the cases considered, both the estimated unstable eigenvalue and the spatial shape of the associated eigenmode are determined by local stability analysis, and they are compared with the reference counterparts obtained at a definitely higher computational cost by a fully three-dimensional global stability analysis. It is shown that the results of local stability analysis, when the higher-order correction term is included, are in excellent agreement with those obtained by global stability analysis. It is also shown that the retained correction term is of crucial importance in this perspective, leading to a remarkable improvement in accuracy with respect to the classical saddle-point estimation.