Search NASA⌕ Search

SEARCH · Search NASA

Results for “calculation method”

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 1,783 records · Page 99

Downward quantum learning from element 118: Automated generation of Fermi–Löwdin orbitals for all atoms

A new algorithm based on a rigorous theorem and quantum data computationally mined from element 118 guarantees automated construction of initial Fermi–Löwdin-Orbital (FLO) starting points for all elements in the Periodic Table. It defines a means for constructing a small library of scalable FLOs for universal use in molecular and solid-state calculations. The method can be systematically improved for greater efficiency and for applications to excited states such as x-ray excitations and optically silent excitations. FLOs were introduced to recast the Perdew–Zunger self-interaction correction (PZSIC) into an explicit unitarily invariant form. The FLOs are generated from a set of N quasi-classical electron positions, referred to as Fermi-Orbital descriptors (FODs), and a set of N-orthonormal single-electron orbitals. FOD positions, when optimized, minimize the PZSIC total energy. However, creating sets of starting FODs that lead to a positive definite Fermi orbital overlap matrix has proven to be challenging for systems composed of open-shell atoms and ions. The proof herein guarantees the existence of a FLOSIC solution and further guarantees that if a solution for N electrons is found, it can be used to generate a minimum of N – 1 and a maximum of 2 N – 2 initial starting points for systems composed of a smaller number of electrons. As a result, applications to heavy and super-heavy atoms are presented. All starting solutions reported here were obtained from a solution for element 118, Oganesson.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

A comparison of QTP functionals against coupled-cluster methods for EAs of small organic molecules

EA-EOM-CCSD electron affinities and LUMO energies of various Kohn–Sham density functional theory (DFT) methods are calculated for an a priori IP benchmark set of 64 small, closed-shell molecules. The purpose of these calculations was to investigate whether the QTP KS-DFT functionals can emulate EA-EOM-CC with only a mean-field approximation. We show that the accuracy of DFT—relative to CCSD—improves significantly when elements of correlated orbital theory are introduced into the parameterization to define the QTP family of functionals. In particular, QTP(02), which has only a single range separation parameter, provides results accurate to a MAD of <0.15 eV for the whole set of 64 molecules compared to EA-EOM-CCSD, far exceeding the results from the non-QTP family of density functionals.

Chemistry↗

Quantum Drude oscillators coupled with Coulomb potential as an efficient model for bonded and non-covalent interactions in atomic dimers

The quantum Drude oscillator (QDO) model has been widely used as an efficient surrogate to describe the electric response properties of matter as well as long-range interactions in molecules and materials. Most commonly, QDOs are coupled within the dipole approximation so that the Hamiltonian can be exactly diagonalized, which forms the basis for the many-body dispersion method [Phys. Rev. Lett. 108, 236402 (2012)]. The dipole coupling is efficient and allows us to study non-covalent many-body effects in systems with thousands of atoms. However, there are two limitations: (i) the need to regularize the interaction at short distances with empirical damping functions and (ii) the lack of multipolar effects in the coupling potential. In this work, we convincingly address both limitations of the dipole-coupled QDO model by presenting a numerically exact solution of the Coulomb-coupled QDO model by means of quantum Monte Carlo methods. We calculate the potential-energy surfaces of homogeneous QDO dimers, analyzing their properties as a function of the three tunable parameters: frequency, reduced mass, and charge. We study the coupled-QDO model behavior at short distances and show how to parameterize this model to enable an effective description of chemical bonds, such as the covalent bond in the H2 molecule.

Chemistry↗

Calculation and experimental verification of a collimator with a Kerr lens for fibre laser mode locking

We present the calculation and experimental verification of the operation of a collimator coupling a light beam into an optical fibre with an As{sub 2}S{sub 3} nonlinear plate. It is shown that a Kerr lens formed in the plate with an increase in the peak pulse power leads to an increase in the output coupling ratio by about 15%, which can be used for the passive mode locking of a fibre laser. (lasers)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Acousto-optic modulators for a controlled frequency shift of light beams in optical and microwave cold-atom frequency standards

A series of specialised acousto-optic modulators/frequency shifters based on paratellurite (TeO{sub 2}) single crystals is developed for systems of laser cooling and frequency stabilisation of laser sources in optical and microwave frequency standards based on cold strontium, rubidium and caesium atoms. The methods of calculation and optimisation of the main parameters of acousto-optic modulators are considered; their basic technical characteristics are experimentally investigated in the spectral range of optical radiation 460 – 900 nm and in the frequency shift range of 5 – 400 MHz. The modulators have a high diffraction efficiency, low power consumption and the ability to work with linearly and randomly polarised optical radiation. (frequency standards)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Verification and validation of detonation-shock-dynamics relations for explosives described by general equation of state and chemical reaction models

Detonation shock dynamics is a powerful method to model the behaviour of High Explosives (HE). However in order to use this method, the underlying relationship between the local radius of curvature and the detonation speed must be known. Previous work has developed methods to calculate this effect using simple, single-step Arrhenius and polytropic gas, models for the chemical reaction and the equation of state, respectively. In recent years, more complex models for both reaction rates and equations of state have been developed which show better agreement with experimental data than these simple models, especially when considering condensed phase explosives.. This work presents the governing equations for solving these problems in a way that is generalised to use arbitrary equations of state as well as reaction models which may have more than a single step and multiple product species. This implementation is verified against exact solutions, demonstrating that the equations were implemented properly. The verified algorithm is then validated against experimental data and high fidelity simulations, showing that it is able to make accurate predictions in a regime where the underlying assumptions of the governing equations are valid. Importantly, this approach has many applications: from creating equivalent detonation shock dynamics models for existing reactive burn calibrations for HE; to developing new functional forms and calibrations of reactive burn models for condensed phase high explosives.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electromagnetic Analysis of ITER Electron Cyclotron Emission Components and Support Structures for the Preliminary Design Review

The ITER electron cyclotron emission (ECE) diagnostic system is located at Design Shield Module (DSM) 2, Equatorial Port (EP) 9 to measure electron temperature profile and electron temperature fluctuations and also to assess nonthermal electron distributions via the oblique view. Therefore, ECE has both radial and oblique views with two couples of mirrors for two different optical views. This ECE diagnostic system shall be exposed to significant power due to unabsorbed electron cyclotron heating power in the plasma. It shall also receive large electromagnetic (EM) loads up to 100 MN/m 3 force density due to the eddy currents generated by the short 16-ms transient plasma disruption. The global EM models of EP12 and EP11 with the worst case of plasma disruption Major Disruption Downward Exponential 16 ms Category III (MD_DW_EXP16MS_CATIII) have been provided by ITER. However, the support structure for these critical components, including four mirrors, is still under the design and development phase, especially during the period of preliminary design reviews (PDRs). For those small in-vessel components that are not modeled explicitly for the PDRs, we can extract the magnetic field (B) data and flux variations (dB/dt) at those locations from the global model and then use the B and dB/dt method to calculate the eddy currents, EM forces, and torques on those components in the local submodels. In conclusion, this paper will present the detailed PDR EM analysis results of the ECE components and support structures.

Electromagnetic↗

Unimodular gravity and the lepton anomalous magnetic moment at one-loop

We work out the one-loop contribution to the lepton anomalous magnetic moment coming from Unimodular Gravity. We use Dimensional Regularization and Dimensional Reduction to carry out the computations. In either case, we find that Unimodular Gravity gives rise to the same one-loop correction as that of General Relativity.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The topology of non-resonant stellarator divertors

We apply topological methods to better understand how the magnetic field in the stellarator edge can be diverted away from the confined region. Our primary method is calculating the winding numbers of closed contours, which gives information on the number and nature of fixed points within a bounded region. We first apply this to the non-resonant divertor Hamiltonian system, and present a simple explanation for the system’s diversion: trajectories are guided away from the confined region by X-points which are ‘unpaired’ i.e. do not have corresponding O-points and therefore do not resemble an island chain. We show how similar phenomena can occur in a similar, axisymmetric Hamiltonian system. Secondly, we find examples of neoclassically optimised stellarators in the quasi-symmetric stellarator repository database which divert the magnetic field via unpaired X-points. We present and discuss three examples, each containing novel phenomena which might be desirable for stellarator divertors. These examples serve as an illustration of new divertor possibilities which exist in realistic stellarators, which may ultimately have application for future experiments and reactors.

X-points, topology↗

Mode coupling and multiquantum vibrational excitations in Feshbach-resonant positron annihilation in molecules

The dominant mechanism of low-energy positron annihilation in polyatomic molecules is through positron capture in vibrational Feshbach resonances (VFR). In this paper, we investigate theoretically the effect of anharmonic terms in the vibrational Hamiltonian on positron annihilation rates. Such interactions enable positron capture in VFRs associated with multiquantum vibrational excitations, leading to enhanced annihilation. Mode coupling can also lead to faster depopulation of VFRs, thereby reducing their contribution to the annihilation rates. To analyze this complex picture, we use coupled-cluster methods to calculate the anharmonic vibrational spectra and dipole transition amplitudes for chloroform, chloroform-d 1 , 1,1-dichloroethylene, and methanol, and use these data to compute positron resonant annihilation rates for these molecules. Theoretical predictions are compared with the annihilation rates measured as a function of incident positron energy. Overall, the results demonstrate the importance of mode coupling in both enhancement and suppression of the VFR. There is also experimental evidence for the direct excitation of multimode VFR. Their contribution is analyzed using a statistical approach, with an outlook towards more accurate treatment of this phenomenon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Perturbative representation of ultrashort nonparaxial elegant Laguerre-Gaussian fields

An analytical method for calculating the electromagnetic fields of a nonparaxial elegant Laguerre-Gaussian (LG) vortex beam is presented for arbitrary pulse duration, spot size, and LG mode. This perturbative approach provides a numerically tractable model for the calculation of arbitrarily high radial and azimuthal LG modes in the nonparaxial regime, without requiring integral representations of the fields. A key feature of this perturbative model is its use of a Poisson-like frequency spectrum, which allows for the proper description of pulses of arbitrarily short duration. Lastly, this model is thus appropriate for simulating laser-matter interactions, including those involving short laser pulses.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Variational Self-Consistent Theory for Beam-Loaded Cavities

A new variational theory is presented for beam loading in microwave cavities. The beam{field interaction is formulated as a dynamical interaction whose stationarity according to Hamilton’s principle will naturally lead to steady-state solutions that indicate how a cavity’s resonant frequency, Q and optimal coupling coefficient will detune as a result of the beam loading. A driven cavity Lagrangian is derived from first principles, including the effects of cavity wall losses, input power and beam interaction. The general formulation is applied to a typical klystron input cavity to predict the appropriate detuning parameters required to maximize the gain (or modulation depth) in the average Lorentz factor boost, (Δ γ). Numerical examples are presented, showing agreement with the general detuning trends previously observed in the literature. The developed formulation carries several advantages for beam-loaded cavity structures. It provides a self-consistent model for the dynamical (nonlinear) beam{field interaction, a procedure for maximizing gain under beam-loading conditions, and a useful set of parameters to guide cavity-shape optimization during the design of beam-loaded systems. Enhanced clarity of the physical picture underlying the problem seems to be gained using this approach, allowing straightforward inclusion or exclusion of different field configurations in the calculation and expressing the final results in terms of measurable quantities. Two field configurations are discussed for the klystron input cavity, using finite magnetic confinement or no confinement at all. Formulating the problem in a language that is directly accessible to the powerful techniques found in Hamiltonian dynamics and canonical transformations may potentially carry an additional advantage in terms of analytical computational gains, under suitable conditions

43 PARTICLE ACCELERATORS↗

Thermal radiation in systems of many dipoles

Systems of many nanoparticles or volume-discretized bodies exhibit collective radiative properties that could be used for enhanced, guided, or tunable thermal radiation. These are commonly treated as assemblies of point dipoles with interactions described by Maxwell's equations and thermal fluctuations correlated by the fluctuation-dissipation theorem. Here, we demonstrate the equivalence of different theories for these systems and provide a complete derivation of many-dipole thermal radiation, showing that the correct use of the fluctuation-dissipation theorem depends on the definitions of fluctuating and induced dipole moments. We formulate a method to calculate the diffusive radiative thermal conductivity of arbitrary collections of nanoparticles; this allows the comparison of thermal radiation to other heat-transfer modes and across different material systems. We calculate the radiative thermal conductivity of ordered and disordered arrays of SiC and SiO 2 nanoparticles and show that thermal radiation can significantly contribute to thermal transport in these systems, because packed nanoparticles have low phonon thermal conductivity. Finally, we demonstrate that the radiative heat transfer strongly depends on the materials and geometrical arrangement of the nanoparticle array, and we verify our calculations by comparison to the exact solution for a one-dimensional particle chain.

42 ENGINEERING↗

Electronic band gaps from quantum Monte Carlo methods

We develop a method for calculating the fundamental electronic gap of semiconductors and insulators using grand canonical quantum Monte Carlo simulations. Further, we discuss the origin of the bias introduced by supercell calculations of finite size and show how to correct the leading and subleading finite size errors either based on observables accessible in the finite-sized simulations or from density-functional theory calculations. Our procedure is applied to carbon, silicon, and molecular hydrogen crystals, and compared to experiment for carbon and silicon.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Landau level of fragile topology

Here, we study the Hofstadter butterfly and Landau levels of the twisted bilayer graphene (TBG). We show that the nontrivial fragile topology of the lowest two bands near the charge neutral point makes their Hofstadter butterfly generically connected with higher bands, closing the gap between the first and second conduction (valence) bands at a certain magnetic flux per unit cell. We also develop a momentum space method for calculating the TBG Hofstadter butterfly, from which we identify three phases where the Hofstadter butterflies of the lowest two bands and the higher bands are connected in different ways. We show this leads to a crossing between the $\textit{ν}$ = 4 Landau fan from the charge neutral point and the zero-field band gap at one flux per moiré unit cell, which corresponds to a magnetic field 25$θ^2$ T (twist angle $\textit{θ}$ in degrees). This provides an experimentally testable feature of the fragile topology. In general, we expect it to be a generic feature that the Hofstadter butterfly of the topological bands are connected with the Hofstadter spectra of other bands. We further show the TBG band theory with Zeeman splitting being the most sizable splitting could result in Landau fans at the charge neutral point and half fillings near the magic angle, and we predict their variations under an in-plane magnetic field.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonlocal and nonadiabatic Pauli potential for time-dependent orbital-free density functional theory

Time-dependent orbital-free density functional theory is an efficient ab initio method for calculating the electronic dynamics of large systems. In comparison to standard time-dependent density functional theory, it computes only a single electronic state regardless of system size, but it requires an additional time-dependent Pauli potential term. Herein we propose a nonadiabatic and nonlocal Pauli potential whose main ingredients are the time-dependent particle and current densities. Our calculations of the optical spectra of metallic and semiconductor clusters indicate that nonlocal and nonadiabatic time-dependent orbital-free density functional theory performs accurately for metallic systems and semiquantitatively for semiconductors. This paper opens the door to wide applicability of time-dependent orbital-free density functional theory for nonequilibrium electron and electron-nuclear dynamics of complex materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Efficient time-dependent orbital-free density functional theory: Semilocal adiabatic response

Orbital-free density functional theory and its time-dependent extension are efficient ab initio methods for calculating the electronic structure and dynamics of large systems. Through the calculation of the optical spectra of selected clusters, we reach three important conclusions: (1) The quality of the spectra is strongly affected by the quality of the corresponding ground-state electron density; (2) the adiabatic part of the electronic response to external perturbations can be safely evaluated at the semilocal level; and (3) the nonadiabatic, current-dependent part of the time-dependent Pauli potential is key to recover correct spectral envelopes.

36 MATERIALS SCIENCE↗