An approximate determination of the lift of slender cylindrical bodies and wing-body combinations at very high supersonic speeds
Explore the source record for details and available documents.
SEARCH · Search NASA
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.
Explore the source record for details and available documents.
Within many-body perturbation theory, Hedin's formalism offers a systematic way to iteratively compute the self-energy Σ of any dynamically correlated interacting system, provided one can evaluate the interaction vertex Γ exactly. This is, however, impossible, in general, for it involves the functional derivative of Σ with respect to the Green's function. Here, we analyze the structure of this derivative, splitting it into four contributions and outlining the type of quasiparticle interactions that each of them generate. Moreover, we show how, in the implementation of self-consistency, the action of these contributions can be classified into two: A quantitative renormalization of previously included interaction terms and the inclusion of qualitatively distinct interaction terms through successive functional derivatives of Γ itself. Implementing this latter type of self-consistency can extend the validity of perturbative approximations based on Hedin's equations toward the high interaction limit, as we show in the example of the Hubbard dimer. Furthermore, our analysis also provides a unifying perspective on the perturbation theory landscape, showing how the T-matrix approach is completely contained in Hedin's formalism.
Abstract Far from equilibrium phenomenon is a central theme of contemporary material research. Such phenomenon can exhibit itself in atomic structure and dynamics, but very often it also happens as non‐equilibrium phenomenon in the electronic structure. In ab initio material simulation, density functional theory (DFT) has played an essential role in studying electronic ground state problems. For excited states, besides many‐body perturbation theory, another powerful tool is the time dependent DFT (TDDFT) method. In particular, the real‐time TDDFT (rt‐TDDFT) method can be used to simulate many non‐equilibrium phenomena directly. Here we introduce our works on some algorithm advances based on our recently rt‐TDDFT method. This method uses the plane‐wave basis set, and significantly accelerates its efficiency by increasing the time step from 0.1–1 as in traditional methods to 0.2–0.5 fs. The noncollinear magnetic moments and spin–orbit coupling have also been included in our rt‐TDDFT method. Furthermore, a Boltzmann‐TDDFT algorithm has been developed to solve the hot carrier overheating problem in Ehrenfest dynamics, and a natural orbital branching algorithm has been developed to overcome the mean‐field approximation in Ehrenfest dynamics nuclear trajectory, thus allows stochastic multiple paths in chemical reactions. Utilizing these methods, we have studied the photoinduced ultrafast demagnetization, ultrafast phase transition, energy transfer between plasmon and hot carriers, as well as the high‐energy ion implantation and low‐energy atomic diffusion in semiconductors. We believe the tools as the ones introduced here can enable us to study a wide range of phenomena which are of great interest in modern day material research. This article is categorized under: Structure and Mechanism > Computational Materials Science Electronic Structure Theory > Ab Initio Electronic Structure Methods Electronic Structure Theory > Density Functional Theory
We develop expressions for electron density defined through the linear response for general density functional approximations, demonstrating results for orbital functionals and for many-body perturbation theory, and explore the connections to developments in DFT.
Wave function methods have offered a robust, systematically improvable means to study ground-state properties in quantum many-body systems. Theories like coupled cluster and their derivatives provide highly accurate approximations to the energy landscape at a reasonable computational cost. Analogues of such methods to study thermal properties, though highly desirable, have been lacking because evaluating thermal properties involve a trace over the entire Hilbert space, which is a formidable task. Besides, excited-state theories are generally not as well studied as ground-state ones. In this mini-review, we present an overview of a finite-temperature wave function formalism based on thermofield dynamics to overcome these difficulties. Thermofield dynamics allows us to map the equilibrium thermal density matrix to a pure state, i.e., a single wave function, albeit in an expanded Hilbert space. Ensemble averages become expectation values over this so-called thermal state. Around this thermal state, we have developed a procedure to generalize ground-state wave function theories to finite temperatures. As explicit examples, we highlight formulations of mean-field, configuration interaction, and coupled cluster theories for thermal properties of Fermions in the grand-canonical ensemble. To assess the quality of these approximations, we also show benchmark studies for the one-dimensional Hubbard model, while comparing against exact results. We will see that the thermal methods perform similarly to their ground-state counterparts, while merely adding a prefactor to the asymptotic computational cost. Furthermore, they also inherit all the properties, good or bad, from the ground-state methods, signifying the robustness of our formalism and the scope for future development.
We report a theoretical study of CdSe nanoplatelets aimed at identifying the main factors determining their photophysical properties. Using atomic configurations optimized with density functional theory calculations, we computed quasiparticle and exciton binding energies of nanoplatelets with two to seven monolayers. We employed many body perturbation theory at the GW level and solved the Bethe-Salpeter equation to obtain absorption spectra and excitonic properties. Our results, which agree well with recent experiments, were then used to design a model that allows us to disentangle the effects of quantum confinement, strain induced by passivating ligands, and dielectric environment on the electronic properties of nanoplatelets. We found that, for the model to accurately reproduce our first principle results, it is critical to account for surface stress and consider a finite potential barrier and energy-dependent effective masses when describing quantum confinement. Our findings call into question previous assumptions on the validity of an infinite barrier to describe carrier confinement in nanoplatelets, suggesting that it may be possible to optimize interfacial charge transfer and extraction by appropriately choosing passivating ligands. The model developed here is generalizable to core-shell platelets and enables the description of system sizes not yet directly treatable by first-principles calculations.
Several closely related ab initio thermal mean-field theories for fermions, both well-established and new ones, are compared with one another at the formalism level and numerically. The theories considered are Fermi–Dirac theory; thermal Hartree–Fock (HF) theory; two modifications of the thermal single-determinant and the first-order finite-temperature many-body perturbation theory based on a zero-temperature or thermal HF reference. Furthermore, thermal full-configuration-interaction theory is used as the benchmark.
Recent developments in ab initio nuclear theory demonstrate promising results in medium- to heavy-mass nuclei. A particular challenge for many of the many-body methodologies, however, is an accurate treatment of the electric-quadrupole, E2, strength associated with collectivity. The valence-space in-medium similarity renormalization group (VS-IMSRG) is a particularly powerful method for accessing medium- and high-mass nuclei but has been found to underpredict E2 strengths. The purpose of this work is to evaluate the isospin dependence of this underprediction. We perform a systematic comparison of VS-IMSRG calculations with available literature. We make use of isoscalar and isovector contributions to the E2 matrix elements to assess isoscalar and isovector contributions to the missing strength. It is found that the E2 strength is consistent throughout T z =|12|, T z =|1|, T z =|32|, and T z =2 pairs within the sd shell. Furthermore, no isovector contribution to the deficiency is identified. A comparison with toy-models and coupled-cluster calculations is used to discuss potential origins of the missing strength, which arises from missing many-particle, many-hole excitations out of the model space. The absence of any significant isovector contribution to the missing E2 strength indicates that the E2 strength discrepancy, and therefore any correction, is largely independent of the isospin of the nuclei in question.
As part of an effort to develop computational techniques for design of natural laminar flow fuselages, a computational study was made of the effect of forebody geometry on laminar boundary layer stability on axisymmetric body shapes. The effects of nose radius on the stability of the incompressible laminar boundary layer was computationally investigated using linear stability theory for body length Reynolds numbers representative of small and medium-sized airplanes. The steepness of the pressure gradient and the value of the minimum pressure (both functions of fineness ratio) govern the stability of laminar flow possible on an axisymmetric body at a given Reynolds number. It was found that to keep the laminar boundary layer stable for extended lengths, it is important to have a small nose radius. However, nose shapes with extremely small nose radii produce large pressure peaks at off-design angles of attack and can produce vortices which would adversely affect transition.
The presence of chalcogen vacancies in monolayer transition metal dichalcogenides (TMDs) leads to excitons with mixed localized-delocalized character and to reduced valley selectivity. Recent experimental advances in defect design in TMDs allow for a close examination of such mixed exciton states as a function of their degree of circular polarization under external magnetic fields, revealing strongly varying defect-induced magnetic properties. A theoretical understanding of these observations and their physical origins demands a predictive, structure-sensitive theory. In this work, we study the effect of chalcogen vacancies on the exciton magnetic properties in monolayer MoS 2 . Using many-body perturbation theory, we show how the complex excitonic picture associated with the presence of defects—with reduced valley and spin selectivity due to hybridized electron-hole transitions—leads to a structurally controllable exciton magnetic response. We find a variety of $\textit{g}$-factors with changing magnitudes and sign depending on the exciton energy and character. Further, our findings suggest a pathway to tune the nature of the excitons—and by that their magneto-optical properties—through defect architecture.
In the first part of this paper, we present a study of the symmetry energy (a s ) and its slope parameter (L) for nuclear matter in the framework of the Fermi liquid theory of Landau and Migdal. We derive an exact relation between a s and L, which involves the nucleon effective masses and three-particle Landau-Migdal parameters. We present simple estimates which suggest that there are two main mechanisms to explain the empirical values of L: The proton-neutron effective mass difference in isospin asymmetric matter and the ℓ = 0 moment of the isovector in-medium three-particle scattering amplitude. In the second part of this paper, we discuss the general structure of three-particle interactions in nuclear matter in the framework of the Fermi liquid theory. The connections to the Bethe-Brueckner-Goldstone theory and other approaches are also discussed. Finally, we show explicitly how the first few terms in the Faddeev series, together with medium induced three-particle interactions, emerge naturally in the Fermi liquid theory.
Here, in this work, we introduce a perturbative correction to the unitary coupled cluster method with single and double excitations (UCCSD) that incorporates the effects of missing triple excitations through fifth-order in many-body perturbation theory (MBPT). Referred to as UCCSD[T-5], this method is benchmarked alongside the previously developed UCCSD[T] to lend insight into the behavior of perturbative triples corrections relative to UCCSDT, which inherently provides an infinite-order treatment of triple excitations, as well as full configuration interaction (FCI). Two key findings emerge from this analysis. First, UCCSD[T] consistently yields ground-state energies in closest agreement with FCI, outperforming both UCCSD[T-5] and UCCSDT. Second, UCCSD[T-5] largely emulates the behavior of UCCSDT, suggesting a close relationship between fifth- and infinite-order triple-excitation contributions. Given the growing interest in UCC ansätze for quantum computing and the limitations of current quantum hardware, these results highlight the potential of classically computed perturbative corrections within UCC theory to capture triple-excitation effects without the additional quantum resources required by the UCCSDT ansatz.
We perform statistically rigorous uncertainty quantification (UQ) for chiral effective field theory (χ EFT) applied to infinite nuclear matter up to twice nuclear saturation density. The equation of state (EOS) is based on high-order many-body perturbation theory calculations with nucleon-nucleon and three-nucleon interactions up to fourth order in the χ EFT expansion. From these calculations our newly developed Bayesian machine-learning approach extracts the size and smoothness properties of the correlated EFT truncation error. Furthermore, we then propose a novel extension that uses multitask machine learning to reveal correlations between the EOS at different proton fractions. The inferred in-medium χ EFT breakdown scale in pure neutron matter and symmetric nuclear matter is consistent with that from free-space nucleon-nucleon scattering. These significant advances allow us to provide posterior distributions for the nuclear saturation point and propagate theoretical uncertainties to derived quantities: the pressure and incompressibility of symmetric nuclear matter, the nuclear symmetry energy, and its derivative. Our results, which are validated by statistical diagnostics, demonstrate that an understanding of truncation-error correlations between different densities and different observables is crucial for reliable UQ. The methods developed here are publicly available as annotated Jupyter notebooks.
We developed and used quantum simulation techniques to predict structural, electronic, optical and thermal properties affecting the scientific design of semiconductor nanostructures. We focused on the study of surfaces and interfaces at the nanoscale, in particular embedding and solvation of semiconductor nanoparticles, and provided a microscopic description of nanostructured materials in realistic environments, directly comparable with experimental conditions. We addressed three major problems: (i) the design of realistic structural models of integrated nanostructures, e.g. quantum dots embedded in solid matrices and nanostructures immersed in fluids; (ii) the development of efficient and accurate calculations to describe the excited state properties of nanostructures, beyond standard Density Functional Theory, so as to understand and predict electronic properties involved in renewable energy applications, for example solar conversion devices. In particular we developed and used methods based on a many body perturbation theory. (iii) the development of calculations of heat transport properties of nanostructured materials, for thermoelectric applications, including predictive calculations of thermal and electrical conductivity and of Seebeck coefficients. Specifically, we developed and used methods based on the Boltzmann transport equation and molecular dynamics.
Quantum embedding theories are promising approaches to investigate strongly correlated electronic states of active regions of large-scale molecular or condensed systems. Notable examples are spin defects in semiconductors and insulators. We present a detailed derivation of a quantum embedding theory recently introduced, which is based on the definition of effective Hamiltonians. The effect of the environment on a chosen active space is accounted for through screened Coulomb interactions evaluated using density functional theory. Importantly, the random phase approximation is not required, and the evaluation of virtual electronic orbitals is circumvented with algorithms previously developed in the context of calculations based on many-body perturbation theory. In addition, we generalize the quantum embedding theory to active spaces composed of orbitals that are not eigenstates of Kohn–Sham Hamiltonians. Finally, we report results for spin defects in semiconductors.
In this work we study ordinary muon capture (OMC) on 24 Mg from a first-principles perspective. Starting from a particular two- and three-nucleon interaction derived from chiral effective field theory, we use the valence-space in-medium similarity renormalization group (VS-IMSRG) framework to construct effective Hamiltonians and muon-capture operators, which nonperturbatively account for many-body physics outside the valence space. Here, the obtained nuclear matrix elements are compared against those from the phenomenological shell model. The impact of including the correlations from the nuclear shell model (NSM) as well as including the induced twobody part is studied in detail. Furthermore, the effects of realistic bound-muon wave function on the operators is studied. Finally, predictions for capture rates to the lowest excited states in 24 Na are given and compared with available data. It is found that the spectroscopic properties of 24 Mg and its OMC daughter 24 Na are fairly well described by both the NSM and VS-IMSRG, and that the effect of the hadronic two-body currents significantly reduces the OMC rates. Both models have some difficulties in matching the measured OMC rates, especially for the 2 + final states. This calls for further studies in other light nuclei with available OMC data.
Herein we demonstrate the prowess of the 3D electron diffraction approach by unveiling the structure of terrylene, the third member in the series of peri-condensed naphthalene analogues, which has eluded structure determination for 65 years. The structure was determined by direct methods using electron diffraction data and corroborated by dispersion-inclusive density functional theory optimizations. Terrylene crystalizes in the monoclinic space group P2 1 /a, arranging in a sandwich-herringbone packing motif, similar to analogous compounds. Having solved the crystal structure, we use many-body perturbation theory to evaluate the excited-state properties of terrylene in the solid-state. We find that terrylene is a promising candidate for intermolecular singlet fission, comparable to tetracene and rubrene.
No known body in the Solar System has a spectrum redder than that of object 5145 Pholus. We use Hapke scattering theory and optical constants measured in this laboratory to examine the ability of mixtures of a number of organic solids and ices to reproduce the observed spectrum and phase variation. The primary materials considered are poly-HCN, kerogen, Murchison organic extract, Titan tholin, ice tholin, and water ice. In a computer grid search of over 10 million models, we find an intraparticle mixture of 15% Titan tholin, 10% poly-HCN, and 75% water ice with 10-micrometers particles to provide an excellent fit. Replacing water ice with ammonia ice improves the fits significantly while using a pure hydrocarbon tholin, Tholin alpha, instead of Titan tholin makes only modest improvements. All acceptable fits require Titan tholin or some comparable material to provide the steep slope in the visible, and poly-HCN or some comparable material to provide strong absorption in the near-infrared. A pure Titan tholin surface with 16-micrometers particles, as well as all acceptable Pholus models, fit the present spectrophotometric data for the transplutonian object 1992 QB(sub 1). The feasibility of gas-phase chemistry to generate material like Titan tholin on such small objects is examined. An irradiated transient atmosphere arising from sublimating ices may generate at most a few centimeters of tholin over the lifetime of the Solar System, but this is insignificant compared to the expected lag deposit of primordial contaminants left behind by the sublimating ice. Irradiation of subsurface N2/CH4 or NH3/CH4 ice by cosmic rays may generate approximately 20 cm of tholin in the upper 10 m of regolith in the same time scale but the identity of this tholin to its gas-phase equivalent has not been demonstrated.