A study of the optimization method used in the NAVY/NASA gas turbine engine computer code

Sources of numerical noise affecting the convergence properties of the Powell's Principal Axis Method of Optimization in the NAVY/NASA gas turbine engine computer code were investigated. The principal noise source discovered resulted from loose input tolerances used in terminating iterations performed in subroutine CALCFX to satisfy specified control functions. A minor source of noise was found to be introduced by an insufficient number of digits in stored coefficients used by subroutine THERM in polynomial expressions of thermodynamic properties. Tabular results of several computer runs are presented to show the effects on program performance of selective corrective actions taken to reduce noise

An efficient method for the Quantum Monte Carlo evaluation of the static density-response function of a many-electron system

In a recent Letter we introduced Hellmann-Feynman operator sampling in diffusion Monte Carlo calculations. Here we derive, by evaluating the second derivative of the total energy, an efficient method for the calculation of the static density-response function of a many-electron system. Our analysis of the effect of the nodes suggests that correlation is described correctly and we find that the effect of the nodes can be dealt with

Study to document low thrust trajectory optimization programs HILTOP and ASTOP

Detailed documentation of the HILTOP and ASTOP computer programs is presented along with results of the analyses of the possible extension of the HILTOP program and results of an extra-ecliptic mission study performed with HILTOP

Approximation for discrete Fourier transform and application in study of three-dimensional interacting electron gas

The discrete Fourier transform is approximated by summing over part of the terms with corresponding weights. The approximation reduces significantly the requirement for computer memory storage and enhances the numerical computation efficiency with several orders without loosing accuracy. As an example, we apply the algorithm to study the three-dimensional interacting electron gas under the renormalized-ring-diagram approximation where the Green's function needs to be self-consistently solved. We present the results for the chemical potential, compressibility, free energy, entropy, and specific heat of the system. The ground-state energy obtained by the present calculation is compared with the existing results of Monte Carlo simulation and random-phase approximation.Comment: 11 pages, 13 figure

Isospin and density dependences of nuclear matter symmetry energy coefficients II

Symmetry energy coefficients of explicitly isospin asymmetric nuclear matter at variable densities (from .5$\rho_0$ up to 2 $\rho_0$) are studied as generalized screening functions. An extended stability condition for asymmetric nuclear matter is proposed. We find the possibility of obtaining stable asymmetric nuclear matter even in some cases for which the symmetric nuclear matter limit is unstable. Skyrme-type forces are extensively used in analytical expressions of the symmetry energy coefficients derived as generalized screening functions in the four channels of the particle hole interaction producing alternative behaviors at different $\rho$ and $b$ (respectively the density and the asymmetry coefficient). The spin and spin-isospin coefficients, with corrections to the usual Landau Migdal parameters, indicate the possibility of occurring instabilities with common features depending on the nuclear density and n-p asymmetry. Possible relevance for high energy heavy ions collisions and astrophysical objects is discussed.Comment: 16 pages (latex) plus twelve figures in four eps files, to be published in I.J.M.P.

Phase diagram of a Bose gas near a wide Feshbach resonance

In this paper, we study the phase diagram of a homogeneous Bose gas with a repulsive interaction near a wide Feshbach resonance at zero temperature. The Bose-Einstein-condensation (BEC) state of atoms is a metastable state. When the scattering length $a$ exceeds a critical value depending on the atom density $n$, $na^3>0.035$, the molecular excitation energy is imaginary and the atomic BEC state is dynamically unstable against molecule formation. The BEC state of diatomic molecules has lower energy, where the atomic excitation is gapped and the molecular excitation is gapless. However when the scattering length is above another critical value, $na^3>0.0164$, the molecular BEC state becomes a unstable coherent mixture of atoms and molecules. In both BEC states, the binding energy of diatomic molecules is reduced due to the many-body effect.Comment: 5 pages, 4 figure

Electric-field correlations in quantum charged fluids coupled to the radiation field

In a recent paper [S.El Boustani, P.R.Buenzli, and Ph.A.Martin, Phys.Rev. E 73, 036113 (2006) cond-mat/0511537], about quantum charges in equilibrium with radiation, among other things the asymptotic form of the electric-field correlation has been obtained by a microscopic calculation. It has been found that this correlation has a long-range algebraic decay (except in the classical limit). The macroscopic approach, in the Course of Theoretical Physics of Landau and Lifshitz, gives no such long-range algebraic decay. In this Brief Report, we revisit and complete the macroscopic approach of Landau and Lifshitz, we confirm their result, and suggest that, perhaps, the use of a classical electromagnetic field by El Boustani et al. was not justified.Comment: 10 pages. Title changed. Minor modifications, including a better justification of eq.(8

Zero sound in triplet-correlated superfluid neutron matter

The linear response of a superfluid neutron liquid onto external vector field is studied for the case of ^{3}P_{2}-\,^{3}F_{2} pairing. The consideration is limited to the case when the wave-length of the perturbation is large as compared to the coherence length in the superfluid matter and the transferred energy is small in comparison with the gap amplitude. The obtained results are used to analyse the collisionless phonon-like excitations of the condensate of superfluid neutrons. In particular, we analyze the case of neutron condensation into the state with $m_{j}=0$ which is conventionally considered as the preferable one in the bulk matter of neutron stars. Zero sound (if it exists) is found to be anisotropic and undergoes strong decrement below some temperature threshold depending substantially on the intensity of Fermi-liquid interactions.Comment: 16 pages, 2 figure

Existence of Long-Range Order for Trapped Interacting Bosons

We derive an inequality governing ``long range'' order for a localized Bose-condensed state, relating the condensate fraction at a given temperature with effective curvature radius of the condensate and total particle number. For the specific example of a one-dimensional, harmonically trapped dilute Bose condensate, it is shown that the inequality gives an explicit upper bound for the Thomas-Fermi condensate size which may be tested in current experiments.Comment: 4 pages, 1 figure, RevTex4. Title changed at the request of editors; to appear in Phys. Rev. Letter

Particle linear theory on a self-gravitating perturbed cubic Bravais lattice

Discreteness effects are a source of uncontrolled systematic errors of N-body simulations, which are used to compute the evolution of a self-gravitating fluid. We have already developed the so-called "Particle Linear Theory" (PLT), which describes the evolution of the position of self-gravitating particles located on a perturbed simple cubic lattice. It is the discrete analogue of the well-known (Lagrangian) linear theory of a self-gravitating fluid. Comparing both theories permits to quantify precisely discreteness effects in the linear regime. It is useful to develop the PLT also for other perturbed lattices because they represent different discretizations of the same continuous system. In this paper we detail how to implement the PLT for perturbed cubic Bravais lattices (simple, body and face-centered) in a cubic simulation box. As an application, we will study the discreteness effects -- in the linear regime -- of N-body simulations for which initial conditions have been set-up using these different lattices.Comment: 9 pages, 4 figures and 4 tables. Minor corrections to match published versio