Variational description of continuum states in terms of integral relations

Two integral relations derived from the Kohn Variational Principle (KVP) are used for describing scattering states. In usual applications the KVP requires the explicit form of the asymptotic behavior of the scattering wave function. This is not the case when the integral relations are applied since, due to their short range nature, the only condition for the scattering wave function $\Psi$ is that it be the solution of $(H-E)\Psi=0$ in the internal region. Several examples are analyzed for the computation of phase-shifts from bound state type wave functions or, in the case of the scattering of charged particles, it is possible to obtain phase-shifts using free asymptotic conditions. As a final example we discuss the use of the integral relations in the case of the Hyperspherical Adiabatic method.Comment: 34 pages, 7 figures, accepted in Phys. Rev.

General integral relations for the description of scattering states using the hyperspherical adiabatic basis

In this work we investigate 1+2 reactions within the framework of the hyperspherical adiabatic expansion method. To this aim two integral relations, derived from the Kohn variational principle, are used. A detailed derivation of these relations is shown. The expressions derived are general, not restricted to relative $s$ partial waves, and with applicability in multichannel reactions. The convergence of the ${\cal K}$-matrix in terms of the adiabatic potentials is investigated. Together with a simple model case used as a test for the method, we show results for the collision of a $^4$He atom on a \dimer dimer (only the elastic channel open), and for collisions involving a $^6$Li and two $^4$He atoms (two channels open).Comment: Accepted for publication in Physical Review

Integral relations and the adiabatic expansion method for 1+2 reactions above the breakup threshold: Helium trimers with soft-core potentials

The integral relations formalism introduced in \cite{bar09,rom11}, and designed to describe 1+$N$ reactions, is extended here to collision energies above the threshold for the target breakup. These two relations are completely general, and in this work they are used together with the adiabatic expansion method for the description of 1+2 reactions. The neutron-deuteron breakup, for which benchmark calculations are available, is taken as a test of the method. The s-wave collision between the $^4$He atom and $^4$He$_2$ dimer above the breakup threshold and the possibility of using soft-core two-body potentials plus a short-range three-body force will be investigated. Comparisons to previous calculations for the three-body recombination and collision dissociation rates will be shown.Comment: To be published in Physical Review

The Three-Nucleon System Near the N-d Threshold

The three-nucleon system is studied at energies a few hundred keV above the N-d threshold. Measurements of the tensor analyzing powers $T_{20}$ and $T_{21}$ for p-d elastic scattering at $E_{c.m.}=432$ keV are presented together with the corresponding theoretical predictions. The calculations are extended to very low energies since they are useful for extracting the p-d scattering lengths from the experimental data. The interaction considered here is the Argonne V18 potential plus the Urbana three-nucleon potential. The calculation of the asymptotic D- to S-state ratio for $^3$H and $^3$He, for which recent experimental results are available, is also presented.Comment: Latex, 11 pages, 2 figures, to be published in Phy.Lett.

Realistic Calculation of the hep Astrophysical Factor

The astrophysical factor for the proton weak capture on 3He is calculated with correlated-hyperspherical-harmonics bound and continuum wave functions corresponding to a realistic Hamiltonian consisting of the Argonne v18 two-nucleon and Urbana-IX three-nucleon interactions. The nuclear weak charge and current operators have vector and axial-vector components, that include one- and many-body terms. All possible multipole transitions connecting any of the p-3He S- and P-wave channels to the 4He bound state are considered. The S-factor at a p-3He center-of-mass energy of 10 keV, close to the Gamow-peak energy, is predicted to be 10.1 10^{-20} keV b, a factor of five larger than the standard-solar-model value. The P-wave transitions are found to be important, contributing about 40 % of the calculated S-factor.Comment: 8 pages RevTex file, submitted to Phys. Rev. Let

Benchmark calculation of p-3H and n-3He scattering

p-3H and n-3He scattering in the energy range above the n-3He but below the d-d thresholds is studied by solving the 4-nucleon problem with a realistic nucleon-nucleon interaction. Three different methods -- Alt, Grassberger and Sandhas, Hyperspherical Harmonics, and Faddeev-Yakubovsky -- have been employed and their results for both elastic and charge-exchange processes are compared. We observe a good agreement between the three different methods, thus the obtained results may serve as a benchmark. A comparison with the available experimental data is also reported and discussed.Comment: 13 pages, 6 figures. arXiv admin note: text overlap with arXiv:1109.362

Total cross section for p-d breakup below 30 MeV

The total cross section for p-d breakup is studied in terms of the elastic S-matrix through the unitary condition. Calculations using the complex Kohn variational method along with the Pair Correlated Hyperspherical Harmonic basis are presented. The results have been restricted to energies below Ep=30 MeV where Coulomb effects are expected to be sizable and are compared to the existing data. Two different measurements have been found in the literature: 40 years ago, Gibbons and Macklin (1959); and 26 years ago, Carlosn et al. (1973). The calculations are found to be in reasonable agreement with these old data, though a discrepancy is observed near the deuteron breakup threshold. Moreover, a detailed analysis of the contributions to the observable from different partial waves has been presented. Unexpectedly, the main contribution for a wide range of energies has been detected in the J=3/2- state

Benchmark calculation for proton-deuteron elastic scattering observables including Coulomb

Two independent calculations of proton-deuteron elastic scattering observables including Coulomb repulsion between the two protons are compared in the proton lab energy region between 3 MeV and 65 MeV. The hadron dynamics is based on the purely nucleonic charge-dependent AV18 potential. Calculations are done both in coordinate space and momentum space. The coordinate-space calculations are based on a variational solution of the three-body Schr\"odinger equation using a correlated hyperspherical expansion for the wave function. The momentum-space calculations proceed via the solution of the Alt-Grassberger-Sandhas equation using the screened Coulomb potential and the renormalization approach. Both methods agree within 1% on all observables, showing the reliability of both numerical techniques in that energy domain. At energies below three-body breakup threshold the coordinate-space method remains favored whereas at energies higher than 65 MeV the momentum-space approach seems to be more efficient.Comment: Submitted to Phys. Rev.

Low energy n-\nuc{3}{H} scattering : a novel testground for nuclear interaction

The low energy n-\nuc{3}{H} elastic cross sections near the resonance peak are calculated by solving the 4-nucleon problem with realistic NN interactions. Three different methods -- Alt, Grassberger and Shandas (AGS), Hyperspherical Harmonics and Faddeev-Yakubovsky -- have been used and their respective results are compared. We conclude on a failure of the existing NN forces to reproduce the n-\nuc{3}{H} total cross section.Comment: To be published in Phys. Rev.

Robust multi-fidelity design of a micro re-entry unmanned space vehicle

This article addresses the preliminary robust design of a small-scale re-entry unmanned space vehicle by means of a hybrid optimization technique. The approach, developed in this article, closely couples an evolutionary multi-objective algorithm with a direct transcription method for optimal control problems. The evolutionary part handles the shape parameters of the vehicle and the uncertain objective functions, while the direct transcription method generates an optimal control profile for the re-entry trajectory. Uncertainties on the aerodynamic forces and characteristics of the thermal protection material are incorporated into the vehicle model, and a Monte-Carlo sampling procedure is used to compute relevant statistical characteristics of the maximum heat flux and internal temperature. Then, the hybrid algorithm searches for geometries that minimize the mean value of the maximum heat flux, the mean value of the maximum internal temperature, and the weighted sum of their variance: the evolutionary part handles the shape parameters of the vehicle and the uncertain functions, while the direct transcription method generates the optimal control profile for the re-entry trajectory of each individual of the population. During the optimization process, artificial neural networks are utilized to approximate the aerodynamic forces required by the optimal control solver. The artificial neural networks are trained and updated by means of a multi-fidelity approach: initially a low-fidelity analytical model, fitted on a waverider type of vehicle, is used to train the neural networks, and through the evolution a mix of analytical and computational fluid dynamic, high-fidelity computations are used to update it. The data obtained by the high-fidelity model progressively become the main source of updates for the neural networks till, near the end of the optimization process, the influence of the data obtained by the analytical model is practically nullified. On the basis of preliminary results, the adopted technique is able to predict achievable performance of the small spacecraft and the requirements in terms of thermal protection materials