A survey of impulsive trajectories Final report

Literature survey of astrodynamics problems on intercept, transfer, and rendezvous trajectorie

Analytical Hartree-Fock gradients for periodic systems

We present the theory of analytical Hartree-Fock gradients for periodic systems as implemented in the code CRYSTAL. We demonstrate how derivatives of the integrals can be computed with the McMurchie-Davidson algorithm. Highly accurate gradients with respect to nuclear coordinates are obtained for systems periodic in 0,1,2 or 3 dimensions.Comment: accepted by International Journal of Quantum Chemistr

Precision orbit computations for an operational environment

Taking advantage of the improvements to the Earth's gravitation field and tracking station coordinates, an orbital computational consistency of the order of 5 meters was achieved for total position differences between orbital solutions for the Seasat and GEOS-3. The main source of error in these solutions was in the mathematical models that are required to generate these results, i.e., gravitation, atmospheric drag, etc. Different Earth gravitation fields and tracking coordinates were analyzed and evaluated in obtaining these computational results. Comparisons and evaluations of the Seasat results were obtained in terms of different solution types such as the Doppler only, Laser only, Doppler and Laser, etc. Other investigation using the Seasat data were made in order to determine their effect on the computational results at this particular level of consistency

Preliminary Orbit Determination System (PODS) for Tracking and Data Relay Satellite System (TDRSS)-tracked target Spacecraft using the homotopy continuation method

The Preliminary Orbit Determination System (PODS) provides early orbit determination capability in the Trajectory Computation and Orbital Products System (TCOPS) for a Tracking and Data Relay Satellite System (TDRSS)-tracked spacecraft. PODS computes a set of orbit states from an a priori estimate and six tracking measurements, consisting of any combination of TDRSS range and Doppler tracking measurements. PODS uses the homotopy continuation method to solve a set of nonlinear equations, and it is particularly effective for the case when the a priori estimate is not well known. Since range and Doppler measurements produce multiple states in PODS, a screening technique selects the desired state. PODS is executed in the TCOPS environment and can directly access all operational data sets. At the completion of the preliminary orbit determination, the PODS-generated state, along with additional tracking measurements, can be directly input to the differential correction (DC) process to generate an improved state. To validate the computational and operational capabilities of PODS, tests were performed using simulated TDRSS tracking measurements for the Cosmic Background Explorer (COBE) satellite and using real TDRSS measurements for the Earth Radiation Budget Satellite (ERBS) and the Solar Mesosphere Explorer (SME) spacecraft. The effects of various measurement combinations, varying arc lengths, and levels of degradation of the a priori state vector on the PODS solutions were considered

Limitation of entanglement due to spatial qubit separation

We consider spatially separated qubits coupled to a thermal bosonic field that causes pure dephasing. Our focus is on the entanglement of two Bell states which for vanishing separation are known as robust and fragile entangled states. The reduced two-qubit dynamics is solved exactly and explicitly. Our results allow us to gain information about the robustness of two-qubit decoherence-free subspaces with respect to physical parameters such as temperature, qubit-bath coupling strength and spatial separation of the qubits. Moreover, we clarify the relation between single-qubit coherence and two-qubit entanglement and identify parameter regimes in which the terms robust and fragile are no longer appropriate.Comment: 7 pages, 3 figures; revised version, accepted for publication in Europhys. Let

Local Variational Principle

A generalization of the Gibbs-Bogoliubov-Feynman inequality for spinless particles is proven and then illustrated for the simple model of a symmetric double-well quartic potential. The method gives a pointwise lower bound for the finite-temperature density matrix and it can be systematically improved by the Trotter composition rule. It is also shown to produce groundstate energies better than the ones given by the Rayleigh-Ritz principle as applied to the groundstate eigenfunctions of the reference potentials. Based on this observation, it is argued that the Local Variational Principle performs better than the equivalent methods based on the centroid path idea and on the Gibbs-Bogoliubov-Feynman variational principle, especially in the range of low temperatures.Comment: 15 pages, 5 figures, one more section adde

Farm plans for a 200-acre central Missouri farm : a comparative analysis of the economic potential for alternative farming systems

Missouri Agricultural Experiment Station in cooperation with Farm Production Economics Division, Economic Research Service, Department of Agriculture.Digitized 2007 AES

Moments of spectral functions: Monte Carlo evaluation and verification

The subject of the present study is the Monte Carlo path-integral evaluation of the moments of spectral functions. Such moments can be computed by formal differentiation of certain estimating functionals that are infinitely-differentiable against time whenever the potential function is arbitrarily smooth. Here, I demonstrate that the numerical differentiation of the estimating functionals can be more successfully implemented by means of pseudospectral methods (e.g., exact differentiation of a Chebyshev polynomial interpolant), which utilize information from the entire interval $(-\beta \hbar / 2, \beta \hbar/2)$. The algorithmic detail that leads to robust numerical approximations is the fact that the path integral action and not the actual estimating functional are interpolated. Although the resulting approximation to the estimating functional is non-linear, the derivatives can be computed from it in a fast and stable way by contour integration in the complex plane, with the help of the Cauchy integral formula (e.g., by Lyness' method). An interesting aspect of the present development is that Hamburger's conditions for a finite sequence of numbers to be a moment sequence provide the necessary and sufficient criteria for the computed data to be compatible with the existence of an inversion algorithm. Finally, the issue of appearance of the sign problem in the computation of moments, albeit in a milder form than for other quantities, is addressed.Comment: 13 pages, 2 figure