Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Schrodinger equations”

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 37 records · Page 2

Numerical Simulations of Light Bullets, Using the Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that we currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Karr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗

Numerical Simulations of Light Bullets, Using The Full Vector, Time Dependent, Nonlinear Maxwell Equations

This paper will present results in computational nonlinear optics. An algorithm will be described that solves the full vector nonlinear Maxwell's equations exactly without the approximations that are currently made. Present methods solve a reduced scalar wave equation, namely the nonlinear Schrodinger equation, and neglect the optical carrier. Also, results will be shown of calculations of 2-D electromagnetic nonlinear waves computed by directly integrating in time the nonlinear vector Maxwell's equations. The results will include simulations of 'light bullet' like pulses. Here diffraction and dispersion will be counteracted by nonlinear effects. The time integration efficiently implements linear and nonlinear convolutions for the electric polarization, and can take into account such quantum effects as Kerr and Raman interactions. The present approach is robust and should permit modeling 2-D and 3-D optical soliton propagation, scattering, and switching directly from the full-vector Maxwell's equations.

Goorjian, Peter M.↗

A method of solving simple harmonic oscillator Schroedinger equation

A usual step in solving totally Schrodinger equation is to try first the case when dimensionless position independent variable w is large. In this case the Harmonic Oscillator equation takes the form (d(exp 2)/dw(exp 2) - w(exp 2))F = 0, and following W.K.B. method, it gives the intermediate corresponding solution F = exp(-w(exp 2)/2), which actually satisfies exactly another equation, (d(exp 2)/dw(exp 2) + 1 - w(exp 2))F = 0. We apply a different method, useful in anharmonic oscillator equations, similar to that of Rampal and Datta, and although it is slightly more complicated however it is also more general and systematic.

Maury, Juan Carlos F.↗

A study of shock-associated magnetohydrodynamic waves in the solar wind

Three major topics were addressed, one theoretical and two observational. The topics were: (1) an attempt to understand the evolution of the large-amplitude magnetohydrodynamic (MHD) waves in the foreshock, using a nonlinear wave equation called the Derivative Nonlinear Schrodinger equation (henceforth DNLS) as a model, (2) using the extensive set of ISE data to test for the presence of various nonlinear wave processes which might be present, and (3) a study of plasma turbulence in the interstellar medium which might be physically similar to that in the solar wind. For these investigations we used radioastronomical techniques. Good progress was made in each of these areas and a separate discussion of each is given.

Spangler, Steven R.↗

Berry phase in Heisenberg representation

We define the Berry phase for the Heisenberg operators. This definition is motivated by the calculation of the phase shifts by different techniques. These techniques are: the solution of the Heisenberg equations of motion, the solution of the Schrodinger equation in coherent-state representation, and the direct computation of the evolution operator. Our definition of the Berry phase in the Heisenberg representation is consistent with the underlying supersymmetry of the model in the following sense. The structural blocks of the Hamiltonians of supersymmetrical quantum mechanics ('superpairs') are connected by transformations which conserve the similarity in structure of the energy levels of superpairs. These transformations include transformation of phase of the creation-annihilation operators, which are generated by adiabatic cyclic evolution of the parameters of the system.

Andreev, V. A.↗

Resonant double-core excitations with ultrafast, intense X-ray pulses

Intense few-to-sub-femtosecond soft X-ray pulses can produce neutral, two-site excited double-core-hole states by promoting two core electrons to the same unoccupied molecular orbital. We theoretically investigate double nitrogen K-edge excitations of nitrous oxide (N$_2$O) with multiconfigurational electronic structure calculations. We show that the second core-excitation energy is reduced with respect to its ground state value. A site-selective double core-excitation mechanism using intense few-femtosecond x-rays is investigated using time-dependent Schrodinger equation (TDSE) simulations. The subsequent two-step Auger-Meitner and two-electrons-one-electron decay spectra of the double core-excited states are analysed using a Mulliken population analysis of the multiconfirational wavefunctions. The change in the electron emission lineshape between the absorption of 1 or 2 photons in the resonant core-excitation is predicted by combining this approach with the TDSE simulations. Furthermore, we examine the possibility of resolving the double core-excited states with X-ray pump-probe techniques by calculating the chemical shifts of the core-electron binding energy of the core-excited states and decay products.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Reduction of the equation for lower hybrid waves in a plasma to a nonlinear Schroedinger equation

Equations describing the nonlinear propagation of waves in an anisotropic plasma are rarely exactly soluble. However it is often possible to make approximations that reduce the exact equations into a simpler equation. The use of MACSYMA to make such approximations, and so reduce the equation describing lower hybrid waves into the nonlinear Schrodinger equation which is soluble by the inverse scattering method is demonstrated. MACSYMA is used at several stages in the calculation only because there is a natural division between calculations that are easiest done by hand, and those that are easiest done by machine.

Karney, C. F. F.↗

Self-focusing of whistler waves

The theory of axially symmetric self-focusing of whistler waves, based on the full system of Maxwell equations, is developed. The plasma is described by the magnetohydrodynamic equations including the ponderomotive force from RF field. The nonlinear Schrodinger equations (NSE) for arbitrary azimuthal modes of whistler waves are derived. It is shown that they differ from the NSE for a scalar field; this is connected with an intrinsic angular momentum due to the rotating polarization of whistlers. It is shown that the self-focusing, as described by the NSE, differs in its final stage from the results following the full set of Maxwell equations. The latter gives defocusing after sufficient narrowing of the initial wave beam, due to transformation of the trapped wave into a nontrapped branch which is not contained in the NSE description. The oscillatory character of the defocusing is demonstrated.

Karpman, V. I.↗

Coherence Phenomena in Coupled Optical Resonators

We predict a variety of photonic coherence phenomena in passive and active coupled ring resonators. Specifically, the effective dispersive and absorptive steady-state response of coupled resonators is derived, and used to determine the conditions for coupled-resonator-induced transparency and absorption, lasing without gain, and cooperative cavity emission. These effects rely on coherent photon trapping, in direct analogy with coherent population trapping phenomena in atomic systems. We also demonstrate that the coupled-mode equations are formally identical to the two-level atom Schrodinger equation in the rotating-wave approximation, and use this result for the analysis of coupled-resonator photon dynamics. Notably, because these effects are predicted directly from coupled-mode theory, they are not unique to atoms, but rather are fundamental to systems of coherently coupled resonators.

Smith, D. D.↗

FD-TD modeling of 2-D dielectric waveguides for propagation and scattering of femtosecond optical solitons

Experimentalists have produced all-optical switches capable of 100-fs responses. To adequately model such switches, nonlinear effects in optical materials (both instantaneous and dispersive) must be included. In principle, the behavior of electromagnetic fields in nonlinear dielectrics can be determined by solving Maxwell's equations subject to the assumption that the electric polarization has a nonlinear relation to the electric field. However, until our previous work, the resulting nonlinear Maxwell's equations have not been solved directly. Rather, approximations have been made that result in a class of generalized nonlinear Schrodinger equations (GNLSE) that solve only for the envelope of the optical pulses. In this paper, we present first-time calculations from the vector nonlinear Maxwell's equations of femtosecond soliton propagation and scattering, including carrier waves, in two-dimensional systems of dielectric waveguides exhibiting the Kerr and Raman quantum effects. We use the finite-difference time-domain (FD-TD) method in an extension of our 1-D work. There, in a fundamental innovation, we treated the linear and nonlinear convolutions for the electric polarization as new dependent variables. By differentiating these convolutions in the time domain, we derived an equivalent system of coupled, nonlinear second-order ODE's. These equations together with Maxwell's equations form the system that is solved to determine the electromagnetic fields in inhomogeneous nonlinear dispersive media. Backstorage in time is limited to only that needed by the time-integration algorithm for the ODE's, rather than that needed to store the time-history of the kernel functions of the convolutions (1000-10,000 time steps). Thus, a 2-D nonlinear optics model from Maxwell's equations is now feasible.

Joseph, Rose↗

IMEX Runge-Kutta Parareal for Non-diffusive Equations

Parareal is a widely studied parallel-in-time method that can achieve meaningful speedup on certain problems. However, it is well known that the method typically performs poorly on non-diffusive equations. This paper analyzes linear stability and convergence for IMEX Runge-Kutta Parareal methods on non-diffusive equations. By combining standard linear stability analysis with a simple convergence analysis, we find that certain Parareal configurations can achieve parallel speedup on non-diffusive equations. These stable configurations all posses low iteration counts, large block sizes, and a large number of processors. Numerical examples using the nonlinear Schrodinger equation demonstrate the analytical conclusions.

97 MATHEMATICS AND COMPUTING↗

The equilibrium and stability of the gaseous component of the galaxy, 2

A time-independent, linear, plane and axially-symmetric stability analysis was performed on a self-gravitating, plane-parallel, isothermal layer of nonmagnetic, nonrotating gas. The gas layer was immersed in a plane-stratified field isothermal layer of stars which supply a self-consistent gravitational field. Only the gaseous component was perturbed. Expressions were derived for the perturbed gas potential and perturbed gas density that satisfied both the Poisson and hydrostatic equilibrium equations. The equation governing the size of the perturbations in the mid-plane was found to be analogous to the one-dimensional time-independent Schrodinger equation for a particle bound by a potential well, and with similar boundary conditions. The radius of the neutral state was computed numerically and compared with the Jeans' and Ledoux radius. The inclusion of a rigid stellar component increased the Ledoux radius, though only slightly. Isodensity contours of the neutrual or marginally unstable state were constructed.

Kellman, S. A.↗

Application of the Finite Element Method in Atomic and Molecular Physics

The finite element method (FEM) is a numerical algorithm for solving second order differential equations. It has been successfully used to solve many problems in atomic and molecular physics, including bound state and scattering calculations. To illustrate the diversity of the method, we present here details of two applications. First, we calculate the non-adiabatic dipole polarizability of Hi by directly solving the first and second order equations of perturbation theory with FEM. In the second application, we calculate the scattering amplitude for e-H scattering (without partial wave analysis) by reducing the Schrodinger equation to set of integro-differential equations, which are then solved with FEM.

Shertzer, Janine↗

Multiple soliton production and the Korteweg-de Vries equation.

Compressive square-wave pulses are launched in a double-plasma device. Their evolution is interpreted according to the Korteweg-de Vries description of Washimi and Taniuti. Square-wave pulses are an excitation for which an explicit solution of the Schrodinger equation permits an analytical prediction of the number and amplitude of emergent solitons. Bursts of energetic particles (pseudowaves) appear above excitation voltages greater than an electron thermal energy, and may be mistaken for solitons.

Hershkowitz, N.↗

Boost-phase discrimination research activities

Theoretical research in two areas was performed. The aerothermodynamics research focused on the hard-body and rocket plume flows. Analytical real gas models to describe finite rate chemistry were developed and incorporated into the three-dimensional flow codes. New numerical algorithms capable of treating multi-species reacting gas equations and treating flows with large gradients were also developed. The computational chemistry research focused on the determination of spectral radiative intensity factors, transport properties and reaction rates. Ab initio solutions to the Schrodinger equation provided potential energy curves transition moments (radiative probabilities and strengths) and potential energy surfaces. These surfaces were then coupled with classical particle reactive trajectories to compute reaction cross-sections and rates.

Cooper, David M.↗

Static and dynamic properties of atomic nuclei with high-resolution potentials

Here, we compute ground-state and dynamical properties of 4 He and 16 O nuclei using as input high-resolution, phenomenological nucleon-nucleon and three-nucleon forces that are local in coordinate space. The nuclear Schrodinger equation for both nuclei is accurately solved employing the auxiliary-field diffusion Monte Carlo approach. For the 4 He nucleus, detailed benchmarks are carried out with the hyperspherical harmonics method. In addition to presenting results for the binding energies and radii, we also analyze the momentum distributions of these nuclei and their Euclidean response function corresponding to the isoscalar density transition. The latter quantity is particularly relevant for lepton-nucleus scattering experiments, as it paves the way to quantum Monte Carlo calculations of electroweak response functions of 16 O.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗