Search NASA⌕ Search

SEARCH · Search NASA

Results for “Perturbation theory”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 253 records · Page 14

Semiclassical theory of unimolecular dissociation induced by a laser field

A semiclassical nonperturbative theory of direct photodissociation in a laser field is developed in which photon absorption and dissociation are treated in a unified fashion. This is achieved by visualizing nuclear dynamics as a representative particle moving on electronic-field surfaces. Methods are described for calculating dissociation rates and probabilities by Monte Carlo selection of initial conditions and integration of classical trajectories on these surfaces. This unified theory reduces to the golden rule expression in the weak-field and short-time limits, and predicts nonlinear behavior, i.e., breakdown of the golden rule expression in intense fields. Field strengths above which lowest-order perturbation theory fails to work have been estimated for some systems. Useful physical insights provided by the electronic-field representation have been illustrated. Intense field effects are discussed which are amenable to experimental observation. The semiclassical methods used here are also applicable to multiple-surface dynamics in fieldfree unimolecular and bimolecular reactions.

Yuan, J.-M.↗

Highly Accurate Quartic Force Field and Rovibrational Spectroscopic Constants for the Azirinyl Cation (c-C2NH2+) and Its Isomers

Azirinyl cation is an aromatic cyclic molecule that is isoelectronic with cyclopropenylidene, c-C3H2, and c-C3H3+. Cyclopropenylidene has been shown to be ubiquitous, existing in many different astrophysical environments. Given the similar chemistry between C and N, and the relative abundances between C and N in astrophysical environments, it is expected that there should be aromatic ringed molecules that incorporate N in the ring, but as yet, no such molecule has been identified. To address this issue, the present study uses high levels of electronic structure theory to compute a highly accurate quartic force field (QFF) for azirinyl cation and its two lowest lying isomers, the cyanomethyl and isocyanomethyl cations. The theoretical approach uses the singles and doubles coupled-cluster method that includes a perturbative correction for connected triple excitations, CCSD(T), together with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core-correlation. The QFF is then used in a second order vibrational perturbation theory analysis (VPT2) to compute the fundamental vibrational frequencies and rovibrational spectroscopic constants for all three C2NH2+ isomers. The reliability of the VPT2 vibrational frequencies is tested by comparison to vibrational configuration interaction (VCI) calculations and excellent agreement is found between the two approaches. Fundamental vibrational frequencies and rovibrational spectroscopic constants for all singly substituted 13C, 15N, and D isotopologues are also reported. It is expected that the highly accurate spectroscopic data reported herein will be useful in the identification of these cations in high resolution experimental or astronomical observations.

Partha P Bera↗

Semi-analytical theory for long-term rotational motion of uncontrolled satellites

Analytical solutions/integrals of the free rotational motion of uncontrolled satellites which contain spinning rotors are employed in conjunction with a suitable analytical/numerical averaging approach to predict long-term attitude motion. This semi-analytical method of averaging is basically similar to the one used in purely analytical 'general perturbation' theories. However, in view of the complexities involved in analytical treatment of the problem, its application is primarily numerical in character. An example is presented illustrating the use of the theory to predict the attitude motion of an axisymmetric satellite containing symmetrical constant-speed rotors under gravity-gradient torque. The use of the theory is considerably more efficient than straightforward numerical integration of the complete set of equations of motion.

Lahr, B. S.↗

Cosmic backgrounds from primeval dust

A general formalism is developed for describing the spectrum and the anisotropy of radiation from primeval dust as a nearly isotropic cosmic background with fluctuations. The theory incorporates spatial and temporal variations in the density of both luminosity and dust and is sufficiently robust to treat emission by point sources, as well as smoother emission, in high and low optical depth situations. The fluctuations are calculated using linear perturbation theory for arbitrary dust and luminosity clumping, and are calculated nonperturbatively for compact IR sources using shot-noise models. For many realistic models, a large accompanying NIR background is predicted. It is shown that, while spectal information alone provides only coarse information about the dust abundance and type, the redshift, and the other parameters, it may provide a powerful probe of pregalactic history when used in combination with the anisotropy signal.

Bond, J. R.↗

Stationary state model for normal metal tunnel junction phenomena

A model applicable to normal metal tunnel junctions is presented. This model, referred to herein as the stationary state model, is an extension of the extended basis function theory of Kleinman and Duke (1972). Under the assumption that elastic tunneling is the dominant transport mechanism under static bias, the theory has been extended to include the case where there is an ac component of bias potential, and the fluctuation spectrum has been derived. In this approach structure eigenstates are used as the basis, allowing observables to be evaluated without recourse to the perturbation theory inherent in the transfer Hamiltonian model. Comparison is made to the appropriate results of the first order transfer Hamiltonian model, and it is found that there is close but not always exact agreement to lowest order in the tunneling exponential. The stationary state model should be accurate with large barrier transmission, as it includes all orders of the tunneling exponential. The model as presented here should be applicable to normal metal tunnel junctions, where elastic tunneling is the dominant transport mechanism (under static bias), up to infrared excitation frequencies.

Whiteley, S. R.↗

The partially averaged field approach to cosmic ray diffusion

The kinetic equation for particles interacting with turbulent fluctuations is derived by a new nonlinear technique which successfully corrects the difficulties associated with quasilinear theory. In this new method the effects of the fluctuations are evaluated along particle orbits which themselves include the effects of a statistically averaged subset of the possible configurations of the turbulence. The new method is illustrated by calculating the pitch angle diffusion coefficient D sub Mu Mu for particles interacting with slab model magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. Results are compared with those of quasilinear theory and also with those of Monte Carlo calculations. The major effect of the nonlinear treatment in this illustration is the determination of D sub Mu Mu in the vicinity of 90 deg pitch angles where quasilinear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated using D sub Mu Mu as calculated by this technique. It is argued that the partially averaged field method is not limited to small amplitude fluctuating fields and is hence not a perturbation theory.

Jones, F. C.↗

Partially averaged field approach to cosmic ray diffusion

A new nonlinear technique is used to derive the kinetic equation for particles interacting with turbulent fluctuations. Difficulties associated with quasi-linear theory are avoided. The new method evaluates the effects of the fluctuations along particle orbits which themselves include the effects of a statistically averaged subset of the possible configuration of the turbulence. As an illustration, the pitch-angle diffusion coefficient is calculated for particles interacting with 'slab model' magnetic turbulence, i.e., magnetic fluctuations linearly polarized transverse to a mean magnetic field. The pitch-angle diffusion coefficient is determined in the vicinity of 90-deg pitch angles where quasi-linear theory breaks down. The spatial diffusion coefficient parallel to a mean magnetic field is evaluated by use of the calculated pitch-angle diffusion coefficient. It is suggested that the partially averaged field method is not limited to small amplitude fluctuating fields and hence is not a perturbation theory.

Jones, F. C.↗

Highly Accurate Quartic Force Fields, Vibrational Frequencies, and Spectroscopic Constants for Cyclic and Linear C3H3(+)

High levels of theory have been used to compute quartic force fields (QFFs) for the cyclic and linear forms of the C H + molecular cation, referred to as c-C H + and I-C H +. Specifically the 33 3333 singles and doubles coupled-cluster method that includes a perturbational estimate of connected triple excitations, CCSD(T), has been used in conjunction with extrapolation to the one-particle basis set limit and corrections for scalar relativity and core correlation have been included. The QFFs have been used to compute highly accurate fundamental vibrational frequencies and other spectroscopic constants using both vibrational 2nd-order perturbation theory and variational methods to solve the nuclear Schroedinger equation. Agreement between our best computed fundamental vibrational frequencies and recent infrared photodissociation experiments is reasonable for most bands, but there are a few exceptions. Possible sources for the discrepancies are discussed. We determine the energy difference between the cyclic and linear forms of C H +, 33 obtaining 27.9 kcal/mol at 0 K, which should be the most reliable available. It is expected that the fundamental vibrational frequencies and spectroscopic constants presented here for c-C H + 33 and I-C H + are the most reliable available for the free gas-phase species and it is hoped that 33 these will be useful in the assignment of future high-resolution laboratory experiments or astronomical observations.

Huang, Xinchuan↗

Finite-field Calculations of Molecular Polarizabilities Using Field-Induced Polarization Functions: Second- and Third-order Perturbation Correlation Corrections to the Coupled Hartree-Fock Polarizability of H2O

Ordinary Rayleigh-Schroudinger perturbation theory with Moller-Plesset (RSMP) partitioning is used to calculate second- and third-order correlation corrections to the CHF polarizability and dipole moment of the water molecule by a finite-field procedure. Pade approximants are found to be useful in accelerating the convergence of the property perturbation expansions. Field-induced polarization functions suitable for polarizability calculations are determined. The average polarizability calculated, neglecting vibrational averaging, with Dunning's (9s5p/4s-4s2p/2s) contracted GTO basis set augmented by field-induced lslp2d/lp polarization functions is within 3 per cent of the experimental result. Correlation corrections to the dipole moment and polarizability of the water molecule calculated by the finite-field RSMP and single + double excitation CI(SDCI) methods for the same basis set are found to be in close agreement. The RSMP approach has the advantages of being size-consistent and of being capable of greater efficiency than the SCDI method. Comparative calculations carried out using Epstein-Nesbet partitioning show that through third order RSEN correlation perturbation expansions for the dipole moment and polarizability are less rapidly convergent than RSMP expansions. However, reasonable accord with RSMP results can be achieved by using Pade approximants to accelerate the convergence of RSEN energy perturbation expansions. The convergence of RSMP property correlation expansions based on the zeroth-order uncoupled-Hartree-Fock (UCHF) and coupled-Hartree-Fock (CHF) approximations are compared through third order. Whereas the CHF + RSMP expansions are for practical purposes fully converged, the UCHF + RSMP expansions are not adequately converged.

Langhoff, S. R.↗

A high level Ab initio study of the anionic hydrogen-bonded complexes FH-CN-, FH-NC-, H2O-CN- and H2O-NC-

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

A high-level ab initio study of the anionic hydrogen-bonded complexes FH-CN(-), FH-NC(-), H2O-CN(-), and H2O-NC(-)

HF, H2O, CN- and their hydrogen-bonded complexes were studied using state-of-the-art ab initio quantum mechanical methods. A large Gaussian one particle basis set consisting of triple zeta plus double polarization plus diffuse s and p functions (TZ2P + diffuse) was used. The theoretical methods employed include self consistent field, second order Moller-Plesset perturbation theory, singles and doubles configuration interaction theory and the singles and doubles coupled cluster approach. The FH-CN- and FH-NC- and H2O-CN-, H2O-NC- pairs of complexes are found to be essentially isoenergetic. The first pair of complexes are predicted to be bound by approx. 24 kcal/mole and the latter pair bound by approximately 15 kcal/mole. The ab initio binding energies are in good agreement with the experimental values. The two being shorter than the analogous C-N hydrogen bond. The infrared (IR) spectra of the two pairs of complexes are also very similar, though a severe perturbation of the potential energy surface by proton exchange means that the accurate prediction of the band center of the most intense IR mode requires a high level of electronic structure theory as well as a complete treatment of anharmonic effects. The bonding of anionic hydrogen-bonded complexes is discussed and contrasted with that of neutral hydrogen-bonded complexes.

Lee, Timothy J.↗

Cosmic ray sources, acceleration and propagation

A review is given of selected papers on the theory of cosmic ray (CR) propagation and acceleration. The high isotropy and a comparatively large age of galactic CR are explained by the effective interaction of relativistic particles with random and regular electromagnetic fields in interstellar medium. The kinetic theory of CR propagation in the Galaxy is formulated similarly to the elaborate theory of CR propagation in heliosphere. The substantial difference between these theories is explained by the necessity to take into account in some cases the collective effects due to a rather high density of relativisitc particles. In particular, the kinetic CR stream instability and the hydrodynamic Parker instability is studied. The interaction of relativistic particles with an ensemble of given weak random magnetic fields is calculated by perturbation theory. The theory of CR transfer is considered to be basically completed for this case. The main problem consists in poor information about the structure of the regular and the random galactic magnetic fields. An account is given of CR transfer in a turbulent medium.

Ptuskin, V. S.↗

An ab initio benchmark study of the H + CO --> HCO reaction

The H + CO --> HCO reaction has been characterized with correlation consistent basis sets at five levels of theory in order to benchmark the sensitivities of the barrier height and reaction ergicity to the one-electron and n-electron expansions of the electronic wave function. Single and multireference methods are compared and contrasted. The coupled cluster method RCCSD(T) was found to be in very good agreement with Davidson-corrected internally-contracted multireference configuration interaction (MRCI+Q). Second-order Moller-Plesset perturbation theory (MP2) was also employed. The estimated complete basis set (CBS) limits for the barrier height (in kcal/mol) for the five methods, including harmonic zero-point energy corrections, are MP2, 4.66; RCCSD, 4.78; RCCSD(T), 4.15; MRCI, 5.10; and MRCI+Q, 4.07. Similarly, the estimated CBS limits for the ergicity of the reaction are: MP2, -17.99; RCCSD, -13.34; RCCSD(T), -13.79; MRCI, -11.46; and MRCI+Q, -13.70. Additional basis set explorations for the RCCSD(T) method demonstrate that aug-cc-pVTZ sets, even with some functions removed, are sufficient to reproduce the CBS limits to within 0.1-0.3 kcal/mol.

NASA Discipline Exobiology↗

A review of averaging techniques and their application to orbit determination systems

The theory of numerical averaging and analytical averaging will be reviewed and the application of these techniques to orbit and parameter estimation problems will be presented. Comparisons will be made between utilizing mean elements versus tracking data as the observation types. Results will be presented comparing the accuracy and efficiency of the combined orbit estimation and orbit prediction problem using averaged equations of motion, the Cowell equations of motion and the Brouwer general perturbation theory. The problem of converting the averaged element space back to osculating element space for orbit operations will also be discussed.

Velez, C. E.↗

Instabilities in astrophysical jets. II - Numerical simulations of slab jets

Numerical simulations of an unstable supersonic slab-symmetric jet are described. The instabilities within the jet are characterized by growing internal body waves and their coupled surface waves that are also predicted in linear perturbation theory. The characteristic theory of fluid dynamics is used to help interpret the wave morphologies. It is demonstrated that these waves can be excited by imposing an arbitrary disturbance. From the numerical simulations, it is found that the sound waves propagating against the flow slow down as they propagate outward, and they grow in amplitude. These waves eventually disrupt the jet at a certain length. This disruption length is related to the jet Mach number and the perturbation intensity. Thus, the Mach number of a jet observed with a radio telescope can be estimated by measuring the disruption length and estimating the perturbation intensity. The jet Mach numbers in radio tailed sources determined in this way agree quite well with estimates from ram pressure bending arguments. The wiggles and flares observed in many extragalactic jets, especially in tailed radio sources, appear to be intimately related to instabilities and the jet disruption process.

Zhao, Jun-Hui↗

Modified SUMT for structural synthesis

The concepts of the singular perturbation theory are employed to modify the SUMT (Sequential Unconstrained Minimization Technique) algorithm of Fiacco-McCormick (1968) in order to make it more robust and reliable. The algorithm, which uses a sequence of penalty parameters to convert a constrained problem into a sequence of unconstrained problems, suffers from the need to minimize an ill-conditioned penalty function. By using the modified SUMT algorithm on two different structural optimization problems, it is shown that the singular perturbation SUMT easily converges to accurate solutions.

Kamat, M. P.↗

Modeling issues in precision orbit determination for Mars orbiter

This paper examines the accuracy of recent Mars gravity models and the importance of perturbations due to the Mars radiation pressure and the Martian moons, Phobos and Deimos, on the trajectories of Mars orbiters. A linear orbit perturbation theory is used to characterize the patterns of gravity field near resonances for the Viking and Mariner 9 spacecraft. These resonances are shown to have considerable power and their potential for contributing to Mars gravity solutions is emphasized. It is shown that some of the same resonance orders which appear in the Viking orbits, dominate the radial orbit error spectrum for Mars Observer. Results of orbit determination simulations at the Goddard Space Flight Center show that the perturbations caused by the Martian moons and the Mars radiation pressure are larger than 0.1 mm/s, the expected precision of the Mars Observer Doppler tracking data. Tests with the Viking Doppler data indicate that best analysis of these data mandates the inclusion of the Phobos gravitational perturbation in the modeling of Viking spacecraft trajectories.

Lemoine, Frank G.↗

Non-stationary and non-linear dispersive medium as external field which generates the squeezed states

The theory of quantum effects in nonlinear dielectric media is developed. The nonlinear dielectric media is influenced by an external pumping field. The diagonalization of the Hamiltonian of a quantized field is obtained by the canonical Bogoliubov transformations. The transformations allow us to obtain the general expressions for the number of created photons and for the degree of squeezing. In the case of a plane pumping wave, for example, the results are calculated by using the zero order of the secular perturbation theory, with small parameters characterizing the medium nonlinearity. The Heisenberg equations of motion are obtained for non-stationary case and a commonly used Hamiltonian is derived from the first principles of quantum electrodynamics.

Lobashov, A. A.↗