Search NASA⌕ Search

SEARCH · Search NASA

Results for “splitting 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 1,387 records · Page 77

Coherent-state quantum cryptography using pseudorandom number generators

Quantum key distribution plays an important role in modern cryptography, since the security of the transmitted keys is guaranteed by fundamental laws of nature. A method using pseudorandom number generators well known from classical cryptography is considered. It is shown that their use in quantum cryptography makes it possible to increase the key generation rate under very weak assumptions about the capabilities of the eavesdropper. A practical scheme of a coherent-state quantum key distribution protocol using pseudorandom sequences is proposed. The cryptographic strength of the proposed protocol against a beam-splitting attack is considered. (quantum cryptography)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Intertwined charge, spin, and pairing orders in doped iron ladders

Motivated by recent experimental progress on iron-based ladder compounds, we study the doped two-orbital Hubbard model for the two-leg ladder BaFe2S3. The model is constructed by using ab initio hopping parameters and the ground state properties are investigated using the density matrix renormalization group method. We show that the (π,0) magnetic ordering at half filling, with ferromagnetic rungs and antiferromagnetic legs, becomes incommensurate upon hole doping. Moreover, depending on the strength of the Hubbard U coupling, other magnetic patterns, such as (0,π), are also stabilized. We found that the binding energy for two holes becomes negative for intermediate Hubbard interaction strength, indicating hole pairing. Due to the crystal-field split among orbitals, the holes primarily reside in one orbital, with the other one remaining half filled. This resembles orbital selective Mott states. The formation of tight hole pairs continues with increasing hole density, as long as the magnetic order remains antiferromagnetic in one direction. The study of pair-pair correlations indicates the dominance of the intraorbital spin-singlet channel, as opposed to other pairing channels. Although in a range of hole doping pairing correlations decay slowly, our results can also be interpreted as corresponding to a charge density wave made of pairs, a precursor of eventual superconductivity after interladder couplings are included. Such a scenario of intertwined orders has been extensively discussed before in the cuprates, and our results suggest a similar physics could exist in ladder iron-based superconductors. Finally, we also show that a robust Hund's coupling is needed for pairing to occur.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗

Rythmos: Solution and Analysis Package for Differential-Algebraic and Ordinary-Differential Equations

Time integration is a central component for most transient simulations. It coordinates many of the major parts of a simulation together, e.g., a residual calculation with a transient solver, solution with the output, various operator-split physics, and forward and adjoint solutions for inversion. Even though there is this variety in these transient simulations, there is still a common set of algorithms and procedures to progress transient solutions for ordinary-differential equations (ODEs) and differential-alegbraic equations (DAEs). Rythmos is a collection of these algorithms that can be used for the solution of transient simulations. It provides common time-integration methods, such as Backward and Forward Euler, Explicit and Implicit Runge-Kutta, and Backward-Difference Formulas. It can also provide sensitivities, and adjoint components for transient simulations. Rythmos is a package within Trilinos, and requires some other packages (e.g., Teuchos and Thrya) to provide basic time-integration capabilities. It also can be coupled with several other Trilinos packages to provide additional capabilities (e.g., AztecOO and Belos for linear solutions, and NOX for non-linear solutions). The documentation is broken down into three parts: Theory Manual, User's Manual, and Developer's Guide. The Theory Manual contains the basic theory of the time integrators, the nomenclature and mathematical structure utilized within Rythmos, and verification results demonstrating that the designed order of accuracy is achieved. The User's Manual provides information on how to use the Rythmos, description of input parameters through Teuchos Parameter Lists, and description of convergence test examples. The Developer's Guide is a high-level discussion of the design and structure of Rythmos to provide information to developers for the continued development of capabilities. Details of individual components can be found in the Doxygen webpages.

97 MATHEMATICS AND COMPUTING↗

Hybrid Architecture Active Wavefront Sensing and Control

A method was developed for performing relatively high-speed wavefront sensing and control to overcome thermal instabilities in a segmented primary mirror telescope [e.g., James Webb Space Telescope (JWST) at L2], by using the onboard fine guidance sensor (FGS) to minimize expense and complexity. This FGS performs centroiding on a bright star to feed the information to the pointing and control system. The proposed concept is to beam split the image of the guide star (or use a single defocused guide star image) to perform wavefront sensing using phase retrieval techniques. Using the fine guidance sensor star image for guiding and fine phasing eliminates the need for other, more complex ways of achieving very accurate sensing and control that is needed for UV-optical applications. The phase retrieval occurs nearly constantly, so passive thermal stability over fourteen days is not required. Using the FGS as the sensor, one can feed segment update information to actuators on the primary mirror that can update the primary mirror segment fine phasing with this frequency. Because the thermal time constants of the primary mirror are very slow compared to this duration, the mirror will appear extremely stable during observations (to the level of accuracy of the sensing and control). The sensing can use the same phase retrieval techniques as the JWST by employing an additional beam splitter, and having each channel go through a weak lens (one positive and one negative). The channels can use common or separate detectors. Phase retrieval can be performed onboard. The actuation scheme would include a coarse stage able to achieve initial alignment of several millimeters of range (similar to JWST and can use a JWST heritage sensing approach in the science camera) and a fine stage capable of continual updates.

Feinberg, Lee↗

Metrology Optical Power Budgeting in SIM Using Statistical Analysis Techniques

The Space Interferometry Mission (SIM) is a space-based stellar interferometry instrument, consisting of up to three interferometers, which will be capable of micro-arc second resolution. Alignment knowledge of the three interferometer baselines requires a three-dimensional, 14-leg truss with each leg being monitored by an external metrology gauge. In addition, each of the three interferometers requires an internal metrology gauge to monitor the optical path length differences between the two sides. Both external and internal metrology gauges are interferometry based, operating at a wavelength of 1319 nanometers. Each gauge has fiber inputs delivering measurement and local oscillator (LO) power, split into probe-LO and reference-LO beam pairs. These beams experience power loss due to a variety of mechanisms including, but not restricted to, design efficiency, material attenuation, element misalignment, diffraction, and coupling efficiency. Since the attenuation due to these sources may degrade over time, an accounting of the range of expected attenuation is needed so an optical power margin can be book kept. A method of statistical optical power analysis and budgeting, based on a technique developed for deep space RF telecommunications, is described in this paper and provides a numerical confidence level for having sufficient optical power relative to mission metrology performance requirements.

stellar interferometry↗

Subharmonic Oscillations and Chaos in Dynamic Atomic Force Microscopy

The increasing use of dynamic atomic force microscopy (d-AFM) for nanoscale materials characterization calls for a deeper understanding of the cantilever dynamics influencing scan stability, predictability, and image quality. Model development is critical to such understanding. Renormalization of the equations governing d- AFM provides a simple interpretation of cantilever dynamics as a single spring and mass system with frequency dependent cantilever stiffness and damping parameters. The renormalized model is sufficiently robust to predict the experimentally observed splitting of the free-space cantilever resonance into multiple resonances upon cantilever-sample contact. Central to the model is the representation of the cantilever sample interaction force as a polynomial expansion with coefficients F(sub ij) (i,j = 0, 1, 2) that account for the effective interaction stiffness parameter, the cantilever-to-sample energy transfer, and the amplitude of cantilever oscillation. Application of the Melnikov method to the model equation is shown to predict a homoclinic bifurcation of the Smale horseshoe type leading to a cascade of period doublings with increasing drive displacement amplitude culminating in chaos and loss of image quality. The threshold value of the drive displacement amplitude necessary to initiate subharmonic generation depends on the acoustic drive frequency, the effective damping coefficient, and the nonlinearity of the cantilever-sample interaction force. For parameter values leading to displacement amplitudes below threshold for homoclinic bifurcation other bifurcation scenarios can occur, some of which lead to chaos.

Cantrell, John H.↗

Detecting Thermally-Induced Spinodal Decomposition with Picosecond Ultrasonics in Cast Austenitic Stainless Steels

Destructive techniques to monitor nuclear reactor component health may not always be available during service, as they are time-consuming and often require pre-installed inspection coupons. Non-destructive evaluation (NDE) techniques can bridge this gap by rapidly identifying the state of mission-critical reactor components, via inference between NDE-measurable material properties and those of ultimate interest, such as ductility and toughness. Here, we demonstrate one such inference about the health of thermally aged cast austenitic stainless steels. Observations of surface acoustic wave peak (SAW) splitting correlate with spinodal decomposition-induced embrittlement as destructively measured by Charpy impact energy. Elastodynamic calculations and molecular dynamics simulations of the effects of spinodal decomposition on elastic moduli support that the new acoustic modes present are due to stiffening in the δ-ferrite domains. Finally, this discovery enables one to probe structure-property relationships in materials in a greatly accelerated manner, suggesting that similar inference methods can be used to determine material fitness-for-service, or to quickly uncover new structure-property relationships.

304-type stainless steel↗

Electrical Control of the Nuclear Spin States of Rare-Earth Adatoms

Rare-earth adatoms on surfaces have been studied for potential atomic-scale magnetic storage, quantum sensing, and quantum computing applications. Despite accumulating experimental efforts, a comprehensive description of the electronic configurations of the adatoms remains elusive. Here, we investigate two charge states and several electronic configurations, including 5d and 6s valence shells, for a Sm adatom on a MgO substrate using multiconfigurational ab initio methods, for the possibility of using the Sm nuclear spin levels as qubits. For the configurations in a neutral charge state, we find that the electronic ground state is a singlet, and thus the hyperfine interaction associated with the 147 Sm nucleus is absent, which may greatly enhance nuclear spin coherence time. The degeneracy of the nuclear levels is lifted by the nuclear quadrupole interaction. We show that the splitting of the nuclear levels can be controlled by a static electric field, and that Rabi oscillations between the nuclear levels can be induced by a time-dependent electric field. For the configurations in a singly charged state, electronic Kramers doublets are formed. The electronic configurations including an unpaired 6s orbital exhibit a strong hyperfine Stark effect due to a large Fermi contact contribution to the hyperfine interaction. In these configurations, electric-field-induced Rabi oscillations between the electronic-nuclear levels can occur at frequencies up to 3 orders of magnitude higher than those for the neutral charge state. The proposed system may be experimentally observed within scanning tunneling microscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field methods

The fully correlated frequency-independent Dirac–Coulomb–Breit Hamiltonian provides the most accurate description of electron–electron interaction before going to a genuine relativistic quantum electrodynamics theory of many-electron systems. In this work, we introduce a correlated Dirac–Coulomb–Breit multiconfigurational self-consistent-field method within the frameworks of complete active space and density matrix renormalization group. In this approach, the Dirac–Coulomb–Breit Hamiltonian is included variationally in both the mean-field and correlated electron treatment. Here, we also analyze the importance of the Breit operator in electron correlation and the rotation between the positive- and negative-orbital space in the no-virtual-pair approximation. Atomic fine-structure splittings and lanthanide contraction in diatomic fluorides are used as benchmark studies to understand the contribution from the Breit correlation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electrical Tuning of Tin-Vacancy Centers in Diamond

Group-IV color centers in diamond have attracted significant attention as solid-state spin qubits because of their excellent optical and spin properties. Among these color centers, the tin-vacancy (Sn-V - ) center is of particular interest because its large ground-state splitting enables long spin coherence times at temperatures above 1 K. However, color centers typically suffer from inhomogeneous broadening, which can be exacerbated by nanofabrication-induced strain, hindering the implementation of quantum nodes emitting indistinguishable photons. Although strain and Raman tuning have been investigated as promising tuning techniques to overcome the spectral mismatch between distinct group-IV color centers, other approaches need to be explored to find methods that can offer more localized control without sacrificing emission intensity. Here, we study the electrical tuning of Sn-V - centers in diamond via the direct-current Stark effect. We demonstrate a tuning range beyond 1.7 GHz. We observe both quadratic and linear dependence on the applied electric field. Further, we also confirm that the tuning effect we observe is a result of the applied electric field and is distinct from thermal tuning due to Joule heating. Stark tuning is a promising avenue toward overcoming detunings between emitters and enabling the realization of multiple identical quantum nodes.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Anomalous random correlations of force constants on the lattice dynamical properties of disordered Au-Fe alloys

Here, gold iron (Au-Fe) alloys are of immense interest due to their biocompatibility, anomalous Hall conductivity, and applications in various medical treatments. However, irrespective of the method of preparation, they often exhibit a high level of disorder with properties sensitive to the thermal or magnetic annealing temperatures. We calculate the lattice dynamical properties of Au 1–x Fe x alloys using density functional theory methods where, being multisite properties, reliable interatomic force constant (IFC) calculations in disordered alloys remain a challenge. We follow a twofold approach: an accurate IFC calculation in an environment with nominally zero chemical pair correlations to mimic the homogeneously disordered alloy and a configurational averaging for the desired phonon properties (e.g., dispersion, density of states, and entropy). We find an anomalous change in the IFC's and phonon dispersion (split bands) near x=0.19, which is attributed to the local stiffening of the Au-Au bonds when Au is in the vicinity of Fe. Other results based on mechanical and thermophysical properties reflect a similar anomaly: Phonon entropy, e.g., becomes negative below x=0.19, suggesting a tendency for chemical unmixing, reflecting the onset of a miscibility gap in the phase diagram. Our results match fairly well with reported data wherever available.

36 MATERIALS SCIENCE↗

Atomic and electron structure of reconstructed (111) surface in ZnSe and CdSe crystals

The atomic and electron structure of four variants of polar (111)-(2 × 2) surfaces in ZnSe and CdSe terminated by a cation, namely, the ideal, relaxed, reconstructed, and relaxed after reconstruction surfaces, are calculated for the first time from the first principles. The surface is simulated by a film with a thickness of 12 atomic layers and a vacuum gap of ~16 Å in the layered superlattice approximation. Four fictitious hydrogen atoms with a charge of 0.5 electrons each are added for closing dangling Se bonds on the opposite side of the film. Ab initio calculations are performed using the QUANTUM ESPRESSO software based on the density functional theory. It is shown that relaxation results in splitting of atomic layers. We calculate and analyze the band structures and total and layer-wise densities of electron states for four variants of the surface.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

High pO2 Flux Growth and Characterization of NdNiO3 Crystals

Single crystals of the perovskite nickelate NdNiO3 with dimensions of up to 50 μm on edge have been successfully grown using the flux method at a temperature of 400 °C and oxygen pressure of 200 bar. The crystals were investigated by a combination of techniques, including high-resolution synchrotron X-ray single-crystal and powder diffraction and physical property measurements such as magnetic susceptibility and resistivity. Resistivity measurements revealed a metal-insulator transition (MIT) at TMIT~180 K with apparent thermal hysteresis; however, no superlattice peaks or peak splitting below TMIT, which corresponds to a structural transition from Pbnm to P21/n, was observed. The successful growth of NdNiO3 crystals at relatively low temperatures and oxygen pressure provides an alternative approach for preparing single crystals of interesting perovskites such as RNiO3 (R = Sm-Lu) and parent phases of superconducting square planar nickelates.

Wang, Xiaoli↗

Improving Chirped Fiber Bragg Grating Resolution for Position-Sensitive Sensors in Shock- and Detonation-Driven Experiments

Chirped fiber Bragg gratings (CFBGs) are robust diagnostic sensors that are widely used to track detonation-driven and shock wave propagation. CFBGs are inscribed with a linearly chirped periodic index of refraction changes that alter the Bragg wavelength along the length of the probe. The light return of each individual Bragg element is captured by a detector at a unique time to map the full reflected spectrum. The CFBG spectrum is measured with a dispersive Fourier transform of the reflected light that temporally stretches the spectrum to increase spatial resolution and make a one-to-one map of the wavelength on a time axis. Here, we propose an improvement of CFBG temporal resolution by incorporating two co-linear laser pulses with orthogonal polarization states and a 5 ns time offset. The two separate signals were split and tracked by two separate detectors. An oscilloscope captured good separation in the signals, and two separate spectrograms were generated and interleaved in the post-processing of the data. This novel technique doubled the CFBG temporal resolution and led to a doubled location resolution. As a proof-of-concept of this technique, the resolution improvement was compared between standard CFBG measurements and the two polarization states method on a position-sensitive CFBG sensor. CFBG resolution doubling will advance sensor capabilities and will have a direct impact on improving capture and analysis in dynamic, high-explosive experiments.

42 ENGINEERING↗

The shape of cooling filaments in old supernova remnants

The evolution of a filament in an old supernova remnant is investigated in two dimensions using analytical methods to determine its structure and stability and to show whether it collapses into a sheet or a cylinder. All parameters are assumed constant along the length of the feature. The cross section of the filament is found to flatten as it cools, and to degenerate into a thin sheet which is unstable to asymmetric warping disturbances. This probably leads to a number of parallel striations as the filament splits up, consistent with observations in the Cygnus Loop. The ratio of depth to width is typically 10:1, and so the feature may appear to be tubelike or a sheet, depending on perspective. The application to Vela and Puppis, IC 443, S 147, and CTB 1 is pointed out.

Smith, M. D.↗

A comparison of internal energy calculation methods for diatomic molecules

Various methods of calculating the internal energy of diatomic molecules are studied. An accurate and efficient method for computing the eigenvalues of the vibrational Schroedinger equation for an arbitrary potential is developed. The method is based on a finite-element discretization using the cubic Lobatto element. A combination of spectrum slicing and the Laguerre algorithm is used to solve for the eigenvalues. A simple method to compute the quasi-bound states is presented. For N2 molecules, all vibrational-rotational states of eleven available electronic potentials are computed, and summed to obtain the exact internal energy function with temperature. The total computation required 314 seconds of CPU-time on NASA's Cray 2 computer. Various approximate models are discussed and compared with the exact numerical simulation. It is shown that the splitting of the macroscopic internal energy into separate electronic, rotational, and vibrational energies is not justified at high temperatures.

Liu, Yen↗

Characterization of the rotation of cometary nuclei

The primary methods used to determine the spin state of cometary nuclei and the pitfalls and successes experienced in their use are reviewed. Attention is given to the theory of rigid body rotation from the standpoint of remote (astronomical) observers, with emphasis on what is known of the effects of nutation of lightcurves, the influence of torques induced by jet activity, the effects of internal energy and mass dissipation, and nuclear splitting. The available knowledge on rotation for eight comets, including Comet P/Halley, is reviewed. The following questions are addressed: achievement of a consensus on the rotational state of P/Halley, the possibility of accurately determining the amplitude of the transverse nongravitational force associated with rotation, and whether the orientations of fanlike comas are a valid indicator of the orientation of the spin vector. The prerequisites for answering these questions are discussed.

Belton, Michael J. S.↗