SEARCH · Engineering Papers

Results for “topological derivative”

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.

SourceResearch connectionFrom yearThrough year

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

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Search for non-resonant production of semi-visible jets using Run 2 data in ATLAS

Semi-visible jets, with a significant contribution to the event's missing transverse momentum, can arise in strongly interacting dark sectors. This results in an event topology where one of the jets can be aligned with the direction of the missing transverse momentum. The first search for semi-visible jets produced via a t-channel mediator exchange is presented. The analysis uses proton-proton collisions with an integrated luminosity of 139 fb —1 and a centre-of-mass energy of 13 TeV, collected with the ATLAS detector during the Run 2 of the LHC. No excess over Standard Model predictions is observed. Assuming a coupling strength of unity between the mediator, a Standard Model quark and a dark quark, mediator masses up to 2.7 TeV are excluded at the 95% confidence level. Upper limits on the coupling strength are also derived.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2020

Electronic correlations and flattened band in magnetic Weyl semimetal candidate Co 3 Sn 2 S 2

The interplay between electronic correlations and topological protection may offer a rich avenue for discovering emergent quantum phenomena in condensed matter. However, electronic correlations have so far been little investigated in Weyl semimetals (WSMs) by experiments. Here, we report a combined optical spectroscopy and theoretical calculation study on the strength and effect of electronic correlations in a magnet Co 3 Sn 2 S 2 . The electronic kinetic energy estimated from our optical data is about half of that obtained from single-particle ab initio calculations in the ferromagnetic ground state, which indicates intermediate-strength electronic correlations in this system. Furthermore, comparing the energy and side-slope ratios between the interband-transition peaks at high energies in the experimental and single-particle-calculation-derived optical conductivity spectra with the bandwidth-renormalization factors obtained by many-body calculations enables us to estimate the Coulomb-interaction strength (U ~ 4 eV) in Co 3 Sn 2 S 2 . Besides, a sharp experimental optical conductivity peak at low energy, which is absent in the single-particle-calculation-derived spectrum but is consistent with the optical conductivity peaks obtained by many-body calculations with U ~ 4 eV, indicates that an electronic band connecting the two Weyl cones is flattened by electronic correlations and emerges near the Fermi energy in Co 3 Sn 2 S 2 . Our work paves the way for exploring flat-band-generated quantum phenomena in WSMs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

SMARTER Rules-Based Distributed Deconfliction of ADMS Applications

A conceptual numerical methodology derived from Grid Architecture principles is introduced for deconflicting setpoints issued by multiple advanced distribution management system applications. The methodology applies technical, economic, environmental, and social rules to eliminate non-viable combinations. The concept of temporal equipment controls budgets is introduced to preserve the health of physical assets and avoid equipment damage through repeated controls cycling. The rules are combined with a multi-criteria decision-making framework to select a near-optimal set of deconflicted setpoints using a set of qualitative and quantitative decision criteria selected by the distribution system operator. Numerical results are demonstrated on the IEEE 123-bus test feeder for three competing applications. Three alternative distributed schemes are used to decompose the problem: by topological area, by phase, and fully decentralized. The fully decentralized implementation is shown to yield near-optimal deconfliction results with significantly reduced computational time.

Anderson, Alexander A.↗

DOE OSTI2025

Uncertainty Visualization of Critical Points of 2D Scalar Fields for Parametric and Nonparametric Probabilistic Models

This paper presents a novel end-to-end framework for closed-form computation and visualization of critical point uncertainty in 2D uncertain scalar fields. Critical points are fundamental topological descriptors used in the visualization and analysis of scalar fields. The uncertainty inherent in data (e.g., observational and experimental data, approximations in simulations, and compression), however, creates uncertainty regarding critical point positions. Uncertainty in critical point positions, therefore, cannot be ignored, given their impact on downstream data analysis tasks. Here, in this work, we study uncertainty in critical points as a function of uncertainty in data modeled with probability distributions. Although Monte Carlo (MC) sampling techniques have been used in prior studies to quantify critical point uncertainty, they are often expensive and are infrequently used in production-quality visualization software. We, therefore, propose a new end-to-end framework to address these challenges that comprises a threefold contribution. First, we derive the critical point uncertainty in closed form, which is more accurate and efficient than the conventional MC sampling methods. Specifically, we provide the closed-form and semianalytical (a mix of closed-form and MC methods) solutions for parametric (e.g., uniform, Epanechnikov) and nonparametric models (e.g., histograms) with finite support. Second, we accelerate critical point probability computations using a parallel implementation with the VTK-m library, which is platform portable. Finally, we demonstrate the integration of our implementation with the ParaView software system to demonstrate near-real-time results for real datasets.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1995

Topological structure of the magnetotail as a function of interplanetary magnetic field direction

Magnetic reconnection between the interplanetary magnetic field (IMF) and the geomagnetic field is thought to play a major role in the transfer of solar wind momentum and energy to the magnetosphere. As the angle between the IMF and the geomagnetic field is changed at the bow of the magnetosphere, the topological record of the location of the reconnection region should be recorded in the magnetosheath and on the magnetopause along the flanks of the tail, because the super fast flow freezes strong magnetic gradients formed in the bow reconnection regions into the plasma downstream. In this report, we present results from a three-dimensional, magnetohydrodynamic (MHD), global numerical simulation code for the location of the separatrix between unconnected IMF magnetosheath field lines and reconnected field lines which penetrate the magnetopause and connect to the polar ionosphere. The angle between the IMF direction and the line where the separatrix crosses the magnetopause is shown to be a sensitive function of the IMF clock angle. We also explain how this behavior can be used to derive an approximate relation for the dependence of the cross-polar voltage on the IMF clock angle. We conclude with a note of caution concerning the importance of physical boundary conditions in magnetoplasma simulations.

Fedder, J. A.↗

DOE OSTI2025

A plant-specific cytochrome b 5 –like protein is essential for phytosterol biosynthesis

Sterols are essential isoprenoid derivatives that contribute to membrane structure and function. In plants, they also serve as precursors to phytohormones and specialized metabolites important for development, defense, and health. Although the sterol biosynthetic pathway is considered well-characterized, we report the discovery of a plant-specific cytochrome b 5 –like protein, CB5LP, as a critical component of phytosterol biosynthesis. Loss of CB5LP in Arabidopsis causes embryonic defects, seedling lethality, and accumulation of 14α-methyl-sterols, with reduced levels of sitosterol and stigmasterol—indicating a defect in sterol 14α-demethylation. TurboID-based proximity labeling and in vitro assays show that CB5LP physically and functionally interacts with CYP51, a cytochrome P450 enzyme catalyzing this demethylation step. Unlike canonical cytochrome b 5 proteins, CB5LP has a reversed topology and is exclusive to plants, acting as an evolutionarily distinct electron donor. This discovery reveals an uncharacterized redox partnership essential for sterol biosynthesis and highlights a promising target for the development of selective herbicide.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS2012

Linked-List-Based Multibody Dynamics (MBDyn) Engine

This new release of MBDyn is a software engine that calculates the dynamics states of kinematic, rigid, or flexible multibody systems. An MBDyn multibody system may consist of multiple groups of articulated chains, trees, or closed-loop topologies. Transient topologies are handled through conservation of energy and momentum. The solution for rigid-body systems is exact, and several configurable levels of nonlinear term fidelity are available for flexible dynamics systems. The algorithms have been optimized for efficiency and can be used for both non-real-time (NRT) and real-time (RT) simulations. Interfaces are currently compatible with NASA's Trick Simulation Environment. This new release represents a significant advance in capability and ease of use. The two most significant new additions are an application programming interface (API) that clarifies and simplifies use of MBDyn, and a link-list infrastructure that allows a single MBDyn instance to propagate an arbitrary number of interacting groups of multibody top ologies. MBDyn calculates state and state derivative vectors for integration using an external integration routine. A Trickcompatible interface is provided for initialization, data logging, integration, and input/output.

MacLean, John↗

NASA NTRS2015

The Time-Dependent Chemistry of Cometary Debris in the Solar Corona

Recent improvements in solar observations have greatly progressed the study of sungrazing comets. They can now be imaged along the entirety of their perihelion passage through the solar atmosphere, revealing details of their composition and structure not measurable through previous observations in the less volatile region of the orbit further from the solar surface. Such comets are also unique probes of the solar atmosphere. The debris deposited by sungrazers is rapidly ionized and subsequently influenced by the ambient magnetic field. Measuring the spectral signature of the deposited material highlights the topology of the magnetic field and can reveal plasma parameters such as the electron temperature and density. Recovering these variables from the observable data requires a model of the interaction of the cometary species with the atmosphere through which they pass. The present paper offers such a model by considering the time-dependent chemistry of sublimated cometary species as they interact with the solar radiation field and coronal plasma. We expand on a previous simplified model by considering the fully time-dependent solutions of the emitting species' densities. To compare with observations, we consider a spherically symmetric expansion of the sublimated material into the corona and convert the time-dependent ion densities to radial profiles. Using emissivities from the CHIANTI database and plasma parameters derived from a magnetohydrodynamic simulation leads to a spatially dependent emission spectrum that can be directly compared with observations. We find our simulated spectra to be consistent with observation.

sun: corona↗

DOE OSTI2021

Magnetic photocurrents in multifold Weyl fermions

Here, we examine the magneto-optical response of chiral multifold fermions. Specifically, we show that they are ideal candidates for observing the helical magnetic effect (HME) previously predicted for simple Weyl fermions. Unlike Weyl fermions, the HME is present in multifold fermions even in the simplest case where the low-energy dispersion is linear and spherically symmetric. In this ideal case, we derive an analytical expression for the HME and find it is proportional to the circular photogalvanic effect; for realistic parameters and accounting for the geometry of the setup, the HME photocurrent could be roughly the same order of magnitude as the circular photogalvanic effect observed in multifold fermions. Additional nonlinear and symmetry-breaking terms will ruin the quantization but not hurt the observation of the HME.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS2006

CAPRI: A Geometric Foundation for Computational Analysis and Design

CAPRI is a software building tool-kit that refers to two ideas; (1) A simplified, object-oriented, hierarchical view of a solid part integrating both geometry and topology definitions, and (2) programming access to this part or assembly and any attached data. A complete definition of the geometry and application programming interface can be found in the document CAPRI: Computational Analysis PRogramming Interface appended to this report. In summary the interface is subdivided into the following functional components: 1. Utility routines -- These routines include the initialization of CAPRI, loading CAD parts and querying the operational status as well as closing the system down. 2. Geometry data-base queries -- This group of functions allow all top level applications to figure out and get detailed information on any geometric component in the Volume definition. 3. Point queries -- These calls allow grid generators, or solvers doing node adaptation, to snap points directly onto geometric entities. 4. Calculated or geometrically derived queries -- These entry points calculate data from the geometry to aid in grid generation. 5. Boundary data routines -- This part of CAPRI allows general data to be attached to Boundaries so that the boundary conditions can be specified and stored within CAPRI s data-base. 6. Tag based routines -- This part of the API allows the specification of properties associated with either the Volume (material properties) or Boundary (surface properties) entities. 7. Geometry based interpolation routines -- This part of the API facilitates Multi-disciplinary coupling and allows zooming through Boundary Attachments. 8. Geometric creation and manipulation -- These calls facilitate constructing simple solid entities and perform the Boolean solid operations. Geometry constructed in this manner has the advantage that if the data is kept consistent with the CAD package, therefore a new design can be incorporated directly and is manufacturable. 9. Master Model access This addition to the API allows for the querying of the parameters and dimensions of the model. The feature tree is also exposed so it is easy to see where the parameters are applied. Calls exist to allow for the modification of the parameters and the suppression/unsuppression of nodes in the tree. Part regeneration is performed by a single API call and a new part becomes available within CAPRI (if the regeneration was successful). This is described in a separate document. Components 1-7 are considered the CAPRI base level reader.

Haimes, Robert↗

DOE OSTI2020

The changing mouse embryo transcriptome at whole tissue and single-cell resolution

During mammalian embryogenesis, differential gene expression gradually builds the identity and complexity of each tissue and organ system 1 . Here we systematically quantified mouse polyA-RNA from day 10.5 of embryonic development to birth, sampling 17 tissues and organs. The resulting developmental transcriptome is globally structured by dynamic cytodifferentiation, body-axis and cell-proliferation gene sets that were further characterized by the transcription factor motif codes of their promoters. We decomposed the tissue-level transcriptome using single-cell RNA-seq (sequencing of RNA reverse transcribed into cDNA) and found that neurogenesis and haematopoiesis dominate at both the gene and cellular levels, jointly accounting for one-third of differential gene expression and more than 40% of identified cell types. By integrating promoter sequence motifs with companion ENCODE epigenomic profiles, we identified a prominent promoter de-repression mechanism in neuronal expression clusters that was attributable to known and novel repressors. Focusing on the developing limb, single-cell RNA data identified 25 candidate cell types that included progenitor and differentiating states with computationally inferred lineage relationships. We extracted cell-type transcription factor networks and complementary sets of candidate enhancer elements by using single-cell RNA-seq to decompose integrative cis -element (IDEAS) models that were derived from whole-tissue epigenome chromatin data. These ENCODE reference data, computed network components and IDEAS chromatin segmentations are companion resources to the matching epigenomic developmental matrix, and are available for researchers to further mine and integrate.

59 BASIC BIOLOGICAL SCIENCES↗

DOE OSTI2022

Elucidating correlated defects in metal organic frameworks using theory-guided inelastic neutron scattering spectroscopy

Metal organic frameworks (MOFs) that incorporate metal oxide cluster nodes, exemplified by UiO-66, have been widely studied, especially in terms of their deviations from the ideal, defect-free crystalline structures. Although defects such as missing linkers, missing nodes, and the presence of adventitious synthesis-derived node ligands (such as acetates and formates) have been proposed, their exact structures remain unknown. Previously, it was demonstrated that defects are correlated and span multiple unit cells. The highly specialized techniques used in these studies are not easily applicable to other MOFs. Thus, there is a need to develop new experimental and computational approaches to understand the structure and properties of defects in a wider variety of MOFs. Here, we show how low-frequency phonon modes measured by inelastic neutron scattering (INS) spectroscopy can be combined with density functional theory (DFT) simulations to provide unprecedented insights into the defect structure of UiO-66. We are able to identify and assign peaks in the fingerprint region (<100 cm –1 ) which correspond to phonon modes only present in certain defective topologies. Specifically, this analysis suggests that our sample of UiO-66 consists of predominantly defect-free fcu regions with smaller domains corresponding to a defective bcu topology with 4 and 2 acetate ligands bound to the Zr 6 O 8 nodes. Importantly, the INS/DFT approach provides detailed structural insights (e.g., relative positions and numbers of acetate ligands) that are not accessible with microscopy-based techniques. Furthermore, the quantitative agreement between DFT simulations and the experimental INS spectrum combined with the relative simplicity of sample preparation, suggests that this methodology may become part of the standard and preferred protocol for the characterization of MOFs, and, in particular, for elucidating the structure defects in these materials.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Near-axis expansion of stellarator equilibrium at arbitrary order in the distance to the axis

A direct construction of equilibrium magnetic fields with toroidal topology at arbitrary order in the distance from the magnetic axis is carried out, yielding an analytical framework able to explore the landscape of possible magnetic flux surfaces in the vicinity of the axis. This framework can provide meaningful analytical insight into the character of high-aspect-ratio stellarator shapes, such as the dependence of the rotational transform and the plasma beta limit on geometrical properties of the resulting flux surfaces. The approach developed here is based on an asymptotic expansion on the inverse aspect ratio of the ideal magnetohydrodynamics equation. The analysis is simplified by using an orthogonal coordinate system relative to the Frenet–Serret frame at the magnetic axis. The magnetic field vector, the toroidal magnetic flux, the current density, the field line label and the rotational transform are derived at arbitrary order in the expansion parameter. Moreover, a comparison with a near-axis expansion formalism employing an inverse coordinate method based on Boozer coordinates (the so-called Garren–Boozer construction) is made, where both methods are shown to agree at lowest order. Finally, as a practical example, a numerical solution using a W7-X equilibrium is presented, and a comparison between the lowest-order solution and the W7-X magnetic field is performed.

Physics↗

DOE OSTI2021

Solitons in lattice field theories via tight-binding supersymmetry

Reflectionless potentials play an important role in constructing exact solutions to classical dynamical systems (such as the Korteweg-de Vries equation), non-perturbative solutions of various large-N field theories (such as the Gross-Neveu model), and closely related solitonic solutions to the Bogoliubov-de Gennes equations in the theory of superconductivity. These solutions rely on the inverse scattering method, which reduces these seemingly unrelated problems to identifying reflectionless potentials of an auxiliary one-dimensional quantum scattering problem. There are several ways of constructing these potentials, one of which is quantum mechanical supersymmetry (SUSY). In this paper, motivated by recent experimental platforms, we generalize this framework to develop a theory of lattice solitons. We first briefly review the classical inverse scattering method in the continuum limit, focusing on the Korteweg-de Vries (KdV) equation and SU(N) Gross-Neveu model in the large N limit. We then generalize this methodology to lattice versions of interacting field theories. Our analysis hinges on the use of trace identities, which are relations connecting the potential of an equation of motion to the scattering data. For a discrete Schrödinger operator, such trace identities had been known as far back as Toda; however, we derive a new set of identities for the discrete Dirac operator. We then use these identities in a lattice Gross-Neveu and chiral Gross-Neveu (Nambu-Jona-Lasinio) model to show that lattice solitons correspond to reflectionless potentials associated with the discrete scattering problem. These models are of significance as they are equivalent to a mean-field theory of a lattice superconductor. To explicitly construct these solitons, we generalize supersymmetric quantum mechanics to tight-binding models. We show that a matrix transformation exists that maps a tight-binding model to an isospectral one which shares the same structure and scattering properties. The corresponding soliton solutions have both modulated hopping and onsite potential, the former of which has no analogue in the continuum limit. We explicitly compute both topological and non-topological soliton solutions as well as bound state spectra in the aforementioned models.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Macroscopic Traffic Modeling Using Probe Vehicle Data: A Machine Learning Approach

Abstract The macroscopic fundamental diagram (MFD) captures an orderly relationship among traffic flow, density, and speed at the network level. It is a simple yet powerful tool for modeling traffic dynamics in large urban networks with broad application in traffic control and management. However, empirically derived MFDs in urban regions require high-resolution traffic data from the network. Having the network flow and vehicular density estimated at the (granular) census tract level using vehicle probe data, we apply machine learning methods to predict the MFDs across U.S. urban areas and capture the impacts of location-specific input features on the network flow–density relationships at a large scale. The results show that, among the four tested machine learning approaches (Random Forest, XGBoost, Support Vector Machine, and Neural Network), XGBoost delivers the best performance in predicting network traffic flow based on vehicular density and location attributes. Using interaction Shapley Additive explanation (SHAP) values and partial correlation analysis, we examine the factors influencing MFD shapes across different locations. Our empirical findings reveal that across U.S. urban areas, network topology, transportation infrastructure, and land use are primary factors shaping MFD curves, while demand and trip-related factors play a lesser role. Specifically, higher ranking roads, centrality, and development levels correlate positively with network capacity and critical density, whereas negative associations are observed for network connectivity, mixed-use development, and road roughness levels.

Jin, Ling↗

NASA NTRS1995

Rapid Airplane Parametric Input Design (RAPID)

RAPID is a methodology and software system to define a class of airplane configurations and directly evaluate surface grids, volume grids, and grid sensitivity on and about the configurations. A distinguishing characteristic which separates RAPID from other airplane surface modellers is that the output grids and grid sensitivity are directly applicable in CFD analysis. A small set of design parameters and grid control parameters govern the process which is incorporated into interactive software for 'real time' visual analysis and into batch software for the application of optimization technology. The computed surface grids and volume grids are suitable for a wide range of Computational Fluid Dynamics (CFD) simulation. The general airplane configuration has wing, fuselage, horizontal tail, and vertical tail components. The double-delta wing and tail components are manifested by solving a fourth order partial differential equation (PDE) subject to Dirichlet and Neumann boundary conditions. The design parameters are incorporated into the boundary conditions and therefore govern the shapes of the surfaces. The PDE solution yields a smooth transition between boundaries. Surface grids suitable for CFD calculation are created by establishing an H-type topology about the configuration and incorporating grid spacing functions in the PDE equation for the lifting components and the fuselage definition equations. User specified grid parameters govern the location and degree of grid concentration. A two-block volume grid about a configuration is calculated using the Control Point Form (CPF) technique. The interactive software, which runs on Silicon Graphics IRIS workstations, allows design parameters to be continuously varied and the resulting surface grid to be observed in real time. The batch software computes both the surface and volume grids and also computes the sensitivity of the output grid with respect to the input design parameters by applying the precompiler tool ADIFOR to the grid generation program. The output of ADIFOR is a new source code containing the old code plus expressions for derivatives of specified dependent variables (grid coordinates) with respect to specified independent variables (design parameters). The RAPID methodology and software provide a means of rapidly defining numerical prototypes, grids, and grid sensitivity of a class of airplane configurations. This technology and software is highly useful for CFD research for preliminary design and optimization processes.

Smith, Robert E.↗

NASA NTRS1991

Power electronics for low power arcjets

In anticipation of the needs of future light-weight, low-power spacecraft, arcjet power electronics in the 100 to 400 W operating range were developed. Limited spacecraft power and thermal control capacity of these small spacecraft emphasized the need for high efficiency. Power topologies similar to those in the higher 2 kW and 5 to 30 kW power range were implemented, including a four transistor bridge switching circuit, current mode pulse-width modulated control, and an output current averaging inductor with an integral pulse generation winding. Reduction of switching transients was accomplished using a low inductance power distribution network, and no passive snubber circuits were necessary for power switch protection. Phase shift control of the power bridge was accomplished using an improved pulse width modulation to phase shift converter circuit. These features, along with conservative magnetics designs allowed power conversion efficiencies of greater than 92.5 percent to be achieved into resistive loads over the entire operating range of the converter. Electromagnetic compatibility requirements were not considered in this work, and control power for the converter was derived from AC mains. Addition of input filters and control power converters would result in an efficiency of on the order of 90 percent for a flight unit. Due to the developmental nature of arcjet systems at this power level, the exact nature of the thruster/power processor interface was not quantified. Output regulation and current ripple requirements of 1 and 20 percent respectively, as well as starting techniques, were derived from the characteristics of the 2 kW system but an open circuit voltage in excess of 175 V was specified. Arcjet integration tests were performed, resulting in successful starts and stable arcjet operation at power levels as low as 240 W with simulated hydrazine propellants.

Hamley, John A.↗

DOE OSTI2025

Conformal geometry from entanglement

In a physical system with conformal symmetry, observables depend on cross-ratios, measures of distance invariant under global conformal transformations (conformal geometry for short). We identify a quantum information-theoretic mechanism by which the conformal geometry emerges at the gapless edge of a 2+1D quantum many-body system with a bulk energy gap. We introduce a novel pair of information-theoretic quantities (\mathfrak{c}_{\textrm{tot}}, \eta) ( 𝔠 tot , η ) that can be defined locally on the edge from the wavefunction of the many-body system, without prior knowledge of any distance measure. We posit that, for a topological groundstate, the quantity \mathfrak{c}_{\textrm{tot}} 𝔠 tot is stationary under arbitrary variations of the quantum state, and study the logical consequences. We show that stationarity, modulo an entanglement-based assumption about the bulk, implies (i) \mathfrak{c}_{\textrm{tot}} 𝔠 tot is a non-negative constant that can be interpreted as the total central charge of the edge theory. (ii) \eta η is a cross-ratio, obeying the full set of mathematical consistency rules, which further indicates the existence of a distance measure of the edge with global conformal invariance. Thus, the conformal geometry emerges from a simple assumption on groundstate entanglement. We show that stationarity of \mathfrak{c}_{\textrm{tot}} 𝔠 tot is equivalent to a vector fixed-point equation involving \eta η , making our assumption locally checkable. We also derive similar results for 1+1D systems under a suitable set of assumptions.

Kim, Isaac H.↗