Search NASA⌕ Search

SEARCH · Search NASA

Results for “semiempirical methods”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
NASA NTRS1991

Disordering and H-embrittlement of Pb(110) surface using embedded atom method and molecular dynamics

The dynamics and structure of the crystal melt interface are of great importance in crystal growth, melting, and other properties of the solid-liquid interface. The Embedded Atom Method (EAM), a semiempirical method developed by Daw and Baskes for determining the energetics of atoms in a bulk environment, is a useful method in the computer simulation. The EAM functions of PB in conjunction with the MD code are employed to study melting and disordering of the Pb(100) surface. The EAM functions are fitted to the bulk experimental data at zero temperature, and it is not clear if it will behave correctly at higher temperatures. In fact, it is shown that use of a similar type of potential at higher temperatures might result in a negative coefficient of thermal expansion. The primary purpose here is to perform MD simulation combined with the EAM functions of PB to study melting of the Pb(100) surface. This provides an accurate test of the EAM functions at higher temperatures. In particular, we calculate the two dimensional structure factor, the two dimensional distribution function, and density profile as a function of temperature up to the melting point of lead.

Karimi, Majid↗
NASA NTRS1994

Modelling of Surfaces. Part 2: Metallic Alloy Surfaces Using the BFS Method

Using BFS, a new semiempirical method for alloys, we study the surface structure of fcc ordered binary alloys. We concentrate on the calculation of surface energies and surface relaxations for the L1(sub 0) and L1(sub 2) ordered structures. Different terminations of the low-index faces are studied. Also, we present results for the interlayer relaxations for planes close to the surface, revealing different relaxations for atoms of different species producing a rippled surface layer.

Bozzolo, Guillermo↗
NASA NTRS1994

Composition Dependence of Bulk Alloy Properties

We introduce an approximate expression for the direct calculation of properties of alloys in terms of the pure components. This rule can be obtained as a particular case from the Bozzolo-Ferrante-Smith (BFS) method, a recently developed semiempirical method for alloys. In particular, we examine the application of this rule to several examples of the concentration dependence of the lattice parameter of binary and ternary alloys.

Bozzolo, Guillermo↗
NASA NTRS1992

Method for calculating alloy energetics

A semiempirical method for the computation of alloy energies is introduced. It is based on the equivalent-crystal theory of defect-formation energies in elemental solids. The method is both simple and accurate. Heats of formation as a function of composition are computed for some binary alloys of Cu, Ni, Al, Ag, Pd, Pt, and Au using the heats of solution in the dilute limit as experimental input. The separation of heats into strain and chemical components helps in understanding the energetics. In addition, lattice-parameter contractions seen in solid solutions of Ag and Au are accurately predicted. Good agreement with experiment is obtained in all cases.

Bozzolo, Guillermo↗
NASA NTRS1992

Determination of parameters of a new method for predicting alloy properties

Recently, a semiempirical method for alloys based on equivalent crystal theory was introduced. The method successfully predicts the concentration dependence of the heat of formation and lattice parameter of binary alloys. A study of the parameters of the method is presented, along with new results for (gamma)Fe-Pd and (gamma)Fe-Ni alloys.

Bozzolo, Guillermo↗
NASA NTRS1986

Aircraft Noise Prediction Program theoretical manual: Propeller aerodynamics and noise

The prediction sequence used in the aircraft noise prediction program (ANOPP) is described. The elements of the sequence are called program modules. The first group of modules analyzes the propeller geometry, the aerodynamics, including both potential and boundary-layer flow, the propeller performance, and the surface loading distribution. This group of modules is based entirely on aerodynamic strip theory. The next group of modules deals with the first group. Predictions of periodic thickness and loading noise are determined with time-domain methods. Broadband noise is predicted by a semiempirical method. Near-field predictions of fuselage surface pressrues include the effects of boundary layer refraction and scattering. Far-field predictions include atmospheric and ground effects.

Zorumski, W. E.↗
NASA NTRS1957

Base Pressure at Supersonic Speeds on Two-dimensional Airfoils and on Bodies of Revolution with and Without Fins Having Turbulent Boundary Layers

An analysis has been made of available experimental data to show the effects of most of the variables that are more predominant in determining base pressure at supersonic speeds. The analysis covers base pressures for two-dimensional airfoils and for bodies of revolution with and without stabilizing fins and is restricted to turbulent boundary layers. The present status of available experimental information is summarized as are the existing methods for predicting base pressure. A simple semiempirical method is presented for estimating base pressure. For two-dimensional bases, this method stems from an analogy established between the base-pressure phenomena and the peak pressure rise associated with the separation of the boundary layer. An analysis made for axially symmetric flow indicates that the base pressure for bodies of revolution is subject to the same analogy. Based upon the methods presented, estimations are made of such effects as Mach number, angle of attack, boattailing, fineness ratio, and fins. These estimations give fair predictions of experimental results. (author)

LOVE EUGENE S↗
NASA NTRS1951

A semi-empirical method for calculating the pitching moment of bodies of revolution at low Mach numbers

A semiempirical method, in which potential theory is arbitrarily combined with an approximate viscous theory, for calculating the aerodynamic pitching moments for bodies of revolution is presented. The method can also be used for calculating the lift and drag forces. The calculated and experimental force and moment characteristics of 15 bodies of revolution are compared.

FLOW, VISCOUS↗
NASA NTRS1993

Applicability of PM3 to transphosphorylation reaction path: Toward designing a minimal ribozyme

A growing body of evidence shows that RNA can catalyze many of the reactions necessary both for replication of genetic material and the possible transition into the modern protein-based world. However, contemporary ribozymes are too large to have self-assembled from a prebiotic oligonucleotide pool. Still, it is likely that the major features of the earliest ribozymes have been preserved as molecular fossils in the catalytic RNA of today. Therefore, the search for a minimal ribozyme has been aimed at finding the necessary structural features of a modern ribozyme (Beaudry and Joyce, 1990). Both a three-dimensional model and quantum chemical calculations are required to quantitatively determine the effects of structural features of the ribozyme on the reaction it catalyzes. Using this model, quantum chemical calculations must be performed to determine quantitatively the effects of structural features on catalysis. Previous studies of the reaction path have been conducted at the ab initio level, but these methods are limited to small models due to enormous computational requirements. Semiempirical methods have been applied to large systems in the past; however, the accuracy of these methods depends largely on a simple model of the ribozyme-catalyzed reaction, or hydrolysis of phosphoric acid. We find that the results are qualitatively similar to ab initio results using large basis sets. Therefore, PM3 is suitable for studying the reaction path of the ribozyme-catalyzed reaction.

Manchester, John I.↗
NASA NTRS1981

VORCAM: A computer program for calculating vortex lift effect of cambered wings by the suction analogy

A user's guide to an improved version of Woodward's chord plane aerodynamic panel computer code is presumed. The guide can be applied to cambered wings exhibiting edge separated flow, including those with leading edge vortex flow at subsonic and supersonic speeds. New orientations for the rotated suction force are employed based on the momentum principal. The supersonic suction analogy method is improved by using an effective angle of attack defined through a semiempirical method.

Lan, C. E.↗
NASA NTRS1990

Application of panel methods to wind-tunnel wall interference corrections

Wind-tunnel wall interference effects were estimated using numerical solutions capable of computing the flow field over complex three-dimensional bodies. This approach was then tested by applying it to two typical wind tunnel experiments involving a high-lift wing and a fighter airplane model. In both cases, the static pressure signature along the test section walls was measured and computed for additional validation. The wind-tunnel wall correction method described here, which inherently includes effects of lift and blockage, provides more details than previous semiempirical methods.

Browne, Lindsey↗
DOE OSTI2025

Extending quantum-mechanical benchmark accuracy to biological ligand-pocket interactions

Predicting the binding affinity of ligands to protein pockets is key in the drug design pipeline. The flexibility of ligand-pocket motifs arises from a range of attractive and repulsive electronic interactions during binding. Accurately accounting for all interactions requires robust quantum-mechanical (QM) benchmarks, which are scarce for ligand-pocket systems. Additionally, disagreement between “gold standard” Coupled Cluster (CC) and Quantum Monte Carlo (QMC) methods casts doubt on many benchmarks for larger non-covalent systems. We introduce the “QUantum Interacting Dimer” (QUID) benchmark framework containing 170 non-covalent (non-)equilibrium systems modeling chemically and structurally diverse ligand-pocket motifs. Symmetry-adapted perturbation theory shows that QUID broadly covers non-covalent binding motifs and energetic contributions. Robust binding energies are obtained using complementary CC and QMC methods, achieving agreement of 0.5 kcal/mol. The benchmark data analysis reveals that several dispersion-inclusive density functional approximations provide accurate energy predictions, though their atomic van der Waals forces differ in magnitude and orientation. Contrarily, semiempirical methods and empirical force fields require improvements in capturing non-covalent interactions (NCIs) for out-of-equilibrium geometries. The wide span of NCIs, highly accurate interaction energies, and analysis of molecular properties take QUID beyond the “gold standard” for QM benchmarks of ligand-protein systems.

Puleva, Mirela [University of Luxembourg, Luxembou↗
DOE OSTI2025

Convergent Protocols for Computing Protein–Ligand Interaction Energies Using Fragment-Based Quantum Chemistry

Fragment-based quantum chemistry methods offer a way to sidestep the steep nonlinear scaling of electronic structure calculations so that large molecular systems can be investigated using high-level methods. Here, we use fragmentation to compute protein–ligand interaction energies in systems with several thousand atoms, using a new software platform for managing fragment-based calculations that implements a screened many-body expansion. Convergence tests using a minimal-basis semiempirical method (HF-3c) indicate that two-body calculations, with single-residue fragments and simple hydrogen caps, are sufficient to reproduce interaction energies obtained using conventional supramolecular electronic structure calculations, to within 1 kcal/mol at about 1% of the computational cost. We also demonstrate that the HF-3c results are illustrative of trends obtained with density functional theory in basis sets up to augmented quadruple-ζ quality. Strategic deployment of fragmentation facilitates the use of converged biomolecular model systems alongside high-quality electronic structure methods and basis sets, bringing ab initio quantum chemistry to systems of hitherto unimaginable size. This will be useful for generation of high-quality training data for machine learning applications.

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

Quick-and-Easy Validation of Protein–Ligand Binding Models Using Fragment-Based Semiempirical Quantum Chemistry

Electronic structure calculations in enzymes converge very slowly with respect to the size of the model region that is described using quantum mechanics (QM), requiring hundreds of atoms to obtain converged results and exhibiting substantial sensitivity (at least in smaller models) to which amino acids are included in the QM region. As such, there is considerable interest in developing automated procedures to construct a QM model region based on well-defined criteria. However, testing such procedures is burdensome due to the cost of large-scale electronic structure calculations. Here, we show that semiempirical methods can be used as alternatives to density functional theory (DFT) to assess convergence in sequences of models generated by various automated protocols. The cost of these convergence tests is reduced even further by means of a many-body expansion. We use this approach to examine convergence (with respect to model size) of protein–ligand binding energies. Fragment-based semiempirical calculations afford well-converged interaction energies in a tiny fraction of the cost required for DFT calculations. Two-body interactions between the ligand and single-residue amino acid fragments afford a low-cost way to construct a “QM-informed” enzyme model of reduced size, furnishing an automatable active-site model-building procedure. This provides a streamlined, user-friendly approach for constructing ligand binding-site models that needs neither a priori information nor manual adjustments. Extension to model-building for thermochemical calculations should be straightforward.

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

Generalization of Equivalent Crystal Theory to Include Angular Dependence

In the original Equivalent Crystal Theory, each atomic site in the real crystal is assigned an equivalent lattice constant, in general different from the ground state one. This parameter corresponds to a local compression or expansion of the lattice. The basic method considers these volumetric transformations and, in addition, introduces the possibility that the reference lattice is anisotropically distorted. These distortions however, were introduced ad-hoc. In this work, we generalize the original Equivalent Crystal Theory by systematically introducing site-dependent directional distortions of the lattice, whose corresponding distortions account for the dependence of the energy on anisotropic local density variations. This is done in the spirit of the original framework, but including a gradient term in the density. This approach is introduced to correct a deficiency in the original Equivalent Crystal Theory and other semiempirical methods in quantitatively obtaining the correct ratios of the surface energies of low index planes of cubic metals (100), (110), and (111). We develop here the basic framework, and apply it to the calculation of Fe (110) and Fe (111) surface energy formation. The results, compared with first principles calculations, show an improvement over previous semiempirical approaches.

Ferrante, John↗
NASA NTRS2014

Acoustics Research of Propulsion Systems

The liftoff phase induces some of the highest acoustic loading over a broad frequency for a launch vehicle. These external acoustic environments are used in the prediction of the internal vibration responses of the vehicle and components. Thus, predicting these liftoff acoustic environments is critical to the design requirements of any launch vehicle but there are challenges. Present liftoff vehicle acoustic environment prediction methods utilize stationary data from previously conducted hold-down tests; i.e. static firings conducted in the 1960's, to generate 1/3 octave band Sound Pressure Level (SPL) spectra. These data sets are used to predict the liftoff acoustic environments for launch vehicles. To facilitate the accuracy and quality of acoustic loading, predictions at liftoff for future launch vehicles such as the Space Launch System (SLS), non-stationary flight data from the Ares I-X were processed in PC-Signal in two forms which included a simulated hold-down phase and the entire launch phase. In conjunction, the Prediction of Acoustic Vehicle Environments (PAVE) program was developed in MATLAB to allow for efficient predictions of sound pressure levels (SPLs) as a function of station number along the vehicle using semiempirical methods. This consisted, initially, of generating the Dimensionless Spectrum Function (DSF) and Dimensionless Source Location (DSL) curves from the Ares I-X flight data. These are then used in the MATLAB program to generate the 1/3 octave band SPL spectra. Concluding results show major differences in SPLs between the hold-down test data and the processed Ares IX flight data making the Ares I-X flight data more practical for future vehicle acoustic environment predictions.

Gao, Ximing↗
DOE OSTI2024

Density Functional Tight-Binding Models for Band Structures of Transition-Metal Alloys and Surfaces across the d -Block

First-principles electronic structure simulations are an invaluable tool for understanding chemical bonding and reactions. While machine-learning models such as interatomic potentials significantly accelerate the exploration of potential energy surfaces, electronic structure information is generally lost. Particularly in the field of heterogeneous catalysis, simulated electron band structures provide fundamental insights into catalytic reactivity. This ab initio knowledge is preserved in semiempirical methods such as density functional tight binding (DFTB), which extend the accessible computational length and time scales beyond first-principles approaches. In this paper here we present Shell-Optimized Atomic Confinement (SOAC) DFTB electronic-part-only parametrizations for bulk and surface band structures of all d-block transition metals that enable efficient predictions of electronic descriptors for large structures or high-throughput studies on complex systems outside the computational reach of density functional theory.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗