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.

234 records · Page 13
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Semiclassical Calculation of Reaction Rate Constants for Homolytical Dissociations

There is growing interest in extending organometallic chemical vapor deposition (OMCVD) to III-V materials that exhibit large thermal decomposition at their optimum growth temperature, such as indium nitride. The group III nitrides are candidate materials for light-emitting diodes and semiconductor lasers operating into the blue and ultraviolet regions. To overcome decomposition of the deposited compound, the reaction must be conducted at high pressures, which causes problems of uniformity. Microgravity may provide the venue for maintaining conditions of laminar flow under high pressure. Since the selection of optimized parameters becomes crucial when performing experiments in microgravity, efforts are presently geared to the development of computational OMCVD models that will couple the reactor fluid dynamics with its chemical kinetics. In the present study, we developed a method to calculate reaction rate constants for the homolytic dissociation of III-V compounds for modeling OMCVD. The method is validated by comparing calculations with experimental reaction rate constants.

Beatriz H. Cardelino↗
DOE OSTI2024

Semiclassical treatment of photon cascades in nuclei

Here, we present a simple semiclassical treatment of photon cascades, suitable for use in nuclear fission simulation codes. The approximation is here developed for E⁢1 and E⁢2 transitions and its quality is illustrated for a variety of two-photon cascades. Implementation of the treatment into Monte Carlo simulations would make it possible to address photon correlation observables quantitatively.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2012

Electron Stark Broadening Database for Atomic N, O, and C Lines

A database for efficiently computing the electron Stark broadening line widths for atomic N, O, and C lines is constructed. The line width is expressed in terms of the electron number density and electronatom scattering cross sections based on the Baranger impact theory. The state-to-state cross sections are computed using the semiclassical approximation, in which the atom is treated quantum mechanically whereas the motion of the free electron follows a classical trajectory. These state-to-state cross sections are calculated based on newly compiled line lists. Each atomic line list consists of a careful merger of NIST, Vanderbilt, and TOPbase line datasets from wavelength 50 nm to 50 micrometers covering the VUV to IR spectral regions. There are over 10,000 lines in each atomic line list. The widths for each line are computed at 13 electron temperatures between 1,000 K 50,000 K. A linear least squares method using a four-term fractional power series is then employed to obtain an analytical fit for each line-width variation as a function of the electron temperature. The maximum L2 error of the analytic fits for all lines in our line lists is about 5%.

Liu, Yen↗
NASA NTRS1995

Rototranslational collision-induced absorption by H2-H2 pairs at temperatures from 600 to 7000 K

The computation of the far-infrared, rototranslational (RT) collision-induced absorption (CIA) spectra of H2-H2 pairs is presented at temperatures from 600 to 7000 K for the first time. Theoretical results are based on the quantum mechanical and semiclassical, three lowest translational spectral moments obtained for H2 pairs. The effective, isotropic H2-H2 interaction potential, suitable for the high-temperature computations, and the ab initio induced dipoles, have been used as input. Special effort has been made to account for the rotational and vibrational states dependence of the dipoles, since it was found to be relevant at the high temperatures employed. The computations of the entire RT band account for all populated vibrational states of hydrogen molecule and include vibrational transitions v tends towards v-prime = v, with v = 0, 1, 2 and 3. The described method makes use of the adequately selected model line shapes with the temperature-dependent parameters. The presented model is useful for the 'model atmospheres' of zero- and low-metallicity, cool and dense stellar atmospheres, where CIA is known to be imporatnt.

Zheng, Chunguang↗
DOE OSTI2018

Cluster truncated Wigner approximation in strongly interacting systems

We present a general method by which linear quantum Hamiltonian dynamics with exponentially many degrees of freedom is replaced by approximate classical nonlinear dynamics with the number of degrees of freedom (phase space dimensionality) scaling polynomially in the system size. This method is based on generalization of the truncated Wigner approximation (TWA) to a higher dimensional phase space, where phase space variables are associated with a complete set of quantum operators spanning finite size clusters. The method becomes asymptotically exact with increasing cluster size. The crucial feature of TWA is fluctuating initial conditions, which we approximate by a Gaussian distribution. We show that such fluctuations dramatically increase accuracy of TWA over traditional cluster mean-field approximations. In this way we can treat on equal footing quantum and thermal fluctuations as well as compute entanglement and various equal and non-equal time correlation functions. The main limitation of the method is exponential scaling of the phase space dimensionality with the cluster size, which can be significantly reduced by using the language of Schwinger bosons and can likely be further reduced by truncating the local Hilbert space variables. We demonstrate the power of this method analyzing dynamics in various spin chains with and without disorder and show that we can capture such phenomena as long time hydrodynamic relaxation, many-body localization and the ballistic spread of entanglement.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1995

A FFT Method for the Quasiclassical Selection of Initial Ro-Vibrational States of Triatomic Molecules

This paper describes the use of an exact fast Fourier transform (FFT) method to prepare specified vibrational-rotational states of triatomic molecules. The method determines the Fourier coefficients needed to describe the coordinates and momenta of a vibrating-rotating triatomic molecule. Once the Fourier coefficients of a particular state are determined, it is possible to easily generate as many random sets of initial cartesian coordinates and momenta as desired. All the members of each set will correspond to the particular vibrational-rotational state selected. For example, in the case of the ground vibrational state of a non-rotating water molecule, the calculated actions of 100 sets of initial conditions produced actions within 0.001 h(bar) of the specified quantization values and energies within 5 cm(sup -1) of the semiclassical eigenvalue. The numerical procedure is straightforward for states in which all the fundamental frequencies are independent. However for states for which the fundamental frequencies become commensurate (resonance states), there are additional complications. In these cases it is necessary to determine a new set of "fundamental" frequencies and to modify the quantization conditions. Once these adjustments are made, good results are obtained for resonance states. The major problems are in labelling the large number of Fourier coefficients and the presence of regions of chaotic motion. Results are presented for the vibrational states of H2O and HCN and the ro-vibrational states of H2O.

Eaker, Charles W.↗
NASA NTRS1994

Computed potential energy surfaces for chemical reactions

Quantum mechanical methods have been used to compute potential energy surfaces for chemical reactions. The reactions studied were among those believed to be important to the NASP and HSR programs and included the recombination of two H atoms with several different third bodies; the reactions in the thermal Zeldovich mechanism; the reactions of H atom with O2, N2, and NO; reactions involved in the thermal De-NO(x) process; and the reaction of CH(squared Pi) with N2 (leading to 'prompt NO'). These potential energy surfaces have been used to compute reaction rate constants and rates of unimolecular decomposition. An additional application was the calculation of transport properties of gases using a semiclassical approximation (and in the case of interactions involving hydrogen inclusion of quantum mechanical effects).

Walch, Stephen P.↗
DOE OSTI2018

From Hartree Dynamics to the Relativistic Vlasov Equation

We derive the relativistic Vlasov equation from quantum Hartree dynamics for fermions with relativistic dispersion in the mean-field scaling, which is naturally linked with an effective semiclassic limit. Similar results in the non-relativistic setting have been recently obtained in Benedikter et al. (Arch Rat Mech Anal 221(1): 273–334, 2016). The new challenge that we have to face here, in the relativistic setting, consists in controlling the difference between the quantum kinetic energy and the relativistic transport term appearing in the Vlasov equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2018

Kinetic theory for strongly coupled Coulomb systems

The calculation of dynamical properties for matter under extreme conditions is a challenging task. The popular Kubo-Greenwood model exploits elements from equilibrium density-functional theory (DFT) that allow a detailed treatment of electron correlations, but its origin is largely phenomenological; traditional kinetic theories have a more secure foundation but are limited to weak ion-electron interactions. The objective here is to show how a combination of the two evolves naturally from the short-time limit for the generator of the effective single-electron dynamics governing time correlation functions without such limitations. This provides a theoretical context for the current DFT-related approach, the Kubo-Greenwood model, while showing the nature of its corrections. The method is to calculate the short-time dynamics in the single-electron subspace for a given configuration of the ions. This differs from the usual kinetic theory approach in which an average over the ions is performed as well. In this way the effective ion-electron interaction includes strong Coulomb coupling and is shown to be determined from DFT. The correlation functions have the form of the random-phase approximation for an inhomogeneous system but with renormalized ion-electron and electron-electron potentials. The dynamic structure function, density response function, and electrical conductivity are calculated as examples. The static local field corrections in the dielectric function are identified in this way. The current analysis is limited to semiclassical electrons (quantum statistical potentials), so important quantum conditions are excluded. However, a quantization of the kinetic theory is identified for broader application while awaiting its detailed derivation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Off-shell Partition Functions in 3d Gravity

We explore three-dimensional gravity with negative cosmological constant via canonical quantization. We focus on chiral gravity which is related to a single copy of PSL(2,R) Chern-Simons theory and is simpler to treat in canonical quantization. Its phase space for an initial value surface Σ is given by the appropriate moduli space of Riemann surfaces. We use geometric quantization to compute partition functions of chiral gravity on three-manifolds of the form Σ×S 1 , where Σ can have asymptotic boundaries. Most of these topologies do not admit a classical solution and are thus not amenable to a direct semiclassical path integral computation. We use an index theorem that expresses the partition function as an integral of characteristic classes over phase space. In the presence of n asymptotic boundaries, we use techniques from equivariant cohomology to localize the integral to a finite-dimensional integral over $\overline{M}$ g,n , which we evaluate in low genus cases. Higher genus partition functions quickly become complicated since they depend in an oscillatory way on Newton's constant. There is a precise sense in which one can isolate the non-oscillatory part which we call the fake partition function. We establish that there is a topological recursion that computes the fake partition functions for arbitrary Riemann surfaces Σ. As a result, there is a scaling limit in which the model reduces to JT gravity and our methods give a novel way to compute JT partition functions via equivariant localization.

Classical and Quantum Gravity↗
DOE OSTI2013

Computer simulation of liquid-vapor coexistence of confined quantum fluids

The liquid-vapor coexistence (LV) of bulk and confined quantum fluids has been studied by Monte Carlo computer simulation for particles interacting via a semiclassical effective pair potential V{sub eff}(r) = V{sub LJ} + V{sub Q}, where V{sub LJ} is the Lennard-Jones 12-6 potential (LJ) and V{sub Q} is the first-order Wigner-Kirkwood (WK-1) quantum potential, that depends on β = 1/kT and de Boer's quantumness parameter Λ=h/σ√(mε), where k and h are the Boltzmann's and Planck's constants, respectively, m is the particle's mass, T is the temperature of the system, and σ and ε are the LJ potential parameters. The non-conformal properties of the system of particles interacting via the effective pair potential V{sub eff}(r) are due to Λ, since the LV phase diagram is modified by varying Λ. We found that the WK-1 system gives an accurate description of the LV coexistence for bulk phases of several quantum fluids, obtained by the Gibbs Ensemble Monte Carlo method (GEMC). Confinement effects were introduced using the Canonical Ensemble (NVT) to simulate quantum fluids contained within parallel hard walls separated by a distance L{sub p}, within the range 2σ ⩽ L{sub p} ⩽ 6σ. The critical temperature of the system is reduced by decreasing L{sub p} and increasing Λ, and the liquid-vapor transition is not longer observed for L{sub p}/σ < 2, in contrast to what has been observed for the classical system.

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

Fourfold Clusters of Rovibrational Energies in H2Te Studied With an Ab Initio Potential Energy Function

We report an ab initio investigation of the cluster effect (i.e., the formation of nearly degenerate, four member groups of rotation-vibration energy levels at higher J and K(sub a). values) in the H2Te molecule. The potential energy function has been calculated ab initio at a total of 334 molecular geometries by means of the CCSD(T) method where the (1s-4f) core electrons of Te were described by an effective core potential. The values of the potential energy function obtained cover the region up to around 10,000/cm above the equilibrium energy. On the basis of the ab initio potential, the rotation-vibration energy spectra of H2Te-130 and its deuterated isotopomers have been calculated with the MORBID (Morse Oscillator Rigid Bender Internal Dynamics) Hamiltonian and computer program. In particular, we have calculated the rotational energy manifolds for J less than or = 40 in the vibrational ground state, the upsilon(sub 2) state, the "first triad" (the upsilon(sub l)/upsilon(sub 3)/2upsilon(sub 2) interacting vibrational states), and the "second triad" (the upsilon(sub 1) + upsilon(sub 2/upsilon(sub 2) + upsilon(sub 3)/3upsilon(sub 2) states) of H2Te-130. We find that the cluster formation in H2Te is very similar to those of of H2Se and H2S, which we have studied previously. However, contrary to semiclassical predictions, we do not determine any significant displacement of the clusters towards lower J values relative to H2Se. Hence the experimental observation of the cluster states in H2Te will be at least as difficult as in H2Se.

Jensen, Per↗
DOE OSTI2022

Magnetism and magnetotransport in the kagome antiferromagnet Mn 3 Ge

We perform classical Monte Carlo and stochastic Landau-Lifshitz-Gilbert simulations to study the temperature-dependent magnetism of the kagome antiferromagnet Weyl metal Mn3Ge, and we find that a long-range chiral order sets in at a transition temperature well below the Neel temperature (T N ). Based on the crystalline symmetries imposed by the chiral magnetic order, we argue for the presence of multiple isoenergetic Weyl nodes (nodes that are at the same energy and with a congruent Fermi surface around them) near the chemical potential. Using the semiclassical Boltzmann equations, we show that the combined contribution to the net longitudinal magnetoconductance (LMC) and the planar Hall conductance (PHC) from tilted Weyl nodes can lead to signatures that are qualitatively distinct from that of a single pair of Weyl nodes. In particular, we show that magnetic orders with different chiralities can give rise to different periods in LMC and PHC as a function of the in-plane magnetic field direction. As a result, this is ultimately related to differences in the symmetry-imposed constraints on the Weyl nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Thermal conductivity of III-V semiconductor superlattices

This paper presents a semiclassical model for the anisotropic thermal transport in III-V semiconductor superlattices (SLs). An effective interface rms roughness is the only adjustable parameter. Thermal transport inside a layer is described by the Boltzmann transport equation in the relaxation time approximation and is affected by the relevant scattering mechanisms (three-phonon, mass-difference, and dopant and electron scattering of phonons), as well as by diffuse scattering from the interfaces captured via an effective interface scattering rate. The in-plane thermal conductivity is obtained from the layer conductivities connected in parallel. The cross-plane thermal conductivity is calculated from the layer thermal conductivities in series with one another and with thermal boundary resistances (TBRs) associated with each interface; the TBRs dominate cross-plane transport. The TBR of each interface is calculated from the transmission coefficient obtained by interpolating between the acoustic mismatch model (AMM) and the diffuse mismatch model (DMM), where the weight of the AMM transmission coefficient is the same wavelength-dependent specularity parameter related to the effective interface rms roughness that is commonly used to describe diffuse interface scattering. The model is applied to multiple III-arsenide superlattices, and the results are in very good agreement with experimental findings. The method is both simple and accurate, easy to implement, and applicable to complicated SL systems, such as the active regions of quantum cascade lasers. It is also valid for other SL material systems with high-quality interfaces and predominantly incoherent phonon transport.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2019

Universal collective modes in two-dimensional chiral superfluids

In this work, we utilize semiclassical kinetic equations to investigate the order parameter collective modes of a class of two-dimensional superfluids. Extending the known results for p-wave superfluids, we show that for any chiral ground state of angular momentum L ≥ 1, there exists a subgap mode with a mass √2Δ in the BCS limit, where Δ is the magnitude of the ground-state gap. We determine the most significant Landau parameter that contributes to the mass renormalization and show explicitly that the renormalized modes become massless at the Pomeranchuk instability of the fermion vacuum. Particularly for L = 1, we propose a continuous field theory to include the Fermi liquid effect in a quadrupolar channel and produce the same result under consistent approximations. Furthermore, they provide potential diagnostics for distinguishing two-dimensional chiral ground states of different angular momenta with the order parameter collective modes and reveal another low-energy degree of freedom near the nematic transition.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2014

Improving long time behavior of Poisson bracket mapping equation: A non-Hamiltonian approach

Understanding nonadiabatic dynamics in complex systems is a challenging subject. A series of semiclassical approaches have been proposed to tackle the problem in various settings. The Poisson bracket mapping equation (PBME) utilizes a partial Wigner transform and a mapping representation for its formulation, and has been developed to describe nonadiabatic processes in an efficient manner. Operationally, it is expressed as a set of Hamilton's equations of motion, similar to more conventional classical molecular dynamics. However, this original Hamiltonian PBME sometimes suffers from a large deviation in accuracy especially in the long time limit. Here, we propose a non-Hamiltonian variant of PBME to improve its behavior especially in that limit. As a benchmark, we simulate spin-boson and photosynthetic model systems and find that it consistently outperforms the original PBME and its Ehrenfest style variant. We explain the source of this improvement by decomposing the components of the mapping Hamiltonian and by assessing the energy flow between the system and the bath. We discuss strengths and weaknesses of our scheme with a viewpoint of offering future prospects.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Correlating isothermal compressibility to nucleon fluctuations in the inner crust of neutron stars

The question of how and which physical observables or thermodynamic parameters can best predict the onset of a possible phase transition in the inner crust of neutron stars remains largely unresolved. Here, using semiclassical Monte Carlo simulations, we investigate the isothermal compressibility and density fluctuations in a region of relevance to the dynamics of the inner crust. We show that the isothermal compressibility serves as a robust observable to characterize the transition from the nonuniform crust to the uniform core for proton fractions over 0.2. Moreover, we show explicitly how the two-component isothermal compressibility, computed using the Kirkwood-Buff theory, is directly connected to the fluctuations in the number density, recorded in the grand canonical ensemble by monitoring the number of particles in a small volume located at the center of the simulation box. That is, we compute mean-square particle fluctuations and compare them against the isothermal compressibility for different proton fractions. Although our results show that the mean-square particle fluctuations are proportional to the isothermal compressibility, the lack of a perfect correlation is attributed to the relatively small number of particles included in the simulations. The nonunity slope observed in the dimensionless isothermal compressibility—total nucleon fluctuation variance relationship suggests that the inner crust of neutron stars is composed of anisotropic and inhomogeneous matter.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2017

Comprehensive study of the dynamics of a classical Kitaev spin liquid

In this work, we study the spin-S Kitaev model in the classical (S→∞) limit using Monte Carlo simulations combined with semiclassical spin dynamics. We discuss differences and similarities in the dynamical structure factors of the spin-1/2 and the classical Kitaev liquids. The quantum behavior is restricted to low temperatures where a gap protects visons from decohering the system. Once this quantum gap is breached, at low temperatures compared to the coupling constant, significant entropic disorder decoheres the Majorana fermions and the system is described quantitatively by classical dynamics. The low-temperature and low-energy spectrum of the classical model exhibits a finite-energy peak, which is the precursor of the one produced by the Majorana modes of the S=1/2 model. The classical peak is spectrally narrowed compared to the quantum result and can be explained by magnon excitations within fluctuating one-dimensional manifolds (loops). Hence the difference from the classical limit to the quantum limit can be understood by the fractionalization of magnons propagating in one-dimensional manifolds. Moreover, we show that the momentum-space distribution of the low-energy spectral weight of the S=1/2 model follows the momentum-space distribution of zero modes of the classical model.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗