Variational supersymmetric approach to evaluate Fokker-Planck probability

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

Stochastic variational approach applied to 2, 3 and 4 alpha particles

A stochastic variational method is applied to calculate the binding energies and root-mean-square radii of 2, 3 and 4 alpha particles using an S-wave Ali-Bodmer potential. The results agree with other calculations. We discuss the application of the present method to study the universality in weakly-bound three and four-body systems in the context of ultracold atomic traps.

COMPLEX KOHN VARIATIONAL PRINCIPLE FOR 2-NUCLEON BOUND-STATE AND SCATTERING WITH THE TENSOR POTENTIAL

Complex Kohn variational principle is applied to the numerical solution of the fully off-shell Lippmann-Schwinger equation for nucleon-nucleon scattering for various partial waves including the coupled S-3(1), D-3(1), channel. Analytic expressions are obtained for all the integrals in the method for a suitable choice of expansion functions. Calculations with the partial waves S-1(0), P-1(1), D-1(2), and S-3(1)-D-3(1) of the Reid soft core potential show that the method converges faster than other solution schemes not only for the phase shift but also for the off-shell t matrix elements. We also show that it is trivial to modify this variational principle in order to make it suitable for bound-state calculation. The bound-state approach is illustrated for the S-3(1)-D-3(1) channel of the Reid soft-core potential for calculating the deuteron binding, wave function, and the D state asymptotic parameters. (c) 1995 Academic Press, Inc.

DIRAC CONFINEMENT BY A VARIATIONAL APPROACH

An analytical approximate method for the Dirac equation with confining power law scalar plus vector potentials, applicable to the problem of the relativistic quark confinement, is presented. The method consists in an improved version of a saddle-point variational approach and it is applied to the fundamental state of massless single quarks for some especial cases of physical interest. Our treatment emphasizes aspects such as the quantum-mechanical relativistic Virial theorem, the saddle-point character of the critical point of the expectation value of the total energy, as well as the Klein paradox and the behaviour of the saddle-point variational energies and wave functions.

Supersymmetric variational energies of 3d confined potentials

Within the approach of supersymmetric quantum mechanics associated with the variational method a recipe to construct the superpotential of three-dimensional confined potentials in general is proposed. To illustrate the construction, the energies of the harmonic oscillator and the Hulthen potential, both confined in three dimensions are evaluated. Comparison with the corresponding results of other approximative and exact numerical results is presented. (C) 2003 Elsevier B.V. All rights reserved.

Supersymmetric variational energies for the confined Coulomb system

The methodology based on the association of the variational method with supersymmetric quantum mechanics is used to evaluate the energy states of the confined hydrogen atom. (C) 2002 Elsevier B.V. B.V. All rights reserved.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Variational studies and replica symmetry breaking in the generalization problem of the binary perceptron

We analyze the average performance of a general class of learning algorithms for the nondeterministic polynomial time complete problem of rule extraction by a binary perceptron. The examples are generated by a rule implemented by a teacher network of similar architecture. A variational approach is used in trying to identify the potential energy that leads to the largest generalization in the thermodynamic limit. We restrict our search to algorithms that always satisfy the binary constraints. A replica symmetric ansatz leads to a learning algorithm which presents a phase transition in violation of an information theoretical bound. Stability analysis shows that this is due to a failure of the replica symmetric ansatz and the first step of replica symmetry breaking (RSB) is studied. The variational method does not determine a unique potential but it allows construction of a class with a unique minimum within each first order valley. Members of this class improve on the performance of Gibbs algorithm but fail to reach the Bayesian limit in the low generalization phase. They even fail to reach the performance of the best binary, an optimal clipping of the barycenter of version space. We find a trade-off between a good low performance and early onset of perfect generalization. Although the RSB may be locally stable we discuss the possibility that it fails to be the correct saddle point globally. ©2000 The American Physical Society.

Variational calculation of positronium-helium-atom scattering length

A basis-set calculation scheme for S-waves Ps-He elastic scattering below the lowest inelastic threshold was formulated using a variational expression for the transition matrix. The scheme was illustrated numerically by calculating the scattering length in the electronic doublet state: a=1.0±0.1 a.u.

Confined systems through variational supersymmetric quantum mechanics

The energy states of the confined harmonic oscillator and the Hulthén potentials are evaluated using the Variational Method associated to Supersymmetric Quantum Mechanics.

An unproblematic method to solve economic and emission dispatch

This paper proposes a method to determine the output of all online units with minimum total cost when the amount of emission is reasonable. A joint economic and emission dispatch is proposed in order to get a significant compromise between costs and emission such that real power supply-demand equilibrium is satisfied. In order to have a meaningful compromise between costs and emission in the problem formulation, two variables are used, weighting factor and price penalty factor. A case study comprising of a 3-unit power system is employed, where various demand is used. Results for the test system indicate the fastness and effectiveness of proposed method. © 2011 IEEE.

Supervised variational relevance learning, an analytic geometric feature selection with applications to omic datasets

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Augmented Lagrangian method with nonmonotone penalty parameters for constrained optimization

At each outer iteration of standard Augmented Lagrangian methods one tries to solve a box-constrained optimization problem with some prescribed tolerance. In the continuous world, using exact arithmetic, this subproblem is always solvable. Therefore, the possibility of finishing the subproblem resolution without satisfying the theoretical stopping conditions is not contemplated in usual convergence theories. However, in practice, one might not be able to solve the subproblem up to the required precision. This may be due to different reasons. One of them is that the presence of an excessively large penalty parameter could impair the performance of the box-constraint optimization solver. In this paper a practical strategy for decreasing the penalty parameter in situations like the one mentioned above is proposed. More generally, the different decisions that may be taken when, in practice, one is not able to solve the Augmented Lagrangian subproblem will be discussed. As a result, an improved Augmented Lagrangian method is presented, which takes into account numerical difficulties in a satisfactory way, preserving suitable convergence theory. Numerical experiments are presented involving all the CUTEr collection test problems.

The boundedness of penalty parameters in an augmented Lagrangian method with constrained subproblems

Augmented Lagrangian methods are effective tools for solving large-scale nonlinear programming problems. At each outer iteration, a minimization subproblem with simple constraints, whose objective function depends on updated Lagrange multipliers and penalty parameters, is approximately solved. When the penalty parameter becomes very large, solving the subproblem becomes difficult; therefore, the effectiveness of this approach is associated with the boundedness of the penalty parameters. In this paper, it is proved that under more natural assumptions than the ones employed until now, penalty parameters are bounded. For proving the new boundedness result, the original algorithm has been slightly modified. Numerical consequences of the modifications are discussed and computational experiments are presented.

Variational second order flow estimation for PIV sequences

[EN] We present in this paper a variational approach to accurately estimate simultaneously the velocity field and its derivatives directly from PIV image sequences. Our method differs from other techniques that have been presented in the literature in the fact that the energy minimization used to estimate the particles motion depends on a second order Taylor development of the flow. In this way, we are not only able to compute the motion vector field, but we also obtain an accurate estimation of their derivatives. Hence, we avoid the use of numerical schemes to compute the derivatives from the estimated flow that usually yield to numerical amplification of the inherent uncertainty on the estimated flow. The performance of our approach is illustrated with the estimation of the motion vector field and the vorticity on both synthetic and real PIV datasets.