High-order harmonic generation of ZnO crystals in chirped and static electric fields

High harmonic generation in ZnO crystals under chirped single-color field and static electric field are investigated by solving the semiconductor Bloch equation(SBE). It is found that when the chirp pulse is introduced, the interference structure becomes obvious while the harmonic cutoff is not extended. Furthermore, the harmonic efficiency is improved when the static electric field is included. These phenomena are demonstrated by the classical recollision model in real space affected by the waveform of laser field and inversion symmetry. Specifically, the electron motion in k-space shows that the change of waveform and the destruction of the symmetry of the laser field lead to the incomplete X-structure of the crystal-momentum-resolved(k-resolved) inter-band harmonic spectrum. Furthermore, a pre-acceleration process in the solid four-step model is confirmed.

UNIQUENESS FOR SOLUTIONS OF NONHOMOGENEOUS A -HARMONIC EQUATIONS WITH VERY WEAK BOUNDARY VALUES

A uniqueness result for nonhomogeneous quasilinear elliptic partial differential equations with very weak boundary values was proved.

Fully discrete Jacobi-spherical harmonic spectral method for Navier-Stokes equations

A fully discrete Jacobi-spherical harmonic spectral method is provided for the Navier-Stokes equations in a ball. Its stability and convergence are proved. Numerical results show efficiency of this approach. The proposed method is also applicable to other problems in spherical geometry.

Interference of harmonics emitted by different tunneling momentum channels in laser fields

By numerically solving the semiconductor Bloch equation(SBEs),we theoretically study the high-harmonic generation of ZnO crystals driven by one-color and two-color intense laser pulses.The results show the enhancement of harmonics and the cut-off remains the same in the two-color field,which can be explained by the recollision trajectories and electron excitation from multi-channels.Based on the quantum path analysis,we investigate contribution of different ranges of the crystal momentum k of ZnO to the harmonic yield,and find that in two-color laser fields,the intensity of the harmonic yield of different ranges from the crystal momentum makes a big difference and the harmonic intensity is depressed from all k channels,which is related to the interferences between harmonics from symmetric k channels.

Numerical Study of the Vibrations of Beams with Variable Stiffness under Impulsive or Harmonic Loading

The behavior of beams with variable stiffness subjected to the action of variable loadings (impulse or harmonic) is analyzed in this paper using the successive approximation method. This successive approximation method is a technique for numerical integration of partial differential equations involving both the space and time, with well-known initial conditions on time and boundary conditions on the space. This technique, although having been applied to beams with constant stiffness, is new for the case of beams with variable stiffness, and it aims to use a quadratic parabola (in time) to approximate the solutions of the differential equations of dynamics. The spatial part is studied using the successive approximation method of the partial differential equations obtained, in order to transform them into a system of time-dependent ordinary differential equations. Thus, the integration algorithm using this technique is established and applied to examples of beams with variable stiffness, under variable loading, and with the different cases of supports chosen in the literature. We have thus calculated the cases of beams with constant or variable rigidity with articulated or embedded supports, subjected to the action of an instantaneous impulse and harmonic loads distributed over its entire length. In order to justify the robustness of the successive approximation method considered in this work, an example of an articulated beam with constant stiffness subjected to a distributed harmonic load was calculated analytically, and the results obtained compared to those found numerically for various steps (spatial h and temporal τ ¯ ) of calculus, and the difference between the values obtained by the two methods was small. For example for ( h=1/8 , τ ¯ =1/ 64 ), the difference between these values is 17%.

On numerical stationary distribution of overdamped Langevin equation in harmonic system

Efficient numerical algorithm for stochastic differential equation has been an important object in the research of statistical physics and mathematics for a long time.In this work we study the highly accurate numerical algorithm for the overdamped Langevin equation.In particular,our interest is in the behaviour of the numerical schemes for solving the overdamped Langevin equation in the harmonic system.Based on the large friction limit of the underdamped Langevin dynamic scheme,three algorithms for overdamped Langevin equation are obtained.We derive the explicit expression of the stationary distribution of each algorithm by analysing the discrete time trajectory for both one-dimensional case and multi-dimensional case.The accuracy of the stationary distribution of each algorithm is illustrated by comparing with the exact Boltzmann distribution.Our results demonstrate that the“BAOA-limit”algorithm generates an accurate distribution of the harmonic system in a canonical ensemble,within a stable range of time interval.The other algorithms do not produce the exact distribution of the harmonic system.

Phase-coherence dynamics of frequency-comb emission via high-order harmonic generation in few-cycle pulse trains

Frequency-comb emission via high-order harmonic generation(HHG)provides an alternative method for the coherent vacuum ultraviolet(VUV)and extreme ultraviolet(XUV)radiation at ultrahigh repetition rates.In particular,the temporal and spectral features of the HHG were shown to carry profound insight into frequency-comb emission dynamics.Here we present an ab initio investigation of the temporal and spectral coherence of the frequency comb emitted in HHG of He atom driven by few-cycle pulse trains.We find that the emission of frequency combs features a destructive and constructive coherences caused by the phase interference of HHG,leading to suppression and enhancement of frequency-comb emission.The results reveal intriguing and substantially different nonlinear optical response behaviors for frequency-comb emission via HHG.The dynamical origin of frequency-comb emission is clarified by analyzing the phase coherence in HHG processes in detail.Our results provide fresh insight into the experimental realization of selective enhancement of frequency comb in the VUV–XUV regimes.

Calibration of quantitative rescattering model for simulating vortex high-order harmonic generation driven by Laguerre–Gaussian beam with nonzero orbital angular momentum

We calibrate the macroscopic vortex high-order harmonic generation(HHG)obtained by the quantitative rescattering(QRS)model to compute single-atom induced dipoles against that by solving the time-dependent Schr?dinger equation(TDSE).We show that the QRS perfectly agrees with the TDSE under the favorable phase-matching condition,and the QRS can accurately predict the main features in the spatial profiles of vortex HHG if the phase-matching condition is not good.We uncover that harmonic emissions from short and long trajectories are adjusted by the phase-matching condition through the time-frequency analysis and the QRS can simulate the vortex HHG accurately only when the interference between two trajectories is absent.This work confirms that it is an efficient way to employ the QRS model in the single-atom response for precisely simulating the macroscopic vortex HHG.

Spectral Method for Three-Dimensional Nonlinear Klein-Gordon Equation by Using Generalized Laguerre and Spherical Harmonic Functions

In this paper,a generalized Laguerre-spherical harmonic spectral method is proposed for the Cauchy problem of three-dimensional nonlinear Klein-Gordon equation. The goal is to make the numerical solutions to preserve the same conservation as that for the exact solution.The stability and convergence of the proposed scheme are proved.Numerical results demonstrate the efficiency of this approach.We also establish some basic results on the generalized Laguerre-spherical harmonic orthogonal approximation,which play an important role in spectral methods for various problems defined on the whole space and unbounded domains with spherical geometry.

Spherical harmonics method for neutron transport equation based on unstructured-meshes

Based on a new second-order neutron transport equation, self-adjoint angular flux (SAAF) equation, the spherical harmonics (PN) method for neutron transport equation on unstructured-meshes is derived. The spherical harmonics function is used to expand the angular flux. A set of differential equations about the spatial variable, which are coupled with each other, can be obtained. They are solved iteratively by using the finite element method on un- structured-meshes. A two-dimension transport calculation program is coded according to the model. The numerical results of some benchmark problems demonstrate that this method can give high precision results and avoid the ray effect very well.

ON VERY WEAK SOLUTIONS OF A-HARMONICEQUATION WITH VERY WEAK BOUNDARYVALUES

In this paper,the following result is given by using Hodge decomposition: There exists r(0) = r(0)(n,p,a,b), such that if u is an element of W-loc(1,r)(Omega) is a very weak solution of (1.1),with C max{1,p - 1} < r < p and u is an element of W-0(1,r)(Omega;partial derivativeOmega\E) where E subset of partial derivativeOmega is a closed set and small in an appropriate capacity sense, then u = 0, a.e. in Omega provided that r(0) < r < p.

Exact Radial Solution of the Non-relativistic Schrdinger Equation for the Helium Atom with the Potential Harmonic Method

We proposed a simple potential harmonic(PH) scheme for calculating the non\|relativistic radial correlation energies of atomic systems. The scheme was applied to the low\|lying \%n\%\+1\%S\%(\%n\%=1,2) and \%n\%\+3\%S\%(\%n\%=2,3) states of the helium atom. The results exhibit a very stable convergence characterization in both the angular and radial directions with PH and generalized Laguerre functions(GLF) respectively, even though the method is non\|variational one. The ninth significant figure of the non\|relativistic radial energy(NRE) calculated for the ground state exactly agrees with that of the most accurate literature data from the modified configuration interaction method. The convergent NRE′s for the excited states 2\+1\%S\%, 2\+3\%S\% and 3\+3\%S\% with the similar accuracy were also obtained.

High-gain adaptive regulator for a string equation with uncertain harmonic disturbance under boundary output feedback control

This paper considers the boundary stabilization and parameter estimation of a one-dimensional wave equation in the case when one end is fixed and control and harmonic disturbance with uncertain amplitude are input at another end. A high-gain adaptive regulator is designed in terms of measured collocated end velocity. The existence and uniqueness of the classical solution of the closed-loop system is proven. It is shown that the state of the system approaches the standstill as time goes to infinity and meanwhile , the estimated parameter converges to the unknown parameter.

Generalized Harmonic Oscillator and the Schrdinger Equation with Position-Dependent Mass

We study the generalized harmonic oscillator that has both the position-dependent mass and the potential depending on the form of mass function in a more general framework. The explicit expressions of the eigenvalue and eigenfunction for such a system are given, they have the same forms as those for the usual harmonic oscillator with constant mass. The coherent state and its properties corresponding effective potentials for several mass functions, for the system with PDM are also discussed. We give the the systems with such potentials are isospectral to the usual harmonic oscillator.

1/3 PURE SUB-HARMONIC SOLUTION AND FRACTAL CHARACTERISTIC OF TRANSIENT PROCESS FOR DUFFING'S EQUATION

The 1/3 sub-harmonic solution for the Duffing＇s with damping equation was investigated by using the methods of harmonic balance and numerical integration. The assumed solution is introduced, and the domain of sub-harmonic frequencies was found. The asymptotical stability of the subharmonic resonances and the sensitivity of the amplitude responses to the variation of damping coefficient were examined. Then, the subharmonic resonances were analyzed by using the techniques from the general fractal theory. The analysis indicates that the sensitive dimensions of the system time-field responses show sensitivity to the conditions of changed initial perturbation, changed damping coefficient or the amplitude of excitation, thus the sensitive dimension can clearly describe the characteristic of the transient process of the subharmonic resonances.

Extremum Principle for Very Weak Solutions of A-Harmonic Equation with Weight

Extremum principle for very weak solutions of A-harmonic equation div A(x,▽u)=0 is obtained, where the operator A:Ω × Rn→Rnsatisfies some coercivity and controllable growth conditions with Mucken-houpt weight.

REGULARITY FOR VERY WEAK SOLUTIONS TO A-HARMONIC EQUATION

In this paper, the following result is given by using Hodge decomposition and weak reverse Holder inequality： For every r1 with P-（2^n＋1 100^n^2 p（2^3＋n/（P-1）＋1）b/a）^-1〈r1〈p,there exists the exponent r2 = r2（n, r1,p） 〉 p, such that for every very weak solution u∈W^1r1,loc（Ω） to A-harmonic equation, u also belongs to W^1r2,loc（Ω） . In particular, u is the weak solution to A-harmonic equation in the usual sense.

The spin-one Duffin Kemmer Petiau equation in the presence of pseudo-harmonic oscillatory ring-shaped potential

The Duffin-Kemmer-Petiau equation （DKP） is studied in the presence of a pseudo-harmonic oscillatory ring-shaped potential in （1 ＋ 3）-dimensional space-time for spin-one particles. The exact energy eigenvalues and the eigenfunctions are obtained using the Nikiforov-Uvarov method.

EVOLUTION ANALYSIS OF TS WAVE AND HIGH-ORDER HARMONIC WAVES IN BOUNDARY LAYERS

Linear and nonlinear evolutions of TS wave and high-order harmonic waves in boundary layers are studied based on the parabolic stability equation （PSE）. Initial conditions are derived by the local method with the Landau expansion. The evolution process and characteristics of the disturbance amplitude and the velocity profile, etc. , especially stronger nonlinear effects, are computed by an efficient numerical method. Effects and regulations of different initial amplitudes, frequencies and pressure gradients on the evolution of disturbances are explored, which are directly relative to the stability and the transition in boundary layers. Simulation results are in good agreement with the data of the accuracy direct numerical simulation （DNS） using full Navier-Stokes equations.

DYNAMIC BEHAVIOR OF TWO PARALLEL SYMMETRY CRACKS IN MAGNETO-ELECTRO-ELASTIC COMPOSITES UNDER HARMONIC ANTI-PLANE WAVES

The dynamic behavior of two parallel symmetry cracks in magneto-electro-elastic composites under harmonic anti-plane shear waves is studied by Schmidt method. By using the Fourier transform, the problem can be solved with a pair of dual integral equations in which the unknown variable is the jumps of the displacements across the crack surfaces. To solve the dual integral equations, the jumps of the displacements across the crack surface were expanded in a series of Jacobi polynomials. The relations among the electric filed, the magnetic flux and the stress field were obtained. From the results, it can be obtained that the singular stresses in piezoelectric/piezomagnetic materials carry the same forms as those in a general elastic material for the dynamic anti-plane shear fracture problem. The shielding effect of two parallel cracks was also discussed.