Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function approximation”

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,207 records · Page 67

Quantum interference in polyenes

The explicit form of the zeroth Green's function in the Hückel model, approximated by the negative of the inverse of the Hückel matrix, has direct quantum interference consequences for molecular conductance. We derive a set of rules for transmission between two electrodes attached to a polyene, when the molecule is extended by an even number of carbons at either end (transmission unchanged) or by an odd number of carbons at both ends (transmission turned on or annihilated). These prescriptions for the occurrence of quantum interference lead to an unexpected consequence for switches which realize such extension through electrocyclic reactions: for some specific attachment modes the chemically closed ring will be the ON position of the switch. Normally the signs of the entries of the Green's function matrix are assumed to have no physical significance; however, we show that the signs may have observable consequences. In particular, in the case of multiple probe attachments – if coherence in probe connections can be arranged – in some cases new destructive interference results, while in others one may have constructive interference. One such case may already exist in the literature.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Concentration dependence of the wings of a dipole-broadened magnetic resonance line in magnetically diluted lattices

The singularities of the time autocorrelation functions (ACFs) of magnetically diluted spin systems with dipole–dipole interaction (DDI), which determine the high-frequency asymptotics of autocorrelation functions and the wings of a magnetic resonance line, are studied. Using the self-consistent fluctuating local field approximation, nonlinear equations are derived for autocorrelation functions averaged over the independent random arrangement of spins (magnetic atoms) in a diamagnetic lattice with different spin concentrations. The equations take into account the specificity of the dipole–dipole interaction. First, due to its axial symmetry in a strong static magnetic field, the autocorrelation functions of longitudinal and transverse spin components are described by different equations. Second, the long-range type of the dipole–dipole interaction is taken into account by separating contributions into the local field from distant and near spins. The recurrent equations are obtained for the expansion coefficients of autocorrelation functions in power series in time. From them, the numerical value of the coordinate of the nearest singularity of the autocorrelation function is found on the imaginary time axis, which is equal to the radius of convergence of these expansions. It is shown that in the strong dilution case, the logarithmic concentration dependence of the coordinate of the singularity is observed, which is caused by the presence of a cluster of near spins whose fraction is small but contribution to the modulation frequency is large. As an example a silicon crystal with different {sup 29}Si concentrations in magnetic fields directed along three crystallographic axes is considered.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Nonlinear ionization mechanism dependence of energy absorption in diamond under femtosecond laser irradiation

We present first-principles calculations for nonlinear photoionization of diamond induced by the intense femtosecond laser field. A real-time and real-space time-dependent density functional theory with the adiabatic local-density approximation is applied to describe the laser-material interactions in the Kohn-Sham formalism with the self-interaction correction. For a certain laser wavelength, the intensity dependence of energy absorption on multiphoton and/or tunnel ionization mechanisms is investigated, where laser intensity regions vary from 10{sup 12} W/cm{sup 2} to 10{sup 16} W/cm{sup 2}. In addition, the effect of laser wavelength on energy absorption at certain ionization mechanism is discussed when the Keldysh parameter is fixed. Theoretical results show that: (1) at the fixed laser wavelength, the relationship between the energy absorption and laser intensity shows a good fit of E = c{sub M}I{sup N} (N is the number of photons absorbed to free from the valence band) when multiphoton ionization dominates; (2) while when tunnel ionization becomes significant, the relationship coincides with the expression of E = c{sub T}I{sup n} (n < N).

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

The bispectrum of cosmic string temperature fluctuations including recombination effects

We calculate the cosmic microwave background temperature bispectrum from cosmic strings, including the contributions from the last scattering surface, using a well-established Gaussian model for the string energy-momentum correlation functions, and a simplified model for the cosmic fluid. We check our approximation for the integrated Sachs-Wolfe (ISW) contribution against the bispectrum obtained from the full sky map of the cosmic string ISW signal used by the Planck team, obtaining good agreement. We validate our model for the last scattering surface contribution by comparing the predicted temperature power spectrum with that obtained from a full Boltzmann code treatment applied to the Unconnected Segment Model of a string network. We find that including the last scattering contribution has only a small impact on the upper limit on the string tension resulting from the bispectrum at Planck resolutions, and argue that the bispectrum is unlikely to be competitive with the power spectrum at any resolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Structural relaxation of water during rapid cooling from ambient temperatures

Experiments investigating the properties of deeply supercooled liquid water are needed to develop a comprehensive understanding of water’s anomalous properties. One approach involves transiently heating nanoscale water films into the supercooled region for several nanoseconds at a time and then interrogating the water films after they have quenched to cryogenic temperatures. To relate the results obtained with this approach to other experiments and simulations on supercooled water, it is important to understand how closely the quenched structure tracks the (metastable) equilibrium structure of water as a function of the transient heating temperature. A key step involves quantifying the extent to which water that is transiently heated to ambient temperatures [hyperquenched water (HQW)] subsequently relaxes toward the structure of low-density amorphous (LDA) ice as it cools. We analyzed the infrared reflection–absorption spectra of LDA, HQW, and crystalline ice films to determine their complex indices of refraction. With this information, we estimate that HQW retains ~50%–60% of a structural motif characteristic of water at high temperatures with the balance comprised of a low-temperature motif. This result, along with results from x-ray diffraction experiments on water and amorphous ices, allows one to quantify the fraction of the high-temperature motif at approximately zero pressure as a function of temperature from 150 to 350 K.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Vibrational coherence and energy transfer in two-dimensional spectra with the optimized mean-trajectory approximation

The optimized mean-trajectory (OMT) approximation is a semiclassical method for computing vibrational response functions from action-quantized classical trajectories connected by discrete transitions that represent radiation-matter interactions. Here, we extend the OMT to include additional vibrational coherence and energy transfer processes. This generalized approximation is applied to a pair of anharmonic chromophores coupled to a bath. The resulting 2D spectra are shown to reflect coherence transfer between normal modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Orbital contributions in the element-resolved valence electronic structure of Bi 2 Se 3

In this work, we studied the bulk band structure of a topological insulator (TI) Bi 2 Se 3 and determined the contributions of the Bi and Se orbital states to the valence bands using standing-wave excited hard x-ray photoemission spectroscopy (SW-HAXPES). This SW technique can provide the element-resolved information and extract individual Bi and Se contributions to the Bi2Se3 valence band. Comparisons with density-functional theory calculations (local density approximation and GW) reveal that the Bi 6s, Bi 6p, and Se 4p states are dominant in the Bi 2 Se 3 HAXPES valence band. Furthermore, these findings pave a way for studying the element-resolved band structure and orbital contributions of this class of TIs.

36 MATERIALS SCIENCE↗

Lattice dynamics and electronic structure of mixed halofluoride scintillators under high pressure

We report the structural, lattice dynamics and electronic structure of mixed halo-fluoride scintillators MClF (M = Ca, Sr) based on density functional theory within generalized gradient approximation. The pressure dependent structural parameters and zone centered phonon frequencies are in reasonable agreement with the available experimental data. Both the structures are dynamically stable up to 50 GPa. MClF compounds are found to have a direct band gap along Γ−Γ direction at ambient as well as under pressure. It is found that these materials can serve as good storage phosphors in the low pressure region from the calculated electronic structure and optical spectra.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Explicitly correlated ring-coupled-cluster-doubles theory

The connection between the random-phase approximation and the ring-coupled-cluster-doubles method bridges the gap between density-functional and wave-function theories and the importance of the random-phase approximation lies in both its broad applicability and this linking role in electronic-structure theory. In this contribution, we present an explicitly correlated approach to the random-phase approximation, based on the direct ring-coupled-cluster-doubles ansatz, which overcomes the problem of slow basis-set convergence, inherent to the random-phase approximation. Benchmark results for a test set of 106 molecules and a selection of 10 organic complexes from the S22 test set demonstrate that convergence to within 99% of the basis-set limit is reached for triple-zeta basis sets for atomisation energies, while quadruple-zeta basis sets are required for interaction energies. Corrections due to single excitations into the complementary auxiliary space reduce the basis-set incompleteness error by one order of magnitude, while contributions due to the coupling of conventional and geminal amplitudes are in general negligible. We find that a non-iterative explicitly correlated correction to first order in perturbation theory exhibits the best ratio of accuracy to computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Diffusion coefficients of transition metals in fcc cobalt

Using first-principles density functional theory (DFT), we calculate the diffusivities of 32 different solute elements—all transition metals, together with Al and Si—in fcc cobalt within the formalism of the five-frequency model. For self-diffusion in fcc cobalt, we compare the accuracy of various approximations to the exchange-correlation energy functional of DFT in estimating the activation energy, and find that only the Perdew-Burke-Ernzerhof (PBE) approximation agrees well with experimental reports and all other functionals largely overestimate it. Our calculations also show that an accurate estimation of the self-diffusion coefficient requires explicit calculation of the effective jump frequency and vacancy formation entropy via phonons. Using accurate self-diffusion data and scaling all solute-related attempt frequencies with respect to the attempt frequency for self-diffusion using a simple relation involving the atomic mass and melting temperature of the solute yields solute diffusivities in excellent agreement with experiments, where such data is available. Additionally, we find that large solutes spontaneously relax toward the nearest neighbor vacancy to relieve the misfit strain, and the extent of this relaxation correlates negatively with the migration energy. Thus, in general, larger solutes have lower migration energies and diffuse faster than smaller solutes in fcc cobalt. However, extremely large solutes, e.g., group III elements Sc, Y, Lu, tend to be trapped in an energy valley located halfway toward the vacancy, and monovacancy mediated diffusion may no longer be valid in such cases. Finally, for all the solutes considered, we systematically tabulate the diffusion-related quantities calculated—diffusion prefactors, migration and activation energies—constructing an extensive and accurate first-principles database for solute diffusion in fcc cobalt.

36 MATERIALS SCIENCE↗

Magnus method for electronic structure calculations at extreme conditions

We present the application of Magnus based methods to the solution of first order coupled ordinary differential equations in High Energy Density (HED) physics applications. Our focus is on the application to quantum mechanical methods, specifically on the solution of the radial Dirac equation for real and complex energies. HED applications require accurate solutions across a wide range of spatial and energy domains, including regimes where the solutions exhibit pronounced oscillatory behavior. Such cases pose significant computational challenges. We demonstrate that Magnus-based integrators can efficiently and accurately address these challenges. We discuss the implementation of the Magnus method for the solution of the radial Dirac equation, including practical considerations such as the evaluation of matrix exponentials, numerical integration, error estimation, and adaptive step size control. We also discuss the application of these methods to complex energy Green’s function techniques and the efficient approximation of integrals of the solutions relevant to HED electronic structure calculations. Here, we demonstrate the accuracy and robustness of the resulting method in applications to the free-particle case, for which analytic solutions are available for comparison, as well as the challenging case of gold at HED conditions.

general physics↗

Electronic structural and magnetic properties of Mn{sub 5}Ge{sub 3} clusters

Theoretical understanding of the stability, ferromagnetism, and spin polarization of Mn{sub 5}Ge{sub 3} clusters has been performed by using the density functional theory with generalized gradient approximation for exchange and correlation. The magnetic moments and magnetic anisotropy energy (MAE) have been calculated for both bulk and clusters, and the enhanced magnetic moment as well as the enlarged MAE have been identified in clusters. The most attractive achievement is that Mn{sub 5}Ge{sub 3} clusters show a fine half-metallic character with large energy scales. The present results may have important implications for potential applications of small Mn{sub 5}Ge{sub 3} clusters as both emerging spintronics and next-generation data-storage technologies.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Real-space formulation of the stress tensor for O (N) density functional theory: Application to high temperature calculations

We present an accurate and efficient real-space formulation of the Hellmann–Feynman stress tensor for O (N) Kohn–Sham density functional theory (DFT). While applicable at any temperature, the formulation is most efficient at high temperature where the Fermi–Dirac distribution becomes smoother and the density matrix becomes correspondingly more localized. We first rewrite the orbital-dependent stress tensor for real-space DFT in terms of the density matrix, thereby making it amenable to O (N) methods. We then describe its evaluation within the O (N) infinite-cell Clenshaw–Curtis Spectral Quadrature (SQ) method, a technique that is applicable to metallic and insulating systems, is highly parallelizable, becomes increasingly efficient with increasing temperature, and provides results corresponding to the infinite crystal without the need of Brillouin zone integration. We demonstrate systematic convergence of the resulting formulation with respect to SQ parameters to exact diagonalization results and show convergence with respect to mesh size to the established plane wave results. Here, we employ the new formulation to compute the viscosity of hydrogen at 10 6 K from Kohn–Sham quantum molecular dynamics, where we find agreement with previous more approximate orbital-free density functional methods.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Calculation of nuclear spin-spin coupling constants using frozen density embedding

We present a method for a subsystem-based calculation of indirect nuclear spin-spin coupling tensors within the framework of current-spin-density-functional theory. Our approach is based on the frozen-density embedding scheme within density-functional theory and extends a previously reported subsystem-based approach for the calculation of nuclear magnetic resonance shielding tensors to magnetic fields which couple not only to orbital but also spin degrees of freedom. This leads to a formulation in which the electron density, the induced paramagnetic current, and the induced spin-magnetization density are calculated separately for the individual subsystems. This is particularly useful for the inclusion of environmental effects in the calculation of nuclear spin-spin coupling constants. Neglecting the induced paramagnetic current and spin-magnetization density in the environment due to the magnetic moments of the coupled nuclei leads to a very efficient method in which the computationally expensive response calculation has to be performed only for the subsystem of interest. We show that this approach leads to very good results for the calculation of solvent-induced shifts of nuclear spin-spin coupling constants in hydrogen-bonded systems. Also for systems with stronger interactions, frozen-density embedding performs remarkably well, given the approximate nature of currently available functionals for the non-additive kinetic energy. As an example we show results for methylmercury halides which exhibit an exceptionally large shift of the one-bond coupling constants between {sup 199}Hg and {sup 13}C upon coordination of dimethylsulfoxide solvent molecules.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accelerating defect predictions in semiconductors using graph neural networks

First-principles computations reliably predict the energetics of point defects in semiconductors but are constrained by the expense of using large supercells and advanced levels of theory. Machine learning models trained on computational data, especially ones that sufficiently encode defect coordination environments, can be used to accelerate defect predictions. Here, we develop a framework for the prediction and screening of native defects and functional impurities in a chemical space of group IV, III–V, and II–VI zinc blende semiconductors, powered by crystal Graph-based Neural Networks (GNNs) trained on high-throughput density functional theory (DFT) data. Using an innovative approach of sampling partially optimized defect configurations from DFT calculations, we generate one of the largest computational defect datasets to date, containing many types of vacancies, self-interstitials, anti-site substitutions, impurity interstitials and substitutions, as well as some defect complexes. We applied three types of established GNN techniques, namely crystal graph convolutional neural network, materials graph network, and Atomistic Line Graph Neural Network (ALIGNN), to rigorously train models for predicting defect formation energy (DFE) in multiple charge states and chemical potential conditions. We find that ALIGNN yields the best DFE predictions with root mean square errors around 0.3 eV, which represents a prediction accuracy of 98% given the range of values within the dataset, improving significantly on the state-of-the-art. We further show that GNN-based defective structure optimization can take us close to DFT-optimized geometries at a fraction of the cost of full DFT. The current models are based on the semi-local generalized gradient approximation-Perdew–Burke–Ernzerhof (PBE) functional but are highly promising because of the correlation of computed energetics and defect levels with higher levels of theory and experimental data, the accuracy and necessity of discovering novel metastable and low energy defect structures at the PBE level of theory before advanced methods could be applied, and the ability to train multi-fidelity models in the future with new data from non-local functionals. The DFT-GNN models enable prediction and screening across thousands of hypothetical defects based on both unoptimized and partially optimized defective structures, helping identify electronically active defects in technologically important semiconductors.

Rahman, Md Habibur (ORCID:000000027705984X)↗

Theoretical estimate of the effect of thermal agitation on ribosome motion generated by stochastic microswimming

The effect of thermal agitation on ribosome motion is evaluated through the Péclet number, assuming that the ribosome is self-propelled along the mRNA during protein synthesis by a swimming stroke consisting of a cycle of stochastically-generated ribosome configurations involving its two subunits. The ribosome velocity probability distribution function is obtained, giving an approximately normal distribution. Its mean and variance together with an estimate of the in vivo free diffusion coefficient of the ribosome and using only configuration changes of small size, give a Péclet number similar to motor proteins and microorganisms. These results suggest the feasibility of the stochastic microswimming hypothesis to explain ribosome motion. - Highlights: • Hypothesis: the ribosome swims along the mRNA during protein synthesis. • How can withstand the thermal agitation of the particles in the surrounding fluid? • The Péclet number Pe is calculated estimating the in vivo diffusion coefficient. • With small size deformations the ribosome has speeds similar to motor proteins. • Similar Pe values for the ribosome and the motor proteins mean feasible hypothesis.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Excitation energies with linear response density matrix functional theory along the dissociation coordinate of an electron-pair bond in N-electron systems

Time dependent density matrix functional theory in its adiabatic linear response formulation delivers exact excitation energies ω{sub α} and oscillator strengths f{sub α} for two-electron systems if extended to the so-called phase including natural orbital (PINO) theory. The Löwdin-Shull expression for the energy of two-electron systems in terms of the natural orbitals and their phases affords in this case an exact phase-including natural orbital functional (PILS), which is non-primitive (contains other than just J and K integrals). In this paper, the extension of the PILS functional to N-electron systems is investigated. With the example of an elementary primitive NO functional (BBC1) it is shown that current density matrix functional theory ground state functionals, which were designed to produce decent approximations to the total energy, fail to deliver a qualitatively correct structure of the (inverse) response function, due to essential deficiencies in the reconstruction of the two-body reduced density matrix (2RDM). We now deduce essential features of an N-electron functional from a wavefunction Ansatz: The extension of the two-electron Löwdin-Shull wavefunction to the N-electron case informs about the phase information. In this paper, applications of this extended Löwdin-Shull (ELS) functional are considered for the simplest case, ELS(1): one (dissociating) two-electron bond in the field of occupied (including core) orbitals. ELS(1) produces high quality ω{sub α}(R) curves along the bond dissociation coordinate R for the molecules LiH, Li{sub 2}, and BH with the two outer valence electrons correlated. All of these results indicate that response properties are much more sensitive to deficiencies in the reconstruction of the 2RDM than the ground state energy, since derivatives of the functional with respect to both the NOs and the occupation numbers need to be accurate.

74 ATOMIC AND MOLECULAR PHYSICS↗

Dynamically consistent method for mixed quantum-classical simulations: A semiclassical approach

We introduce a new semiclassical (SC) framework, the Mixed Quantum-Classical Initial Value Representation (MQC-IVR), that can be tuned to reproduce existing quantum-limit and classical-limit SC approximations to quantum real-time correlation functions. Applying a modified Filinov transformation to a quantum-limit SC formulation leads to the association of a Filinov parameter with each degree of freedom in the system; varying this parameter from zero to infinity controls the extent of quantization of the corresponding mode. The resulting MQC-IVR expression provides a consistent dynamic framework for mixed quantum-classical simulations and we demonstrate its numerical accuracy in the calculation of real-time correlation functions for a model 1D system and a model 2D system over the full range of quantum- to classical-limit behaviors.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗