Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Semiclassical methods”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Analytical study of bound states in graphene nanoribbons and carbon nanotubes: The variable phase method and the relativistic Levinson theorem

The problem of localized states in 1D systems with a relativistic spectrum, namely, graphene stripes and carbon nanotubes, is studied analytically. The bound state as a superposition of two chiral states is completely described by their relative phase, which is the foundation of the variable phase method (VPM) developed herein. Based on our VPM, we formulate and prove the relativistic Levinson theorem. The problem of bound states can be reduced to the analysis of closed trajectories of some vector field. Remarkably, the Levinson theorem appears as the Poincaré index theorem for these closed trajectories. The VPM equation is also reduced to the nonrelativistic and semiclassical limits. The limit of a small momentum p{sub y} of transverse quantization is applicable to an arbitrary integrable potential. In this case, a single confined mode is predicted.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2001

Nano-Transistor Modeling: Two Dimensional Green's Function Method

Two quantum mechanical effects that impact the operation of nanoscale transistors are inversion layer energy quantization and ballistic transport. While the qualitative effects of these features are reasonably understood, a comprehensive study of device physics in two dimensions is lacking. Our work addresses this shortcoming and provides: (a) a framework to quantitatively explore device physics issues such as the source-drain and gate leakage currents, DIBL (Drain Induced Barrier Lowering), and threshold voltage shift due to quantization, and b) a means of benchmarking quantum corrections to semiclassical models (such as density-gradient and quantum-corrected MEDICI).

Svizhenko, Alexei↗
DOE OSTI2014

Cross Sections for Inner-Shell Ionization by Electron Impact

An analysis is presented of measured and calculated cross sections for inner-shell ionization by electron impact. We describe the essentials of classical and semiclassical models and of quantum approximations for computing ionization cross sections. The emphasis is on the recent formulation of the distorted-wave Born approximation by Bote and Salvat [Phys. Rev. A 77, 042701 (2008)] that has been used to generate an extensive database of cross sections for the ionization of the K shell and the L and M subshells of all elements from hydrogen to einsteinium (Z = 1 to Z = 99) by electrons and positrons with kinetic energies up to 1 GeV. We describe a systematic method for evaluating cross sections for emission of x rays and Auger electrons based on atomic transition probabilities from the Evaluated Atomic Data Library of Perkins et al. [Lawrence Livermore National Laboratory, UCRL-ID-50400, 1991]. We made an extensive comparison of measured K-shell, L-subshell, and M-subshell ionization cross sections and of Lα x-ray production cross sections with the corresponding calculated cross sections. We identified elements for which there were at least three (for K shells) or two (for L and M subshells) mutually consistent sets of cross-section measurements and for which the cross sections varied with energy as expected by theory. The overall average root-mean-square deviation between the measured and calculated cross sections was 10.9% and the overall average deviation was −2.5%. This degree of agreement between measured and calculated ionization and x-ray production cross sections was considered to be very satisfactory given the difficulties of these measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Level density within a micro-macroscopic approach

Statistical level density $ρ$($E, A$) is derived for nucleonic system with a given energy $E$, particle number $A$ and other integrals of motion in the micro-macroscopic approximation beyond the standard saddle-point method of the Fermi gas model. This level density reaches the two limits; the well-known Fermi gas grand-canonical ensemble limit for a large entropy $S$ related to large excitation energies, and the finite micro-canonical limit for a small combinatorical entropy $S$ at low excitation energies. In conclusion, the inverse level density parameter $K$ as function of the particle number $A$ in the semiclassical periodic orbit theory, taking into account the extended Thomas-Fermi and Strutinsky shell corrections, is calculated and compared with experimental data.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Gate-Tunable Multiband Transport in ZrTe 5 Thin Devices

Interest in ZrTe 5 has been reinvigorated in recent years owing to its potential for hosting versatile topological electronic states and intriguing experimental discoveries. However, the mechanism of many of its unusual transport behaviors remains controversial: for example, the characteristic peak in the temperature-dependent resistivity and the anomalous Hall effect. Here, through employing a clean dry-transfer fabrication method in an inert environment, we successfully obtain high-quality ZrTe 5 thin devices that exhibit clear dual-gate tunability and ambipolar field effects. Such devices allow us to systematically study the resistance peak as well as the Hall effect at various doping densities and temperatures, revealing the contribution from electron–hole asymmetry and multiple-carrier transport. By comparing with theoretical calculations, we suggest a simplified semiclassical two-band model to explain the experimental observations. Finally, our work helps to resolve the longstanding puzzles on ZrTe 5 and could potentially pave the way for realizing novel topological states in the two-dimensional limit.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1979

The adiabatic semiclassical perturbation theory for vibrationally inelastic scattering. I - Collinear calculations. II - Three-dimensional treatment

A semiclassical approximation to treat vibrationally inelastic scattering is developed. The vibrational basis set used is adiabatic with respect to a reference potential which is chosen to be as close as possible to the true potential and also gives easily obtainable solutions to the vibrational wave equation. The radial wave functions are obtained using the WKB approximation, and the coupled Schroedinger equations are solved by a first-order perturbation method to yield a phase shift matrix which is exponentiated to give the full scattering matrix. Results were obtained for all the cases computed by Secrest and Johnson and by Clark and Dickinson, and the agreement is better than 10% for half of the cross-sections and rarely off by more than a factor of 2.

Cross, R. J., Jr.↗
DOE OSTI2018

Extended Hamiltonian molecular dynamics: semiclassical trajectories with improved maintenance of zero point energy

It is well known that classical trajectories, even if they are initiated with zero point energy (ZPE) in each mode (trajectories initiated this way are commonly called quasiclassical trajectories), do not maintain ZPE in the final states. The energy of high-frequency modes will typically leak into low-frequency modes or relative translation of subsystems during the time evolution. This can lead to severe problems such as unphysical dissociation of a molecule, production of energetically disallowed reaction products, and unphysical product energy distributions. Here a new molecular dynamics method called extended Hamiltonian molecular dynamics (EHMD) is developed to improve the ZPE problem in classical molecular dynamics. In EHMD, two images of a trajectory are connected by one or more springs. The EHMD method is tested with the Henon–Heiles Hamiltonian in reduced and real units and with a Hamiltonian with quartic anharmonicity in real units, and the method is found to improve zero-point maintenance as intended.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Understanding the nature of mean-field semiclassical light-matter dynamics: An investigation of energy transfer, electron-electron correlations, external driving, and long-time detailed balance

Semiclassical electrodynamics (with quantum matter plus classical electrodynamics fields) is an appealing approach for studying light-matter interactions, especially for realistic molecular systems. However, there is no unique semiclassical scheme. On the one hand, intermolecular interactions can be described instantaneously by static two-body interactions connecting two different molecules, while a classical transverse E field acts as a spectator at short distance; we will call this Hamiltonian no. I. On the other hand, intermolecular interactions can also be described as effects that are mediated exclusively through a classical one-body Ε field without any quantum effects at all (assuming we ignore electronic exchange); we will call this Hamiltonian no. II. Moreover, one can also mix these two different Hamiltonians into a third, hybrid Hamiltonian, which preserves quantum electron-electron correlations for lower excitations but describes higher excitations in a mean-field way. To investigate which semiclassical scheme is most reliable for practical use, here we study the real-time dynamics of a minimalistic many-site model—a pair of identical two-level systems (TLSs)—undergoing either resonance energy transfer (RET) or collectively driven dynamics. While both approaches (no. 1 and no. 2) perform reasonably well when there is no strong external excitation, we find that no single approach is perfect for all conditions (and all methods fail when a strong external field is applied). Each method has its own distinct problems: Hamiltonian no. I performs best for RET but behaves in a complicated manner for driven dynamics; Hamiltonian no. II is always stable, but obviously fails for RET at short distances. One key finding is that, for externally driven dynamics, a full configuration-interaction description of Hamiltonian no. I strongly overestimates the long-time electronic energy, highlighting the not obvious fact that, if one plans to merge quantum molecules with classical light, a full, exact treatment of electron-electron correlations can actually lead to worse results than a simple mean-field electronic structure treatment. Future work will need to investigate (i) how these algorithms behave in the context of more than a pair of TLSs and (ii) whether or not these algorithms can be improved in general by including crucial aspects of spontaneous emission.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1979

Proton collisional excitation in the lowest lying 3P terms of ions in the Be and Mg isoelectronic sequences

Proton collisional excitation cross sections and rate constants are presented for transitions between the 3P(J) fine-structure levels of the lowest-lying sp configurations in a number of astrophysically important ions belonging to the Be and Mg isoelectronic sequences. The calculations were made by direct integration of the Schroedinger equation resulting from semiclassical Coulomb excitation theory. The cross sections and rate constants for the 3P(J) transitions in the lowest-lying P(2) configurations are expected to be similar to those for the corresponding sp configuration transitions, and this is illustrated for C III. For the high-temperature ion Ca XVII alpha particle excitation is shown to be unimportant for situations involving ordinary values of the He/H abundance ratio. A simple, but apparently accurate method for determining certain radial integrals for low-lying excited configurations is proposed.

Landman, D. A.↗
DOE OSTI2013

Semiclassical Monte-Carlo approach for modelling non-adiabatic dynamics in extended molecules

Modelling of non-adiabatic dynamics in extended molecular systems and solids is a next frontier of atomistic electronic structure theory. The underlying numerical algorithms should operate only with a few quantities (that can be efficiently obtained from quantum chemistry), provide a controlled approximation (which can be systematically improved) and capture important phenomena such as branching (multiple products), detailed balance and evolution of electronic coherences. Here we propose a new algorithm based on Monte-Carlo sampling of classical trajectories, which satisfies the above requirements and provides a general framework for existing surface hopping methods for non-adiabatic dynamics simulations. In particular, our algorithm can be viewed as a post-processing technique for analysing numerical results obtained from the conventional surface hopping approaches. Presented numerical tests for several model problems demonstrate efficiency and accuracy of the new method.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2015

Multidimensional effects in nonadiabatic statistical theories of spin- forbidden kinetics. A case study of 3 O + CO → CO 2

The appropriateness of treating crossing seams of electronic states of different spins as nonadiabatic transition states in statistical calculations of spin-forbidden reaction rates is considered. We show that the spin-forbidden reaction coordinate, the nuclear coordinate perpendicular to the crossing seam, is coupled to the remaining nuclear degrees of freedom. We found that this coupling gives rise to multidimensional effects that are not typically included in statistical treatments of spin-forbidden kinetics. Three qualitative categories of multidimensional effects may be identified: static multidimensional effects due to the geometry-dependence of the local shape of the crossing seam and of the spin–orbit coupling, dynamical multidimensional effects due to energy exchange with the reaction coordinate during the seam crossing, and nonlocal(history-dependent) multidimensional effects due to interference of the electronic variables at second, third, and later seam crossings. Nonlocal multidimensional effects are intimately related to electronic decoherence, where electronic dephasing acts to erase the history of the system. A semiclassical model based on short-time full-dimensional trajectories that includes all three multidimensional effects as well as a model for electronic decoherence is presented. The results of this multidimensional nonadiabatic statistical theory (MNST) for the 3 O + CO → CO 2 reaction are compared with the results of statistical theories employing one-dimensional (Landau–Zener and weak coupling) models for the transition probability and with those calculated previously using multistate trajectories. The MNST method is shown to accurately reproduce the multistate decay-of-mixing trajectory results, so long as consistent thresholds are used. Furthermore, the MNST approach has several advantages over multistate trajectory approaches and is more suitable in chemical kinetics calculations at low temperatures and for complex systems. The error in statistical calculations that neglect multidimensional effects is shown to be as large as a factor of 2 for this system, with static multidimensional effects identified as the largest source of error.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Benchmarking semiclassical approaches to strong-field QED: Nonlinear Compton scattering in intense laser pulses

The recoil associated with photon emission is key to the dynamics of ultrarelativistic electrons in strong electromagnetic fields, as found in high-intensity laser-matter interactions and astrophysical environments such as neutron star magnetospheres. When the energy of the photon becomes comparable to that of the electron, it is necessary to use quantum electrodynamics (QED) to describe the dynamics accurately. However, computing the appropriate scattering matrix element from strong-field QED is not generally possible due to multiparticle effects and the complex structure of the electromagnetic fields. Therefore, these interactions are treated semiclassically, coupling probabilistic emission events to classical electrodynamics using rates calculated in the locally constant field approximation. Here, we provide comprehensive benchmarking of this approach against the exact QED calculation for nonlinear Compton scattering of electrons in an intense laser pulse. We find agreement at the percentage level between the photon spectra, as well as between the models' predictions of absorption from the background field, for normalized amplitudes a0 > 5. We discuss possible routes towards improved numerical methods and the implications of our results for the study of QED cascades.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Metaplectic geometrical optics for ray-based modeling of caustics: Theory and algorithms

The optimization of radio frequency-wave (RF) systems for fusion experiments is often performed using ray-tracing codes, which rely on the geometrical-optics (GO) approximation. However, GO fails at caustics such as cutoffs and focal points, erroneously predicting the wave intensity to be infinite. This is a critical shortcoming of GO, since the caustic wave intensity is often the quantity of interest, e.g., RF heating. Full-wave modeling can be used instead, but the computational cost limits the speed at which such optimizations can be performed. Here, we have developed a less expensive alternative called metaplectic geometrical optics (MGO). Instead of evolving waves in the usual x (coordinate) or k (spectral) representation, MGO uses a mixed X$\equiv$Ax+Bk representation. By continuously adjusting the matrix coefficients A and B along the rays, one can ensure that GO remains valid in the X coordinates without caustic singularities. The caustic-free result is then mapped back onto the original x space using metaplectic transforms. Here, we overview the MGO theory and review algorithms that will aid the development of an MGO-based ray-tracing code. We show how using orthosymplectic transformations leads to considerable simplifications compared to previously published MGO formulas. We also prove explicitly that MGO exactly reproduces standard GO when evaluated far from caustics (an important property that until now has only been inferred from numerical simulations), and we relate MGO to other semiclassical caustic-removal schemes published in the literature. Finally this discussion is then augmented by an explicit comparison of the computed spectrum for a wave bounded between two cutoffs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

Electron transmission through a class of anthracene aldehyde molecules

Transmission of electrons via metal-molecule-metal junctions, involving rotor-stator anthracene aldehyde molecules is investigated. Two model barriers having input parameters evaluated from accurate ab initio calculations are proposed and the transmission coefficients are obtained by using the quasiclassical approximation. Transmission coefficients further enter in the integral for the net current, utilizing Simmons’ method. Conformational dependence of the tunneling processes is evident and the presence of the side groups enhances the functionality of the future single-molecule based electronic devices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2016

Classical molecular dynamics simulation of electronically non-adiabatic processes

Both classical and quantum mechanics (as well as hybrids thereof, i.e., semiclassical approaches) find widespread use in simulating dynamical processes in molecular systems. For large chemical systems, however, which involve potential energy surfaces (PES) of general/arbitrary form, it is usually the case that only classical molecular dynamics (MD) approaches are feasible, and their use is thus ubiquitous nowadays, at least for chemical processes involving dynamics on a single PES (i.e., within a single Born-Oppenheimer electronic state). Our paper reviews recent developments in an approach which extends standard classical MD methods to the treatment of electronically non-adiabatic processes, i.e., those that involve transitions between different electronic states. This approach treats nuclear and electronic degrees of freedom (DOF) equivalently (i.e., by classical mechanics, thereby retaining the simplicity of standard MD), and provides "quantization" of the electronic states through a symmetrical quasi-classical (SQC) windowing model. The approach is seen to be capable of treating extreme regimes of strong and weak coupling between the electronic states, as well as accurately describing coherence effects in the electronic DOF (including the de-coherence of such effects caused by coupling to the nuclear DOF). A survey of recent applications is presented to illustrate the performance of the approach. We also described is a newly developed variation on the original SQC model (found universally superior to the original) and a general extension of the SQC model to obtain the full electronic density matrix (at no additional cost/complexity).

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2019

Combining the mapping Hamiltonian linearized semiclassical approach with the generalized quantum master equation to simulate electronically nonadiabatic molecular dynamics

The generalized quantum master equation (GQME) provides a powerful framework for simulating electronically nonadiabatic molecular dynamics. Within this framework, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density matrix is fully captured by a memory kernel superoperator. In this paper, we consider two different procedures for calculating the memory kernel of the GQME from projection-free inputs obtained via the combination of the mapping Hamiltonian (MH) approach and the linearized semiclassical (LSC) approximation. Here, the accuracy and feasibility of the two procedures are demonstrated on the spin-boson model. We find that although simulating the electronic dynamics by direct application of the two LSC-based procedures leads to qualitatively different results that become increasingly less accurate with increasing time, restricting their use to calculating the memory kernel leads to an accurate description of the electronic dynamics. Comparison with a previously proposed procedure for calculating the memory kernel via the Ehrenfest method reveals that MH/LSC methods produce memory kernels that are better behaved at long times and lead to more accurate electronic dynamics

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1974

Quantum graviton creation in a model universe

Consideration of the mechanism of production of gravitons in the empty, anisotropic, spatially inhomogeneous Gowdy three-torus cosmology. The Gowdy cosmology is an exact solution of the vacuum Einstein equations and is obtained as a generalization of the homogeneous empty Bianchi Type I (Kasner) cosmology by permitting the metric components to depend on one of the space variables in addition to time. The Hamiltonian methods of Arnowitt, Deser, and Misner are employed to identify the dynamical variables which are to be quantized. The WKB regime solution is identical to that found by Doroshkevich, Zel'dovich, and Novikov (DZN) for a universe containing collisionless anisotropic radiation. Using a procedure similar to that of Parker (1971) or Zel'dovich and Starobinskii (1971) for defining quantum number, it is found that the DZN large-time radiation consists of quanta (gravitons) created from an initial vacuum. The quantum behavior is much like the semiclassical enhancement of quantum number with the added feature of creation of quanta from vacuum fluctuations.

Berger, B. K.↗
DOE OSTI2018

Entry-level spin distributions and relative γ -neutron branching ratios of samarium isotopes populated by the ( p , t ) reaction

The spin distribution populated in the $^{150}\mathrm{Sm}(p,t)$ reaction with 25 MeV protons has been investigated here using the silicon telescope and high-purity germanium array Hyperion. Angular distributions of outgoing tritons and relative intensities of ${\gamma}$ rays from low-lying levels are used to deduce the initial spin distribution using statistical calculations and are compared with a semiclassical model of the $(p,t)$ reaction. Experimental and theoretical results suggest a peak strength near $4{\hbar}$ when gated on outgoing tritons detected in Hyperion, while the full $4{\pi}$ distribution may peak near $2{\hbar}$. Additionally, the average ${\gamma}$-neutron branching ratio is determined from the data with the aid of one of two new codes presented in this work. Finally, the analysis is extended to other recent studies of $^{152,154}\mathrm{Sm}(p,t)$ to form systematics of the spin-dependence of the surrogate method applied to heavy nuclei in the transitional to deformed region.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗