Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Wave Function”

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,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2014

Simulation of ultrasonic surface waves with multi-Gaussian and point source beam models

In the past decade, multi-Gaussian beam models have been developed to solve many complicated bulk wave propagation problems. However, to date those models have not been extended to simulate the generation of Rayleigh waves. Here we will combine Gaussian beams with an explicit high frequency expression for the Rayleigh wave Green function to produce a three-dimensional multi-Gaussian beam model for the fields radiated from an angle beam transducer mounted on a solid wedge. Simulation results obtained with this model are compared to those of a point source model. It is shown that the multi-Gaussian surface wave beam model agrees well with the point source model while being computationally much more efficient.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Large perpendicular magnetic anisotropy of single Co atom on MgO monolayer: A first-principles study

Realizing the magnetic bit with a single atom is the ultimate goal for magnetic storage. Based on density functional theory, the magnetic anisotropy (MA) of single Co atom on MgO monolayer has been investigated. Results show that this two dimensional system possesses a large perpendicular MA, about 5.8 meV per Co atom. Besides, there exists remarkable unquenched orbital moments for different magnetization directions, which can be attributed to the reduction of coordination number in two dimensional system and is responsible for the enhanced MA. The Bloch pseudo-wavefunction and band structure of Co d-orbitals have been calculated to elucidate the origin of the perpendicular MA.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2015

Dust heating by Alfvén waves using non-Maxwellian distribution function

Quasilinear theory is employed in order to evaluate the resonant heating rate by Alfvén waves, of multiple species dust particles in a hot, collisionless, and magnetized plasma, with the underlying assumption that the dust velocity distribution function can be modeled by a generalized (r, q) distribution function. The kinetic linear dispersion relation for the electromagnetic dust cyclotron Alfvén waves is derived, and the dependence of the heating rate on the magnetic field, mass, and density of the dust species is subsequently investigated. The heating rate and its dependence on the spectral indices r and q of the distribution function are also investigated. It is found that the heating is sensitive to negative value of spectral index r.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2016

On the existence of the optimal order for wavefunction extrapolation in Born-Oppenheimer molecular dynamics

Wavefunction extrapolation greatly reduces the number of self-consistent field (SCF) iterations and thus the overall computational cost of Born-Oppenheimer molecular dynamics (BOMD) that is based on the Kohn–Sham density functional theory. Going against the intuition that the higher order of extrapolation possesses a better accuracy, we demonstrate, from both theoretical and numerical perspectives, that the extrapolation accuracy firstly increases and then decreases with respect to the order, and an optimal extrapolation order in terms of minimal number of SCF iterations always exists. We also prove that the optimal order tends to be larger when using larger MD time steps or more strict SCF convergence criteria. By example BOMD simulations of a solid copper system, we show that the optimal extrapolation order covers a broad range when varying the MD time step or the SCF convergence criterion. Therefore, we suggest the necessity for BOMD simulation packages to open the user interface and to provide more choices on the extrapolation order. Another factor that may influence the extrapolation accuracy is the alignment scheme that eliminates the discontinuity in the wavefunctions with respect to the atomic or cell variables. We prove the equivalence between the two existing schemes, thus the implementation of either of them does not lead to essential difference in the extrapolation accuracy.

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

Towards quantifying the role of exact exchange in the prediction hydrogen bond spin-spin coupling constants involving fluorine

The effect of a fraction of Hartree-Fock exchange on the calculated spin-spin coupling constants involving fluorine through a hydrogen bond is analyzed in detail. Coupling constants calculated using wavefunction methods are revisited in order to get high-level calculations using the same basis set. Accurate MCSCF results are obtained using an additive approach. These constants and their contributions are used as a reference for density functional calculations. Within the density functional theory, the Hartree-Fock exchange functional is split in short- and long-range using a modified version of the Coulomb-attenuating method with the SLYP functional as well as with the original B3LYP. Results support the difficulties for calculating hydrogen bond coupling constants using density functional methods when fluorine nuclei are involved. Coupling constants are very sensitive to the Hartree-Fock exchange and it seems that, contrary to other properties, it is important to include this exchange for short-range interactions. Best functionals are tested in two different groups of complexes: those related with anionic clusters of type [F(HF){sub n}]{sup −} and those formed by difluoroacetylene and either one or two hydrogen fluoride molecules.

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

Geometry of quantum Hall states: Gravitational anomaly and transport coefficients

We show that universal transport coefficients of the fractional quantum Hall effect (FQHE) can be understood as a response to variations of spatial geometry. Some transport properties are essentially governed by the gravitational anomaly. We develop a general method to compute correlation functions of FQH states in a curved space, where local transformation properties of these states are examined through local geometric variations. We introduce the notion of a generating functional and relate it to geometric invariant functionals recently studied in geometry. We develop two complementary methods to study the geometry of the FQHE. One method is based on iterating a Ward identity, while the other is based on a field theoretical formulation of the FQHE through a path integral formalism.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2014

Electronic correlation without double counting via a combination of spin projected Hartree-Fock and density functional theories

Several schemes to avoid the double counting of correlations in methods that merge multireference wavefunctions with density functional theory (DFT) are studied and here adapted to a combination of spin-projected Hartree-Fock (SUHF) and DFT. The advantages and limitations of the new method, denoted SUHF+f{sub c}DFT, are explored through calculations on benchmark sets in which the accounting of correlations is challenging for pure SUHF or DFT. It is shown that SUHF+f{sub c}DFT can greatly improve the description of certain molecular properties (e.g., singlet-triplet energy gaps) which are not improved by simple addition of DFT dynamical correlation to SUHF. However, SUHF+f{sub c}DFT is also shown to have difficulties dissociating certain types of bonds and describing highly charged ions with static correlation. Possible improvements to the current SUHF+f{sub c}DFT scheme are discussed in light of these results.

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

Cost-effective description of strong correlation: Efficient implementations of the perfect quadruples and perfect hextuples models

Novel implementations based on dense tensor storage are presented here for the singlet-reference perfect quadruples (PQ) [J. A. Parkhill et al., J. Chem. Phys. 130, 084101 (2009)] and perfect hextuples (PH) [J. A. Parkhill and M. Head-Gordon, J. Chem. Phys. 133, 024103 (2010)] models. The methods are obtained as block decompositions of conventional coupled-cluster theory that are exact for four electrons in four orbitals (PQ) and six electrons in six orbitals (PH), but that can also be applied to much larger systems. PQ and PH have storage requirements that scale as the square, and as the cube of the number of active electrons, respectively, and exhibit quartic scaling of the computational effort for large systems. Applications of the new implementations are presented for full-valence calculations on linear polyenes (C n H n+2 ), which highlight the excellent computational scaling of the present implementations that can routinely handle active spaces of hundreds of electrons. The accuracy of the models is studied in the π space of the polyenes, in hydrogen chains (H 50 ), and in the π space of polyacene molecules. In all cases, the results compare favorably to density matrix renormalization group values. With the novel implementation of PQ, active spaces of 140 electrons in 140 orbitals can be solved in a matter of minutes on a single core workstation, and the relatively low polynomial scaling means that very large systems are also accessible using parallel computing.

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

An evaluation for geometries, formation enthalpies, and dissociation energies of diatomic and triatomic (C, H, N, O), NO 3 , and HNO 3 molecules from the PAW DFT method with PBE and optB88-vdW functionals

The structural geometries, formation enthalpies, and dissociation energies of all diatomic and triatomic molecules consisting of the four basic elements C, H, N, and/or O are calculated using the projector augmented wave density functional theory (DFT) method with the Perdew–Burke–Ernzerhof and optB88-vdW exchange-correlation functionals. The calculations are also extended to two larger molecules NO 3 and HNO 3 , which consist of four and five atoms, respectively. In total, 82 molecules or isomers are considered in the calculations. The geometric parameters including 42 bond lengths and 15 bond angles of these molecules calculated using the planewave DFT method are highly satisfactory, relative to the available experimental data. The error analysis is also performed for 49 formation enthalpies and 138 dissociation energies (including 51 atomization energies as well as the corresponding bond dissociation energies). The results are also compared with the previous data from various atomic-orbital-based methods for molecules and from similar or different planewave DFT methods for various solids and other molecules. This provides an informative and instructive evaluation especially for calculating the large-size material systems containing these small molecules as well as for developing the DFT methods further.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2020

Wave-by-Wave Forecasting via Assimilation of Marine Radar Data

Abstract This work describes a phase-resolving wave-forecasting algorithm that is based on the assimilation of marine radar image time series. The algorithm is tested against synthetic data and field observations. The algorithm couples X-band marine radar observations with a phase-resolving wave model that uses the linear mild slope equations for reconstruction of water surface elevations over a large domain of O (km) and a prescribed time window of O (min). The reconstruction also enables wave-by-wave forecasting through forward propagation in space and time. Marine radar image time series provide the input wave observations through a previously given relationship between backscatter intensity and the radial component of the sea surface slope. The algorithm assimilates the wave slope information into the model via a best-fit wave source function at the boundary that minimizes the slope reconstruction error over an annular region at the outer ranges of the radar images. The wave model is then able to propagate the waves across a polar domain to a location of interest at nearer ranges. The constraints on the method for achieving real-time forecasting are identified, and the algorithm is verified against synthetic data and tested using field observations.

Engineering↗
DOE OSTI2023

Importance of Dispersion in the Molecular Geometries of Mn(III) Spin-Crossover Complexes

The computational investigation of the molecular geometries of a pair of manganese(III) spin-crossover complexes is reported. For the geometry of the quintet high-spin state, density functionals significantly overestimate Mn–Namine bond distances, although the geometry for the triplet intermediate-spin state is well described. Here, comparisons with several wave function-based methods demonstrate that this error is due to the limited ability of commonly used density functionals to recover dispersion beyond a certain extent. Among the methods employed for geometry optimization, restricted open-shell Møller–Plesset perturbation theory (MP2) appropriately describes the high-spin geometry but results in a slightly shorter Mn–O distance in both spin states. On the other hand, extended multistate complete active space second-order perturbation theory (XMS-CASPT2) provides a good description of the geometry for the intermediate-spin state but also sufficiently recovers dispersion, performing well for the high-spin state. Despite the fact that the electronic structure of both spin states is dominated by one-electron configuration, XMS-CASPT2 offers a balanced approach, leading to molecular geometries with much better agreement with experiment than MP2 and DFT. A scan along the Mn–N amine bond demonstrates that for these complexes coupled cluster methods (i.e., DLPNO-CCSD(T)) also yield bond distances in agreement with experiment while multiconfiguration pair density functional theory (MC-PDFT) is unable to recover dispersion well enough, analogous to single-reference DFT.

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

Multicomponent density functional theory embedding formulation

Multicomponent density functional theory (DFT) methods have been developed to treat two types of particles, such as electrons and nuclei, quantum mechanically at the same level. In the nuclear-electronic orbital (NEO) approach, all electrons and select nuclei, typically key protons, are treated quantum mechanically. For multicomponent DFT methods developed within the NEO framework, electron-proton correlation functionals based on explicitly correlated wavefunctions have been designed and used in conjunction with well-established electronic exchange-correlation functionals. Herein a general theory for multicomponent embedded DFT is developed to enable the accurate treatment of larger systems. In the general theory, the total electronic density is separated into two subsystem densities, denoted as regular and special, and different electron-proton correlation functionals are used for these two electronic densities. In the specific implementation, the special electron density is defined in terms of spatially localized Kohn-Sham electronic orbitals, and electron-proton correlation is included only for the special electron density. The electron-proton correlation functional depends on only the special electron density and the proton density, whereas the electronic exchange-correlation functional depends on the total electronic density. This scheme includes the essential electron-proton correlation, which is a relatively local effect, as well as the electronic exchange-correlation for the entire system. This multicomponent DFT-in-DFT embedding theory is applied to the HCN and FHF{sup −} molecules in conjunction with two different electron-proton correlation functionals and three different electronic exchange-correlation functionals. The results illustrate that this approach provides qualitatively accurate nuclear densities in a computationally tractable manner. The general theory is also easily extended to other types of partitioning schemes for multicomponent systems.

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

Single-Layer Plasmonic Metasurface Half-Wave Plates with Wavelength-Independent Polarization Conversion Angle

Manipulation of polarization state is of great fundamental importance and plays a crucial role in modern photonic applications such as optical communication, imaging, and sensing. Metamaterials and metasurfaces have attracted increasing interest in this area because they facilitate designer optical response through engineering the composite subwavelength structures. In this paper, we propose a general methods of designing half-wave plate and demonstrate in the near-infrared wavelength range an optically thin plasmonic metasurface half-wave plates that rotate the polarization direction of the linearly polarized incident light with a high degree of linear polarization. Finally, the half-wave plate functionality is realized through arranging the orientation of the nanoantennas to form an appropriate spatial distribution profile, which behave exactly as in classical half-wave plates but over in a wavelength-independent way.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2017

The Acceleration of Charged Particles at a Spherical Shock Moving through an Irregular Magnetic Field

We investigate the physics of charged-particle acceleration at spherical shocks moving into a uniform plasma containing a turbulent magnetic field with a uniform mean. This has applications to particle acceleration at astrophysical shocks, most notably, to supernovae blast waves. We numerically integrate the equations of motion of a large number of test protons moving under the influence of electric and magnetic fields determined from a kinematically defined plasma flow associated with a radially propagating blast wave. Distribution functions are determined from the positions and velocities of the protons. The unshocked plasma contains a magnetic field with a uniform mean and an irregular component having a Kolmogorov-like power spectrum. The field inside the blast wave is determined from Maxwell’s equations. The angle between the average magnetic field and unit normal to the shock varies with position along its surface. It is quasi-perpendicular to the unit normal near the sphere’s equator, and quasi-parallel to it near the poles. We find that the highest intensities of particles, accelerated by the shock, are at the poles of the blast wave. The particles “collect” at the poles as they approximately adhere to magnetic field lines that move poleward from their initial encounter with the shock at the equator, as the shock expands. The field lines at the poles have been connected to the shock the longest. We also find that the highest-energy protons are initially accelerated near the equator or near the quasi-perpendicular portion of the shock, where the acceleration is more rapid.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Calculation of room temperature conductivity and mobility in tin-based topological insulator nanoribbons

Monolayers of tin (stannanane) functionalized with halogens have been shown to be topological insulators. Using density functional theory (DFT), we study the electronic properties and room-temperature transport of nanoribbons of iodine-functionalized stannanane showing that the overlap integral between the wavefunctions associated to edge-states at opposite ends of the ribbons decreases with increasing width of the ribbons. Obtaining the phonon spectra and the deformation potentials also from DFT, we calculate the conductivity of the ribbons using the Kubo-Greenwood formalism and show that their mobility is limited by inter-edge phonon backscattering. We show that wide stannanane ribbons have a mobility exceeding 10{sup 6} cm{sup 2}/Vs. Contrary to ordinary semiconductors, two-dimensional topological insulators exhibit a high conductivity at low charge density, decreasing with increasing carrier density. Furthermore, the conductivity of iodine-functionalized stannanane ribbons can be modulated over a range of three orders of magnitude, thus rendering this material extremely interesting for classical computing applications.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Analytic derivative couplings in time-dependent density functional theory: Quadratic response theory versus pseudo-wavefunction approach

We revisit the formalism for analytic derivative couplings between excited states in time-dependent density functional theory (TDDFT). We derive and implement these couplings using quadratic response theory, then numerically compare this response-theory formulation to couplings implemented previously based on a pseudo-wavefunction formalism and direct differentiation of the Kohn-Sham determinant. Numerical results, including comparison to full configuration interaction calculations, suggest that the two approaches perform equally well for many molecular systems, provided that the underlying DFT method affords accurate potential energy surfaces. The response contributions are found to be important for certain systems with high symmetry, but can be calculated with only a moderate increase in computational cost beyond what is required for the pseudo-wavefunction approach. In the case of spin-flip TDDFT, we provide a formal proof that the derivative couplings obtained using response theory are identical to those obtained from the pseudo-wavefunction formulation, which validates our previous implementation based on the latter formalism.

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

Ensemble density functional theory method correctly describes bond dissociation, excited state electron transfer, and double excitations

State-averaged (SA) variants of the spin-restricted ensemble-referenced Kohn-Sham (REKS) method, SA-REKS and state-interaction (SI)-SA-REKS, implement ensemble density functional theory for variationally obtaining excitation energies of molecular systems. In this work, the currently existing version of the SA-REKS method, which included only one excited state into the ensemble averaging, is extended by adding more excited states to the averaged energy functional. A general strategy for extension of the REKS-type methods to larger ensembles of ground and excited states is outlined and implemented in extended versions of the SA-REKS and SI-SA-REKS methods. The newly developed methods are tested in the calculation of several excited states of ground-state multi-reference systems, such as dissociating hydrogen molecule, and excited states of donor–acceptor molecular systems. For hydrogen molecule, the new method correctly reproduces the distance dependence of the lowest excited state energies and describes an avoided crossing between the doubly excited and singly excited states. For bithiophene–perylenediimide stacked complex, the SI-SA-REKS method correctly describes crossing between the locally excited state and the charge transfer excited state and yields vertical excitation energies in good agreement with the ab initio wavefunction methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗