Search NASA⌕ Search

SEARCH · Search NASA

Results for “transit method”

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 19 records

Slater transition methods for core-level electron binding energies

Methods for computing core-level ionization energies using self-consistent field (SCF) calculations are evaluated and benchmarked. These include a “full core hole” (or “ΔSCF”) approach that fully accounts for orbital relaxation upon ionization, but also methods based on Slater’s transition concept in which the binding energy is estimated from an orbital energy level that is obtained from a fractional-occupancy SCF calculation. A generalization that uses two different fractional-occupancy SCF calculations is also considered. The best of the Slater-type methods afford mean errors of 0.3–0.4 eV with respect to experiment for a dataset of K-shell ionization energies, a level of accuracy that is competitive with more expensive many-body techniques. An empirical shifting procedure with one adjustable parameter reduces the average error below 0.2 eV. Here, this shifted Slater transition method is a simple and practical way to compute core-level binding energies using only initial-state Kohn–Sham eigenvalues. It requires no more computational effort than ΔSCF and may be especially useful for simulating transient x-ray experiments where core-level spectroscopy is used to probe an excited electronic state, for which the ΔSCF approach requires a tedious state-by-state calculation of the spectrum. As an example, we use Slater-type methods to model x-ray emission spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

The method of transition boundary for the solution of diffraction problems

A new method, which enables one to solve some diffraction problems, is put forth. The technique is based on a relation between the diffracted and scattered geometric optics waves at the transition boundaries. The concept of initial waves is also used instead of the incident field. The general algorithm of the method is outlined and applied to the diffraction problems of waves by perfect electric conductor half-screen, resistive half-plane and impedance sheet junction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The Opus radiative opacity code, a new implementation of the super-transition array method

This report describes the development of our new radiative opacity code Opus, submitted as part of a 2020 Level 2 milestone requirement. This code aims to incorporate the capabilities in existing but separate opacity codes previously developed at LLNL (Opal, Topaz, and Tycho). The first phase in the development, which was to implement the super-transition (STA) method of Bar-Shalom and co-workers is now completed. In this work, we describe the underlying physics, implementation, usage, and regression testing as well as providing numerical examples for the opacity tables of boron, chromium, iron, and nickel. The latter are benchmarked against the LLNL legacy STA code Tycho, with an average absolute relative difference of ~1.5 and 0.5% for Planck and Rosseland opacity means, respectively. The cases where differences are larger are well understood in terms of numerical instabilities near discontinuities of the average degree of ionization inherent to the assumed opacity/equation-of-state model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

On the potentially transformative role of auxiliary-field quantum Monte Carlo in quantum chemistry: A highly accurate method for transition metals and beyond

Approximate solutions to the ab initio electronic structure problem have been a focus of theoretical and computational chemistry research for much of the past century, with the goal of predicting relevant energy differences to within “chemical accuracy” (1 kcal/mol). For small organic molecules, or in general, for weakly correlated main group chemistry, a hierarchy of single-reference wave function methods has been rigorously established, spanning perturbation theory and the coupled cluster (CC) formalism. For these systems, CC with singles, doubles, and perturbative triples is known to achieve chemical accuracy, albeit at O(N7) computational cost. In addition, a hierarchy of density functional approximations of increasing formal sophistication, known as Jacob’s ladder, has been shown to systematically reduce average errors over large datasets representing weakly correlated chemistry. However, the accuracy of such computational models is less clear in the increasingly important frontiers of chemical space including transition metals and f-block compounds, in which strong correlation can play an important role in reactivity. A stochastic method, phaseless auxiliary-field quantum Monte Carlo (ph-AFQMC), has been shown to be capable of producing chemically accurate predictions even for challenging molecular systems beyond the main group, with relatively low O(N3 − N4) cost and near-perfect parallel efficiency. Herein, we present our perspectives on the past, present, and future of the ph-AFQMC method. We focus on its potential in transition metal quantum chemistry to be a highly accurate, systematically improvable method that can reliably probe strongly correlated systems in biology and chemical catalysis and provide reference thermochemical values (for future development of density functionals or interatomic potentials) when experiments are either noisy or absent. Finally, we discuss the present limitations of the method and where we expect near-term development to be most fruitful.

Chemistry↗

α-Weighted Transition Rate Matrix Method

A more direct approach in calculating the Transition Rate Matrix Method operators and the significance of this method are presented. Consequences of using α-spectrum as the Transition Rate Matrix Method weighting function instead of the κ-spectrum are investigated. A prompt supercritical infinite medium multigroup system is considered. Results on continuous problem using OpenMC will be provided in the full paper.

Variansyah, I.↗

α-weighted Transition Rate Matrix Method

Relative effectiveness of an α-weighted Transition Rate Matrix Method (TRMM) as compared to the normally used κ- weighted, for solving reactor transient problem is investigated. An improved Monte Carlo κ-α iteration with delayed neutron is implemented to the OpenMC to generate the fundamental α-eigenfunction to be used as the weighting function in calculating the TRMM forward and adjoint operators. A prompt supercritical transient problem and a deep subcritical transient problem of a two-region slab system are devised to compare the relative performance of the α- and κ-weighted TRMM.

Variansyah, Ilham↗

Giant Outer Transiting Exoplanet Mass (GOT ‘EM) Survey. I. Confirmation of an Eccentric, Cool Jupiter with an Interior Earth-sized Planet Orbiting Kepler-1514

Despite the severe bias of the transit method of exoplanet discovery toward short orbital periods, a modest sample of transiting exoplanets with orbital periods greater than 100 days is known. Long-term radial velocity (RV) surveys are pivotal to confirming these signals and generating a set of planetary masses and densities for planets receiving moderate to low irradiation from their host stars. Here we conduct RV observations of Kepler-1514 from the Keck I telescope using the High Resolution Echelle Spectrometer. From these data, we measure the mass of the statistically validated giant (1.108 ± 0.023 R {sub J}) exoplanet Kepler-1514 b with a 218-day orbital period as 5.28 ± 0.22 M {sub J}. The bulk density of this cool (∼390 K) giant planet is 4.82{sub −0.25}{sup +0.26}g cm{sup −3}, consistent with a core supported by electron degeneracy pressure. We also infer an orbital eccentricity of 0.401{sub −0.014}{sup +0.013} from the RV and transit observations, which is consistent with planet–planet scattering and disk cavity migration models. The Kepler-1514 system contains an Earth-size, Kepler Object of Interest on a 10.5-day orbit that we statistically validate against false-positive scenarios, including those involving a neighboring star. The combination of the brightness (V = 11.8) of the host star and the long period, low irradiation, and high density of Kepler-1514 b places this system among a rare group of known exoplanetary systems and as one that is amenable to continued study.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

On the Utility of Transmission Color Analysis i: Differentiating Super-Earths and Sub-Neptunes

The majority of exoplanets found to date have been discovered via the transit method, and transmission spectroscopy represents the primary method of studying these distant worlds. Currently, in-depth atmospheric characterization of transiting exoplanets entails the use of spectrographs on large telescopes, requiring significant observing time to study each planet. Previous studies have demonstrated trends for solar system worlds using color–color photometry of reflectance spectra, as well as trends within transmission spectra for hot Jupiters. Building on these concepts, we have investigated the use of transmission color photometric analysis for efficient, coarse categorization of exoplanets and for assessing the nature of these worlds, with a focus on resolving the bulk composition degeneracy to aid in discriminating super-Earths and sub-Neptunes. We present our methodology and first results, including spectrum models, model comparison frameworks, and wave band selection criteria. We present our results for different transmission “color” metrics, filter selection methods, and numbers of filters. Assuming noise-free spectra of isothermal atmospheres in chemical equilibrium, with our pipeline, we are able to constrain atmospheric mean molecular weight in order to distinguish between super-Earth and sub-Neptune atmospheres with >90% overall accuracy using specific low-resolution filter combinations, . We also found that increasing the number of filters does not substantially impact this performance. This method could allow for broad characterization of large numbers of planets much more efficiently than current methods permit, enabling population- and system-level studies. Additionally, data collected via this method could inform follow-up observing time by large telescopes for more detailed studies of worlds of interest.

47 OTHER INSTRUMENTATION↗

Toward Complete Characterization: Prospects for Directly Imaging Transiting Exoplanets

High-contrast direct imaging of exoplanets can provide many important observables, including measurements of the orbit, spectra that probe the lower layers of the atmosphere, and phase variations of the planet, but cannot directly measure planet radius or mass. Our future understanding of directly imaged exoplanets will therefore rely on extrapolated models of planetary atmospheres and bulk composition, which need robust calibration. We estimate the population of extrasolar planets that could serve as calibrators for these models. Critically, this population of “standard planets” must be accessible to both direct imaging and the transit method, allowing for radius measurement. We show that the search volume of a direct imaging mission eventually overcomes the transit probability falloff with semimajor axis, so that as long as cold planets are not exceedingly rare, the population of transiting planets and directly imageable planets overlaps. Using current extrapolations of Kepler occurrence rates, we estimate that ∼8 standard planets could be characterized shortward of 800 nm with an ambitious future direct imaging mission like LUVOIR-A and several dozen could be detected at the V band. We show the design space that would expand the sample size and discuss the extent to which ground- and space-based surveys could detect this small but crucial population of planets.

79 ASTRONOMY AND ASTROPHYSICS↗

Limb Darkening and Planetary Transits: Testing Center-to-limb Intensity Variations and Limb-darkening Directly from Model Stellar Atmospheres

The transit method, employed by Microvariability and Oscillation of Stars ( MOST ), Kepler , and various ground-based surveys has enabled the characterization of extrasolar planets to unprecedented precision. These results are precise enough to begin to measure planet atmosphere composition, planetary oblateness, starspots, and other phenomena at the level of a few hundred parts per million. However, these results depend on our understanding of stellar limb darkening, that is, the intensity distribution across the stellar disk that is sequentially blocked as the planet transits. Typically, stellar limb darkening is assumed to be a simple parameterization with two coefficients that are derived from stellar atmosphere models or fit directly. In this work, we revisit this assumption and compute synthetic planetary-transit light curves directly from model stellar atmosphere center-to-limb intensity variations (CLIVs) using the plane-parallel Atlas and spherically symmetric SAtlas codes. We compare these light curves to those constructed using best-fit limb-darkening parameterizations. We find that adopting parametric stellar limb-darkening laws leads to systematic differences from the more geometrically realistic model stellar atmosphere CLIV of about 50–100 ppm at the transit center and up to 300 ppm at ingress/egress. While these errors are small, they are systematic, and they appear to limit the precision necessary to measure secondary effects. Our results may also have a significant impact on transit spectra.

79 ASTRONOMY AND ASTROPHYSICS↗

First critical repressive H3K27me3 marks in embryonic stem cells identified using designed protein inhibitor

The polycomb repressive complex 2 (PRC2) histone methyltransferase plays a central role in epigenetic regulation in development and in cancer, and hence to interrogate its role in a specific developmental transition, methods are needed for disrupting function of the complex with high temporal and spatial precision. The catalytic and substrate recognition functions of PRC2 are coupled by binding of the N-terminal helix of the Ezh2 methylase to an extended groove on the EED trimethyl lysine binding subunit. Disrupting PRC2 function can in principle be achieved by blocking this single interaction, but there are few approaches for blocking specific protein–protein interactions in living cells and organisms. Here, in this work, we describe the computational design of proteins that bind to the EZH2 interaction site on EED with subnanomolar affinity in vitro and form tight and specific complexes with EED in living cells. Induction of the EED binding proteins abolishes H3K27 methylation in human embryonic stem cells (hESCs) and at all but the earliest stage blocks self-renewal, pinpointing the first critical repressive H3K27me3 marks in development.

60 APPLIED LIFE SCIENCES↗

A Jacobian-free pseudo-arclength continuation method for phase transitions in inhomogeneous thermodynamic systems

Developing phase diagrams for inhomogeneous systems in thermodynamics is difficult, in part, due to the large phase space and the possibility of unstable and metastable solutions arising from first-order phase transitions. Pseudo-arclength continuation (PAC) is a method that allows one to trace out stable and unstable solutions of nonlinear systems. Typically, PAC utilizes the Jacobian in order to implement Newton (or quasi-Newton) steps. In this work, we present a Jacobian-free PAC method that is amenable to the usual workflows in inhomogeneous thermodynamics. We demonstrate our method in systems that have first-order phase transitions, including a novel example of polyelectrolyte complex coacervation in confinement, where multiple surface phase transitions occur and can overlap with one another.

Chemistry↗

Universal Passive Synchronization Method for Grid-Forming Inverters Without Mode Transition

Power systems are transforming with increasing levels of inverter-based resources (IBRs). This transformation requires critical roles of grid-forming (GFM) inverters replacing synchronous generators for bulk power system stabilization and ancillary services, also allowing flexible power system operation, such as microgrid that is operated by multiple GFM IBRs to achieve system resilience against contingencies. To realize the resilient power systems allowing flexible in-and-out operation of GFM IBRs potentially programmed with different primary controls, a synchronization method universally applicable, i.e., independent of control types, would be beneficial to ease the integration process, but it has not been actively studied. To fill the gap, this paper proposes a universal synchronization method that achieves a passive synchronization to enable a smooth transition in a grid with off-nominal system parameters, i.e., voltage and frequency. The logic proposed requires no modification on the primary control, thus applicable to any type of GFMs with a voltage reference input. To validate the concept, a simulation of an IEEE 13-bus benchmark system modified with 3 GFM inverters is presented. It simulates an inverter-driven black start scenario in which GFM inverters autonomously turn on and connect to the grid under heavy loading, using the synchronization logic. The case study demonstrates that GFM inverters can tune their voltage reference to smoothly synchronize without severe transients, and contribute to a seamless black start of the grid under unbalanced load conditions. Two GFM methods - Droop and dispatchable virtual oscillator control - are used for the demo to validate feasibility and interoperability of the passive synchronization.

black start↗

Wannier–Koopmans method calculations for transition metal oxide band gaps

The widely used density functional theory (DFT) has a major drawback of underestimating the band gaps of materials. Wannier–Koopmans method (WKM) was recently developed for band gap calculations with accuracy on a par with more complicated methods. WKM has been tested for main group covalent semiconductors, alkali halides, 2D materials, and organic crystals. Here we apply the WKM to another interesting type of material system: the transition metal (TM) oxides. TM oxides can be classified as either with d 0 or d 10 closed shell occupancy or partially occupied open shell configuration, and the latter is known to be strongly correlated Mott insulators. We found that, while WKM provides adequate band gaps for the d0 and d10 TM oxides, it fails to provide correct band gaps for the group with partially occupied d states. This issue is also found in other mean-field approaches like the GW calculations. We believe that the problem comes from a strong interaction between the occupied and unoccupied d-state Wannier functions in a partially occupied d-state system. We also found that, for pseudopotential calculations including deep core levels, it is necessary to remove the electron densities of these deep core levels in the Hartree and exchange–correlation energy functional when calculating the WKM correction parameters for the d-state Wannier functions.

36 MATERIALS SCIENCE↗

Generator coordinate method for transition-state dynamics in nuclear fission

Since its beginnings, fission theory has assumed that low-energy induced fission takes place through transition-state channels at the barrier tops. Nevertheless, up to now there is no microscopic theory applicable to those conditions. We suggest that modern reaction theory is suitable for this purpose, and propose a methodology based on a configuration-interaction framework using the generator coordinate method (GCM). Simple reaction-theoretic models are constructed with the Gaussian overlap approximation to parametrize both the dynamics within the channels and their incoherent couplings to states outside the barrier. The physical characteristics of the channels examined here are their effective bandwidths and the quality of the coupling to compound-nucleus states as measured by the transmission factor T. We also investigate the spacing of GCM states with respect to their degree of overlap. We find that a rather coarse mesh provides an acceptable accuracy for estimating the bandwidths and transmission factors. The common numerical stability problem in using the GCM is avoided due to the choice of meshes and the finite bandwidths of the channels. Here, the bandwidths of the channels are largely controlled by the zero-point energy with respect to the collective coordinate in the GCM configurations.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗