Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “basis functions”

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 865 records · Page 48

Myosin dilated cardiomyopathy mutation S532P disrupts actomyosin interactions, leading to altered muscle kinetics, reduced locomotion, and cardiac dilation in Drosophila

Dilated cardiomyopathy (DCM), a life-threatening disease characterized by pathological heart enlargement, can be caused by myosin mutations that reduce contractile function. To better define the mechanistic basis of this disease, we employed the powerful genetic and integrative approaches available in Drosophila melanogaster. To this end, we generated and analyzed the first fly model of human myosin–induced DCM. The model reproduces the S532P human β-cardiac myosin heavy chain DCM mutation, which is located within an actin-binding region of the motor domain. In concordance with the mutation’s location at the actomyosin interface, steady-state ATPase and muscle mechanics experiments revealed that the S532P mutation reduces the rates of actin-dependent ATPase activity and actin binding and increases the rate of actin detachment. The depressed function of this myosin form reduces the number of cross-bridges during active wing beating, the power output of indirect flight muscles, and flight ability. Further, S532P mutant hearts exhibit cardiac dilation that is mutant gene dose–dependent. Our study shows that Drosophila can faithfully model various aspects of human DCM phenotypes and suggests that impaired actomyosin interactions in S532P myosin induce contractile deficits that trigger the disease.

59 BASIC BIOLOGICAL SCIENCES↗

Pushing Cu uphill of the volcano curve: Impact of a WC support on the catalytic activity of copper toward the hydrogen evolution reaction

Here, the adsorption of atomic H and H 2 on copper mono- and submonolayers supported on hexagonal WC(0001) surfaces has been investigated using density functional theory with the Perdew–Burke–Ernzerhof exchange correlation functional and D2 van der Waals corrections. Results evidence the impact of the termination of the carbide substrate on fundamental properties of Cu adatoms, and, hence, on the stability of molecular and atomic hydrogen, defining copper's catalytic activity for hydrogen evolution reaction. Using H adsorption energy as a descriptor, catalytic activity of Cu adlayers for hydrogen evolution reaction was estimated using traditional volcano curves and a curve, obtained at low hydrogen coverage. Obtained results evidence that copper adlayers supported on the WC may present a viable low-cost alternative to noble metal-based catalysts, with improved catalytic activity compared to that of copper. This, potentially, can be a useful basis for designing and developing novel functional materials with predetermined catalytic properties.

08 HYDROGEN↗

Excitation of the metastable E(3 Sigma g plus) state of N2 by electron impact

The contribution of the N2(E(3 Sigma g plus)) state to the total metastable excitation function of N2 assessed on the basis of time-of-flight studies of metastable nitrogen molecules. The cross section for electron impact excitation state was determined in the domain of the resonance form threshold (11.87 eV) to an energy of about 13 eV. The maximum value of the cross section was found to be (7.0 + or - 4.0) x 10 to the -18th power sq cm at an energy of 12.2 eV. The measurement was made absolute by using the previously determined yield of the metastable detector, the lifetime of the E state, and by eliminating the energy spread in the electron beam from the raw data. The half-width (FWHM) of the resonance-like excitation function near threshold was found to be about 0.4 eV. No substantial evidence was obtained from the present data for the presence of the nonresonant part of the excitation function for the state studied.

Borst, W. L.↗

Improvements in the orbitalwise scaling down of Perdew–Zunger self-interaction correction in many-electron regions

The Perdew-Zunger (PZ) method offers a way to remove the self-interaction (SI) error from density functional approximations on an orbital by orbital basis. The PZ method provides significant improvements for the properties such as barrier heights or dissociation energies but results in over-correcting the properties well described by SI-uncorrected semi-local functional. One cure to rectify the over-correcting tendency is to scale down the magnitude of SI-correction of each orbital in the many electron region. We have implemented the orbitalwise scaled down SI-correction (OSIC) scheme of Vydrov et al. [J. Chem. Phys. 124, 094108 (2006)] using the Fermi-L owdin SI-correction method. After validating the OSIC implementation with previously reported OSICLSDA results, we examine its performance with the most successful non-empirical SCAN meta-GGA functional. Using different forms of scaling factors to identify one-electron regions, we assess the performance of OSIC-SCAN for a wide range of properties: total energies, ionization potentials and electron affinities for atoms, atomization energies, dissociation and reaction energies, and reaction barrier heights of molecules. Our results show that OSIC-SCAN provides superior results than the previously reported OSIC-LSDA, -PBE, and -TPSS results. Furthermore, we propose selective scaling of OSIC (SOSIC) to remove its major shortcoming that destroys the -1/r asymptotic behavior of the potentials. The SOSIC method gives the highest occupied orbital eigenvalues practically identical to those in PZSIC and unlike OSIC provides bound atomic anions even with larger powers of scaling factors. Furthermore, SOSIC compared to PZSIC or OSIC provides more balanced description of total energies and barrier heights.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Semiempirical Theories of the Affinities of Negative Atomic Ions

The determination of the electron affinities of negative atomic ions by means of direct experimental investigation is limited. To supplement the meager experimental results, several semiempirical theories have been advanced. One commonly used technique involves extrapolating the electron affinities along the isoelectronic sequences, The most recent of these extrapolations Is studied by extending the method to Include one more member of the isoelectronic sequence, When the results show that this extension does not increase the accuracy of the calculations, several possible explanations for this situation are explored. A different approach to the problem is suggested by the regularities appearing in the electron affinities. Noting that the regular linear pattern that exists for the ionization potentials of the p electrons as a function of Z, repeats itself for different degrees of ionization q, the slopes and intercepts of these curves are extrapolated to the case of the negative Ion. The method is placed on a theoretical basis by calculating the Slater parameters as functions of q and n, the number of equivalent p-electrons. These functions are no more than quadratic in q and n. The electron affinities are calculated by extending the linear relations that exist for the neutral atoms and positive ions to the negative ions. The extrapolated. slopes are apparently correct, but the intercepts must be slightly altered to agree with experiment. For this purpose one or two experimental affinities (depending on the extrapolation method) are used in each of the two short periods. The two extrapolation methods used are: (A) an isoelectronic sequence extrapolation of the linear pattern as such; (B) the same extrapolation of a linearization of this pattern (configuration centers) combined with an extrapolation of the other terms of the ground configurations. The latter method Is preferable, since it requires only experimental point for each period. The results agree within experimental error with all data, except with the most recent value of C, which lies 10% lower.

Edie, John W.↗

The Impact of Nonequilibrium Ionization on SDO/AIA and Hinode/EIS Observations of Impulsively Heated Plasmas

Most plasma diagnostics assume the emitting material is in a state of ionization equilibrium. For example, the AIA temperature response functions have been derived on this basis. The assumption is reasonable whenever the plasma is evolving slowly or is very dense, but these are not the conditions that apply during impulsive heating events. It is now widely believed that many coronal loops are bundles of unresolved strands that are heated quasi-randomly by nanoflares. Full blown flares are thought to have similar sub-structure. We have studied the importance of nonequilibrium effects in these circumstances by modeling nanoflare-heated loops and simulating their observation by AIA and the EIS spectrometer on Hinode. We find that the intensities of hot emission lines can be highly suppressed and that the net emission from the loop tends to be dominated by strands that have entered a slow cooling phase, well after the impulsive energy release has ended. The hottest strands are relatively invisible, both because they are tenuous and because they cool rapidly by thermal conduction. Thus, AIA channels that are normally thought of as being sensitive to hot plasma, such 131 and 94, are in fact frequently not able to detect the hot plasma that is present. The magnitude of the effect is case dependent. Great care must be exercised when using the standard temperature response functions in situations where nonequilibrium ionization is likely to be important.

Klimchuk, James A.↗

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗

QMCPACK: Advances in the development, efficiency, and application of auxiliary field and real-space variational and diffusion quantum Monte Carlo

We review recent advances in the capabilities of the open source ab initio Quantum Monte Carlo (QMC) package QMCPACK and the workflow tool Nexus used for greater efficiency and reproducibility. The auxiliary field QMC (AFQMC) implementation has been greatly expanded to include k-point symmetries, tensor-hypercontraction, and accelerated graphical processing unit (GPU) support. These scaling and memory reductions greatly increase the number of orbitals that can practically be included in AFQMC calculations, increasing the accuracy. Advances in real space methods include techniques for accurate computation of bandgaps and for systematically improving the nodal surface of ground state wavefunctions. Results of these calculations can be used to validate application of more approximate electronic structure methods, including GW and density functional based techniques. To provide an improved foundation for these calculations, we utilize a new set of correlation-consistent effective core potentials (pseudopotentials) that are more accurate than previous sets; these can also be applied in quantum-chemical and other many-body applications, not only QMC. These advances increase the efficiency, accuracy, and range of properties that can be studied in both molecules and materials with QMC and QMCPACK.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Near Hartree-Fock quality GTO basis sets for the second-row atoms

Energy optimized, near Hartree-Fock quality Gaussian basis sets ranging in size from (17s12p) to (20s15p) are presented for the ground states of the second-row atoms for Na(2P), Na(+), Na(-), Mg(3P), P(-), S(-), and Cl(-). In addition, optimized supplementary functions are given for the ground state basis sets to describe the negative ions, and the excited Na(2P) and Mg(3P) atomic states. The ratios of successive orbital exponents describing the inner part of the 1s and 2p orbitals are found to be nearly independent of both nuclear charge and basis set size. This provides a method of obtaining good starting estimates for other basis set optimizations.

Partridge, Harry↗

Near Hartree-Fock quality GTO basis sets for the second-row atoms

Energy optimized, near Hartree-Fock quality Gaussian basis sets ranging in size from (17s12p) to (20s15p) are presented for the ground states of the second-row atoms for Na(2P), Na(+), Na(-), Mg(3P), P(-), S(-), and Cl(-). In addition, optimized supplementary functions are given for the ground state basis sets to describe the negative ions, and the excited Na(2P) and Mg(3P) atomic states. The ratios of successive orbital exponents describing the inner part of the 1s and 2p orbitals are found to be nearly independent of both nuclear charge and basis set size. This provides a method of obtaining good starting estimates for other basis set optimizations.

Partridge, Harry↗

TiCl, TiH and TiH+ Bond Energies, a Test of a Correlation Consistent Ti Basis Set

Correlation consistent basis sets are developed for Ti atom. The polarization functions are optimized for the average of the 3F and 5F states. One series of correlation consistent basis sets is for 3d and 4s correlation, while the second series includes 3s and 3p correlation as well as 3d and 4s correlation. These basis sets are tested using the Ti 3F-5F separation and the dissociation energies of TiCl X4Phi, TiH X4Phi, and TiH(+) X3Phi. The CCSD(T) complete basis set limit values are determined by extrapolation. The Douglas Kroll approach is used to compute the scalar relativistic effect. Spin-orbit effects are taken from experiment and/or computed at the CASSCF level. The Ti 3F-5F separation is in excellent agreement with experiment, while the TiCl, TiH, and TiH(+) bond energies are in good agreement with experiment. Extrapolation with the valence basis set is consistent with other atoms, while including 3s and 3p correlation appears to make extrapolation.

Bauschlicher, Charles W., Jr.↗

Finite-field Calculations of Molecular Polarizabilities Using Field-Induced Polarization Functions: Second- and Third-order Perturbation Correlation Corrections to the Coupled Hartree-Fock Polarizability of H2O

Ordinary Rayleigh-Schroudinger perturbation theory with Moller-Plesset (RSMP) partitioning is used to calculate second- and third-order correlation corrections to the CHF polarizability and dipole moment of the water molecule by a finite-field procedure. Pade approximants are found to be useful in accelerating the convergence of the property perturbation expansions. Field-induced polarization functions suitable for polarizability calculations are determined. The average polarizability calculated, neglecting vibrational averaging, with Dunning's (9s5p/4s-4s2p/2s) contracted GTO basis set augmented by field-induced lslp2d/lp polarization functions is within 3 per cent of the experimental result. Correlation corrections to the dipole moment and polarizability of the water molecule calculated by the finite-field RSMP and single + double excitation CI(SDCI) methods for the same basis set are found to be in close agreement. The RSMP approach has the advantages of being size-consistent and of being capable of greater efficiency than the SCDI method. Comparative calculations carried out using Epstein-Nesbet partitioning show that through third order RSEN correlation perturbation expansions for the dipole moment and polarizability are less rapidly convergent than RSMP expansions. However, reasonable accord with RSMP results can be achieved by using Pade approximants to accelerate the convergence of RSEN energy perturbation expansions. The convergence of RSMP property correlation expansions based on the zeroth-order uncoupled-Hartree-Fock (UCHF) and coupled-Hartree-Fock (CHF) approximations are compared through third order. Whereas the CHF + RSMP expansions are for practical purposes fully converged, the UCHF + RSMP expansions are not adequately converged.

Langhoff, S. R.↗

On the light-front wave functions of quarkonia

The light-front wave functions of hadrons allow us to calculate a wide range of physical observables; however, the wave functions themselves cannot be measured. We discuss recent results for quarkonia obtained in basis light-front quantization using an effective Hamiltonian with a confining model in both the transverse and longitudinal directions and with explicit one-gluon exchange. In particular, we focus on the numerical convergence of the basis expansion, as well as the asymptotic behavior of the light-front wave functions. Finally, we also illustrate that, for mesons with unequal quark masses, the maxima of the light-front wave functions depend in a non-trivial way on the valence quark-mass difference.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Magnetic exchange interactions in binuclear and tetranuclear iron(III) complexes described by spin‐flip DFT and Heisenberg effective Hamiltonians

Abstract Low‐energy spectra of single‐molecule magnets (SMMs) are often described by Heisenberg Hamiltonians. Within this formalism, exchange interactions between magnetic centers determine the ground‐state multiplicity and energy separation between the ground and excited states. In this contribution, we extract exchange coupling constants ( J ) for a set of iron (III) binuclear and tetranuclear complexes from all‐electron calculations using non‐collinear spin‐flip time‐dependent density functional theory (NC‐SF‐TDDFT). For 12 binuclear complexes with J ‐values ranging from −6 to −132 cm −1 , our benchmark calculations using the short‐range hybrid ω PBEh functional and 6‐31G(d,p) basis set agree well with the experimentally derived values (mean absolute error of 4.7 cm −1 ). For the tetranuclear SMMs, the computed J constants are within 6 cm −1 from the experimentally derived values. We explore the range of applicability of the Heisenberg model by analyzing bonding patterns in these Fe(III) complexes using natural orbitals (NO), their occupations, and the number of effectively unpaired electrons. The results illustrate the efficiency of the spin‐flip protocol for computing the exchange couplings and the utility of the NO analysis in assessing the validity of effective spin Hamiltonians.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Leading-color two-loop amplitudes for four partons and a $\textit{W}$ boson in QCD

We present the leading-color two-loop QCD corrections for the scattering of four partons and a W boson, including its leptonic decay. The amplitudes are assembled from the planar two-loop helicity amplitudes for four partons and a vector boson decaying to a lepton pair, which are also used to determine the planar two-loop amplitudes for four partons and a $Z/γ^\ast$ boson with a leptonic decay. The analytic expressions are obtained by setting up a dedicated Ansatz and constraining the free parameters from numerical samples obtained within the framework of numerical unitarity. The large linear systems that must be solved to determine the analytic expressions are constructed to be in Vandermonde form. Such systems can be very efficiently solved, bypassing the bottleneck of Gaussian elimination. Our results are expressed in a basis of one-mass pentagon functions, which opens the possibility of their efficient numerical evaluation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A case study of high-temperature polyetherimide film capacitor fabrication

Capacitors, as the basic building block of power electronics and electrical systems, are the major constraint of the increasingly integrated power systems that require new capable polymer dielectric films operating at higher temperatures >125 °C. Despite tremendous research efforts on the lab scale, there remains a considerable barrier and knowledge basis for converting polymeric films to functioning capacitors desired in actual applications. Here, the authors tackled the polyetherimide (PEI) film scale up issues and developed various engineering processes for film de-wrinkling, optimal metallization, static elimination, and capacitor fabrication improvement. The authors discovered various fabrication challenges such as the electrostatic charge, metallization scheme, winding tension, and end electrode disconnects. To fabricate a high yield of PEI capacitor bobbins, it is necessary to utilize static eliminators (radioactive ionizers), moderately thick aluminum metallization (15–30 Ω/sq), and winding tension (60–100 g). It is also effective to evaluate capacitance, dielectric loss, equivalent series resistance, and thermal cycling stability of capacitors. This work sheds some light on converting PEI films and other polymeric films to capacitors on a scale up fabrication effort. The practical learning from film handling to capacitor fabrication in this work provided the necessary knowledge for manufacturing high-temperature polar film capacitors.

36 MATERIALS SCIENCE↗

Energy-Screened Many-Body Expansion for Protein–Ligand Interactions: Examining Convergence for Metalloenzymes Through Seven–Body Interactions

Fragment-based quantum chemistry is a powerful strategy for calculating protein−ligand interaction energies using quantum chemistry methods. Rigorous convergence often requires hundreds of atoms in the protein binding-site model, especially if that model is constructed using distance-based criteria to select amino acid residues, while three- and four-body calculations exhibit instability related to combinatorial proliferation in the number of subsystem calculations. Here, we report an energy-based screening protocol for the many-body expansion applied to protein−ligand interactions, implemented in the open-source FRAGME∩T code. Using a combination of aggressive screening based on semiempirical quantum chemistry, with an improved graph-theoretical algorithm to eliminate unimportant subsystems, we are able to perform n-body calculations up to n = 7 using density functional theory in triple-ζ basis sets. Distance cutoffs further reduce the cost without compromising accuracy. Rapid and stable convergence of the many-body expansion is obtained by n = 4, for a pair of metalloenzymes in which a divalent ion coordinates directly to the ligand. As compared to previous results that relied solely on distance cutoffs, oscillations in the n-body corrections are reduced or eliminated, although residual errors remain in one case. This work demonstrates that benchmark-quality protein−ligand interaction energies can be systematically converged using a method with excellent parallel efficiency and scalability.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Attractive Way to Correct for Missing Singles Excitations in Unitary Coupled Cluster Doubles Theory

Coupled cluster methods based exclusively on double excitations are comparatively “cheap” and interesting model chemistries, as they are typically able to capture the bulk of the dynamic electron correlation effects. The trade-off in such approximations is that the effect of neglected excitations, particularly single excitations, can be considerable. Using standard and electron-pair-restricted T 2 operators to define two flavors of unitary coupled cluster doubles (UCCD) methods, we investigate the extent to which missing single excitations can be recovered from low-order corrections in many-body perturbation theory (MBPT) within the unitary coupled cluster (UCC) formalism. Here, our analysis includes the derivations of finite-order UCC energy functionals, which are used as a basis to define perturbative estimates of missed single excitations. This leads to the novel UCCD[4S] and UCCD[6S] methods, which consider energy corrections for missing single excitations through fourth- and sixth-order in MBPT, respectively. We also apply the same methodology to the electron-pair-restricted ansatz, but the improvements are only marginal. Our findings show that augmenting UCCD with these post hoc perturbative corrections can lead to UCCSD-quality results.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗