Time-dependent variational approach to solve multi-dimensional time-dependent Schr?dinger equation

We present an efficient approach to solve multi-dimensional time-dependent Schr?dinger equation(TDSE)in an intense laser field.In this approach,each spatial degree of freedom is treated as a distinguishable quasi-particle.The non-separable Coulomb potential is regarded as a two-body operator between different quasi-particles.The time-dependent variational principle is used to derive the equations of motion.Then the high-order multi-dimensional problem is broken down into several lower-order coupled equations,which can be efficiently solved.As a demonstration,we apply this method to solve the two-dimensional TDSE.The accuracy is tested by comparing the direct solutions of TDSE using several examples such as the strong-field ionization and the high harmonic generation.The results show that the present method is much more computationally efficient than the conventional one without sacrificing accuracy.The present method can be straightforwardly extended to three-dimensional problems.Our study provides a flexible method to investigate the laser-atom interaction in the nonperturbative regime.

Variational Approach for the Adapted Solution of Backw ard Stochastic Differential Equations with Locally Lipschitz Diffusion Coefficients

One existence integral condition was obtained for the adapted solution of the general backward stochastic differential equations(BSDEs). Then by solving the integral constraint condition, and using a limit procedure, a new approach method is proposed and the existence of the solution was proved for the BSDEs if the diffusion coefficients satisfy the locally Lipschitz condition. In the special case the solution was a Brownian bridge. The uniqueness is also considered in the meaning of "F0-integrable equivalent class" . The new approach method would give us an efficient way to control the main object instead of the "noise".

A NEW PROOF OF GAFFNEY’S INEQUALITY FOR DIFFERENTIAL FORMS ON MANIFOLDS-WITH-BOUNDARY:THE VARIATIONAL APPROACH à LA KOZONO-YANAGISAWA

Let(M,g_(0))be a compact Riemannian manifold-with-boundary.We present a new proof of the classical Gaffney inequality for differential forms in boundary value spaces over M,via a variational approach a la Kozono-Yanagisawa[Lr-variational inequality for vector fields and the Helmholtz-Weyl decomposition in bounded domains,Indiana Univ.Math.J.58(2009),1853-1920],combined with global computations based on the Bochner technique.

Fractional Variational Approach for Dissipative Mechanical Systems

More recently, a variational approach has been proposed by Lin and Wang for damping motion with a Lagrangian holding the energy term dissipated by a friction force. However, the modified Euler-Lagrange equation obtained within their for- malism leads to an incorrect Newtonian equation of motion due to the nonlocality of the Lagrangian. In this communication, we generalize this approach based on the fractional actionlike variational approach and we show that under some simple restric- tions connected to the fractional parameters introduced in the fractional formalism, this problem may be solved.

Improvement of variational approach in an interacting two-fermion system

A more reasonable trial ground state wave function is constructed for the relative motion of an interacting two-fermion system in a one-dimensional(1D) harmonic potential. At the boundaries both the wave function and its first derivative are continuous and the quasi-momentum is determined by a more practical constraint condition which associates two variational parameters. The upper bound of the ground state energy is obtained by applying the variational principle to the expectation value of the Hamiltonian of relative motion on the trial wave function. The resulted energy and wave function show better agreement with the analytical solution than the original proposal.

A VARIATIONAL APPROACH FOR DETECTING FEATURE LINES ON MESHES

Feature lines are fundamental shape descriptors and have been extensively applied to computer graphics, computer-aided design, image processing, and non-photorealistic renderingi This paper introduces a unified variational framework for detecting generic feature lines on polygonal meshes. The classic Mumford-Shah model is extended to surfaces. Using F-convergence method and discrete differential geometry, we discretize the proposed variational model to sequential coupled sparse linear systems. Through quadratic polyno- mials fitting, we develop a method for extracting valleys of functions defined on surfaces. Our approach provides flexible and intuitive control over the detecting procedure, and is easy to implement. Several measure functions are devised for different types of feature lines, and we apply our approach to various polygonal meshes ranging from synthetic to measured models. The experiments demonstrate both the effectiveness of our algorithms and the visual quality of results.

Convergence Analysis of a Numerical Scheme for the Porous Medium Equation by an Energetic Variational Approach

The porous medium equation(PME)is a typical nonlinear degenerate parabolic equation.We have studied numerical methods for PME by an energetic vari-ational approach in[C.Duan et al.,J.Comput.Phys.,385(2019),pp.13–32],where the trajectory equation can be obtained and two numerical schemes have been devel-oped based on different dissipative energy laws.It is also proved that the nonlinear scheme,based on f logf as the total energy form of the dissipative law,is uniquely solv-able on an admissible convex set and preserves the corresponding discrete dissipation law.Moreover,under certain smoothness assumption,we have also obtained the sec-ond order convergence in space and the first order convergence in time for the scheme.In this paper,we provide a rigorous proof of the error estimate by a careful higher or-der asymptotic expansion and two step error estimates.The latter technique contains a rough estimate to control the highly nonlinear term in a discrete W 1,∞norm and a refined estimate is applied to derive the optimal error order.

Assimilated Tidal Results of Tide Gauge and TOPEX/POSEIDON Data over the China Seas Using a Variational Adjoint Approach with a Nonlinear Numerical Model

In order to obtain an accurate tide description in the China Seas, the 2-dimensional nonlinear numerical Princeton Ocean Model （POM） is employed to incorporate in situ tidal measurements both from tide gauges and TOPEX/POSEIDON （T/P） derived datasets by means of the variational adjoint approach in such a way that unknown internal model parameters, bottom topography, friction coefficients and open boundary conditions, for example, are adjusted during the process. The numerical model is used as a forward model. After the along-track T/P data are processed, two classical methods, i.e. harmonic and response analysis, are implemented to estimate the tide from such datasets with a domain covering the model area extending from 0° to 41°N in latitude and from 99°E to 142°E in longitude. And the results of these two methods are compared and interpreted. The numerical simulation is performed for 16 major constituents. In the data assimilation experiments, three types of unknown parameters （water depth, bottom friction and tidal open boundary conditions in the model equations） are chosen as control variables. Among the various types of data assimilation experiments, the calibration of water depth brings the most promising results. By comparing the results with selected tide gauge data, the average absolute errors are decreased from 7.9 cm to 6.8 cm for amplitude and from 13.0° to 9.0° for phase with respect to the semidiurnal tide M2 constituent, which is the largest tidal constituent in the model area. After the data assimilation experiment is performed, the comparison between model results and tide gauge observation for water levels shows that the RMS errors decrease by 9 cm for a total of 14 stations, mostly selected along the coast of China's Mainland, when a one-month period is considered, and the correlation coefficients improve for most tidal stations among these stations.

Monthly gravity field recovery from GRACE orbits and K-band measurements using variational equations approach

The Gravity Recovery and Climate Experiment（GRACE） mission can significantly improve our knowledge of the temporal variability of the Earth＇s gravity field.We obtained monthly gravity field solutions based on variational equations approach from GPS-derived positions of GRACE satellites and K-band range-rate measurements.The impact of different fixed data weighting ratios in temporal gravity field recovery while combining the two types of data was investigated for the purpose of deriving the best combined solution.The monthly gravity field solution obtained through above procedures was named as the Institute of Geodesy and Geophysics（IGG） temporal gravity field models.IGG temporal gravity field models were compared with GRACE Release05（RL05） products in following aspects：（i） the trend of the mass anomaly in China and its nearby regions within 2005-2010; （ii） the root mean squares of the global mass anomaly during 2005-2010; （iii） time-series changes in the mean water storage in the region of the Amazon Basin and the Sahara Desert between 2005 and 2010.The results showed that IGG solutions were almost consistent with GRACE RL05 products in above aspects（i）-（iii）.Changes in the annual amplitude of mean water storage in the Amazon Basin were 14.7 ± 1.2 cm for IGG,17.1 ± 1.3 cm for the Centre for Space Research（CSR）,16.4 ± 0.9 cm for the GeoForschungsZentrum（GFZ） and 16.9 ± 1.2 cm for the Jet Propulsion Laboratory（JPL） in terms of equivalent water height（EWH）,respectively.The root mean squares of the mean mass anomaly in Sahara were 1.2 cm,0.9 cm,0.9 cm and 1.2 cm for temporal gravity field models of IGG,CSR,GFZ and JPL,respectively.Comparison suggested that IGG temporal gravity field solutions were at the same accuracy level with the latest temporal gravity field solutions published by CSR,GFZ and JPL.

Variational solutions for Hermite-Gaussian solitons in nonlocal nonlinear media

The physical features exhibited by Hermite-Gaussian （HC） beams propagating in nonlocal nonlinear media with Gaussian-shaped response are discussed with an approximate variational method.Using direct numerical simulations,we find that the beam properties in the normalized system are different with the change of the degree of nonlocality.It is shown that initial HG profiles break up into several individual beams with propagation when the degree of nonlocality α is small.HG beams can propagate stably when a is large enough.

Rapid determination of Cl^- by light-absorption ratio variation approach using Cl^-eosin Y-Ag^＋ adsorptive precipitation

The adsorptive precipitation between eosin Y （EY） and AgCl colloids at pH 3.73 caused the sensitive color change of the solution. The reaction mechanism between EY and AgCl was analyzed and this reaction was used for determination of Cl^- in trace level by the light-absorption radio variation approach.

Adiabatic evolution of optical beams of arbitrary shapes in nonlocal nonlinear media

We discuss evolution of Hermite–Gaussian beams of different orders in nonlocal nonlinear media whose characteristic length is set as different functions of propagation distance,using the variational approach.It is proved that as long as the characteristic length varies slowly enough,all the Hermite–Gaussian beams can propagate adiabatically.When the characteristic length gradually comes back to its initial value after changes,all the Hermite–Gaussian beams can adiabatically restore to their own original states.The variational results agree well with the numerical simulations.Arbitrary shaped beams synthesized by Hermite–Gaussian modes can realize adiabatic evolution in nonlocal nonlinear media with gradual characteristic length.

Sufficient conditions of Rayleigh-Taylor stability and instability in equatorial ionosphere

Rayleigh-Taylor （R-T） instability is known as the fundamental mechanism of equatorial plasma bubbles （EPBs）. However, the sufficient conditions of R-T instability and stability have not yet been derived. In the present paper, the sufficient conditions of R-T stability and instability are preliminarily^derived. Linear equations for small perturbation are first obtained from the electron/ion continuity equations, momentum equations, and the current continuity equation in the equatorial ionosphere. The linear equations can be casted as an eigenvalue equation using a normal mode method. The eigenvalue equation is a variable coefficient linear equation that can be solved using a variational approach. With this approach, the sufficient conditions can be obtained as follows： if the minimum systematic eigenvalue is greater than one, the ionosphere is R-T unstable; while if the maximum systematic eigenvalue is less than one, the ionosphere is R-T stable. An approximate numerical method for obtaining the systematic eigenvalues is introduced, and the R-T stable/unstable areas are calculated. Numerical experiments axe designed to validate the sufficient conditions. The results agree with the derived suf- ficient conditions.

Effects of localized impurity on a dark soliton in a Bose-Einstein condensate with an external magnetic trap

The dynamics of a dark soliton has been investigated in a Bose Einstein condensate with an external magnetic trap, and the effects of localized impurity on the dynamics are discussed by the variational approach based on the renormalized integrals of motion. The reciprocal movement of the dark soliton is discussed by performing a standard linear analysis, and it is found that the effects of the localized impurity depend strictly on the positive or negative value of the impurity strength corresponding to the repulsive or attractive impurity. The numerical results confirm the theoretical analysis, and show that the effects also depend on the effective nonlinear coefficient and the harmonic frequency.

Collapse arrest in the space-fractional Schrodinger equation with an optical lattice

The soliton solution and collapse arrest are investigated in the one-dimensional space-fractional Schrodinger equation with Kerr nonlinearity and optical lattice.The approximate analytical soliton solutions are obtained based on the variational approach,which provides reasonable accuracy.Linear-stability analysis shows that all the solitons are linearly stable.No collapses are found when the Levy index 1<α≤2.Forα=1,the collapse is arrested by the lattice potential when the amplitude of perturbations is small enough.It is numerically proved that the energy criterion of collapse suppression in the two-dimensional traditional Schrodinger equation still holds in the one-dimensional fractional Schr odinger equation.The physical mechanism for collapse prohibition is also given.

Dynamics of vector solitons in two-component Bose-Einstein condensates with time-dependent interactions and harmonic potential

We present two kinds of exact vector-soliton solutions for coupled nonlinear Schrodinger equations with time- varying interactions and time-varying harmonic potential. Using the variational approach, we investigate the dynamics of the vector solitons. It is found that the two bright sol/tons oscillate about slightly and pass through each other around the equilibration state which means that they are stable under our modeh At the same time, we obtain the opposite situation for dark-dark solitons.

Coherent interaction and action–counteraction theory in small polaron systems,and ground state properties

Based on the coherent interaction and action–counteraction principles,we investigate the ground state properties for small polaron systems,the coherent-squeezed fluctuation correction,and the anomalous lattice quantum fluctuation,with the new variational generator containing correlated squeezed-coherent coupling and quantum entanglement.Noting tha t-2t is the T.B.A.energy,for the coherent interaction effect,we find the ground-state energy E_(0)to be-2.428t,in which the coherent squeezed fluctuation correction-A_(0)t is-0.463t(where t is the hopping integral,ωis the phonon frequency),with the electron–one-phonon coupling constant g=1 and the electron–two-phonon coupling constant g_(1)=-0.1.However,as a result of the action–counteraction effect,E_(0)is-2.788t,but-E_(0)t is-0.735t.As to the polaron binding energy(EP),for the coherent interaction effect,E_(P) is-1.38ω,but for the action–counteraction effect,E_(P) is-1.88ω.In particular,the electron–two-phonon interaction noticeably enlarges the coherent interaction and the coherent squeezed quantum fluctuation correction.By intervening with the quantum entanglement,the evolutions of the squeezed coherent state and the lattice quantum fluctuation begin to take control.At that time,we encounter a new quantum phase coherence phenomenon—the collapse and revival of inversion repeatedly for the coherent state in the entangled evolution.

Common Fixed Point Theorems and Q-property for Quasi-contractive Mappings under c-distance on TVS-valued Cone Metric Spaces without the Normality

In this paper, we derive the stochastic maximum principle for optimal control problems of the forward-backward Markovian regime-switching system. The control system is described by an anticipated forward-backward stochastic pantograph equation and modulated by a continuous-time finite-state Markov chain. By virtue of classical variational approach, duality method, and convex analysis, we obtain a stochastic maximum principle for the optimal control.

Seismic ground motion variations resulting from site conditions

Amplification of seismic ground motions in the territory of Almaty city is evaluated by using different methods. The pattern and probable causes of ground motion variations in different engineering-geological conditions are characterized. An expeditious application of these techniques within a complex methodical approach for Almaty city microzonation is considered.

Existence for a Higher Order Coupled System of Korteweg-de Vries Equations

Consider the following system of coupled Korteweg-de Vries equations, where u, v ⊆ W2,2, 2≤N≤7 and λi,β > 0, β denotes a real coupling parameter. Firstly, we prove the existence of the solutions of a coupled system of Korteweg-de Vries equations using variation approach and minimization techniques on Nehari manifold. Then, we show the multiplicity of the equations by a bifurcation theory which is rare for studying higher order equations.