SEARCH · Engineering Papers

Results for “variational principles”

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 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS2000

A Parallel Genetic Algorithm for Automated Electronic Circuit Design

Parallelized versions of genetic algorithms (GAs) are popular primarily for three reasons: the GA is an inherently parallel algorithm, typical GA applications are very compute intensive, and powerful computing platforms, especially Beowulf-style computing clusters, are becoming more affordable and easier to implement. In addition, the low communication bandwidth required allows the use of inexpensive networking hardware such as standard office ethernet. In this paper we describe a parallel GA and its use in automated high-level circuit design. Genetic algorithms are a type of trial-and-error search technique that are guided by principles of Darwinian evolution. Just as the genetic material of two living organisms can intermix to produce offspring that are better adapted to their environment, GAs expose genetic material, frequently strings of 1s and Os, to the forces of artificial evolution: selection, mutation, recombination, etc. GAs start with a pool of randomly-generated candidate solutions which are then tested and scored with respect to their utility. Solutions are then bred by probabilistically selecting high quality parents and recombining their genetic representations to produce offspring solutions. Offspring are typically subjected to a small amount of random mutation. After a pool of offspring is produced, this process iterates until a satisfactory solution is found or an iteration limit is reached. Genetic algorithms have been applied to a wide variety of problems in many fields, including chemistry, biology, and many engineering disciplines. There are many styles of parallelism used in implementing parallel GAs. One such method is called the master-slave or processor farm approach. In this technique, slave nodes are used solely to compute fitness evaluations (the most time consuming part). The master processor collects fitness scores from the nodes and performs the genetic operators (selection, reproduction, variation, etc.). Because of dependency issues in the GA, it is possible to have idle processors. However, as long as the load at each processing node is similar, the processors are kept busy nearly all of the time. In applying GAs to circuit design, a suitable genetic representation 'is that of a circuit-construction program. We discuss one such circuit-construction programming language and show how evolution can generate useful analog circuit designs. This language has the desirable property that virtually all sets of combinations of primitives result in valid circuit graphs. Our system allows circuit size (number of devices), circuit topology, and device values to be evolved. Using a parallel genetic algorithm and circuit simulation software, we present experimental results as applied to three analog filter and two amplifier design tasks. For example, a figure shows an 85 dB amplifier design evolved by our system, and another figure shows the performance of that circuit (gain and frequency response). In all tasks, our system is able to generate circuits that achieve the target specifications.

Long, Jason D.↗

NASA NTRS2015

TTEthernet for Integrated Spacecraft Networks

Aerospace projects have traditionally employed federated avionics architectures, in which each computer system is designed to perform one specific function (e.g. navigation). There are obvious downsides to this approach, including excessive weight (from so much computing hardware), and inefficient processor utilization (since modern processors are capable of performing multiple tasks). There has therefore been a push for integrated modular avionics (IMA), in which common computing platforms can be leveraged for different purposes. This consolidation of multiple vehicle functions to shared computing platforms can significantly reduce spacecraft cost, weight, and design complexity. However, the application of IMA principles introduces significant challenges, as the data network must accommodate traffic of mixed criticality and performance levels - potentially all related to the same shared computer hardware. Because individual network technologies are rarely so competent, the development of truly integrated network architectures often proves unreasonable. Several different types of networks are utilized - each suited to support a specific vehicle function. Critical functions are typically driven by precise timing loops, requiring networks with strict guarantees regarding message latency (i.e. determinism) and fault-tolerance. Alternatively, non-critical systems generally employ data networks prioritizing flexibility and high performance over reliable operation. Switched Ethernet has seen widespread success filling this role in terrestrial applications. Its high speed, flexibility, and the availability of inexpensive commercial off-the-shelf (COTS) components make it desirable for inclusion in spacecraft platforms. Basic Ethernet configurations have been incorporated into several preexisting aerospace projects, including both the Space Shuttle and International Space Station (ISS). However, classical switched Ethernet cannot provide the high level of network determinism required by real-time spacecraft applications. Even with modern advancements, the uncoordinated (i.e. event-driven) nature of Ethernet communication unavoidably leads to message contention within network switches. The arbitration process used to resolve such conflicts introduces variation in the time it takes for messages to be forwarded. TTEthernet1 introduces decentralized clock synchronization to switched Ethernet, enabling message transmission according to a time-triggered (TT) paradigm. A network planning tool is used to allocate each device a finite amount of time in which it may transmit a frame. Each time slot is repeated sequentially to form a periodic communication schedule that is then loaded onto each TTEthernet device (e.g. switches and end systems). Each network participant references the synchronized time in order to dispatch messages at predetermined instances. This schedule guarantees that no contention exists between time-triggered Ethernet frames in the network switches, therefore eliminating the need for arbitration (and the timing variation it causes). Besides time-triggered messaging, TTEthernet networks may provide two additional traffic classes to support communication of different criticality levels. In the rate-constrained (RC) traffic class, the frame payload size and rate of transmission along each communication channel are limited to predetermined maximums. The network switches can therefore be configured to accommodate the known worst-case traffic pattern, and buffer overflows can be eliminated. The best-effort (BE) traffic class behaves akin to classical Ethernet. No guarantees are provided regarding transmission latency or successful message delivery. TTEthernet coordinates transmission of all three traffic classes over the same physical connections, therefore accommodating the full spectrum of traffic criticality levels required in IMA architectures. Common computing platforms (e.g. LRUs) can share networking resources in such a way that failures in non-critical systems (using BE or RC communication modes) cannot impact flight-critical functions (using TT communication). Furthermore, TTEthernet hardware (e.g. switches, cabling) can be shared by both TTEthernet and classical Ethernet traffic.

Loveless, Andrew↗

DOE OSTI2021

Entanglement rearrangement in self-consistent nuclear structure calculations

Background: Entanglement plays a central role in a diverse array of increasingly important research areas, including quantum computation, simulation, measurement, sensing, and communication. Extensive suites of investigations have been performed to better understand entanglement in atomic and molecular quantum many-body systems, while the exploration of entanglement in the structure of nuclei and their reactions is presently in its infancy. Purpose: The goal of this work is to begin investigating the entanglement properties of nuclei from first-principles nuclear many-body calculations. We attempt to identify common features and emergent structures of entanglement that could ultimately lead to new and natural many-body schemes. With an eye toward quantum accelerators in future hybrid-supercomputers, criteria for partitioning nuclear many-body calculations into quantum and classical components may provide advantages in future large-scale computations. Along the way we look for explanations of the relative success of phenomenological models such as the nuclear shell model, and for better ways to match to low-energy nuclear effective field theories and lattice QCD calculations to nuclear many-body techniques that are based upon entanglement. Method: We explore the entanglement between single-particle states in 4 He and 6 He. The patterns of entanglement emerging from different single-particle bases are compared, and possible links with the convergence of observables are explored, in particular, ground-state energies. The nuclear wave functions are obtained by performing active-space no-core configuration-interaction calculations using a two-body nucleon-nucleon interaction derived from chiral effective field theory. Entanglement measures within single-particle bases exhibiting different degrees of complexity are determined, in particular, harmonic oscillator (HO), Hartree-Fock (HF), natural (NAT) and variational natural (VNAT) bases. Specifically, single-orbital entanglement entropy, two-orbital mutual information, and negativity are studied. Results: The entanglement structures in 4 He and 6 He are found to be more localized within NAT and VNAT bases than within a HO basis for the optimal HO parameters we have worked with. In particular a core-valence structure clearly emerges from the full no-core calculation of 6 He. Here, the two-nucleon mutual information shows that the VNAT basis, which typically exhibits good convergence properties, effectively decouples the active and inactive spaces. Conclusions: Measures of one- and two-nucleon entanglement are found to be useful in analyzing the structure of nuclear wave functions, in particular the efficacy of basis states, and may provide useful metrics toward developing more efficient schemes for ab initio computations of the structure and reactions of nuclei, and quantum many-body systems more generally.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2020

Shell-model studies of the astrophysical r p -process reactions S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35

Background: Dust grains condensed in the outflows of presolar classical novae should have been present in the protosolar nebula. Candidates for such presolar nova grains have been found in primitive meteorites and can in principle be identified by their isotopic ratios, but the ratios predicted by state-of-the-art one-dimensional hydrodynamic models are uncertain due to nuclear-physics uncertainties. Purpose: To theoretically calculate the thermonuclear rates and uncertainties of the S 34 ( p , γ ) Cl 35 and Cl 34 g , m ( p , γ ) Ar 35 reactions and investigate their impacts on the predicted 34 S/ 32 S isotopic ratio for presolar nova grains. Method: A shell-model approach in a (0 +1 ) ℏ ω model space was used to calculate the properties of resonances in the 34 S(p ,γ ) 35 Cl and 34 g,m Cl (p ,γ ) 35 Ar reactions and their thermonuclear rates. Uncertainties were estimated using a Monte Carlo method. The implications of these rates and their uncertainties on sulfur isotopic nova yields were investigated using a postprocessing nucleosynthesis code. The rates for transitions from the ground state of 34 Cl as well as from the isomeric first excited state of 34 Cl were explicitly calculated. Results: At energies in the resonance region near the proton-emission threshold, many negative-parity states appear. Energies, spectroscopic factors, and proton-decay widths are reported. The resulting thermonuclear rates are compared with previous determinations. Conclusions: The shell-model calculations alone are sufficient to constrain the variation of the 34 S/ 32 S ratios to within about 30%. Uncertainties associated with other reactions must also be considered, but in general we find that the 34 S/ 32 S ratios are not a robust diagnostic to clearly identify presolar grains made from nova ejecta.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS1991

System and component design and test of a 10 hp, 18,000 rpm AC dynamometer utilizing a high frequency AC voltage link, part 1

Hard and soft switching test results conducted with one of the samples of first generation MOS-controlled thyristor (MCTs) and similar test results with several different samples of second generation MCT's are reported. A simple chopper circuit is used to investigate the basic switching characteristics of MCT under hard switching and various types of resonant circuits are used to determine soft switching characteristics of MCT under both zero voltage and zero current switching. Next, operation principles of a pulse density modulated converter (PDMC) for three phase (3F) to 3F two-step power conversion via parallel resonant high frequency (HF) AC link are reviewed. The details for the selection of power switches and other power components required for the construction of the power circuit for the second generation 3F to 3F converter system are discussed. The problems encountered in the first generation system are considered. Design and performance of the first generation 3F to 3F power converter system and field oriented induction moter drive based upon a 3 kVA, 20 kHz parallel resonant HF AC link are described. Low harmonic current at the input and output, unity power factor operation of input, and bidirectional flow capability of the system are shown via both computer and experimental results. The work completed on the construction and testing of the second generation converter and field oriented induction motor drive based upon specifications for a 10 hp squirrel cage dynamometer and a 20 kHz parallel resonant HF AC link is discussed. The induction machine is designed to deliver 10 hp or 7.46 kW when operated as an AC-dynamo with power fed back to the source through the converter. Results presented reveal that the proposed power level requires additional energy storage elements to overcome difficulties with a peak link voltage variation problem that limits reaching to the desired power level. The power level test of the second generation converter after the addition of extra energy storage elements to the HF link are described. The importance of the source voltage level to achieve a better current regulation for the source side PDMC is also briefly discussed. The power levels achieved in the motoring mode of operation show that the proposed power levels achieved in the generating mode of operation can also be easily achieved provided that no mechanical speed limitation were present to drive the induction machine at the proposed power level.

Lipo, Thomas A.↗

NASA NTRS2006

High-Performance Scanning Acousto-Ultrasonic System

A high-performance scanning acousto-ultrasonic system, now undergoing development, is designed to afford enhanced capabilities for imaging microstructural features, including flaws, inside plate specimens of materials. The system is expected to be especially helpful in analyzing defects that contribute to failures in polymer- and ceramic-matrix composite materials, which are difficult to characterize by conventional scanning ultrasonic techniques and other conventional nondestructive testing techniques. Selected aspects of the acousto-ultrasonic method have been described in several NASA Tech Briefs articles in recent years. Summarizing briefly: The acousto-ultrasonic method involves the use of an apparatus like the one depicted in the figure (or an apparatus of similar functionality). Pulses are excited at one location on a surface of a plate specimen by use of a broadband transmitting ultrasonic transducer. The stress waves associated with these pulses propagate along the specimen to a receiving transducer at a different location on the same surface. Along the way, the stress waves interact with the microstructure and flaws present between the transducers. The received signal is analyzed to evaluate the microstructure and flaws. The specific variant of the acousto-ultrasonic method implemented in the present developmental system goes beyond the basic principle described above to include the following major additional features: Computer-controlled motorized translation stages are used to automatically position the transducers at specified locations. Scanning is performed in the sense that the measurement, data-acquisition, and data-analysis processes are repeated at different specified transducer locations in an array that spans the specimen surface (or a specified portion of the surface). A pneumatic actuator with a load cell is used to apply a controlled contact force. In analyzing the measurement data for each pair of transducer locations in the scan, the total (multimode) acousto-ultrasonic response of the specimen is utilized. The analysis is performed by custom software that extracts parameters of signals in the time and frequency domains. The computer hardware and software provide both real-time and postscan processing and display options. For example, oscilloscope displays of waveforms and power spectral densities are available in real time. Images can be computed while scanning continues. Signals can be digitally preprocessed and/or post-processed by filtering, windowing, time-segmenting, and running-waveform-averaging algorithms. In addition, the software affords options for off-line simulation of the waveform-data-acquisition and scanning processes. In tests, the system has been shown to be capable of characterizing microstructural changes and defects in SiC/SiC and C/SiC ceramic-matrix composites. Delaminations, variations in density, microstructural changes attributable to infiltration by silicon, and crack-space indications (defined in the next sentence) have been revealed in images formed from several time- and frequency-domain parameters of scanning acousto-ultrasonic signals. The crack-space indications were image features that were not revealed by other nondestructive testing methods and are so named because they turned out to mark locations where cracking eventually occurred.

Roth, Don↗

DOE OSTI2021

Design principles of tandem cascade photoelectrochemical devices

Cascade photoelectrocatalysis (PEC) is a possible method to improve the selectivity of solar-driven CO 2 reduction (CO 2 R). This concept can be realized by coupling different CO 2 R catalysts to different subcells in a multijunction photovoltaic (PV) stack. Efficient implementation will require careful tuning of the photocurrents and design of the photovoltages provided by the subcells to the CO 2 R catalysts in such a way as to facilitate the target reaction. Here, we outline the design principles of the tandem PEC approach using two-step conversion of CO 2 to ethylene in aqueous electrolyte, via a CO intermediate, as a model system. To perform this reaction, the first coupled PV-catalyst component should provide 4 electrons to reduce 2 molecules of CO 2 to CO; the second component should provide 8 electrons to reduce 2 CO molecules to C 2 H 4 . Based on known CO 2 R catalysts, the overpotential required to produce CO can be less than that required to reduce it to ethylene, creating the opportunity for improved efficiency. Cascade PEC can be realized in a three-terminal tandem (3TT) configuration using III–V-semiconductor based subcells coupled to Au (produces CO intermediate) and Cu (converts CO to ethylene). The current to each catalyst can be controlled by the area of the subcell exposed to the electrolyte, and the photovoltage is determined by the materials selected and device configuration. Operating conditions are found by simulating the coupled system using the open-source circuit simulator SPICE (Simulation Program with Integrated Circuits Emphasis). We identify conditions under which a 3TT configuration can have a higher solar to chemical conversion efficiency compared to a two-terminal two-junction tandem (2T 2J) with the same absorbers and a Cu catalyst only. We also show that 3TT PEC devices can be less sensitive to variations in catalyst activity compared to 2T devices. Lastly, we discuss the applications of cascade PEC to CO 2 reduction, using different intermediates, and to other chemical networks.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2023

Variation of carrier density in semimetals via short-range correlation: A case study with nickelate NdNiO 2

Carrier density is one of the key controlling factors of material properties, particularly in controlling the essential correlations in strongly correlated materials. Typically, carrier density is externally tuned by doping or gating and remains fixed below room temperature. Strangely, the carrier density in correlated semimetals is often found to vary sensitively against weak external controls such as temperature, magnetic field, and pressure. Here, we develop a realistic simulation scheme that incorporates interatomic noncollinear magnetic correlation without a long-range order. Using the recently discovered nickelate superconductor as an example, we demonstrate a rather generic low-energy mechanism that in semimetals short-range correlation can reversely modulate the carrier density as well. Such a mutual influence between correlation and carrier density provides an extra ingredient for sensitive bifurcating behavior. This special feature of correlated semimetals explains their versatile carrier density at low energy and opens up new possibilities of functionalizing these materials.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Mirror symmetry breaking and lateral stacking shifts in twisted trilayer graphene

Here, we construct a continuum model of twisted trilayer graphene using ab initio density-functional-theory calculations, and we apply it to address twisted trilayer electronic structure. Our model accounts for moiré variation in site energies, hopping between outside layers, and hopping within layers. We focus on the role of a mirror symmetry present in ABA graphene trilayers with a middle layer twist. The mirror symmetry is lost intentionally when a displacement field is applied between layers, and unintentionally when the top layer is shifted laterally relative to the bottom layer. We calculate two band-structure characteristics that are directly relevant to transport measurements, namely the Drude weight and the weak-field Hall conductivity. We relate them via the Hall density to assess the influence of the accidental lateral stacking shifts currently present in all experimental devices on electronic properties. Finally, we comment on the role of the possible importance of accidental lateral stacking shifts for superconductivity in twisted trilayers.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRSDate not supplied

Preliminary Results Cycling GEOS-JEDI with GSI-based Background Errors

The first phase of transitioning the NASA GMAO GEOS atmospheric data assimilation capabilities to JEDI involves the replacement of the Grid-point Statistical Interpolation (GSI) with a corresponding JEDI analysis. This includes taking JEDI's Unified Observation Operator (UFO), its underlying dependencies, and the JEDI solver that enables a hybrid 4DEnVar strategy similar to what is used in the current GEOS-GSI system. Variational analysis involves at least two main components associated with the observation and background cost function terms. The first is directly related to the UFO, which is being carefully validated in a joint collaboration between GMAO and NCEP to demonstrate consistency with corresponding observations usage in GSI. The second component is the background term, which in a hybrid system involves the ability to set up both a climatologically-based term and an ensemble-based term. JEDI provides the means to implement both terms through its BUMP component. Use of BUMP would require a complete re-tune of both climatological and ensemble, which is a non-trivial exercise we would prefer to avoid. As an alternative, the work here studies the results of interfacing the GSI-background error capability (GSIBEC) into JEDI through SABER. With this, the exact same background error covariance formulation used in GSI can be employed in JEDI without need for re-tuning. This brief summary covers the work done to interface GSIBEC into JEDI and shows preliminary results where the background error covariances of the control (GEOS-GSI) and experiment (GEOS-JEDI) are identical in corresponding cycling experiments. The cycling exercise is obviously preliminary and so much can be expected from GEOS-JEDI when compared to GEOS-GSI. There is still a number of features that need closer attention and although in some cases in principle ready to cycle have been intentionally either turned off or not fully exercised (e.g., VarBC is applied but not cycled). Other features are still pending implementation, one such example is the implementation of the Tangent Linear Normal Mode Constraint. Still, results are quite encouraging as hopefully the discussion here illustrates.

Ricardo Todling↗

DOE OSTI2023

Spectral-partitioned Kohn-Sham density functional theory

Here we introduce a general, variational scheme for systematic approximation of a given Kohn-Sham free-energy functional by partitioning the density matrix into distinct spectral domains, each of which may be spanned by an independent diagonal representation without requirement of mutual orthogonality. It is shown that by generalizing the entropic contribution to the free energy to allow for independent representations in each spectral domain, the free energy becomes an upper bound to the exact (unpartitioned) Kohn-Sham free energy, attaining this limit as the representations approach Kohn-Sham eigenfunctions. A numerical procedure is devised for calculation of the generalized entropy associated with spectral partitioning of the density matrix. The result is a powerful framework for Kohn-Sham calculations of systems whose occupied subspaces span multiple energy regimes. As a case in point, we apply the proposed framework to warm- and hot-dense matter described by finite-temperature density functional theory, where at high energies the density matrix is represented by that of the free-electron gas, while at low energies it is variationally optimized. We derive expressions for the spectral-partitioned Kohn-Sham Hamiltonian, atomic forces, and macroscopic stresses within the projector-augmented wave (PAW) and the norm-conserving pseudopotential methods. It is demonstrated that at high temperatures, spectral partitioning facilitates accurate calculations at dramatically reduced computational cost. Moreover, as temperature is increased, fewer exact Kohn-Sham states are required for a given accuracy, leading to further reductions in computational cost. Finally, it is shown that standard multiprojector expansions of electronic orbitals within atomic spheres in the PAW method lack sufficient completeness at high temperatures. Spectral partitioning provides a systematic solution for this fundamental problem.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1996

Algebraic Multigrid by Smoothed Aggregation for Second and Fourth Order Elliptic Problems

Multigrid methods are very efficient iterative solvers for system of algebraic equations arising from finite element and finite difference discretization of elliptic boundary value problems. The main principle of multigrid methods is to complement the local exchange of information in point-wise iterative methods by a global one utilizing several related systems, called coarse levels, with a smaller number of variables. The coarse levels are often obtained as a hierarchy of discretizations with different characteristic meshsizes, but this requires that the discretization is controlled by the iterative method. To solve linear systems produced by existing finite element software, one needs to create an artificial hierarchy of coarse problems. The principal issue is then to obtain computational complexity and approximation properties similar to those for nested meshes, using only information in the matrix of the system and as little extra information as possible. Such algebraic multigrid method that uses the system matrix only was developed by Ruge. The prolongations were based on the matrix of the system by partial solution from given values at selected coarse points. The coarse grid points were selected so that each point would be interpolated to via so-called strong connections. Our approach is based on smoothed aggregation introduced recently by Vanek. First the set of nodes is decomposed into small mutually disjoint subsets. A tentative piecewise constant interpolation (in the discrete sense) is then defined on those subsets as piecewise constant for second order problems, and piecewise linear for fourth order problems. The prolongation operator is then obtained by smoothing the output of the tentative prolongation and coarse level operators are defined variationally.

Vanek, Petr↗

NASA NTRS2013

Hydrometeor Size Distribution Measurements by Imaging the Attenuation of a Laser Spot

The optical extinction of a laser due to scattering of particles is a well-known phenomenon. In a laboratory environment, this physical principle is known as the Beer-Lambert law, and is often used to measure the concentration of scattering particles in a fluid or gas. This method has been experimentally shown to be a usable means to measure the dust density from a rocket plume interaction with the lunar surface. Using the same principles and experimental arrangement, this technique can be applied to hydrometeor size distributions, and for launch-pad operations, specifically as a passive hail detection and measurement system. Calibration of a hail monitoring system is a difficult process. In the past, it has required comparison to another means of measuring hydrometeor size and density. Using a technique recently developed for estimating the density of surface dust dispersed during a rocket landing, measuring the extinction of a laser passing through hail (or dust in the rocket case) yields an estimate of the second moment of the particle cloud, and hydrometeor size distribution in the terrestrial meteorological case. With the exception of disdrometers, instruments that measure rain and hail fall make indirect measurements of the drop-size distribution. Instruments that scatter microwaves off of hydrometeors, such as the WSR-88D (Weather Surveillance Radar 88 Doppler), vertical wind profilers, and microwave disdrometers, measure the sixth moment of the drop size distribution (DSD). By projecting a laser onto a target, changes in brightness of the laser spot against the target background during rain and hail yield a measurement of the DSD's second moment by way of the Beer-Lambert law. In order to detect the laser attenuation within the 8-bit resolution of most camera image arrays, a minimum path length is required. Depending on the intensity of the hail fall rate for moderate to heavy rainfall, a laser path length of 100 m is sufficient to measure variations in optical extinction using a digital camera. For hail fall only, the laser path may be shorter because of greater scattering due to the properties of hailstones versus raindrops. A photodetector may replace the camera in automated installations. Laser-based rain and hail measurement systems are available, but they are based on measuring the interruption of a thin laser beam, thus counting individual hydrometeors. These systems are true disdrometers since they also measure size and velocity. The method reported here is a simple method, requiring far less processing, but it is not a disdrometer.

Lane, John↗

DOE OSTI2020

Starting-point-independent quantum Monte Carlo calculations of iron oxide

Quantum Monte Carlo (QMC) methods are useful for studies of strongly correlated materials because they are many body in nature and use the physical Hamiltonian. Typical calculations assume as a starting point a wave function constructed from single-particle orbitals obtained from one-body methods, e.g., density functional theory. However, mean-field-derived wave functions can sometimes lead to systematic QMC biases if the mean-field result poorly describes the true ground state. In this study, we examine the accuracy and flexibility of QMC trial wave functions using variational and fixed-node diffusion QMC estimates of the total spin density and lattice distortion of antiferromagnetic iron oxide (FeO) in the ground state B1 crystal structure. We found that for relatively simple wave functions the predicted lattice distortion was controlled by the choice of single-particle orbitals used to construct the wave function, rather than by subsequent wave function optimization techniques within QMC. By optimizing the orbitals with QMC, we then demonstrate starting-point independence of the trial wave function with respect to the method by which the orbitals were constructed by demonstrating convergence of the energy, spin density, and predicted lattice distortion for two qualitatively different sets of orbitals. The results suggest that orbital optimization is a promising method for accurate many-body calculations of strongly correlated condensed phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1994

Proposed law of nature linking impacts, plume volcanism, and Milankovitch cycles to terrestrial vertebrate mass extinctions via greenhouse-embryo death coupling

A greenhouse-physiological coupling killing mechanism active among mammals, birds, and reptiles has been identified. Operating via environmental thermal effects upon the maternal core-skin blood flow critical to the survival and development of embryos, it reduces the flow of blood to the uterine tract. Today, during hot summers, this phenomena kills embryos on a vast, global scale. Because of sensitivity of many mammals to modern heat, a major modern greenhouse could reduce population numbers on a global scale, and potentially trigger population collapses in the more vulnerable parts of the world. In the geological past, the killing mechanism has likely been triggered into action by greenhouse warming via impact events, plume volcanism, and Earth orbital variations (Milankovitch cycles). Earth's biosphere is maintained and molded by the flow of energy from the solar energy source to Earth and on to the space energy sink (SES). This SES energy flow maintains Earth's biosphere and its living components, as open, intermediate, dissipative, nonequilibrium systems whose states are dependent upon the rate of energy flowing through them. Greenhouse gases such as CO2 in the atmosphere influence the SES energy flow rate. Steady-state flow is necessary for global ecological stability (autopoiesis). Natural fluctuations of the C cycle such as rapid releases of CO2 from the mantle, or oceans, disrupt steady-state SES flow. These fluctuations constantly challenge the biosphere; slowdown of SES energy flow drives it toward thermodynamical equilibrium and stagnation. Fluctuations induced by impact event, mantle plume volcanism, and Milankovitch cycles can grow into structure-breaking waves triggering major perturbations of Earth's C cycle and mass extinctions. A major C cycle perturbation involves readjustment of the outer physiochemical spheres of the Earth: the atmosphere, hydrosphere, and lithosphere; and by necessity, the biosphere. A greenhouse, one manifestation of a major C cycle perturbation, is the most dangerous natural phenomenon that life on Earth can experience. Greenhouse conditions existed during the KT mass extinctions of 65 m.y. ago, and the Pleistocene-Holocene (P-H) mammalian extinctions of 10,000-12,000 yr ago. Coupling climatology to reproductive physiology via effects of ambient air temperature upon uterine blood flow to developing embryos accounts for the extinctions via established physiological principles.

Mclean, D. M.↗

DOE OSTI2022

Interlayer spacing in pillared and grafted MCM-22 type silicas: density functional theory analysis versus experiment

Pillaring of synthetic layered crystalline silicates and aluminosilicates provides a strategy to enhance their adsorption and separation performance, and can facilitate the understanding of such behavior in more complex natural clays. In this study, we perform the first-principles density functional theory calculations for the pillaring of the pure silica polymorph of an MCM-22 type molecular sieve. Starting with a precursor material MCM-22P with fully hydroxylated layers, a pillaring agent, (EtO) 3 SiR, can react with hydroxyl groups (–OH) on adjacent internal surfaces, 2(–OH) + (EtO) 3 SiR + H 2 O → (–O) 2 SiOHR + 3EtOH, to form a pillar bridging these surfaces, or with a single hydroxyl, –OH + (EtO) 3 SiR + 2H 2 O → (–O)Si(OH) 2 R + 3EtOH, grafting to one surface. For computational efficiency, we replace the experimental organic ligand, R, by a methyl group. We find that the interlayer spacing in MCM-22 is reduced by 2.66 Å relative to weakly bound layers in the precursor MCM-22P. Including (–O) 2 SiR bridges for 50% (100%) of the hydroxyl sites in MCM-22P increases the interlayer spacing relative to MCM-22 by 2.52 Å (2.46 Å). For comparison, we also analyze the system where all –OH groups in MCM-22P are replaced by non-bridging grafted (–O)Si(OH) 2 R which results in a smaller interlayer spacing expansion of 2.17 Å relative to MCM-22. Our results for the interlayer spacing in the pillared materials are compatible with experimental observations for a similar MCM-22 type material with low Al content (Si : Al = 51 : 1) of an expansion relative to MCM-22 of roughly 2.8 Å and 2.5 Å from our x-ray diffraction and scanning transmission electron microscopy analyses, respectively. The latter analysis reveals significant variation in individual layer spacings.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2025

Mechanically and electrically tunable Rashba-Edelstein effect in ferroelectric semiconductors, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl)

The ability of materials to convert charge current into spin current is fundamental to many spintronics applications. One means of realizing this conversion is via Rashba-Edelstein effect (REE). Using density functional theory simulations, we predict that REE can be induced in the recently discovered family of semiconducting ferroelectrics, CsGe⁢𝑋 3 (𝑋 = I, Br, Cl). The effect is quantified through Rashba-Edelstein coefficients, 𝜒 𝑥⁢𝑦 =−𝜒 𝑦⁢𝑥 , which are nonzero in valence and conduction bands. The largest values, obtained for CsGeI 3 , are 3.45 × 10 10 $\frac{ℏ}{Λ cm}$ and 0.97 × 10 10 $\frac{ℏ}{Λ cm}$ in the conduction and valence bands, respectively. The values are comparable to, and sometimes exceed, those of other inorganic materials, although the maximal values occur away from the band edges. The coefficients' sign couples to the direction of spontaneous polarization, offering opportunities for nonvolatile spin current manipulation via external electric field. Furthermore, these coefficients are highly tunable through strain engineering owing to strain-induced variations in spin textures and energy isosurfaces. Specifically, 𝜒 in the valence band of CsGeI 3 is enhanced to 3.61 × 10 10 $\frac{ℏ}{Λ cm}$ under 5% biaxial strain. Given this potential, we synthesized CsGeI 3 to validate our simulation structure and found excellent agreement between experiment and simulation, thereby allowing extrapolation of our findings to practically significant temperatures. Furthermore, our study identifies promising materials for nonvolatile, multifunctional spintronic applications.

Electronic structure↗

NASA NTRS2009

Spatial Modulation Improves Performance in CTIS

Suitably formulated spatial modulation of a scene imaged by a computed-tomography imaging spectrometer (CTIS) has been found to be useful as a means of improving the imaging performance of the CTIS. As used here, "spatial modulation" signifies the imposition of additional, artificial structure on a scene from within the CTIS optics. The basic principles of a CTIS were described in "Improvements in Computed- Tomography Imaging Spectrometry" (NPO-20561) NASA Tech Briefs, Vol. 24, No. 12 (December 2000), page 38 and "All-Reflective Computed-Tomography Imaging Spectrometers" (NPO-20836), NASA Tech Briefs, Vol. 26, No. 11 (November 2002), page 7a. To recapitulate: A CTIS offers capabilities for imaging a scene with spatial, spectral, and temporal resolution. The spectral disperser in a CTIS is a two-dimensional diffraction grating. It is positioned between two relay lenses (or on one of two relay mirrors) in a video imaging system. If the disperser were removed, the system would produce ordinary images of the scene in its field of view. In the presence of the grating, the image on the focal plane of the system contains both spectral and spatial information because the multiple diffraction orders of the grating give rise to multiple, spectrally dispersed images of the scene. By use of algorithms adapted from computed tomography, the image on the focal plane can be processed into an image cube a three-dimensional collection of data on the image intensity as a function of the two spatial dimensions (x and y) in the scene and of wavelength (lambda). Thus, both spectrally and spatially resolved information on the scene at a given instant of time can be obtained, without scanning, from a single snapshot; this is what makes the CTIS such a potentially powerful tool for spatially, spectrally, and temporally resolved imaging. A CTIS performs poorly in imaging some types of scenes in particular, scenes that contain little spatial or spectral variation. The computed spectra of such scenes tend to approximate correct values to within acceptably small errors near the edges of the field of view but to be poor approximations away from the edges. The additional structure imposed on a scene according to the present method enables the CTIS algorithms to reconstruct acceptable approximations of the spectral data throughout the scene.

Bearman, Gregory H.↗